lon = ncread(ncFilePath,'lon')'; %读取经度(范围、精度)
lat = ncread(ncFilePath,'lat')'; %读取纬度
P = ncread(ncFilePath,'precipitation'); %读取降雨量
用以上代码将NC4文件转TIF格式过程中,遇到了读取NC文件经纬度以及变量信息失败的问题:
问题是Variables中明明有lon、lat以及precipitation,为什么说找不到呢?
经排查,发现是出在这个NC文件的结构上,看下图:
所有的变量都存储在了组Grid中,而ncread函数对于层级关系需要明确说明,因此代码需改为:
lon = ncread(ncFilePath,'/Grid/lon')'; %读取经度(范围、精度)
lat = ncread(ncFilePath,'/Grid/lat')'; %读取纬度
P = ncread(ncFilePath,'/Grid/precipitation'); %读取月降雨量
如果想要看nc文件的内部结构,可以用函数ncdisp