<div dir="ltr"><div><div>Done: <a href="https://trac.osgeo.org/gdal/ticket/6225#ticket">https://trac.osgeo.org/gdal/ticket/6225#ticket</a><br><br></div>Thanks once again!!<br></div>Simon<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 19, 2015 at 10:59 AM, 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"><span class="">Le jeudi 19 novembre 2015 10:47:50, Simon Lyngby Kokkendorff a écrit :<br>
> I can confirm also that it works when I set a target srs.<br>
> But must admit that the current behaviour was confusing. I knew how to<br>
> interpret the input data - under the assumption that a geotransform is not<br>
> changed when I do no transformations or interpolations:<br>
> If I know that my data are "points" rather than "pixels" - the first<br>
> "point" is lying half a cellsize downwards and rightwards from the upper<br>
> left corner that GDAL reports. That's what I always assume - not relying on<br>
> other metadata items than the GeoTransform.<br>
><br>
> So I agree that it might be worth a ticket.<br>
<br>
</span>Please file one.<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
> Cheers,<br>
> Simon<br>
><br>
> On Thu, Nov 19, 2015 at 10:29 AM, Ari Jolma <<a href="mailto:ari.jolma@gmail.com">ari.jolma@gmail.com</a>> wrote:<br>
> > On 19.11.2015 11:26, Even Rouault wrote:<br>
> >> GDAL version:<br>
> >> GDAL 1.11.2, released 2015/02/10<br>
> >> Simon,<br>
> >><br>
> >> Regarding the fact that the geotiff file looses the half-pixel shift, it<br>
> >> seems<br>
> >> to be due to the absence of a SRS. If I do "gdal_translate test.asc<br>
> >> test.tif -<br>
> >> a_srs EPSG:32630" for example, then the geotransform I get back is<br>
> >> correct.<br>
> >> This is a defect of the GeoTIFF writer that writes the GeoTIFF Tiepoints<br>
> >> with<br>
> >> the PixelIsPoint convention (since your original data is flagged as<br>
> >> such), but<br>
> >> fails to write GTRasterTypeGeoKey=RasterPixelIsPoint as a geotiff key<br>
> >> when there's no SRS. Could be worth a ticket.<br>
> ><br>
> > My conclusion is the same. I think the key should be set regardless of<br>
> > the existence of a SRS.<br>
> ><br>
> > Ari<br>
> ><br>
> >> Some background:<br>
> >> <a href="https://trac.osgeo.org/gdal/wiki/rfc33_gtiff_pixelispoint" rel="noreferrer" target="_blank">https://trac.osgeo.org/gdal/wiki/rfc33_gtiff_pixelispoint</a><br>
> >><br>
> >> A workaround :<br>
> >> gdal_translate test.asc test.tif -mo AREA_OR_POINT=AREA<br>
> >><br>
> >> This will force the geotiff writer to use the default convention, and<br>
> >> thus avoiding going back to the centerofpixel convention when writing<br>
> >> the GeoTIFF<br>
> >> tiepoints.<br>
> >><br>
> >> Even<br>
> ><br>
> > _______________________________________________<br>
> > gdal-dev mailing list<br>
> > <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
> > <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
</div></div><div class="HOEnZb"><div class="h5">--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</div></div></blockquote></div><br></div>