GDAL下无需下载完整文件提取在线netCDF时间序列方法

基础背景 #

• 本人为Python与GDAL新手，此前主要使用R语言开展相关工作，为完成本次需求才初次接触上述工具，如需补充信息提升问题复现度、优化表述清晰度可随时告知。
• 核心需求：无需下载完整文件，从多个在线netCDF文件中提取指定点位的时间序列数据。目前已知gdallocationinfo工具可对本地存储的.nc文件实现点位数值提取，但GDAL提供的在线文件访问虚拟文件系统/vsicurl/无法兼容.nc格式文件。
• 预期效果：实现和在线GeoTIFF文件取值一致的使用体验，对指定在线netCDF文件完成WGS84坐标系下指定经纬度点位的数值提取。
  • 测试用在线netCDF文件：CHELSAcmip5ts_pr_ACCESS1-3_rcp45_2030-2049_V1.1.nc
  • 参考的在线GeoTIFF有效取值命令如下：

gdallocationinfo -wgs84 -valonly "/vsicurl/https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V2/GLOBAL/monthly/cmi/CHELSA_cmi_01_1981_V.2.1.tif" 14.27585 48.9335

• 可接受方案范围：所有可落地的实现路径均可，包括但不限于通过配置让vsicurl支持.nc文件读取、其他替代工具/代码实现方案。

补充说明 #

最初诉求是修复如下命令行的运行问题，经确认vsicurl与.nc文件搭配存在固有兼容问题，因此调整问题核心为寻找可满足需求的通用可行方案：

gdallocationinfo -wgs84 -valonly "/vsicurl/" 14.27585 48.9335

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