[gdal-dev] gdal_translate via wmts.xml produces an empty tif

Helmut Kudrnovsky hellik at web.de
Wed Jun 14 12:16:06 PDT 2017


hi,

following the example in [1] with the WMTS GetCapabilities of [2], I do:

gdal_translate
"WMTS:https://www.basemap.at/wmts/1.0.0/WMTSCapabilities.xml,layer=bmaporthofoto30cm"
wmts_bm.xml -of WMTS

gives this file:  wmts_bm.xml
<http://osgeo-org.1560.x6.nabble.com/file/n5324225/wmts_bm.xml>  

then gdalinfo of wmts_bm.xml shows:

gdalinfo wmts_bm.xml
Driver: WMTS/OGC Web Mab Tile Service
Files: wmts_bm.xml
Size is 6268828, 2969099
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
+no_defs"],
    AUTHORITY["EPSG","3857"]]
Origin = (977649.895640023052692,6281289.982633858919144)
Pixel Size = (0.149291070871600,-0.149291070871600)
Metadata:
  ABSTRACT=Basemap als farbiges Orthofoto
  TITLE=Geoland Basemap Orthofoto
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  977649.896, 6281289.983) (  8d46'56.56"E, 49d 2'16.34"N)
Lower Left  (  977649.896, 5838030.013) (  8d46'56.56"E, 46d21'31.57"N)
Upper Right ( 1913529.941, 6281289.983) ( 17d11'22.31"E, 49d 2'16.34"N)
Lower Right ( 1913529.941, 5838030.013) ( 17d11'22.31"E, 46d21'31.57"N)
Center      ( 1445589.918, 6059659.998) ( 12d59' 9.44"E, 47d42'55.93"N)
Band 1 Block=128x128 Type=Byte, ColorInterp=Red
  Overviews: 3134414x1484550, 1567207x742275, 783604x371137, 391802x185569,
195901x92784, 97950x46392, 48975x23196, 24488x11598, 12244x5799, 6122x2900,
3061x1450, 1530x725, 765x362, 383x181
  Mask Flags: PER_DATASET ALPHA
  Overviews of mask band: 3134414x1484550, 1567207x742275, 783604x371137,
391802x185569, 195901x92784, 97950x46392, 48975x23196, 24488x11598,
12244x5799, 6122x2900, 3061x1450, 1530x725, 765x362, 383x181
Band 2 Block=128x128 Type=Byte, ColorInterp=Green
  Overviews: 3134414x1484550, 1567207x742275, 783604x371137, 391802x185569,
195901x92784, 97950x46392, 48975x23196, 24488x11598, 12244x5799, 6122x2900,
3061x1450, 1530x725, 765x362, 383x181
  Mask Flags: PER_DATASET ALPHA
  Overviews of mask band: 3134414x1484550, 1567207x742275, 783604x371137,
391802x185569, 195901x92784, 97950x46392, 48975x23196, 24488x11598,
12244x5799, 6122x2900, 3061x1450, 1530x725, 765x362, 383x181
Band 3 Block=128x128 Type=Byte, ColorInterp=Blue
  Overviews: 3134414x1484550, 1567207x742275, 783604x371137, 391802x185569,
195901x92784, 97950x46392, 48975x23196, 24488x11598, 12244x5799, 6122x2900,
3061x1450, 1530x725, 765x362, 383x181
  Mask Flags: PER_DATASET ALPHA
  Overviews of mask band: 3134414x1484550, 1567207x742275, 783604x371137,
391802x185569, 195901x92784, 97950x46392, 48975x23196, 24488x11598,
12244x5799, 6122x2900, 3061x1450, 1530x725, 765x362, 383x181
Band 4 Block=128x128 Type=Byte, ColorInterp=Alpha
  Overviews: 3134414x1484550, 1567207x742275, 783604x371137, 391802x185569,
195901x92784, 97950x46392, 48975x23196, 24488x11598, 12244x5799, 6122x2900,
3061x1450, 1530x725, 765x362, 383x181

when I try to get some smaller area of this WMTS by gdal_translating this 
wmts_bm2.xml
<http://osgeo-org.1560.x6.nabble.com/file/n5324225/wmts_bm2.xml>  :

gdal_translate wmts_bm2.xml test.tif
Input file size is 3363, 3363
0...10...20...30...40...50...60...70...80...90...100 - done.

the result is an empty tif:

D:\basemap_orthos\2017\inn>gdalinfo test.tif -mm
Driver: GTiff/GeoTIFF
Files: test.tif
Size is 3363, 3363
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Mercator_1SP"],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
+no_defs"],
    AUTHORITY["EPSG","3857"]]
Origin = (1168011.985145896673203,5932420.039413984864950)
Pixel Size = (0.149291070871600,-0.149291070871600)
Metadata:
  ABSTRACT=Basemap als farbiges Orthofoto
  AREA_OR_POINT=Area
  TITLE=Geoland Basemap Orthofoto
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 1168011.985, 5932420.039) ( 10d29'32.75"E, 46d56'26.96"N)
Lower Left  ( 1168011.985, 5931917.974) ( 10d29'32.75"E, 46d56'15.87"N)
Upper Right ( 1168514.051, 5932420.039) ( 10d29'48.99"E, 46d56'26.96"N)
Lower Right ( 1168514.051, 5931917.974) ( 10d29'48.99"E, 46d56'15.87"N)
Center      ( 1168263.018, 5932169.006) ( 10d29'40.87"E, 46d56'21.42"N)
Band 1 Block=3363x1 Type=Byte, ColorInterp=Red
    Computed Min/Max=0.000,0.000
  Mask Flags: PER_DATASET ALPHA
Band 2 Block=3363x1 Type=Byte, ColorInterp=Green
    Computed Min/Max=0.000,0.000
  Mask Flags: PER_DATASET ALPHA
Band 3 Block=3363x1 Type=Byte, ColorInterp=Blue
    Computed Min/Max=0.000,0.000
  Mask Flags: PER_DATASET ALPHA
Band 4 Block=3363x1 Type=Byte, ColorInterp=Alpha
    Computed Min/Max=0.000,0.000

displaying this WMTS in QGIS with
https://www.basemap.at/wmts/1.0.0/WMTSCapabilities.xml works. the layers are
all displayed.

any hint to solve this?

[1] http://www.gdal.org/frmt_wmts.html
[2] https://www.basemap.at/






-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/gdal-translate-via-wmts-xml-produces-an-empty-tif-tp5324225.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.


More information about the gdal-dev mailing list