netcdf geotiff java_在python中提取特定的netcdf信息并转换为GeoTIFF

我试图从netCDF文件中提取一组特定的数据，然后将这些数据转换为GeoTIFF。在

到目前为止，我已经成功地使用netCDF4提取了我想要的数据，文件中的所有数据都存储为1d数组(lat、lon、我想要的数据)，并将它们分配给一个2d数组。我正在使用的netcdf文件已子集到特定区域。但从这里我不知所措。在

通过以下链接，我对geotiff转换的工作原理有一点了解：

以下是我目前的情况：import netCDF4

import numpy as np

from osgeo import gdal

from osgeo import osr

#Reading in data from files and extracting said data

ncfile = netCDF4.dataset("data.nc", 'r')

dataw = ncfile.variables["dataw"][:]

lat = ncfile.variables["Latitude"][:]

long = ncfile.variables["Longitude"][:]

n = len(dataw)

x = np.zeros((n,3), float)

x[:,0] = long[:]

x[:,1] = lat[:]

x[:,2] = dataw[:]

nx = len(long)

ny = len(lat)

xmin, ymin, xmax, ymax = [long.min(), lat.min(), long.max(), lat.max()]

xres = (xmax - xmin) / float(nx)

yres = (ymax - ymin) / float(ny)

geotransform = (xmin, xres, 0, ymax, 0, -yres)

#Creates 1 raster band file

dst_ds = gdal.GetDriverByName('GTiff').Create('myGeoTIFF.tif', ny, nx, 1, gdal.GDT_Float32)

dst_ds.SetGeoTransform(geotransform) # specify coords

srs = osr.SpatialReference() # establish encoding

srs.ImportFromEPSG(3857) # WGS84 lat/long

dst_ds.SetProjection(srs.ExportToWkt()) # export coords to file

dst_ds.GetRasterBand(1).WriteArray(x) # write r-band to the raster

dst_ds.FlushCache() # write to disk

dst_ds = None # save, close

以上geotiff的创作过程主要来源于：

我尝试这样做是基于这样一个理解：我要写入光栅的数组是一个包含3列、2个单词和1个数据的2d数组。我在snap中检查的结果是一个黑色页面，左侧有一条白线。在

所以我的问题是：

如何从netcdf文件中提取必要的geotransform数据，适当地调整geotransform参数，然后使用提取的lat long+dataw数组写入geotiff文件？在