Proj4 Definition Frage

SunCobalt · 2011-10-30T08:48:03.000Z

Weiß jemand, wie ich damit eine *.prj Datei bastele?

projection=HABSB3_75e
*proj=+proj=lcc +lat_0=47 +lon_0=19 +lat_1=45 +lat_2=49 +ellps=bessel +towgs84=600.0,205.0,437.0
*par_projection=WGS84
*par_datum=MGI
pixel_x=6.3346036011529652
pixel_y=-6.3346036011529652
*rotation=0’0"

kellerma · 2011-10-30T17:55:00.000Z

Hi,

ich würde mir auf http://spatialreference.org/ einfach ein paar .prj und pro4 Dateien zusammensuchen
und versuchen, etwas zusammenzustöpseln,
…,PROJECTION[“Lambert_Conformal_Conic”],…

Ciao,
Frank

ajoessen · 2011-10-31T06:39:53.000Z

SunCobalt:

Weiß jemand, wie ich damit eine *.prj Datei bastele?

projection=HABSB3_75e
*proj=+proj=lcc +lat_0=47 +lon_0=19 +lat_1=45 +lat_2=49 +ellps=bessel +towgs84=600.0,205.0,437.0
*par_projection=WGS84
*par_datum=MGI
pixel_x=6.3346036011529652
pixel_y=-6.3346036011529652
*rotation=0’0"

Hast du mal Quantum GIS probiert? Das kann auch selbst definierte SRS verwenden, und legt für shapefiles auch .prj-Dateien an.

Gruß,
ajoessen

kellerma · 2011-10-31T08:44:11.000Z

ajoessen:

Hast du mal Quantum GIS probiert? Das kann auch selbst definierte SRS verwenden, und legt für shapefiles auch .prj-Dateien an.

Cooler Trick
Geht schneller als selber suchen (qgis 1.4.0):

PROJCS[“Lambert_Conformal_Conic”,GEOGCS[“Bessel 1841”,DATUM[“D_unknown”,SPHEROID[“bessel”,6377397.155,299.1528128]],PRIMEM[“Greenwich”,0],
UNIT[“Degree”,0.017453292519943295]],PROJECTION[“Lambert_Conformal_Conic”],PARAMETER[“standard_parallel_1”,45],
PARAMETER[“standard_parallel_2”,49],PARAMETER[“latitude_of_origin”,47],PARAMETER[“central_meridian”,19],
PARAMETER[“false_easting”,0],PARAMETER[“false_northing”,0],UNIT[“Meter”,1]]

Die Frage ist, ob’s stimmt

usm78-gis · 2011-10-31T14:47:35.000Z

kellerma:

Die Frage ist, ob’s stimmt

Natürlich stimmt es nicht, hier fehlt der TOWGS84 Tag.
Bei den meisten GDAL-basierten Transformationen (d.h auch mit qgis) kommen die Datum-Shifts unter die Räder:
http://home.gdal.org/projects/opengis/wktproblems.html

kellerma · 2011-11-04T19:33:01.000Z

So, da hätten wir noch das kleine python-script von [1]:

$ proj2wkt.py '+proj=lcc +lat_0=47 +lon_0=19 +lat_1=45 +lat_2=49 +ellps=bessel +towgs84=600.0,205.0,437.0'

PROJCS["unnamed",GEOGCS["Bessel 1841",
DATUM["unknown",SPHEROID["bessel",6377397.155,299.1528128],TOWGS84[600,205,437,0,0,0,0]],
PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",45],PARAMETER["standard_parallel_2",49],PARAMETER["latitude_of_origin",47],
PARAMETER["central_meridian",19],PARAMETER["false_easting",0],PARAMETER["false_northing",0]]

[1] http://spatialnotes.blogspot.com/2010/11/converting-wkt-projection-info-to-proj4.html