Hi

I tried running the following code using the shapefiles for the city of Zurich

require("rgdal") # requires sp, will use proj.4 if installed
require("maptools")

zh = readOGR(".",layer="Zurich")
and got following error:
OGR data source with driver: ESRI Shapefile
Source: ".", layer: "Zurich"
with 34 features and 6 fields
Feature type: wkbPolygon with 2 dimensions
Fehler in CRS(p4s) : lat_0 = 0 or 90 or alpha = 90


ogrInfo gives me this:

ogrInfo(".","Zurich")

Source: ".", layer: "Zurich"

Driver: ESRI Shapefile number of rows 34

Feature type: wkbPolygon with 2 dimensions

Extent: (676227 241585) - (689671 254305.7)

CRS: +proj=omerc +lat_0=46.95240555555556 +lonc=7.439583333333333 +alpha=90 +k=1 +x_0=-9419820.590700001 +y_0=200000 +ellps=bessel +units=m +no_defs

LDID: 87

Number of fields: 6

        name type length typeName

1    GMDEQNR    0      7  Integer

2         NR    0      9  Integer

3       NAME    4     45   String

4       GMDE    0      9  Integer

5 FLAECHE_HA    0      9  Integer

6   Quartier    0      4  Integer


I can see that lat_0 is not equal to zero, but don't know how to resolve this issue.

Thanks in advance for help

Renger


_________________________________________
Renger van Nieuwkoop
Centre for Energy Policy and Economics, ETH
Zürichbergstrasse 18 (ZUE)
CH - 8032 Zürich
+41 44 632 02 63
mailto: rengerv@etzh.ch<mailto:rengerv@etzh.ch>
blog.modelworks.ch


	[[alternative HTML version deleted]]

