将GTiff栅格的时间序列堆栈转换为单个NetCDF
从gdal-dev邮件列表移出: 在2013年9月2日星期一晚上7:09,David Shean写道: 嗨,大家好,我正在尝试将具有相同投影/范围/分辨率的GTiff栅格时间序列打包为单个NetCDF文件进行分发。我花了一个小时来咨询在线文档,并且使用gdal_translate,gdalbuildvrt和gdalwarp都没有成功。 有使用现有的gdal命令行实用程序执行此操作的简单方法吗?在使用NetCDF Python API求助于自定义解决方案之前,我想过要问一下。 谢谢。-大卫 在2013年9月3日,星期二,上午10:15,Etienne Tourigny写道: 您想要的可能超出了gdal的范围。它将需要一些聪明的元数据管理,以便gdal_translate将它们放在一个文件中... 我建议您使用gdal_translate将它们全部转换为netcdf,然后使用python-netcdf4(而不是numpy / scipy中的那个)将它们堆叠在时间维度上。 2013年9月3日,星期二,上午7:55,“ Signell,Richard”写道: 大卫,如果您将问题发布到GIS stackexchange组 /gis//上, 我将提供示例代码,该代码应该会有所帮助。 -丰富 ==================== PDT更新9/3/13 17:04 这是我的输入数据集之一的gdalinfo输出: gdalinfo 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Driver: GTiff/GeoTIFF Files: 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Size is 10666, 13387 Coordinate System is: PROJCS["unnamed", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Polar_Stereographic"], …