Topo50 Karten von Österreich

peter_rv · 30 Juli 2009

Für den Deutschlandteil der Karten sind wir ja eigentlich ganz gut versorgt.
Bei Österreich und der Schweiz schauts dagegen nicht so gut aus.
Aber es gibt doch geoland.at für Österreich bzw die österreichischen Bundesländer mit den topografischen Karten:
http://www.geoland.at/geolandWMS/service.aspx?Name=geoland_at_wms&Service=WMS&Version=1.1.1&Request=GetCapabilities
Ich hab da einigermassen erfolglos versucht, eine at_topo_50.wms zu zaubern, allerdings bin ich nicht weit gekommen.

Hat jemand einen Tip?

pfeffer · 30 Juli 2009

Aus der Capabilities-Antwort geht hervor, dass der WMS nur EPSG-Code 31287 unterstützt.
Man kann hier http://www.epsg-registry.org/ den EPSG-Code 31287 nachgucken und feststellen, dass dies eine Lambert-Projektion ist. Lambert-Projektion kann Cachewolf leider derzeit nicht. Cachewolf kann leider nur Gauß-Krüger-Projektionen (mit unterschiedlichen Parameter, d.h. auch UTM) und direkt WGS84.

Wer will, kann gerne die Lambert-Projektion programmieren ;-)
Ich helfe, wo ich kann.

Gruß,
Pfeffer.

greiol · 30 Juli 2009

könnte http://www.i3s.unice.fr/~johan/gps/ helfen? ist gpl und java.
http://www.jhlabs.com/java/maps/proj/ käme mit apache license

pfeffer · 31 Juli 2009

sehr hilfreich ist immer: http://www.crs-geo.eu/nn_124226/crseu/EN/CRS__Description/crs-national__node.html?__nnn=true

Die nötige Helmert-Transformation kann Cachewolf schon, es fehlt nur das Projizieren / Zurückprojezieren zu / von Lambert.

Gruß,
Pfeffer.

EDIT:
PS: Hatte bei meinem ersten Posting übersehen, dass auch EPSG:4326 (=WGS84) irgendwo angeboten wird. Ob das funktioniert kann man gut mit GoogleEarth testen: Hinzufügen/Bild Overlay/Aktualisieren/WMS Parameter/Hinzufügen... und dann "http://www.geoland.at/geolandWMS/service.aspx?Name=geoland_at_wms" eingeben. GE kann nämlich nur WGS84. Ich hab's schon probiert: Es scheint ein ähnliches Problem zu geben wie bei einigen anderen WMS auch: in WGS84 stimmten die Koordinaten nicht. Es müsste also tatsächlich die Lambertprojeotion programmiert werden. Die Franzosen werden es uns danken, denn die haben auch Lambert.

pfeffer · 31 Juli 2009

in dem von Greiol angegebenen Quellcode findet sich:

Code:

/**
   * initalizes from geographic coordinates
   *
   * @param coord geographic coordinates triplet
   * @param Xs false east (coordinate system origin) in meters
   * @param Ys false north (coordinate system origin) in meters
   * @param c projection constant
   * @param n projection exponent
   * @param e reference ellipsoid excentricity
   * @param lg0 longitude of origin wrt to the Greenwich meridian (in radian)
   */
  public ConicProjection(Geographic coord, double Xs, double Ys,
			 double c, double n, double e, double lg0) {
    double eslt = e * Math.sin(coord.lt());
    double l = Math.log(Math.tan(Math.PI/4.0 + coord.lt()/2.0) *
			Math.pow((1.0 - eslt)/(1.0 + eslt), e/2.0));
    east = Xs + c * Math.exp(-n * l) * Math.sin(n * (coord.lg() - lg0));
    north = Ys - c * Math.exp(-n * l) * Math.cos(n * (coord.lg() - lg0));
  }
}

Wenn das alles sein sollte, dann kann es nicht so schwierig sein.
Hilfreich könnte auch noch dieser Thread sein: http://www.geoclub.de/viewtopic.php?f=54&t=23912

Gruß,
Pfeffer.

