[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