[Mapbender-users] fehlerhafte Darstellung Flurstückssuche

Wed Jun 24 03:47:19 PDT 2020

Hallo Armin,

wenn ich Deinen Test der Transformation durchführe, bekomme ich folgenden Fehler:

HINWEIS:  PostGIS was unable to transform the point because either no grid shift files were found, or the point does not lie within the range for which the grid shift is defined. Refer to the ST_Transform() section of the PostGIS manual for details on how to configure PostGIS to alter this behaviour.

ERROR:  FEHLER:  transform: couldn't project point (2.49e+06 5.652e+06 0): failed to load datum shift file (-38)

SQL state: XX000

Diesen Fehler bekomme ich nur in meiner PostgreSQL 11 Datenbank (mit Postgis 3.0), wenn ich das alles in meiner alten 9.6 DB (mit Postgis 2.3) durchführe, klappt alles.

Die GSB-Datei liegen bei mir wie die anderen Projektionsdateien auch, unter C:\Program Files\PostgreSQL\11\share\contrib\postgis-3.0\proj bzw. C:\Program Files\PostgreSQL\9.6\share\contrib\postgis-2.3\proj

Wie kann es sein, dass das datum shift file in der alten Datenbank angesteuert werden kann und in der neuen nicht? Passen die Versionen nicht zusammen?

Ergebnis SELECT PostGIS_Full_Version();

PostgreSQL 9.6: POSTGIS="2.3.0 r15146" GEOS="3.5.0-CAPI-1.9.0 r4090" SFCGAL="1.3.0" PROJ="Rel. 4.9.1, 04 March 2015" GDAL="GDAL 2.1.1, released 2016/07/07" LIBXML="2.7.8" LIBJSON="0.12" TOPOLOGY RASTER

PostgreSQL 11: POSTGIS="3.0.0 r17983" [EXTENSION] PGSQL="110" GEOS="3.8.0-CAPI-1.13.1 " SFCGAL="1.3.2" PROJ="Rel. 5.2.0, September 15th, 2018" GDAL="GDAL 2.4.3, released 2019/10/28" LIBXML="2.9.9" LIBJSON="0.12" LIBPROTOBUF="1.2.1" WAGYU="0.4.3 (Internal)" TOPOLOGY RASTER

Aber scheinbar wird die Datei nicht gefunden..

Viele Grüße,

Markus

-----Ursprüngliche Nachricht-----
Von: Armin Retterath [mailto:armin.retterath at vermkv.rlp.de] 
Gesendet: Freitag, 19. Juni 2020 07:57
An: Markus Lütke Brintrup <geodaten at wipflerplan.de>; Mapbender_users at lists.osgeo.org
Betreff: Re: AW: [Mapbender-users] fehlerhafte Darstellung Flurstückssuche

hallo markus,

du musst jetzt kontrollieren, ob der beta2007 transformationsansatz sowohl in der epsg datei, als auch in der postgis eingestellt sind.

bei debian kenne ich mich aus - aber mit windows hatte ich 20 jahre nichts mehr zu tun ;-).

der mapserver nutzt die epsg datei für die transformation bei den wms/wfs - da steht ja pro crs ein eintrag drin:

z.b.:

# DHDN / 3-degree Gauss-Kruger zone 2

<31466> +proj=tmerc +lat_0=0 +lon_0=6 +k=1 +x_0=2500000 +y_0=0 

+datum=potsdam +units=m +no_defs  <>

wenn man mit shift werten genau transformieren will, muss man die shift gsb datei angeben d.h. die zeile heißt später:

<31466> +proj=tmerc +lat_0=0 +lon_0=6 +k=1 +x_0=2500000 +y_0=0 

+datum=potsdam +ellps=bessel  <mailto:+nadgrids=@BETA2007.gsb,null> +nadgrids=@BETA2007.gsb,null +units=m 

++no_defs

<> 

die shift datei sollte im ordner liegen, in der auch die epsg datei liegt.

