[Proj] Pixel Shifts When Warping

[Yuta Sato]

Dear GDAL and Proj Developers and Users:

I am really frustrated when the MODIS HDF file converted into GeoTiff with
GCS projection never match with the same MODIS file  converted from MODIS
Reprojection tool.

The file I am working is here:

http://e4ftl01.cr.usgs.gov/MOTA/MCD43A4.005/2015.01.01/MCD43A4.A2015001.h26v05.005.2015027224017.hdf

I used the following parameters for gdalwarp:

-t_srs wgs84 -tr 0.004166666666667 0.004166666666667


I used the following parameters for MRT:

OUTPUT_PIXEL_SIZE = 0.004166666666667
DATUM = WGS84
OUTPUT_PROJECTION_TYPE = GEO

The image produced from MRT has:
size = 6025*2400
point-based
ulx, uly = 92.3799071,  40.0000000

The image produced from gdalwarp has:
size = 6027*2400
area based
ulx, uly = 92.3760431,  40.0000000

Do you have any idea on which is more precise: gdalwarp or MRT? Or there is
any way in gdalwarp to get the same image as MRT gives?

Any idea and experience is highly appreciated. Thank you very much in the
advance!

Sincerely yours,
yuta sato
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.maptools.org/pipermail/proj/attachments/20150517/988d5074/attachment.htm