I am working with the Modis version6 NDVI/EVI 16 day composite dataset (MOD13Q1). I am developing a automated algorithm to download, aggregate, and processes the dataset. Therefore, I do not want to use the MODIS specific tools, but use the gdal library to reproject the dataset to UTM 13N. However, I have found conflicting information regarding the correct projection information. spatialreference.org has several different projects:
And here is the gdalsrsinfo for one of the HDF4 datasets:
gdalsrsinfo 'HDF4_EOS:EOS_GRID:"MOD13Q1.A2017049.h10v04.005.2017066041938.hdf":MODIS_Grid_16DAY_250m_500m_VI:250m 16 days NDVI' PROJ.4 : '+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs ' OGC WKT : 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["Meter",1]]
Can anyone provide some clarification to the correct SRID information?