(man kann sie, fall noch nicht vorhanden, vom bkg runterladen: 

 <http://crs.bkg.bund.de/crseu/crs/descrtrans/BeTA/BETA2007.gsb> http://crs.bkg.bund.de/crseu/crs/descrtrans/BeTA/BETA2007.gsb)

für die postgis muss erfolgt die angabe der transformation in der spatial_ref_sys tabelle

für gk2 kannst du das mal checken:

select * from spatial_ref_sys where srid = 31466;

falls der genaue ansatz da nicht hinterlegt ist, kannst du das per sql

anpassen:

update spatial_ref_sys set proj4text = '+proj=tmerc +lat_0=0 +lon_0=6 +k=1 

+x_0=2500000 +y_0=0 +datum=potsdam +ellps=bessel 

 <mailto:+nadgrids=@BETA2007.gsb,null> +nadgrids=@BETA2007.gsb,null +units=m +no_defs' where srid = 31466;

uns eine transformation testen:

SELECT ST_AsText(ST_Transform(ST_GeomFromText('POINT(2490000.00

5652000.00)',31466),25832)) As wgs_geom;

da müssen ungefähr folgende koordinaten rauskommen:

279488.01 5654871.71

der transformationsansatz über die epsg datei - z.b. für mapserver kann über das kommandozeilentool cs2cs getestet werden:

cs2cs +init=epsg:31466 +to +init=epsg:25832 dann eingabe von:

2490000.00 5652000.00

führt auch zu:

279488.01           5654871.71 0.00

die transformationsansätze müssen natürlich für alle gk epsg angepasst werden, die du nutzt.

falls du weitere infos brauchst, frag einfach nach ;-) .

viele grüße

armin

> Markus Lütke Brintrup < <mailto:geodaten at wipflerplan.de> geodaten at wipflerplan.de> hat am 15. Juni 2020 

> um

> 12:09 geschrieben:

> 

> 

> Hallo Armin,

> 

> Mapbender und postgreSQL laufen jeweils auf einem Windows Server 2012.

> 

> Viele Grüße,

> 

> Markus

> 

> 

> 

> -----Ursprüngliche Nachricht-----

> Von: Armin Retterath [ <mailto:armin.retterath at vermkv.rlp.de> mailto:armin.retterath at vermkv.rlp.de]

> Gesendet: Mittwoch, 10. Juni 2020 15:52

> An: Mapbender User List < <mailto:mapbender_users at lists.osgeo.org> mapbender_users at lists.osgeo.org>; Markus 

> Lütke Brintrup < <mailto:geodaten at wipflerplan.de> geodaten at wipflerplan.de>

> Betreff: Re: [Mapbender-users] fehlerhafte Darstellung Flurstückssuche

> 

> hallo markus,

> wahrscheinlich hast du den genauen beta2007 transformationsansatz nur 

> in der epsg datei konfiguriert - daher klappt die transformation durch 

> den mapserver on the fly. in der spatial_ref_sys tabelle der postgis 

> müssen die epsg defs genauso konfiguriert sein, wie in der epsg datei.

> was hast du für ein betriebsystem?

> viele grüsse

> armin

> > Markus Lütke Brintrup < <mailto:geodaten at wipflerplan.de> geodaten at wipflerplan.de> hat am 10. Juni 2020 

> > um

> > 15:24 geschrieben:

> >

> >

> > Hallo Liste,

> >

> >

> >

> > in unserem Mapbender-Projekt haben wir die ALKIS-Daten aus der 

> > postgreSQL als Mapserverdatei eingebunden. Der Dienst wird 

> > lagerichtig angezeigt, obwohl unser Mapbender noch in GK4 ist und 

> > der Dienst in UTM32.

> >

> > Dann habe ich in der postgreSQL einen View gebastelt, den ich als 

> > Flurstücks-Suche im Mapbender integriert habe. So weit, so gut.

> >

> > Ergebnis ist folgendes:

> >

> >

> >

> >

> >

> >

> >

> > Die Lage ist verschoben, d.h. die Transformation des 

> > Koordinatensystem hat augenscheinlich nicht funktioniert.

> >

> >

> >

> > Mein View sieht so aus (ich habe versucht die Geometriespalte zu

> > transformieren):

> >

> >

> >

> > CREATE OR REPLACE VIEW postnas.flurstueckssuche

> >

> > AS

> >

> > SELECT ax_flurstueck.gid,

> >

> > st_transform(ax_flurstueck.wkb_geometry,

> > 31468)::geometry(Geometry,31468) AS geom,

> >

> > gem_shl.gemname AS gemeindename,

> >

> > gema_shl.gemarkung AS gemarkungsname,

> >

> > ax_flurstueck.zaehler,

> >

> > ax_flurstueck.nenner

> >

> > FROM postnas.ax_flurstueck

> >

> > JOIN postnas.gema_shl ON

> > "substring"(ax_flurstueck.flurstueckskennzeichen::text, 1, 6) = 

> > gema_shl.gemashl::text

> >

> > JOIN postnas.gem_shl ON gem_shl.gemshl = gema_shl.gemshl

> >

> > WHERE gema_shl.gemshl IS NOT NULL

> >

> > ORDER BY gema_shl.gemarkung, (ax_flurstueck.zaehler::integer),

> > (ax_flurstueck.nenner::integer) NULLS FIRST;

> >

> >

> >

> > Freundliche Grüße

> >

> >

> >

> > Markus Lütke Brintrup

> >

> >

> >

> > WipflerPLAN

> >

> > Niederlassung Donauries

> >

> >

> >

> > Telefon: 09081 27509-44

> >

> > Telefax: 09081 27509-50

> >

> > < <mailto:mlb at wipflerplan.de> mailto:mlb at wipflerplan.de>  <mailto:mlb at wipflerplan.de> mlb at wipflerplan.de

> >

> > < <http://www.wipflerplan.de/> http://www.wipflerplan.de/>  <http://www.wipflerplan.de> www.wipflerplan.de

> >

> >

> >

> >

> >

> > WipflerPLAN Planungsgesellschaft mbH

> >

> > Standort Nördlingen

> >

> > An der Lach 11 a

> >

> > 86720 Nördlingen

> >

> > Amtsgericht Ingolstadt HRB 190595

> >

> > Geschäftsführer:

> >

> > Wilhelm Wipfler

> >

> > Ursula Burkart

> >

> > Adelheid Buchard

> >

> >

> >

> >

> >

> >

> >

> > _______________________________________________

> > Mapbender_users mailing list

> >  <mailto:Mapbender_users at lists.osgeo.org> Mapbender_users at lists.osgeo.org

> >  <https://lists.osgeo.org/mailman/listinfo/mapbender_users> https://lists.osgeo.org/mailman/listinfo/mapbender_users

> 

--

Armin Retterath

Zentrale Stelle GDI-RP

LVermGeo Rheinland-Pfalz

Von-Kuhl-Straße 49

56070 Koblenz

Telefon 0261 492-466

Telefax 0261 492-492

 <mailto:armin.retterath at vermkv.rlp.de> armin.retterath at vermkv.rlp.de

 <http://www.geoportal.rlp.de> http://www.geoportal.rlp.de

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapbender_users/attachments/20200624/c97b3041/attachment-0001.html>