[Proj] MODIS' sinusoidal to latlon

Horacio Samaniego horacio.samaniego at gmail.com
Thu Jan 13 06:09:50 EST 2011

Thank you all...

The URL Markus N sent shows that it is possible to convert LL to
Sinusoidal (better to say that parameters for the conversion are
successfully implemented in GRASS)

However the question on how to convert the MODIS file is not
explicitly answered. So, here is my experience with it!

Problem: I had a MODIS cover (for Southamerica) previously processed
using GDAL in GeoTiff format that needes to be reprojected to some UTM
(19S, as I will only be interested in the Chilean portion of the

Solution: Thanks to all the members who provided answers....

 the simplest answer to this simple problem was: just use the correct
-s_srs definition with gdalwarp!

gdalwarp -s_srs '+proj=sinu +R=6371007.181 +nadgrids=@null +wktext'
-t_srs EPSG:32719 -r cubic  ${INFILE} ${OUTFILE}

special thanks to Vincent Schut,

I hope this wil help others


On Tue, Jan 11, 2011 at 4:54 AM, Markus Neteler <neteler at osgeo.org> wrote:
> Horacio,
> on our Web site, at
>  http://gis.cri.fmach.it/modis-sinusoidal-gis-files/
> you can find a ZIP file with a Sinusoidal/LL SHAPE file
> and a README how the reprojection was done.
> Hopefully helpful for you,
> Markus
> On Tue, Jan 11, 2011 at 4:22 AM, Horacio Samaniego
> <horacio.samaniego at gmail.com> wrote:
>> I am trying to figure out how I can transform some MODIS data.
>> Unfortunately, I have performed some operations on a set of tiles
>> (merge, sum, etc...) using GDAL and now have two Gtif files with yearly
>> EVI for South America. However, I'm not successfully putting my result
>> in a common projection with the rest of my dataset. MRT is only good for
>> hdf files.
>> Any suggestion on how to automate the process of moving from MODIS
>> sinusoidal to old plain lat/lon (wgs84)?
>> An little enhancement on this topic in the forums will help future
>> users... and any help will be greatly appreciate.
>> thanks in advance!
>> Horacio
>> **here's the original MODIS projection definition:
>> PROJCS["unnamed",
>>    GEOGCS["Unknown datum based upon the custom spheroid",
>> DATUM["Not_specified_based_on_custom_spheroid",
>>            SPHEROID["Custom
>> spheroid",6371007.181,0]],
>> PRIMEM["Greenwich",0],
>> UNIT["degree",0.0174532925199433]],
>> PROJECTION["Sinusoidal"],
>> PARAMETER["longitude_of_center",0],
>> PARAMETER["false_easting",0],
>> PARAMETER["false_northing",0],
>> UNIT["metre",1,
>>        AUTHORITY["EPSG","9001"]]]
>> --
>> Horacio Samaniego <horacio.samaniego at gmail.com>
>> _______________________________________________
>> Proj mailing list
>> Proj at lists.maptools.org
>> http://lists.maptools.org/mailman/listinfo/proj
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj

Instituto de Silvicultura
Universidad Austral de Chile
Isla Teja, Valdivia, Chile

+56 (63) 222480 [of.]

More information about the Proj mailing list