<div dir="ltr">Dear Even:<br><br>Thanks for your response. I'm sorry to specify the MODIS data I'm using. It's not the swath product the thread is talking that GDAL is not appropriate. Yes, in case of swath product which requires gridding, thanks for your improvement of the gdal_grid, I gonna use of it for swath data.<br><br><br>My question is about the MODIS data which was already gridded. The HDF data has 2400 by 2400 pixels in sinusoidal projection. My concern is during conversion into GCS, what parameters gdalwarp uses?<br>PIXEL_STEP = ?<br>LINE_STEP = ?<br>ULX and ULY = ? (HDF file's topleft pixel's topleft corner or center corner)<br><br>I am using nearest interpolation.<div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jul 19, 2015 at 9:02 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Yuta,<br>
<span class="">><br>
> When I convert a HDF file in Sinusoidal Projection into Geographic<br>
> Coordinate System using gdalwarp.exe, how can I specify the following<br>
> parameters?<br>
><br>
> PIXEL_OFFSET = 0<br>
> LINE_OFFSET = 0<br>
> PIXEL_STEP = 1<br>
> LINE_STEP = 1<br>
<br>
</span>You could override the values in a VRT produced by gdal_translate -of VRT<br>
source_dataset out.vrt<br>
But it is dubious that changing the original values reported by the driver is<br>
going to produce something correct.<br>
<span class=""><br>
><br>
> I want to make sure that there are not any pixel shifts during<br>
> re-projection.<br>
> If gdalwarp.exe is not a good option for it, what are other methods for it?<br>
> I want to reproject a MODIS file.<br>
<br>
</span>There have been recent discussions about why GDAL warping is not appropriate<br>
for datasets that have bow-tie : <a href="http://osgeo-org.1560.x6.nabble.com/gdal-dev-
gdalwarp-problems-with-MODIS-L1B-td5206583.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/gdal-dev-<br>
gdalwarp-problems-with-MODIS-L1B-td5206583.html</a><br>
<br>
There are a few hints in <a href="http://stackoverflow.com/questions/27968890/reproject-
and-mosaic-modis-level2-using-gdal-python" rel="noreferrer" target="_blank">http://stackoverflow.com/questions/27968890/reproject-<br>
and-mosaic-modis-level2-using-gdal-python</a><br>
<br>
At that time, I see you were looking for a fast(er) gdal_grid. Recently I've<br>
added in trunk (post 2.0) a linear interpolation method which runs several<br>
magnitude of order faster than the other methods.<br>
<span class="HOEnZb"><font color="#888888"><br>
Even<br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br></div>