Dear users...
i'm trying to use ogr2ogr to transform a UTM-wgs84 shp file in
Gauss-Boaga Roma40 (italian system) shp file...
I use this:
ogr2ogr -t_srs '+proj=tmerc +lat_0=0.0 +lon_0=15.0 +x_0=2520000 +y_0=0.0
+ellps=intl +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
+pm=rome' prova_GB1.shp prova.shp
but the map I obtain is quite far from the right position (some
kilometers)...
where I'm wrong???
thank you
--
ivan marchesini <marchesini@unipg.it>
On Fri, Jan 21, 2005 at 02:10:52AM +0100, ivan marchesini wrote:
Dear users...
i'm trying to use ogr2ogr to transform a UTM-wgs84 shp file in
Gauss-Boaga Roma40 (italian system) shp file...I use this:
ogr2ogr -t_srs '+proj=tmerc +lat_0=0.0 +lon_0=15.0 +x_0=2520000 +y_0=0.0
+ellps=intl +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68
+pm=rome' prova_GB1.shp prova.shp
Not sure, but isn't
+k=0.999600
missing?
but the map I obtain is quite far from the right position (some
kilometers)...where I'm wrong???
thank you
I have tried it with VMAP0 data:
# VMAP0 LatLong
# cut out a piece: -spat xmin ymin xmax ymax: spatial query extents
ogr2ogr builtupa_bolzano_wgs84.shp -spat 10 45 13 47 builtupa_eurnasia_wgs84.shp
#verify:
ogrinfo -summary builtupa_bolzano_wgs84.shp builtupa_bolzano_wgs84
...
Layer name: builtupa_bolzano_wgs84
Geometry: Polygon
Feature Count: 34
Extent: (10.004833, 45.063251) - (12.678473, 46.682194)
Layer SRS WKT:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]]
id: Integer (10.0)
f_code: String (5.0)
nam: String (64.0)
tile_id: Integer (5.0)
fac_id: Integer (10.0)
# Reprojection from LL to GaussBoaga:
# From source projection (as defined in .prj file):
ogr2ogr builtupa_bolzano_GB1.shp \
-t_srs "+init=epsg:26591 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68" \
builtupa_bolzano_wgs84.shp
#Import into Gauss Boaga Fuso Ovest (GB1):
v.in.ogr builtupa_bolzano_GB1.shp out=vmap0_builtupa_BZ
#doesn't look so bad
Please try again adding +k=0.999600 (and report...).
Markus