Thank you, using gdal.reprojectimage(..) brought the expected result.
> Am 02.03.2015 um 18:43 schrieb David Strip <[email protected]>: > > This result is not entirely unexpected, though it depends on the two > projections. A square box in one projection will map onto something entirely > different in many cases, often introducing skew which will change the aspect > ratio and size of the resulting raster. If you view your reprojected dataset > you will almost certainly find a skewed region surrounded by no data values. > > On 3/2/2015 9:01 AM, Robert Oehler wrote: >> >> Hi, >> >> I’m developing a raster library using GDAL. I use the Java bindings. >> >> I would like to do a reprojection of an image which I have in memory. >> >> I’m using an in-memory driver to create a dataset. I fill it with the >> properties of the image and its data. >> >> Now I’m doing >> >> Dataset warped = gdal.AutoCreateWarpedVRT(ds, src_proj.ExportToWkt(), >> dst_proj.ExportToWkt()); >> >> but the resulting raster is very different from the original. (720,720 --> >> 635,810) while the resulting bounding box is correct (compared to my own >> reprojection implementation) >> >> What am I doing wrong? >> >> I’ve seen gdal.reprojectImage(..) What is the difference? Is it preferable >> in my case ? >> >> Thank you >> >> Robert Oehler >> >> >> >> _______________________________________________ >> gdal-dev mailing list >> [email protected] <mailto:[email protected]> >> http://lists.osgeo.org/mailman/listinfo/gdal-dev >> <http://lists.osgeo.org/mailman/listinfo/gdal-dev> > Robert Oehler Tel ( D ) : +49 179 4759609 Tel ( I ) : +39 370 3161263 -------------------------------------------- CEO RoProducts www.roproducts.de
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
