利用MATLAB读取.nc文件单像元数值并转为Excel格式（以中国日降雨量月均数据为例）

lemonade723

已于 2022-08-02 11:22:43 修改

阅读量1.1w

点赞数 15

文章标签： matlab

于 2022-07-22 13:08:02 首次发布

本文链接：https://blog.csdn.net/lemonade723/article/details/125929645

版权

该文介绍如何在MATLAB中处理NetCDF文件，通过经纬度定位特定像元，并批量提取多个月份的降雨量数据。首先确定经纬度对应的行列号，然后利用for循环提取指定位置的每月数据，并将其保存为Excel文件。文中提出两个问题：一是如何不依赖经纬度范围，而是直接用find函数批量提取多经纬度的行列号；二是如何处理大范围数据的批量提取。

摘要由CSDN通过智能技术生成

以中国日降雨量月均数据（nc文件包含12月）为例，提取某经纬度下的多月份像元值。

(【数据分享】1960-2020年中国1公里分辨率月降水数据集)

一、确定经纬度所在行列号

以92.18E，30.475N为例，首先在Matlab中输入以下代码：

ncdisp('CN_Prec_MonthlyMean_1km_2019.nc');		%读取nc文件信息

lat = ncread('CN_Prec_MonthlyMean_1km_2019.nc','latitude');
lon = ncread('CN_Prec_MonthlyMean_1km_2019.nc','longitude');
pre = ncread('CN_Prec_MonthlyMean_1km_2019.nc','pre') ;

工作区获取到lat、lon和pre的信息，打开lat和lon文件

92.18E定位至92.1754~92.1838，确定的单元格为92.1754所在的2248（选择数值小的单元格）；

30.475N定位至30.4718~30.4634，确定的单元格为30.4718所在的2771（选择数值小的单元格）；

二、提取多月单像元数值

将得到的行列号输入至以下代码，得到多月单像元数值——2019.xls：

ncdisp('CN_Prec_MonthlyMean_1km_2019.nc');		%读取nc文件信息

lat = ncread('CN_Prec_MonthlyMean_1km_2019.nc','latitude');
lon = ncread('CN_Prec_MonthlyMean_1km_2019.nc','longitude');
pre = ncread('CN_Prec_MonthlyMean_1km_2019.nc','pre') ;

pre1= [];
for i = [1:12]
    pre1 = [pre1;pre(2771,2248,i)];
end

xlswrite('2019',pre1)