[Mapserver-DE] Projektion

[Dirk Jesko]

Hallo Andreas,

das sieht doch gut aus. Also, das ist eine Lambertprojektion mit zwei 
Standardparallelen auf Bessel-Ellipsoid. Ob es dafür einen epsg-Code 
gibt, kann ich nicht sagen. Daher ist es wahrscheinlich am besten, die 
Parameter explizit anzugeben. Lambert hatte ich im Mapserver noch nicht. 
Daher kenn ich die Parameter nicht genau. Ich empfehle mal in Google 
nach "proj.4" und/oder "mapserver" und "proj=lcc" (Lambert Conformal 
Conic) zu suchen. Dann finden sich bestimmt Beispiele. Das müsste 
irgendwie so aussehen:

    PROJECTION
        "proj=lcc"
        "lon_0=..."
        "x_0=..."
        "lat_0=..."
        "lat_1=..."
        "lat_2=..."
        "y_0=..."
        "ellps=bessel"
        "datum=..."
        "units=m"
        "no_defs"
    END

Evtl. muss noch ein "towgs84"-Parameter angegeben werden, falls die 
Datenquellen mit unterschiedlichen Ellipsoiden arbeiten. Hier gibt es 
aber unterschiedliche Parametersätze je nach Gebiet. Da hilft ggf. nur 
suchen und ausprobieren.

Den PROJCS-String zu interpretieren ist einfacher, als es auf den ersten 
Blick aussieht (die Kenntnis der Projektion vorausgesetzt). Hier mal ein 
paar Tipps:

>PROJCS["Custom",
>
Name der Projektion. In diesem Fall anscheinend etwas selbst definiertes.

>GEOGCS["GCS_Bessel_1841",DATUM["D_Bessel_1841",
>
Ellipsoid und Datum. In diesem Fall Bessel 1841.

>SPHEROID["Bessel_1841",6377397.155,299.1528128]],
>
Das sind nochmal die Parameter des Ellipsoiden. Die sind unerheblich, da 
proj.4 Bessel kennt.

>PRIMEM["Greenwich",0],
>
Der Nullmeridian. Sollte auch unerheblich sein, kann aber bei Bedarf

>UNIT["Degree",0.0174532925199433]],
>
Hier bin ich mir nicht ganz sicher.

>PROJECTION["Lambert_Conformal_Conic"],
>
Name der Projektion, "Lambert Conformal Conic" -> "lcc"

>PARAMETER["False_Easting",400000],
>
Falscher Rechtswert, vermutlich x_0

>PARAMETER["False_Northing",400000],
>
Falscher Nordwert, vermutlich y_0

>PARAMETER["Central_Meridian",13.3333333333333],
>
Bezugsmeridian, vermutlich lon_0

>PARAMETER["Standard_Parallel_1",46],
>
1. Standardparallele, vermutlich lat_1

>PARAMETER["Standard_Parallel_2",49],
>
2. Standardparallele, vermutlich lat_2

>PARAMETER["Central_Parallel",48],
>
Zentraler Breitenkreis, vermutlich lat_0

>UNIT["Meter",1]]
>  
>
Einheit für die Koordinaten. Das sollte der "units"-Eintrag sein.

Ich hoffe das hilft etwas weiter. Evtl. hilft auch eine Suche bei ESRI. 
Wenn ich mich nicht täusche, haben die dieses Format spezifiziert.


Dirk


>lg andreas
>
>----- Original von:  Dirk Jesko <ml.dje at geocontent.de>:
>
>  
>
>>Hallo,
>>    
>>
>
>ich müsste mal eine der Projektionsdateien sehen.
>
>Gruß,
>Dirk
>
>
>aon.912021231.a at aon.at wrote:
>
>  
>
>>Guten morgen,
>>
>>danke für deine mail... also die shape files besitzen eine projektionsdatei. aus dieser ist es mir aber unmöglich den epsg code oder andere angaben für das map file abzuleiten!
>>vielleicht hast du da noch einen tipp. die daten sind aus österreich.
>>Danke für deine hilfe!
>>Lg Andreas
>>
>>    
>>
>