SunCobalt
(Thomas)
October 30, 2011, 8:48am
1
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
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
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
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