[GRASS-user] Import CH1903 Referenced Geodata

[Massimiliano Cannata]

Hi Adrian,
the data are in CH1903/LV03 that correspond to the EPSG:21781.
The .prj produced with ArcMap is somehow not compatible with the right EPSG.
I usually create my location with that EPSG code, and then import with 
the override projection flag.

See below the WTK differences from ESRI AND OGC formats....

ESRI WTK:

PROJCS["CH1903 / LV03",
	GEOGCS["CH1903",
		DATUM["D_CH1903",
			SPHEROID["Bessel_1841",6377397.155,299.1528128]],
			PRIMEM["Greenwich",0],
			UNIT["Degree",0.017453292519943295]],
		PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],
		PARAMETER["latitude_of_center",46.95240555555556],
		PARAMETER["longitude_of_center",7.439583333333333],
		PARAMETER["azimuth",90],PARAMETER["rectified_grid_angle",90],
		PARAMETER["scale_factor",1],PARAMETER["false_easting",600000],
		PARAMETER["false_northing",200000],
		UNIT["Meter",1]]

OGC WTK:

PROJCS["CH1903 / LV03",
    GEOGCS["CH1903",
        DATUM["CH1903",
            SPHEROID["Bessel 1841",6377397.155,299.1528128,
                AUTHORITY["EPSG","7004"]],
            TOWGS84[674.374,15.056,405.346,0,0,0,0],
            AUTHORITY["EPSG","6149"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4149"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    PROJECTION["Hotine_Oblique_Mercator"],
    PARAMETER["latitude_of_center",46.95240555555556],
    PARAMETER["longitude_of_center",7.439583333333333],
    PARAMETER["azimuth",90],
    PARAMETER["rectified_grid_angle",90],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",200000],
    AUTHORITY["EPSG","21781"],
    AXIS["Y",EAST],
    AXIS["X",NORTH]]


Another option is to modify the .prj file with the correct one.


Adrian Waddell wrote:
> Hello There,
>
> I'm trying to import the sample shapefile from
>
> http://www.swisstopo.admin.ch/internet/swisstopo/en/home/products/downloads/landscape/vector25.html 
>
>
> (ArcView Shapefile) into GRASS 6.3.
>
> I guess the Coordinate Reference System for the data is CH1903+ /LV95, 
> EPSG code 2056. However if I create a new location with these CRS 
> settings, I get the following error while trying to import the data:
>
> > v.in.ogr geb_25_l.shp out=geb25
> ERROR: Projection of dataset does not appear to match current location.
>
> I get the same error message if I use the EPSG 4149,4150, 21781 settings.
>
> However if I import the data by creating a new location, i.e.
>
> > v.in.ogr geb_25_l.shp out=geb25 location=CHTopo
>
> the data import works fine, however I get the following error if I try 
> to change the projection settings:
>
> XY-location cannot be projected
>
> Does anybody has an idea how to fix this problem?
>
> Sincerely,
>
> Adrian Waddell
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>


-- 
-----------------------------------------------------
Dr. Massimiliano Cannata
Environmental & Geomatic Engineer

Institute of Earth Sciences - SUPSI
Trevano, C.P. 72, CH-6952 Canobbio, SWITZERLAND

Tel:    +41 (0)58 / 666 62 14  
Fax:    +41 (0)58 / 666 62 09
E-mail: massimiliano.cannata at supsi.ch

Web: 
http://www.ist.supsi.ch
http://istgis.ist.supsi.ch:8001/geomatica/
-------------------------------------------------------