python - 在 Data.VRT 中为 GDALTranslaste 参考纬度/经度 VRT 文件

Question

作为这个问题的后续行动

我有三个栅格：

• 带有纬度值的 lat.tif

• 具有经度值的 lon.tif

• 带有数据值的 data.tif（对应于纬度/经度值）

从那里我创建了 lat.vrt、lon.vrt 和 data.vrt 文件。

如何引用 data.vrt 中的两个 lat.vrt 和 lon.vrt 文件？

参考上一个问题：我将如何“在 data.vrt 元数据中的 geoarrays 中”编写两个 lat/lon .vrt 文件？

标签： pythongdalgeotiff

Answer 1

弄清楚了！

我将以下内容添加到 data.vrt 文件以引用 lat.vrt 和 lon.vrt。如果找不到 lat.vrt/lon.vrt 添加路径：

<Metadata domain="GEOLOCATION">
      <MDI key="X_DATASET">lon_1.vrt</MDI>
      <MDI key="X_BAND">1</MDI>
      <MDI key="Y_DATASET">lat_1.vrt</MDI>
      <MDI key="Y_BAND">1</MDI>
      <MDI key="PIXEL_OFFSET">0</MDI>
      <MDI key="LINE_OFFSET">0</MDI>
      <MDI key="PIXEL_STEP">1</MDI>
      <MDI key="LINE_STEP">1</MDI>    
 </Metadata>