[gdal-dev] MODIS reprojection

[laura0]

Hi,
I'm working with GDAL tool in order to use MODIS (MOD14.A2 and MYD14.A2)
data but I have problems with reprojection in UTM 32. The lat/lon pixel
coordinates I get are not correct (similar but a bit shifted). 
This is what I have done:
In order to geolocate swaths, there is a separate product that contains two
arrays, Latitude and Longitude. Each row, column in these arrays contains
the lat and lon for each corresponding data pixel. I have used gdal to get
from this product, the following information:
Coordinate System is `'
GCP[  0]: Id=, Info=
          (0.5,0.5) -> (36.5436744689941,25.0500869750977,0)
GCP[  1]: Id=, Info=
          (123.5,0.5) -> (32.4521675109863,24.7417182922363,0)
GCP[  2]: Id=, Info=
....
....
GCP[143]

and then I added ALL these rows as option to the command gdal_translate -of
Gtiff -co "tfw=yes" -gcp 0.5 0.5 36.5436744689941 25.0500869750977 -gcp
123.5 0.5 32.4521675109863 24.7417182922363 ...
MOD14.A2013078.2035.005.NRT.hdf_001_0.gaiag.tif temp.tif

with gdalwarp -t_srs "+proj=utm +zone=32 +datum=WGS84" temp.tif output.tif I
get an image that looks ok but the lat/lon coordinates are still not precise
(particularly for longitude).

Can someone help me?
Thanks

Laura



--
View this message in context: http://osgeo-org.1560.x6.nabble.com/MODIS-reprojection-tp5047700.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.