[gdal-dev] MOD13C2 assign coord. sist.

Eloi eloi at openmailbox.org
Fri Sep 6 03:33:33 PDT 2013


-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

On 09/06/2013 11:53 AM, Jose Gomez-Dans wrote:
> Hi Eloi,
> 
> I think that the projection for the global CMG MOD13 products is
> lat/long, and is stored in the dataset itself (I deleted the
> MOD13's last weeek, so can't check ;-D). You shouldn't need to
> specify it, GDAL will pick it up. In case it doesn't pick it up,
> you want "EPSG:4326", but I'm sure it will do.
> 
> Hope that helps, Jose
> 
> 
> On 6 September 2013 10:36, Eloi <eloi at openmailbox.org> wrote:
> 
> Hi all,
> 
> I am trying so assign the coordinate system to MODIS product
> MOD13C2 but do not know what string use to do this.
> 
> gdal_translate -a_srs 'do not know what to put here' -of GTiff
> 
> 'HDF4_EOS:EOS_GRID:"MOD13C2.A2000032.005.2006272104028.hdf":MOD_Grid_monthly_CMG_VI:CMG
>
> 
0.05 Deg Monthly NDVI' ndvi_sinu.tif
> 
> Thanks in advance,
> 
> Eloi
>> _______________________________________________ gdal-dev mailing
>> list gdal-dev at lists.osgeo.org 
>> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>> 
> 

Hi again,

Jose, thanks for your quick reply.

Assigning EPSG:4326 nothing comes up when loading the image with
others data in EPSG:4326.

gdal_translate -a_srs EPSG:4326 -of GTiff
'HDF4_EOS:EOS_GRID:"MOD13C2.A2000032.005.2006272104028.hdf":MOD_Grid_monthly_CMG_VI:CMG
0.05 Deg Monthly NDVI' ndvi_4326.tif

If I try to use gdalwarp and put EPSG:4326 as target I get following
error.
gdalwarp -t_srs EPSG:4326 -of GTiff -r average
'HDF4_EOS:EOS_GRID:"MOD13C2.A2000032.005.2006272104028.hdf":MOD_Grid_monthly_CMG_VI:CMG
0.05 Deg Monthly NDVI' ndvi_4326.tif
ERROR 1: No PROJ.4 translation for source SRS, coordinate
transformation initialization has failed.
Creating output file that is -2147483648P x -2147483648L.
ERROR 1: Attempt to create -2147483648x-2147483648 dataset is
illegal,sizes must be larger than zero.

The only way it works is using a_ullr -180 90 180 -90 but I am not
sure if this is the correct way of doing it.

gdal_translate -a_ullr -180 90 180 -90 -a_srs EPSG:4326 -of GTiff
'HDF4_EOS:EOS_GRID:"MOD13C2.A2000032.005.2006272104028.hdf":MOD_Grid_monthly_CMG_VI:CMG
0.05 Deg Monthly NDVI' ndvi_4326.tif

Cheers,

Eloi
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.11 (GNU/Linux)
Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/

iQEcBAEBAgAGBQJSKa99AAoJEP9qulIE6g+yE+kH/0dHqNokOfnSBtFMRV9gpy9n
N6PYI1qDVv+p08qC2kx6bvNa6Ae6PepPJmWebKCB8IepAsc3YBWYBrjrWVFkWslO
1cj4FthxhMR6KrPzLzYXCjmgy3SUBnYB5aVFu0SJR/+KMLkg/kWoW5SC4RQJQ8Tj
kicfT0fAm/aVvvX+USX8PdXTjEyGtR9uqwIOqmp6xgzwPM1A/QPqvUKKIZL5ujkZ
sta7BFmlpqHtsh3gkNxIfpy/gkaYUlpn/BJMy3OywjvB8BWO4j+Vvut9OM7wSQUS
PibSNY7/6CtWoUxjuOwDYfFHKMisxQebRyVD9uqXehzn1CKzjjkojWhJGdttc/8=
=l5A9
-----END PGP SIGNATURE-----


More information about the gdal-dev mailing list