Proj4 Definition Frage

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"

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

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

Gruß,
ajoessen

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

Die Frage ist, ob’s stimmt :wink:

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

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