QMapShack: WMS durch die Hintertür

Hi,

da QMS leider (noch) kein WMS unterstützt, habe ich nach einem Workaround gesucht. Und gefunden, bzw. gebastelt mit JavaScript: Dieser funktioniert für alle WMS-Karten, die die Web-Mercator-Projektion unterstützen (EPSG-Codes 3857 oder 900913).

Hier als Beispiel die Topokarte von NRW:

<TMS>
<Layer idx="0">
<Title>NRW DTK Sammeldienst</Title>
<Script><![CDATA[(
function convert(z1,x1,y1)
{ function M(n){return 256*n*(156543.03392804062/(1<<z1))-20037508.342789244}
  y=(1<<z1)-1-y1;
  bbox=M(x1)+','+M(y)+','+M(x1+1)+','+M(y+1);
  return 'http://www.wms.nrw.de/geobasis/wms_nw_dtk?LAYERS=nw_dtk_col&FORMAT=image/png&SRS=EPSG:3857&EXCEPTIONS=application/vnd.ogc.se_inimage&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&BBOX=' + bbox + '&WIDTH=256&HEIGHT=256';
}
)]]></Script>
</Layer>
<Copyright>Geobasis NRW 2015</Copyright>
</TMS>

Den Code einfach z.B. als nrw.tms im QMS-Kartenverzeichnis speichern und loslegen.

Viel Spaß
Helmut

Zitat von kiozen

Manche Server benutzen noch die alte EPSG Nummer für die TMS Projektion. Die wäre: EPSG:900913

Hallo Oliver,

ja, ausgerechnet bei dieser verbreiteten Projektion herrscht so viel Chaos:
http://gis.stackexchange.com/q…-mercator-projection-code

Übrigens kann man wohl auch die meisten ArcGIS-Kartenquellen als TMS abbilden. Der bbox-Parameter ist identisch mit dem beim WMS.

Beispiel Brasilien 1:1Mio:

<TMS>
<Layer idx="0">
<Title>Brasilien 1:1.000.000</Title>
<Script><![CDATA[(
function convert(z1,x1,y1)
{ function M(n){return 256*n*(156543.03392804062/(1<<z1))-20037508.342789244}
  y=(1<<z1)-1-y1;
  bbox=M(x1)+','+M(y)+','+M(x1+1)+','+M(y+1);
  return 'http://sigel.aneel.gov.br/arcgis/rest/services/CartasTopograficas/MapServer/export?dpi=96&transparent=true&format=png8&bbox=' + bbox + ' &bboxSR=3857&imageSR=3857&size=256,256&layers=show:17&f=image';
}
)]]></Script>
</Layer>
<Copyright>IBGE</Copyright>
</TMS>

Info-Quelle dazu wäre:
http://sigel.aneel.gov.br/arcg…tasTopograficas/MapServer

ArcGIS von ESRI ist eins der führenden GIS-Systeme. Mit geschickter Suche findet man so noch sehr viele Online-Karten!

Natürlich kannst du das alles ins Wiki übernehmen. Evtl. kann ich dabei noch helfen.

Grüße
Helmut

Zitat von kiozen

Wie sieht das eigentlich aus: Ich habe Rückmeldungen bekommen, dass die Windows Version bei TMS und WMTS Karten nicht immer stabil läuft. Kommt das Problem öfters vor oder ist das nur ein verkorkster Computer?

Ich werde das mal genauer beobachten. Ich glaube, ich habe da auch schon Abstürze erlebt. Ich melde mich.

Zitat von limex

Kartenfreak: Ich bekomme die Tiles, aber alle haben ein Schachbrettmuster.
Hast du eine Idee woran es liegt?
Und: ja, der Kartenausschnitt war in Brasilien!

Auf der Website gibt's das gleiche Phänomen:
http://sigel.aneel.gov.br/sigel.html

Keine Ahnung, warum. Wollen die uns ärgern?

Edit: Habe stattdessen ein Beispiel aus Italien ins Wiki gestellt.

Zitat von limex

Ich hab auch brav gespripted und was für Neuseeland bebastelt.
Hintergrund: Die Kiwis sind besonders schlau, die negativen Lat Koordinaten behandeln sie wie positive
Nicht mit mir
Ich habs auch ins Wiki gestellt. Hoffe es passt. War mein erster Wiki Eintrag überhaupt.

Code

<TMS>
    <Title>NZTopo</Title>
    <MinZoomLevel>6</MinZoomLevel>
    <MaxZoomLevel>15</MaxZoomLevel>
    <Layer idx="0">
        <Title>NZTopo</Title>
        <Script><![CDATA[(
        function convert(z1,x1,y1)
        {    
            function lat2tile(lat,zoom) { return (Math.floor((1-Math.log(Math.tan(lat*Math.PI/180)+1/Math.cos(lat*Math.PI/180))/Math.PI)/2*Math.pow(2,zoom))); }
            function tile2lat(y,z) { var n=Math.PI-2*Math.PI*y/Math.pow(2,z); return (180/Math.PI*Math.atan(0.5*(Math.exp(n)-Math.exp(-n)))); }
            lat = tile2lat(y1,z1); neg_lat = -lat; y2 = lat2tile(neglat,z1);
            return 'http://nz1.nztopomaps.com/' + z1 + '/' + x1 + '/' + y2 + '.png'; 
        }
        )]]></Script>
    </Layer>
</TMS>

Alles anzeigen

Hallo limex,

das kannst du aber auch einfacher haben:

function convert(z1,x1,y1)
 {  return "http://nz1.nztopomaps.com/" + z1 + "/" + x1 + "/" + ((1<<z1) - y1 - 1) + ".png";
 }

Wobei 1<<z1 das Gleiche ist wie Math.pow(2,z1).

Was die Neuseeländer benutzen, entspricht der Indexierung vom OSGEO TMS:
http://www.maptiler.org/google…s-tile-bounds-projection/

Man muss also genau aufpassen, was jeweils mit TMS gemeint ist.