I was doing some tests on Gvsig for raster reprojection, where gdal warp is used. Gvsig uses a default error threshold of 0.125, which AFAIK makes gdal warp fall back to an approximate transformation. Anyway the reprojection is quite slow, even if significantly faster then GDALGenImgProjTransform.<div>
<br></div><div>My test is against a 21698P x 24647L ArcInfo Binary Grid, with a single band with block dimensions 256x4. Bot SRSs are projected.</div><div>The command is: &gt; gdalwarp -s_srs EPSG:3003 -t_srs EPSG:23032 -et 0.125 hdr.adf test.tiff</div>
<div>It took about 9 minutes to complete.</div><div><br></div><div>Loading the same raster on QGis, which implements its own approximate transformation for on-the-fly reprojection [1], is istantaneous.</div><div>I was wondering if a similar algorithm could be use within Gdal...</div>
<div><br></div><div>Giovanni</div><div><br></div><div>[1] <a href="http://trac.osgeo.org/qgis/browser/trunk/qgis/src/core/qgsrasterprojector.cpp">http://trac.osgeo.org/qgis/browser/trunk/qgis/src/core/qgsrasterprojector.cpp</a></div>
<div><br><div><br></div></div>