pfeffer · 31 Juli 2009

Könnte mal jemand überprüfen, ob die im zitierten Code-Block verwendeten Formeln zu dem gleichen Ergebnis fürhren, wie die, die im Thread zitiert sind?

Danke,
Pfeffer.

arbor95 · 1 Aug. 2009

Für den Server gibt es doch eine fertige wms-Datei.

Habt ihr die schon mal getestet?

pfeffer · 1 Aug. 2009

erst ärgere ich mich, dass der Erstposter es nicht gesehen hat.
Dann sehe ich zwar die .wms, aber nicht in der Liste in Cachewolf. Grund: Die Sortierung der WMS ändert sich nicht, wenn man das Zentrum um das Karten geladen werden sollen, neu festlegt, in dem man beim Runterladen-Dialog es ändert. --> Könnte in CW verbessert werden.
Dann lese ich den Kommentar, den wahrscheinlich Du beim erstellen der .wms in die .wms geschrieben hast:

Code:

# nicht überall Luftbilder
# nicht überall gut kalibriert
# bis Skalierung 4,3 gibts noch Fotos, ab 4,5 nur noch Strassen (schlecht kalibriert)

Die Frage ist
a) Sind die Karten selbst schlecht oder liegt es daran, dass die die Umrechnung auf WGS84 nicht richtig hinbekommen (wie beispielsweise die Server von NRW)?

b) Sollte man vielleicht andere Layer nehmen (Topo)?

zu a):
Ich weiß nicht, wie man a) herausbekommen soll, außer durch Lambertprojektion programmieren. - hmm - oder doch: vielleicht weiß ja ein Österreicher Lambert-Koordinaten und kann die in der Online-Karte überprüfen?
Was ich schon probiert habe: man kann bei geoland.at andere WMS mit einbinden. Ich habe mal die EU-OSM-WMS eingebunden. Die liegen recht deutlich daneben. Das kann aber wiederum daran liegen, dass der geoland-Server die OSM-Karten in WGS84 abrufen muss und dass er diese Umrechnung nicht besopnders präzise macht, war ja schon meine Vermutung. Dieser Test bringt also nicht viel.

zu b):
Das möge mal bitte jemand ausprobieren (kann jeder auf geoland.at) und uns sagen, welches die Layer sind, die wir in Cachewolf einstellen sollen.

Gruß,
Pfeffer.

pfeffer · 2 Aug. 2009

ha puhh.
das ist ja alles eine Bastellei, furchtbar.
In Österreich verwenden sie 2 Standardparallelen, in Frankreich nur 1. - Wie unterscheiden sie die Projektionsformeln? - Ich habe hier ( http://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=Lambert_Conformal_Conic ) die Information gefunden, dass
a) 2 standard-paralle das gleich wie eine standard-parallele mit Skalierungsfaktor < 1 ergibt
b) man für 1 Standardparallele die gleiche Formel verwenden kann wir für 2, wenn man einfach die 2. Standardparallele auf die erste setzt.
Das hört sich doch schon mal ganz gut an.

