Luftbilder von Aerowest

Communities Deutschland (Germany)
94

Hier die versprochene weitere Möglichkeit viele verschiedene *.tif mittels Shapefile-Index zu benutzen. Man sieht hier bspw auch, dass die Daten nicht in WGS84 vorliegen (PROJECTION “init=epsg:31277”) , aber vom Mapserver als WGS84 (“wms_srs” “EPSG:4326”) ausgeliefert werden.
Das ist vielleicht einfacher als die Dateien zu verschmelzen. Wie es bei Überlappungen aussieht, weiß ich nicht. Müsste man probieren

gdaltindex senta.shp /osm/fileserver/SentaDOP/*.tif

Das Layer im Mapfile sieht dann so aus


  LAYER
    NAME "SentaDOP"
    STATUS ON
    PROJECTION
          "init=epsg:31277"
    END
    METADATA
     "wms_title"  "SentaDOP"
     "wms_srs"    "EPSG:4326"
     "copyright"    ""
    END


    TYPE RASTER
    TILEINDEX "/osm/fileserver/SentaDOP/senta.shp"
    TILEITEM "LOCATION"
  END

Heute morgen habe ich eine sauber in WGS84 georefernzierte Katasterkarte für jemanden bereitgestellt. Ich habe mal mitgeschrieben. An den ESPG Codes kann ja jeder selber rumändern



thomas@suncobalt:/osm/mapserver/bajmok$ gdalinfo Bajmok_wgs84.tif | more
Driver: GTiff/GeoTIFF
Files: Bajmok_wgs84.tif
       Bajmok_wgs84.tif.ovr
       Bajmok_wgs84.tif.aux.xml
Size is 70810, 52531
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.2572235630016,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
..............



thomas@suncobalt:/osm/mapserver/bajmok$ python /osm/mapserver/print_wms_extent.py Bajmok_wgs84.tif
"wms_extent" "19.396992 45.956686 19.446389 45.993332"
thomas@suncobalt:/osm/mapserver/bajmok$ gdal_translate -of GTiff -co "TILED=YES" Bajmok_wgs84.tif Bajmok_wgs84_tiled.tif
Input file size is 70810, 52531
0...10...20...30...40...50...60...70...80...90...100 - done.
thomas@suncobalt:/osm/mapserver/bajmok$ gdaladdo -ro --config USE_RRD YES -r gauss Bajmok_wgs84_tiled.tif 2 4 8 16 32 64 128 256
0...10...20...30...40...50...60...70...80...90...100 - done.

Abhängig von der Größe kann das u.U. Stunden bis Tage dauern.

Layer in Mapfile definieren


 LAYER
    NAME "Bajmok"
    STATUS ON
    TYPE RASTER

    DATA "/osm/mapserver/bajmok/Bajmok_wgs84_tiled.tif"

    PROJECTION
          "init=epsg:4326"
    END

    METADATA
     "wms_title"  "Bajmok"
     "wms_extent" "19.396992 45.956686 19.446389 45.993332"
     "wms_srs"    "EPSG:4326"
     "copyright"  "OSM use only"
    END
  END

Und dann in JOSM einbinden