Converts MOD04 'swath' data to gridded data, using internal ground control point data.

swath2grid.sh

#!/bin/bash
#
# swath to grid conversion for
# MOD04 reflectance data
#

# get the reprojection information
gdal_translate -of VRT HDF4_EOS:EOS_SWATH:"$1":mod04:Mean_Reflectance_Land land.vrt
gdal_translate -of VRT HDF4_EOS:EOS_SWATH:"$1":mod04:Mean_Reflectance_Ocean ocean.vrt

# delete the bad ground control points
sed '/X=\"-9.990000000000E+02\"/d' land.vrt
sed '/X=\"-9.990000000000E+02\"/d' ocean.vrt

filename=$(basename "$1")
filename="${filename%.*}"

# reproject the data using the cleaned up VRT file
gdalwarp -overwrite -of GTIFF -tps land.vrt $filename.land.reflectance.tif
gdalwarp -overwrite -of GTIFF -tps ocean.vrt $filename.ocean.reflectance.tif

#gdalbuildvrt -separate stack.vrt lsat1.tif lsat2.tif ...
#gdal_translate stack.vrt stack.tif
#gdal_calc.py -A input1.tif -B input2.tif --outfile=result.tif --calc="A+B"