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. _______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
