吉奥蒂夫
您的文件不是常规的tiff文件,而是需要特殊库的geoTiff文件。在
对于python,有一个georasters库来读取这些文件。然后可以用matplotlib显示它们。在
在我看来,使用requests比urllib有更好的接口:import requests
from PIL import Image
import georasters as gr
import matplotlib.pyplot as plt
url = 'http://wdc.dlr.de/wdcservices/wcs.php'
query = {
'COVERAGE': '17e72d93-76d9-4af6-9899-b7b04e2763c8',
'service': 'wcs',
'version': '1.0.0',
'crs': 'epsg:4326',
'bbox': '-25,30,45,70',
'RESX': '0.1',
'RESY': '0.1',
'request': 'getcoverage',
'format': 'application/x-tiff-32f',
'TIME': '2015-12-13T00',
'elevation': '0',
'OUTPUTFILENAME': '17e72d93-76d9-4af6-9899-b7b04e2763c8_2015-12-13T00_0'
}
with open('test.tiff', 'wb') as f:
ret = requests.get(url, stream=True, params=query)
for data in ret.iter_content(1024):
f.write(data)
data = gr.from_file('test.tiff')
plt.imshow(data.raster, cmap='gray')
plt.show()
结果: