[GRASSLIST:5464] Re: ogr2ogr

[Markus Neteler]

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