Question

我正在尝试使用GDAL通过python从NOAA NETCDF文件（云的最高高度）中打开并读取一个数组层（必须通过GDAL完成，以便以后可以重新投影数据）。到目前为止，我一直很成功，直到数据格式从“字节”更改为“短”。现在，将gdal.open命令与.ReadAsArray结合使用，将仅返回填充值填充的数据。

from osgeo import gdal # Add the GDAL library
from osgeo import osr
import numpy as np # Import the Numpy package

# Load the Data =======================================================================================
# Path to the GOES-16 image file
file = "\\OR_ABI-L2-ACHAC-M3_G16_s20182342037325_e20182342040098_c20182342041146.nc"
file_dir = 'PATH OF THE NC FILE''
path = file_dir + file


# Open netcdf file.nc with gdal
ds = gdal.Open('NETCDF:"'+path+'":HT') # HT is the layer name

# Read full data from netcdf
data = ds.ReadAsArray()

OUTPUT:
data
array([[-1, -1, -1, ..., -1, -1, -1],
       [-1, -1, -1, ..., -1, -1, -1],
       [-1, -1, -1, ..., -1, -1, -1],
       ...,
       [-1, -1, -1, ..., -1, -1, -1],
       [-1, -1, -1, ..., -1, -1, -1],
       [-1, -1, -1, ..., -1, -1, -1]], dtype=int16)

我似乎无法将NETCFD文件从短格式转换为字节格式，希望对解决此问题有任何帮助。 **我尝试在命令行中使用gdal_translate，但我想我没有正确的语法。

谢谢！

gdal.open使用填充值填充数据

0 个答案: