从odo.resource(源)获取数据到odo.resource(target)

Question

我正在尝试使用功能扩展ODO库，以将GDAL数据集（带空间信息的栅格）转换为NetCDF文件。

在GDAL数据集中读取良好。但是在NETCDF的创建阶段，我需要GDAL数据集的一些元数据（在调用odo.odo（源，目标）时尚不知道的元数据）。我怎么能实现这一目标？

到目前为止，我的代码的简短版本：

import odo
from odo import resource, append
import gdal
import netCDF4 as nc4
import numpy as np
@resource.register('.+.tif')
def resource_gdal(uri, **kwargs):
    ds = gdal.Open(uri)
    # metadata I need to transfer to netcdf
    b = ds.GetGeoTransform() #bbox, interval
   return ds

@resource.register('.+.nc')
def resource_netcdf(uri, dshape=None, **kwargs):
    ds = nc4.Dataset(uri,'w')
    # create lat lon dimensions and variables
    ds.createDimension(lat, dshape[0].val)
    ds.createDimension(lon, dshape[1].val)
    lat = ds.createVariable('lat','f4', ('lat',))
    lon = ds.createVariable('lon','f4', ('lon',))
    # create a range from the **gdal metadata**
    lat_array = np.arange(dshape[0].val)*b[1]+b[0]
    lon_array = np.arange(dshape[1].val)*b[5]+b[3]
    # assign the range to the netcdf variable
    lat[:] = lat_array
    lon[:] = lon_array
    # create the variable which will hold the gdal data
    data = ds.createVariable('data', 'f4', ('lat', 'lon',))
   return data

@append.register(nc4.Variable, gdal.Dataset)
def append_gdal_to_nc4(tgt, src, **kwargs):
    arr = src.ReadAsArray()
    tgt[:] = arr
    return tgt

谢谢！

Answer 1

我对odo没有太多经验，但是从浏览源代码和文档来看，看来resource_netcdf()不应参与将GDAL数据转换为NetCDF。翻译应是convert.register装饰的gdal_to_netcdf()功能的工作。在这种情况下，由resource_gdal返回的gdal.Dataset对象将具有所有足够的信息（地理标准，像素大小）来制作NetCDF。