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@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: [EMAIL PROTECTED]
Web:
http://www.ist.supsi.ch
http://istgis.ist.supsi.ch:8001/geomatica/
-------------------------------------------------------
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
