如何用Python从WCS的multipart响应中提取TIFF高程数据？

当WCS服务返回multipart/related类型响应时，需要先解析多部分内容，分离出XML元数据和TIFF二进制数据。以下是可行的实现方法：

关键步骤说明 #

1. 处理压缩内容：响应头包含Content-Encoding: gzip，必须先解压原始响应字节，否则无法正确解析多部分结构。
2. 解析多部分响应：使用Python标准库email的BytesParser处理二进制格式的multipart内容，避免编码问题。
3. 提取并保存TIFF：遍历所有响应部分，根据Content-Type筛选出image/tiff类型的内容，写入文件。

完整代码示例 #

import gzip
from email.parser import BytesParser
from email.policy import default
from owslib.wcs import WebCoverageService

wcs_url = "https://geo.api.vlaanderen.be/DHMV/wcs?"
bbox_vla = (102200, 192500, 102700, 193000)

# 连接WCS服务并获取响应
wcs = WebCoverageService(wcs_url)
response = wcs.getCoverage(
    identifier=['DHMVII_DTM_1m'],
    bbox=bbox_vla,
    format='image/tiff',
    crs='EPSG:31370',
    subsets=[('x', 102200, 102700), ('y', 192500, 193000)],
    scalefactor=2
)

# 读取原始响应内容并处理gzip压缩
raw_content = response.read()
if response.headers.get('Content-Encoding') == 'gzip':
    raw_content = gzip.decompress(raw_content)

# 解析multipart/related结构
msg = BytesParser(policy=default).parsebytes(raw_content)

# 遍历所有部分，提取TIFF文件
for part in msg.walk():
    if part.get_content_type() == 'image/tiff':
        tiff_data = part.get_payload(decode=True)
        with open('elevation_dtm.tif', 'wb') as tiff_file:
            tiff_file.write(tiff_data)
        print("高程TIFF文件已成功保存")
        break
else:
    print("响应中未找到image/tiff类型的内容")

常见问题排查 #

• 如果仍无法提取TIFF，检查响应头的boundary值是否与解析时一致：email库会自动从Content-Type头中提取边界值，无需手动指定。
• 确认getCoverage参数正确：format='image/tiff'必须明确指定，避免服务返回默认格式。

内容的提问来源于stack exchange，提问作者stray_dog