Sidebar

spatial-ecology.org

Trainings:

Learn:

Data:

Community:
teachers
students
projects
Matera 2015
Vancouver 2015
Santa Barbara 2015
Site design
Hands on:
Installations

Donations
USD

GBP

wiki:exercise5a

Exercise 5a: Find thresholds for cloud mask

In this exercise, the following tools will be used:

• pktools: pkascii2ogr, pkextract, pkstatogr, pkdumpogr
• bash: cat, redirection

Create a directory for this exercise

mkdir -p ${OUTPUTDIR}/exercise5 Landsat bands: • band 1: blue (0.45 - 0.515 \mu m) • band 2: green (0.525 - 0.605 \mu m) • band 3: red (0.63 - 0.69 \mu m) • band 4: nir (0.75 - 0.90 \mu m) • band 5: swir 1 (1.55 - 1.75 \mu m) • band 6 (6a): tir 1 (10.40 - 12.5 \mu m) • band 7 (6b): tir 2 • band 8 (7): swir 2 (2.09 - 2.35 \mu m) • band9: pan-chromatic Find the geographic coordinates for cloud pixels cat >${OUTPUTDIR}/exercise5/clouds.txt <<EOF
cloud 3222350 3123150
cloud 3123150 3123635
cloud 3221195 3123750
cloud 3219120 3122590
EOF

Find the geographic coordinates for non cloud pixels

cat > ${OUTPUTDIR}/exercise5/clouds.txt <<EOF nocloud 3269100 3135160 nocloud 3269150 3134935 nocloud 3272665 3135925 nocloud 3272460 3135730 nocloud 3273560 3136620 EOF Create vector file for clouds pkascii2ogr -i${OUTPUTDIR}/exercise5/clouds.txt -x 1 -y 2 -n label --ot String -a_srs epsg:3035 -o ${OUTPUTDIR}/exercise5/clouds.shp Create vector file for cloud free pixels pkascii2ogr -i${OUTPUTDIR}/exercise5/noclouds.txt -x 1 -y 2 -n label --ot String -a_srs epsg:3035 -o ${OUTPUTDIR}/exercise5/noclouds.shp Extract LANDSAT cloud pixels pkextract -i${LANDSATDIR}/${LANDSATIMG} -s${OUTPUTDIR}/exercise5/clouds.shp -o ${OUTPUTDIR}/exercise5/clouds_landsat.shp  Extract LANDSAT cloud free pixels pkextract -i${LANDSATDIR}/${LANDSATIMG} -s${OUTPUTDIR}/exercise5/noclouds.shp -o ${OUTPUTDIR}/exercise5/noclouds_landsat.shp  Calculate mean and stdev cloud pixels pkstatogr -i${OUTPUTDIR}/exercise5/clouds_landsat.shp -n B0 -n B1 -n B2 -n B3 -n B4 -n B5 -n B6 -n B7 -mean -stdev| awk '{print $4,$6}' > ${OUTPUTDIR}/exercise5/clouds_landsat_mean_stdev.txt Calculate mean and stdev nocloud pixels pkstatogr -i${OUTPUTDIR}/exercise5/noclouds_landsat.shp -n B0 -n B1 -n B2 -n B3 -n B4 -n B5 -n B6 -n B7 -mean -stdev| awk '{print $4,$6}' > ${OUTPUTDIR}/exercise5/noclouds_landsat_mean_stdev.txt Dump ogr to ascii text file pkdumpogr -i${OUTPUTDIR}/exercise5/clouds_landsat.shp -n label -n B0 -n B1 -n B2 -n B3 -n B4 -n B5 -n B6 -n B7 > ${OUTPUTDIR}/exercise5/clouds_landsat.txt pkdumpogr -i${OUTPUTDIR}/exercise5/noclouds_landsat.shp -n label -n B0 -n B1 -n B2 -n B3 -n B4 -n B5 -n B6 -n B7 > \${OUTPUTDIR}/exercise5/noclouds_landsat.txt