Frankreich ist in mehrere Zonen aufgeteilt, die Ostwerte haben eine _meist_ 200.000 addiert (false Easting), aber man kann nicht (ohne weiteres) erkennen, zu welcher Zone eine Koordinate gehört (wenn ich das richtig recherchiert habe). (siehe http://www.crs-geo.eu/nn_124226/crseu/EN/CRS__Description/crs-national__node.html?__nnn=true )

Ich weiß echt nicht, wie ich da nen vernünftiges Konzept draus basteln soll.

Kennt sich jemand mit Lambert-Projektionen aus?

Gruß,
Pfeffer.

pfeffer · 5 Aug. 2009

ich wollte mal sagen:
Ich bastele gerade intensiv daran, dass Cachewolf den entsprechenden EPSG unterstützt. Jetzt, wo der Code dafür schon fast fertig ist, muss ich leider feststellen, dass die von geoland angebotenen Karten nicht so dolle sind. Die Abdeckung scheint sehr schlecht zu sein. Ich habe noch nicht herausbekommen, welchen Bereich in welcher Auflösung in welchem Layer der Server bietet.
Das das mal jemand für mich herausfinden?

Oder alternativ: Kann mir jemand vielleicht einen WMS (zum Testen) empfehlen, der gute Karten liefert und auf eine Lambert-Projektion angewiesen ist?

Gruß,
Pfeffer.

arbor95 · 5 Aug. 2009

Hi,
die beiden wms-Dateien von Niederösterreich liefern
<SRS>EPSG:31287</SRS>
<SRS>EPSG:31466</SRS>
<SRS>EPSG:31467</SRS>
<SRS>EPSG:4326</SRS>
Bei den layern steht zwar nur EPSG:4326.

Dann gibt es noch 2 wms
http://www.centropemap.org/cgi-bin/mapserv?map=/data/wms/centrope/ce_topogr_background_2008.map&
und
http://www.centropemap.org/cgi-bin/mapserv?map=/data/wms/centrope/ce_sat_image.map&

deren wms-Dateien hab ich mangels Zeit noch nicht zur Veröffentlichung fertig. Für Ö liefern Sie Daten für Niederösterreich, Wien und das Burgenland. Wie das bei den anderen abgedeckten Ländern ist konnt ich noch nicht prüfen.

pfeffer · 5 Aug. 2009

Meine Güte, sind in Österreich viele Koordinaten-Referenzsysteme in Gebrauch!
WGS84 funktioniert bei denen ja offenbar allen nicht richtig.
Bei Niederösterreich vermute ich, dass das Kartenmaterial ursprünglich in EPSG:31259 (="MGI Austria GK M34" auch "BMN - M34") vorliegt und vermutlich dann richtig kalibriert geliefert wird, wenn in diesem System die Karten abgerufen werden.
Das ist aber eine Gauß-Krüger-Projektion, die lässt sich ohnehin bereits jetzt leicht in Cachewolf integrieren.

Warum wird in der .wms 4326 (WGS84) verwendet und nicht 31466/7? - Das ist das deutsche Gauß-Krüger, vielleicht funktioniert ja bereits das besser als WGS84.

Gruß,
Pfeffer.

arbor95 · 5 Aug. 2009

pfeffer schrieb:
Warum wird in der .wms 4326 (WGS84) verwendet und nicht 31466/7? - Das ist das deutsche Gauß-Krüger, vielleicht funktioniert ja bereits das besser als WGS84.

Gruß,
Pfeffer.

Weil ich mich mit Ö nicht wirklich beschäftigt habe! Das war damals einfach schnell was übernommen und hochgeladen. (sorry)

pfeffer · 7 Aug. 2009

so.
Habe die österreichische Lambertprojektion praktisch fertig.
Allerdings scheinen mir die meisten WMS, die damit zusätzlich genutzt werden können auf den ersten Blick recht untauglich.
http://www.centropemap.org liefert keine Hintergrundkarte, nur Grenzen und Naturschutzgebiete o.ä., jedenfalls laut http://www.centropemap.org/cgi-bin/mapserv?map=/data/wms/centrope/ce_administrative_units_bordersonly.map&VERSION=1.1.1&REQUEST=GetCapabilities&SERVICE=WMS

Evtl. sind bei centropemap.org noch weitere WMS-Server verwendet, von denen wir aber (bislang) nicht die URL kennen.

Ich würde es eigentlich gerne an einem ordentlichen Server testen, bevor ich es einchecke...
Kennt jemand noch einen?

Vielleicht taugt ja auch der geoland.at was - nur ich bekomme keine hochaufgelöste Karte / Orthophotos von dem. Könnte vielleicht mal jemand mit deejump / openjump versuchen, entsprechende Layer und Region ausfindig zu machen?

Gruß,
Pfeffer.

arbor95 · 7 Aug. 2009

check mal hiermit:

pfeffer · 7 Aug. 2009

ahh - super, danke!
Da gibt es ja doch vernünftige Karten von Centrop. Damit kann ich jetzt ordnetlich testen und debuggen.

Gruß,
Pfeffer.

pfeffer · 8 Aug. 2009

sooo - ich glaub, Cachewolf funktioniert jetzt richtig, aber centrop liefer in lambert noch schlecher kalibirierte Karten (50km falsch). Aber vielleicht bin ich auch zu müde.
Könnte noch mal jemand folgendes testen:

hiermit eine Karte von centromap abrufen:
http://www.centropemap.org/cgi-bin/mapserv?map=/data/wms/centrope/ce_topogr_background_2008.map&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&SRS=EPSG:31287&BBOX=647425.84,369772.89,657425.84,379772.89&WIDTH=1000&HEIGHT=1000&LAYERS=m50k_01,m50k_02,m50k_03,m50k_04,m50k_05,m50k_06,m50k_07,m50k_08,m50k_09,m50k_10,m50k_11,m50k_12&STYLES=&FORMAT=image/png

dann die Lambert-Koordinaten (x: 647425.84,y:369772.89 (links unten),x:657425.84,y:379772.89 = rechts oben)) in WGS84 umrechnen (lassen), z.B. mit Transdat oder mit http://www.geoland.at/geoland2/(slnvbpml3noq2cbb3niwi355)/coordtrans.aspx - falls dieser Link nicht funktioniert: www.geoland.at, viewer starten und dann oben auf den Button x/y.
als Quellsystem Lambert(alt) [= EPSG:31287] wählen und in WGS84 umrechnen.
Die so gewonnen Koordinaten in maps.google.de eingeben und prüfen, ob man den gleichen Kartenausschnitt gewinnt.
Bei mir (Cachewolf-Lambert) lag er ca. 50km südlich... wenn das auch mit anderen Umrechnungstools so passiert, dann macht der WMS-Server wohl Mist.

