Issues Downloading WMTS Data as Raster in ArcGIS Pro

PerAndersen — Fri, 25 Oct 2024 12:08:04 GMT

Hello GIS community,

I’m trying to create a Python script to download a specific area of imagery from a WMTS service and save it as a raster image (e.g., PNG or TIFF). This script is intended to be used within an ArcGIS Pro toolbox, where the user defines the area they want to download. The downloaded image will be used for Deep Learning in ArcGIS Pro.

Here’s what I’ve tried so far:

I have used GDAL in my script to connect to the WMTS service and attempted to download the data as a raster.
I’m working with the following WMTS URL: https://api.dataforsyningen.dk/orto_foraar_wmts_DAF?service=WMTS&request=GetCapabilities&token=577b8b556967c27dd029d1acbd544650
The process seems to complete without error, but the output image file is empty (0 KB).

I’m not using a layer_name parameter because the WMTS service should provide that based on the GetCapabilities response. The user also defines the extent within ArcGIS Pro.

Below is the script I’ve been working on:

import os
import logging
from osgeo import gdal

# Set up logging
logging.basicConfig(filename=os.path.join(output_folder, 'wmts_download.log'), level=logging.INFO, format='%(asctime)s - %(levelname)s - %(message)s')
logging.info('Log file created.')

# Parameters
wmts_url = 'https://api.dataforsyningen.dk/orto_foraar_wmts_DAF?service=WMTS&request=GetCapabilities&token=577b8b556967c27dd029d1acbd544650'
output_folder = r'D:\Ortofoto'
workspace = '713203,415311581 6173517,82703254 713491,018988423 6173662,51763386'
output_file = os.path.join(output_folder, 'geodanmark_2023_10cm.png')

# Log parameters
logging.info(f'URL: {wmts_url}')
logging.info(f'Output Folder: {output_folder}')

# Open the WMTS dataset
logging.info('Attempting to open WMTS dataset.')
ds = gdal.Open(wmts_url)
if ds is None:
logging.error('Failed to open WMTS dataset.')
raise RuntimeError('Failed to open WMTS dataset.')

# Set spatial reference
logging.info('Setting spatial reference.')
spatial_ref = gdal.osr.SpatialReference()
spatial_ref.ImportFromWkt(workspace)

# Attempt to create output file
logging.info(f'Attempting to create output file: {output_file}')
driver = gdal.GetDriverByName('PNG')
out_ds = driver.CreateCopy(output_file, ds)
if out_ds is None:
logging.error('Failed to create output file.')
raise RuntimeError('Failed to create output file.')

logging.info(f'Downloaded image saved as: {output_file}')

Has anyone faced a similar issue when using a WMTS service to download data as a raster image in ArcGIS Pro?

I would appreciate any advice or examples of scripts that successfully download and save WMTS data, or any alternative approaches for accomplishing this task.

Thanks in advance for your help!

//Per

topic Issues Downloading WMTS Data as Raster in ArcGIS Pro in Python Questions

Issues Downloading WMTS Data as Raster in ArcGIS Pro