The goal is to perform automatically a SDS extraction from a MODIS product (here a LEVEL1B image) for several defined areas (here AERONET photometer locations)
verbose
clear_data
set_home_directory /users/louis
set_tmp_directory /tmp
set_map_directory /users/louis/MAPS/
set_output_directory /users/louis/OUTPUTS
set_input_directory /users/louis/RIEDI/
set_input_hdf_file MOD021KM.A2003040.1810.004.2003041082614.hdf
# OUTPUT RECORD DESCRIPTION
pixel_dump_set linefeed \
SDSName="Height" format="%.0f" \
SDSName="SensorZenith" format="%.2f" \
SDSName="SensorAzimuth" format="%.2f" \
SDSName="SolarZenith" format="%.2f" \
SDSName="SolarAzimuth" format="%.2f" \
linefeed \
SDSName="EV_1KM_RefSB" scaling=Reflectance format="%.2f" DataNb=yes \
linefeed \
SDSName="EV_1KM_Emissive" scaling=Radiance format="%.2f" DataNb=yes \
linefeed \
SDSName="EV_250_Aggr1km_RefSB" scaling=Reflectance format="%.2f" DataNb=yes \
linefeed \
SDSName="EV_500_Aggr1km_RefSB" scaling=Reflectance format="%.2f" DataNb=yes
# PERFORM EXTRACTIONS
# EXTRACTION FOR A SINGLE POINT
pixel_dump Annotation="Maricopa pixel" FileName= "AERONET" Point=33.071,-111.972
pixel_dump Annotation="La Jolla" Point=32.5,-117.16
# EXTRACTION FOR A SMALL AREA
pixel_dump Annotation="Maricopa area" LatitudeMin=33.06 LatitudeMax=33.08 LongitudeMin=-112.0 LongitudeMax=-111.9
# SHOW OUTPUT RESULT USER SELECTED FORMAT
scaling: None : "Box number"
scaling: None : "Raw"
scaling: None : "Column"
scaling: None : "Latitude"
scaling: None : "Longitude"
LineFeed
scaling: Default : "Height" format=" %.0f"
scaling: Default : "SensorZenith" format=" %.2f"
scaling: Default : "SensorAzimuth" format=" %.2f"
scaling: Default : "SolarZenith" format=" %.2f"
scaling: Default : "SolarAzimuth" format=" %.2f"
LineFeed
scaling: Reflectance : "EV_1KM_RefSB" format=" %.2f" DataNb=Yes
LineFeed
scaling: Radiance : "EV_1KM_Emissive" format=" %.2f" DataNb=Yes
LineFeed
scaling: Reflectance : "EV_250_Aggr1km_RefSB" format=" %.2f" DataNb=Yes
LineFeed
scaling: Reflectance : "EV_500_Aggr1km_RefSB" format=" %.2f" DataNb=Yes
# OUTPUT FOR SINGLE POINTS
Maricopa pixel 2003/40 18:10 1 121 460 33.0761 -111.9826
338 19.51 99.76 52.49 151.61
15 0.13 0.12 0.10 0.10 0.10 0.15 0.11 0.15 0.08 0.15 0.19 0.15 0.09 0.12 0.01
16 0.47 0.53 0.52 0.49 0.14 0.40 1.09 2.60 6.94 4.73 7.55 7.09 4.65 3.62 3.16 2.41
2 0.11 0.17
5 0.11 0.10 0.20 0.18 0.14
La Jolla 2003/40 18:10 2 244 149 32.5095 -117.1636
-33 49.34 96.26 54.30 145.85
15 0.20 0.18 0.15 0.13 0.24 0.15 0.11 0.15 0.08 0.15 0.19 0.09 0.08 0.09 0.06
16 0.20 0.22 0.23 0.22 0.07 0.20 0.83 1.78 4.24 3.08 4.78 4.52 3.21 2.72 2.49 2.07
2 0.10 0.10
5 0.16 0.12 0.09 0.07 0.05
# OUTPUT FOR THE SMALL AREA
Maricopa area 2003/40 18:10 3 119 463 33.0794 -111.9335
332 19.50 98.41 52.52 151.64
15 0.13 0.12 0.11 0.10 0.10 0.15 0.11 0.15 0.08 0.15 0.19 0.13 0.08 0.10 0.01
16 0.44 0.51 0.47 0.44 0.13 0.38 1.15 2.60 6.56 4.41 7.27 6.92 4.56 3.65 3.12 2.34
2 0.11 0.15
5 0.11 0.10 0.17 0.16 0.13
Maricopa area 2003/40 18:10 3 119 464 33.0779 -111.9216
332 19.50 98.41 52.52 151.64
15 0.13 0.12 0.10 0.10 0.10 0.15 0.11 0.15 0.08 0.15 0.19 0.13 0.08 0.09 0.01
16 0.42 0.45 0.45 0.43 0.13 0.38 1.14 2.57 6.48 4.37 7.18 6.82 4.51 3.63 3.10 2.36
2 0.11 0.14
5 0.11 0.10 0.16 0.15 0.13
Maricopa area 2003/40 18:10 3 119 465 33.0763 -111.9098
354 19.05 98.44 52.49 151.70
15 0.14 0.13 0.12 0.11 0.11 0.15 0.11 0.15 0.08 0.15 0.19 0.14 0.09 0.11 0.02
16 0.41 0.41 0.44 0.41 0.12 0.35 1.09 2.41 6.09 4.19 6.72 6.33 4.26 3.46 3.00 2.33
2 0.12 0.16
5 0.12 0.11 0.18 0.17 0.14
..............................
..............................
..............................
..............................
Maricopa area 2003/40 18:10 3 122 463 33.0621 -111.9486
338 19.51 99.76 52.49 151.61
15 0.11 0.09 0.07 0.07 0.06 0.06 0.11 0.06 0.08 0.15 0.19 0.11 0.06 0.08 0.01
16 0.33 0.42 0.42 0.42 0.14 0.41 1.16 2.70 6.92 4.54 7.63 7.21 4.79 3.80 3.27 2.45
2 0.06 0.12
5 0.08 0.06 0.12 0.10 0.06
Maricopa area 2003/40 18:10 3 122 464 33.0606 -111.9367
338 19.51 99.76 52.49 151.61
15 0.11 0.09 0.07 0.06 0.06 0.06 0.11 0.06 0.08 0.15 0.19 0.09 0.05 0.07 0.01
16 0.34 0.44 0.41 0.41 0.14 0.40 1.15 2.67 6.81 4.50 7.53 7.16 4.74 3.77 3.27 2.41
2 0.06 0.10
5 0.08 0.06 0.10 0.09 0.06