Gruß,
Pfeffer.

arbor95 · 8 Aug. 2009

bin jetzt erst mal bis ende august im urlaub , aber wenn ich zeit / internet finde werd ich das mal checken.

Axel2 · 9 Aug. 2009

pfeffer schrieb:
hiermit eine Karte von centromap abrufen:
http://www.centropemap.org/cgi-bin/mapserv?map=/data/wms/centrope/ce_topogr_background_2008.map&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&SRS=EPSG:31287&BBOX=647425.84,369772.89,657425.84,379772.89&WIDTH=1000&HEIGHT=1000&LAYERS=m50k_01,m50k_02,m50k_03,m50k_04,m50k_05,m50k_06,m50k_07,m50k_08,m50k_09,m50k_10,m50k_11,m50k_12&STYLES=&FORMAT=image/png

dann die Lambert-Koordinaten (x: 647425.84,y:369772.89 (links unten),x:657425.84,y:379772.89 = rechts oben)) in WGS84 umrechnen (lassen), z.B. mit Transdat oder mit http://www.geoland.at/geoland2/(slnvbpml3noq2cbb3niwi355)/coordtrans.aspx - falls dieser Link nicht funktioniert: http://www.geoland.at, viewer starten und dann oben auf den Button x/y.
als Quellsystem Lambert(alt) [= EPSG:31287] wählen und in WGS84 umrechnen.

Hallo Pfeffer,
habe es ausprobiert.

Es paßt gut, wenn ich bei der Koordinatentransformation mit geoland.at nur Lambert, _nicht_ Lambert-alt, angebe.

Viele Grüße,
Axel

pfeffer · 9 Aug. 2009

danke!
Dann werde ich damit mal versuchen, ob ich weiter komme, obwohl Lambert(neu) eigentlich falsch ist, denn EPSG:31287 bezieht sich eindeutig auf den Bessel-Ellipsoid.

Gruß,
Pfeffer.

Topo50 Karten von Österreich

Geocacher

Geowizard

Geoguru

Geowizard

Geowizard

Geowizard

Geoguru

Geowizard

Geowizard

Geowizard

Geoguru

Geowizard

Geoguru

Geowizard

Geoguru

Anhänge

Geowizard

Geowizard

Geoguru

Geocacher

Geowizard