pandas - 将 Pandas 数据框转换为 geoTIFF (python)

Question

我有一个带有 X/Y 和纬度/经度坐标的 pandas df。

我想使用纬度/经度列转换数据框并将 TIFF 图像存储在 WGS84 CRS 中。

谢谢

Answer 1

为您推荐几个套餐:xarray和 rioxarray .

xarray 是 pydata 对带标签的 N 维数组的解决方案(想想 pandas，但在 3 维以上，或者想想 numpy ND 数组，但使用 pandas 索引而不仅仅是位置索引)。

rioxarray 是一个结合了 xarray 和 rasterio 的扩展包，提供了读取和写入包括 GeoTIFF 在内的光栅文件的能力。 rioxarray 在 converting xarray DataArrays to rasters 上有文档.另请参阅用于转换的 API 文档 RasterArray和 RasterDataset对象到栅格。

在您的情况下，假设您的正交维度是 (lat, lon) 并且 model_lat 和 model_lon 实际上都由 lat 索引和 lon(例如，它们处于 3D 投影中)，并且 res 是您要编码的波段，您的结果将如下所示:

import xarray as xr
import rioxarray

da = df.set_index(['lat', 'lon']).to_xarray()

# promote the data variables model lat/long to 2d coordinates
da = da.set_coords(['model_lat', 'model_long'])

da.res.rio.to_raster(filepath)