<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<p>gdalinfo for NTv2LSBB_LSA.gsb gave the following information:</p>
<p>Size is 256, 264<br>
Coordinate System is: EPSG 4326<br>
Origin: 10.5, 53.1<br>
SYSTEM_F: S4283<br>
SYSTEM_T: ETRS89<br>
Upper Left 10.5, 53.1<br>
Lower Left 10.5, 50.9<br>
Upper Right 13.3, 53.1<br>
Lower Right 13.3, 50.9<br>
</p>
<p>ogrinfo for punkte_31468.dxf gave the following information:</p>
<p>Geometry: Unknown (any)<br>
Feature Count: 36<br>
Extent: (4523806.5, 5758974.7) - (4524134.7, 5759234.0)<br>
Layer SRS WKT: (unknown)</p>
<p>according to QGIS punkte_31468.dxf should be at 11.6 52.1 which
is inside the extent of the NTv2LSBB_LSA.gsb. after checking the
documentation for the NTv2LSBB_LSA.gsb i noticed i couldn't use it
for transformations from 31468 to 25833, that's why i tried it
with BETA2007.gsb instead, but it didn't work either - it gave the
same errors.</p>
<p>gdalinfo for BETA2007.gsb gave the following information:</p>
<p>Size is 62, 84<br>
Coordinate System is: EPSG 4326<br>
Origin: 5.4, 55.4<br>
SYSTEM_F: DHDN90<br>
SYSTEM_T: ETRS89<br>
Upper Left 5.4, 55.4<br>
Lower Left 5.4, 47.0<br>
Upper Right 15.8, 55.4<br>
Lower Right 15.6, 47.0</p>
<p>my adjusted code for BETA2007.gsb looks like:</p>
<p>ogr2ogr -s_srs "EPSG:31468" -t_srs "EPSG:25833" -dim XYZ -ct
"+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12 +k=1
+x_0=4500000 +y_0=0 +ellps=bessel +step +proj=hgridshift
+grids=BETA2007.gsb +step +proj=utm +zone=33 +ellps=GRS80" -f DXF
punkte_25833.dxf punkte_31468.dxf</p>
<p>the points i used to test this were digitized around Magdeburg,
Saxony-Anhalt in EPSG:31468. according to aerial data and
openstreetmap, the position and coordinate system are correct when
i check it in QGIS and switch to different coordinate systems.<br>
</p>
<div class="moz-cite-prefix">Am 06.02.2021 um 22:24 schrieb Markus
Metz:<br>
</div>
<blockquote type="cite"
cite="mid:CAG+h=FHVrLcbZ3Gtz0jXaV0V8mCEKKOxSjd7Msa0PbToGEDVpQ@mail.gmail.com">
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
<div dir="ltr"><br>
<br>
On Fri, Feb 5, 2021 at 3:40 PM Karsten Tebling <<a
href="mailto:tebling@masuch.de" moz-do-not-send="true">tebling@masuch.de</a>>
wrote:<br>
><br>
> Hello,<br>
><br>
> i'm struggling with how to use the -ct option in ogr2ogr
and could need<br>
> some help.<br>
><br>
> i copied a NTv2.gsb to osgeo4w64/share/proj and edited the
proj.db by<br>
> inserting a new row in the grid_transformation-table. using
this<br>
> transformation in qgis works. however when i try to use it
in ogr2ogr<br>
> with the same data i always get "failed to reproject
feature (geometry<br>
> probably out of source or destination srs). no known way to
write<br>
> feature with geometry 'none'". here is my code:<br>
><br>
> ogr2ogr -s_srs "EPSG:31468" -t_srs "EPSG:25833" -dim XYZ
-ct<br>
> "+proj=pipeline +step +inv +proj=tmerc +lat_0=0 +lon_0=12
+k=1<br>
> +x_0=4500000 +y_0=0 +ellps=bessel +step +proj=hgridshift<br>
> +grids=NTv2LSBB_LSA.gsb +step +proj=utm +zone=33
+ellps=GRS80" -f DXF<br>
> punkte_25833.dxf punkte_31468.dxf<br>
><br>
<div><br>
</div>
<div>Is the grid NTv2LSBB_LSA.gsb used for hgridshift really
covering the full extent of punkte_31468.dxf ?</div>
<div>You can check with "gdalinfo NTv2LSBB_LSA.gsb" and "ogrinfo
-so punkte_31468.dxf". </div>
<div><br>
</div>
<div>Markus M<br>
</div>
<div><br>
</div>
<div>> i'm using the latest qgis-ltr and osgeo4w shell</div>
><br>
> greetings<br>
> _______________________________________________<br>
> gdal-dev mailing list<br>
> <a href="mailto:gdal-dev@lists.osgeo.org"
moz-do-not-send="true">gdal-dev@lists.osgeo.org</a><br>
> <a
href="https://lists.osgeo.org/mailman/listinfo/gdal-dev"
moz-do-not-send="true">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a></div>
</blockquote>
</body>
</html>