Using a command file we will extract several parameters but compare only the .6 band of the 3 sensors. (Msphinx command.com)
The area of projection will be
Latitudes: 3.0 to 43.0 degrees (4400 pixel size over 1km2)
Logitudes: -23.0 to 3.0 degrees (3000 pixel size over 1km2)
At the end of the process Msphinx data planes content:
plane 1: Meteosat .6 reprojected band
plane 2: Meteosat .8 reprojected band
plane 3: Meteosat 1.6 reprojected band
plane 4: MODIS 250m(1km) .620 reprojected band
plane 5: MODIS 250m(1km) .841 reprojected band
plane 6: MODIS 500m(1km) 1.23 reprojected band
plane 7: POLDER param 29 670 nm (P) .670 reprojected
plane 8: POLDER param 32 865 nm (P) .865 reprojected
plane 9: POLDER param 33 910 nm .910 reprojected
plane 10: plane 7 shifted
plane 11: plane 1 shifted
plane 4,10,11 are the final result ( .6 band for the 3 sensors)
Planes 100 to 108 intermediary scratch data plane
verbose
send_set_verbose
set_project ProjectName="MULTI_SENSORS" DataPlanesOnDisk=Yes
set_workpace_geometry PlaneWidth=10000 PlaneHeight=10000 DataType=0 CellWidth=10000 CellHeight=10000
send_destroy_all_data
set_view_planes ViewPlane=0 red=1 green=2 blue=3
# START METLOOK PROCESS AND FILL PLANES 100 to 102
send_METLook <<
verbose
clear_data
set_input_directory /Users/louis/METEOSAT_MODIS/
set_input_met_file MSG1-SEVI-MSG15-0008-NA-20030607111234.950000000Z-1941.nat
set_projection_to_geometry ProjectionTo="SINUSOIDAL"
select_output_bands 1 2 3
export_met_planes plane=100 x=0 y=0 width=3712 height=3712 scaling=yes
send_exit
>>
# SET FINAL PROJECTION PARAMETERS (LINEAR PROJECTION)
set_projection_to_geometry XTo=0 YTo=0 WidthTo=3000 HeightTo=4400 ProjectionTo="LINEAR"\
LatitudeMinTo=3 LatitudeMaxTo=43 \
LongitudeMinTo=-23 LongitudeMaxTo=3
# REPROJECT PLANE 100 to 102 to 1 to 3 (INPUT SINUSOIDAL OUTPUT LINEAR)
send_projection PixelOverlap=No PlaneTo=1 PlaneFrom=100 \
XFrom=0 YFrom=0 WidthFrom=0 HeightFrom=0 ProjectionFrom=""
send_projection PixelOverlap=No PlaneTo=2 PlaneFrom=101 \
XFrom=0 YFrom=0 WidthFrom=0 HeightFrom=0 ProjectionFrom=""
send_projection PixelOverlap=No PlaneTo=3 PlaneFrom=102 \
XFrom=0 YFrom=0 WidthFrom=0 HeightFrom=0 ProjectionFrom=""
# START HDFLOOK PROCESS AND FILL PLANES 4 to 6 IN THE FINAL PROJECTION
send_hdflook <<
verbose
clear_data
set_input_directory /Users/louis/METEOSAT_MODIS/
set_projection_to_geometry ProjectionTo="LINEAR" \
LatitudeMinTo=3 LatitudeMaxTo=43 LongitudeMinTo=-23 LongitudeMaxTo=3
set_input_hdf_file /Users/louis/METEOSAT_MODIS/MYD021KM.A20*
select_SDS SDSName="EV_250_Aggr1km_RefSB" index=1,height,width
create_MODIS_SDS_Mosaic plane=4 WidthTo=3000 HeightTo=4400 scaling="Radiance"
set_input_hdf_file /Users/louis/METEOSAT_MODIS/MYD021KM.A20*
select_SDS SDSName="EV_250_Aggr1km_RefSB" index=2,height,width
create_MODIS_SDS_Mosaic plane=5 WidthTo=3000 HeightTo=4400 scaling="Radiance"
set_input_hdf_file /Users/louis/METEOSAT_MODIS/MYD021KM.A20*
select_SDS SDSName="EV_500_Aggr1km_RefSB" index=3,height,width
create_MODIS_SDS_Mosaic plane=6 WidthTo=3000 HeightTo=4400 scaling="Radiance"
send_exit
>>
# START POLDER PROCESS AND FILL PLANES 103 to 108 IN POLDER GEOMETRY
send_Mpol <<
verbose
clear_data
set_input_directory /Users/louis/METEOSAT_MODIS
set_input_MPOL_file P2L1TBG1045036CD
select_viewing_directions 7
select_parameters_by_name "Normalised Radiance at 670 nm (P)" \
"Normalised Radiance at 865 nm (P)" \
"Normalised Radiance at 910 nm"
export_MPOL_Parameters plane=103 x=0 y=0 scaling=yes
set_input_MPOL_file P2L1TBG1045037CD
select_viewing_directions 7
select_parameters_by_name "Normalised Radiance at 670 nm (P)" \
"Normalised Radiance at 865 nm (P)" \
"Normalised Radiance at 910 nm"
export_MPOL_Parameters plane=106 x=0 y=0 scaling=yes
send_exit
>>
# BUILD A REPROJECTED MOSAIC USING PLANE 103 to 108 TO 4 TO 6 (INPUT POLDER OUTPUT LINEAR)
send_projection PixelOverlap=No PlaneTo=7 PlaneFrom=103 \
XFrom=0 YFrom=0 WidthFrom=6480 HeightFrom=3240 ProjectionFrom="POLDER"
send_projection PixelOverlap=No PlaneTo=8 PlaneFrom=104 \
XFrom=0 YFrom=0 WidthFrom=6480 HeightFrom=3240 ProjectionFrom="POLDER"
send_projection PixelOverlap=No PlaneTo=9 PlaneFrom=105 \
XFrom=0 YFrom=0 WidthFrom=6480 HeightFrom=3240 ProjectionFrom="POLDER"
send_projection PixelOverlap=Yes PlaneTo=7 PlaneFrom=106 \
XFrom=0 YFrom=0 WidthFrom=6480 HeightFrom=3240 ProjectionFrom="POLDER"
send_projection PixelOverlap=Yes PlaneTo=8 PlaneFrom=107 \
XFrom=0 YFrom=0 WidthFrom=6480 HeightFrom=3240 ProjectionFrom="POLDER"
send_projection PixelOverlap=Yes PlaneTo=9 PlaneFrom=108 \
XFrom=0 YFrom=0 WidthFrom=6480 HeightFrom=3240 ProjectionFrom="POLDER"
# SHIFT POLDER AND METEOSAT REPROJECTED PLANES TO MATCH MODIS RESULTS ( A SMALL SHIFT IS REMAINING)
# THE SHIFT IS COMPUTED USING THE MSPHINX COST CONTOUR
send_copy_plane FromMask=no PlaneFrom= 7 XFrom=5 YFrom=5 WidthFrom=3000 HeightFrom=4400 \
ToMask=no PlaneTo=10 XTo=0 YTo=0 WidthTo=3000 HeightTo=4400
send_copy_plane FromMask=no PlaneFrom= 1 XFrom=12 YFrom=12 WidthFrom=3000 HeightFrom=4400 \
ToMask=no PlaneTo=11 XTo=0 YTo=0 WidthTo=3000 HeightTo=4400
# SELECT THE RESULT PLANES BAND .6 of THE SENSORS: PLANE 4 MODIS, PLANE 10 POLDER, PLANE 11 METEOSAT
set_view_planes ViewPlane=0 red=4 green=10 blue=11
# THE IMAGE RESULT MUST BE IN GREY LEVELS ALL REMAINIG COLORS ARE FROM SHIFTED CLOUDS OR PIXEL NOT MATCHING
send_enhance ViewPlane=0 X=0 Y=0 Width=3000 Height=4400 Mask=no A=0 B=0 log=no
THE IMAGE RESULT MUST BE IN GREY LEVELS
ALL REMAINING COLORS ARE FROM SHIFTED CLOUDS OR PIXEL NOT MATCHING
| |