need .prj for i guess SREF98 (4073/4074/4075) in Serbia

Communities Deutschland (Germany)
1

Hi,

ich hoffe jemand kann helfen. Versuche seit gestern vergebens das “.prj” File für GeoTIFFs zu finden.
Ich habe einige die “.tif” und die “.tfw” Datei aber leider keine “.prj” Datei.

Der Inhalt der “.tfw” Datei ist folgender:

0.100000
0.00
0.00
-0.100000
7623000.0500
4899499.9500

Das PicLayer-Plugin kalibriert(ungenau!!!) das Bild folgendermaßen:

JOSM PicLayer-Plugin Kalibrierdatendaten
#Mon Aug 01 10:50:49 CEST 2011
SHEARX=0.0
INITIAL_POS_y=5500828.871677955
SCALEY=0.16260835096521864
SCALEX=0.16260835096521864
POSITION_Y=5500717.37611237
INITIAL_POS_X=2508437.53713302
POSITION_X=2509008.775239358
INITIAL_SCALE=61.08115391578207
SHEARY=0.0
ANGLE=0.0

Ich glaube dass es sich hierbei um SREF98 (http://georepository.com/crs_4075/SREF98.html) bezogene Daten handelt. Sicher bin ich mir allerdings nicht. Es könnten auch 4075 oder 4073 sein.
Da die .prj Datei fehlt weigert sich JOSM das Image zu importieren. Auch der Import in den geoserver ist mir bisher nicht gelungen.
Die Genauigkeit von PicLayer genügt leider nicht.
Das Ziel ist es den GeoServer über PostGIS zu versorgen.

Wer kann helfen?

PS: Gibt es ein Plugin für JOSM welches die Koordinaten von SREF98 nach WGS84 (http://georepository.com/transformation_4077/SREF98-to-WGS-84-1.html) umrechnet? Ich habe leider keines gefunden. Würde ich aber auch brauchen!!!

–
happy mapping
lg aus Negotin
Kotoj

2

Hi,

mit dem Plugin “ImportImagePlugin” solltest Du SREF98 (4073/4074/4075) importieren können,
findet sich auf jendenfalls in der srs-Liste dieses Plugins.

Das Plugin “proj” hat auch die SREF98.

Ciao,
Frank

3

Hallo Frank,

danke für die rasche Antwort.

Das ImportImagePlugin kann ohne der .prj Datei nicht verwendet werden, da es sich hierbei um ein Spatial Reference System handelt. Alle dies-bezogenen Versuche schlugen fehl.

Ein Plugin “proj” habe ich nicht gefunden.
Meinst du das Plugin proj4j? Dieses scheint keine tiff öffnen zu können.

–
LG Kotoj

4

Hallo Kotoj,

ja sorry, hab’ Deine Frage zu schnell überlesen: Du hast ja noch gar keine GeoTiffs, Du willst jene ja erst herstellen.

Richtig, ich meinte “proj4j” (‘Ableger’ des proj); das kann “nur” umprojezieren (wie das commandline-tool auch), aber der Datenebene :frowning:

Ciao,
Frank

5

Hallo Frank,

ja genau. Meine Daten sind unvollständig. Die hiesigen Vermessungstechniker verwenden sie, einfach so…
Die verwenden irgendeine Autocat Spezialversion.

Und ich kann mit den Daten nichts anfangen :frowning:

–
LG Kotoj

6

Vielleicht hilft dir das hier:


--- EPSG 4075 : SREF98
---
INSERT INTO "spatial_ref_sys" ("srid","auth_name","auth_srid","srtext","proj4text") VALUES (4075,'EPSG',4075,'GEOGCS["SREF98",DATUM["Serbian_Reference_Network_1998",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","1034"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4075"]]','+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs ');
---

# SREF98
<4075> +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs  <>

4073,SREF98,3534,geocentric,6500,1034,,,,,Geodesy.,,"Delcev et al, FIG Working Week May 2009; http://www.gig.net/pub",OGP,2009/05/11,2009.015,1,0
4074,SREF98,3534,geographic 3D,6423,1034,4073,15592,,,Geodesy.,,"Delcev et al, FIG Working Week May 2009; http://www.gig.net/pub",OGP,2009/05/11,2009.015,1,0
4075,SREF98,3534,geographic 2D,6422,1034,4074,15593,,,Horizontal component of 3D system.,Replaces MGI 1901 (CRS code 3906) in Serbia.,OGP. See 3D CRS for original information source.,OGP,2009/05/11,2009.015,1,0

EDIT: Demnach müsste EPSG 4075 die gleichen Parameter haben wie 3824 und 3889.

gruß,
ajoessen

7

Nachtrag: Geoserver hat EPSG 4075 im Angebot. Damit kannst du die tif-dateien als lokalen WMS-Dienst ausgeben lassen. Anleitung siehe meinen Beitrag Nr. 164 bei
http://forum.openstreetmap.org/viewtopic.php?id=13207
(nur 31466 gegen 4075 tauschen)

Oder mit den obenstehenden proj-Angaben in Quantum GIS ein benutzerdefiniertes Koordinatensystem anlegen.

Gruß,
ajoessen

8

Hallo ajoessen,

danke für die Tipps. Mittlerweile bin ich (mit der Hilfe eines Freundes) darauf gekommen, dass es sich nicht um SREF98 sondern wahrscheinlich um MGI/Balkan-Zone-7 referenzierte Bilder es sich handelt. Leider weicht das Datum vom Standard ab und damit die Photos von der richtigen Position.

Sobald ich ein passendes .prj File habe, werde ich es hier posten!

EDIT: Ach ja, die Vermessungstechniker hierzulande haben von SREF98 noch nie etwas gehört.
Das scheint wohl eher eine Totgeburt zu sein. - meine persönliche Meinung -

Danke und liebe Grüße
Kotoj

9

Hallo,

dieser Link hier war eigentlich für Thomas (SunCobalt) im “Aerowest-Thread” gedacht (er hat aber nicht darauf reagiert :wink: ),
vielleicht ist er hier nützlich:
http://www.spatial-analyst.net/wiki/index.php?title=MGI_/_Balkans_coordinate_systems

Ciao,
Frank

10

Danke Frank. Ich kenne den Link. Aber Kotojs EPSG weicht leider von dem dortigen Standard ab.

Oh sorry Frank. Bei Deinem Link waren die Params richtig. Bei spatialreference.org, wo man auch prj und die in die epsg File gehörenden Parameter abfragen kann ( http://spatialreference.org/ref/epsg/31277/ ) waren sie falsch. In der originalen epsg von proj ist eine noch eine dritte (auch falsche) Variante drin.

11

Wenn es einen konstanten Versatz gibt, könnte man das ja im josm-WMS korrigieren. Deshalb mein tip mit dem Geoserver-WMS.

Oder in Quantum GIS die dort vorhandenen EPSG durchprobieren, oder ein eigenes KBS definieren.

Gruß,
ajoessen

12

Als WMS ist es kein Problem. Ich habe einfach einen Dummy EPSG Code angelegt

MGI / Balkans zone 7 (deprecated)

<31277> +proj=tmerc +lat_0=0 +lon_0=21 +k=0.9999 +x_0=7500000 +y_0=0 +ellps=bessel +datum=hermannskogel +units=m +no_defs <>

MGI / Balkans zone 7 - 7 Param shift

<900916> +proj=tmerc +lat_0=0 +lon_0=21 +k=0.9999 +x_0=7500000 +y_0=0 +ellps=bessel +towgs84=574.02732,170.17492,401.5453,4.88786,-0.66524,-13.24673,6.88933 +units=m +no_defs

Das mit Quantum wollte ich heute abend mal versuchen

13

In Quantum GIS ist 31277 nicht deprecated, steht da mit:


+proj=tmerc +lat_0=0 +lon_0=21 +k=0.9999 +x_0=7500000 +y_0=0 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs

Daneben gibts noch ein 31267 MGI/3-degree Gauss Zone 7 (deprecated)


+proj=tmerc +lat_0=0 +lon_0=21 +k=1 +x_0=7500000 +y_0=0 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs

Vielleicht passt das ja.

Gruß,
ajoessen

14

nee, leider nicht. Mit diesen Parametern “+towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232” hat man recht hohe Abweichungen. Der 7 Parameter Shift des Datums muss damit erfolgen “+towgs84=574.02732,170.17492,401.5453,4.88786,-0.66524,-13.24673,6.88933”

15

In spatialreference.org steht (für 31267):
Deprecated: ProjCRS related to inappropriate projection (coordinate operation) and CS resulting in incorrect scale factor and axes order and abbreviation. Use EPSG:31277 instead.

Ciao,
Frank