Example reproject data with MRT

In this short example you can understand how to concatenate the scripts to obtain a GeoTIFF file for each band of the chosen product using as backend MODIS Reprojection Tools (MRT).

Warning

This example is based on a Linux based system. Please if you use other OS change the paths where data will be saved

Downloading data

For first you need to obtain data, so you need to use modis_download.py

Warning

Remember to register in the NASA portal following the instructions at User and password session

modis_download.py -I -f 2012-12-05 -O -t h28v05,h29v05,h28v04 /tmp

Warning

In this example we are working on Japan extension, so please change the name of tiles according with your region.

User and password are passed through standard input.

In this example we download data for only one day (2012-12-05) using the option “-O”.

Inside /tmp/ directory you will find a file called listfileMOD11A1.005.txt containing the names of files downloaded. The name of file it is related to the product that you download.

Warning

Every time that you download new files of same product it will be overwrite, so if you need it, you should rename the file

Mosaic data

modis_mosaic.py is the script to use.

modis_mosaic.py -m /path/to/mrt/ -o /tmp/outputfile /tmp/listfileMOD11A1.005.txt

Warning

/path/to/mrt/ is the directory where Modis Reprojection Tools is stored

The output of this command are outputfile.hdf and outputfile.hdf.xml inside the directory /tmp. It’s reading the input files contained in listfileMOD11A1.005.txt

Convert data

The last part of the procedure is to convert the mosaic, from HDF format and sinusoidal projection, to GeoTIFF with several projection. You have to use modis_convert.py

modis_convert.py -s '( 1 1 1 1 1 1 1 1 1 1 1 1 )' -m /path/to/mrt/ -o /tmp/finalfile.tif -g 250 /tmp/outputfile.hdf

Extract quality information

If necessary, you can extract specific quality type from the chosen quality layer. In this particular case, we extract the Mandatory QA flag of the daytime temperature. You have to use modis_quality.py

modis_quality.py -p MOD11A1 -l 1 -t 1 /tmp/outputfile.hdf /tmp/mod11a1_daytime_qaflag.tif