?以中國日降雨量月均數(shù)據(jù)(nc文件包含12月)為例,提取某經(jīng)緯度下的多月份像元值。
(【數(shù)據(jù)分享】1960-2020年中國1公里分辨率月降水?dāng)?shù)據(jù)集)
一、確定經(jīng)緯度所在行列號(hào)
以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') ;
工作區(qū)獲取到lat、lon和pre的信息,打開lat和lon文件
92.18E定位至92.1754~92.1838,確定的單元格為92.1754所在的2248(選擇數(shù)值小的單元格);
30.475N定位至30.4718~30.4634,確定的單元格為30.4718所在的2771(選擇數(shù)值小的單元格);
二、提取多月單像元數(shù)值
將得到的行列號(hào)輸入至以下代碼,得到多月單像元數(shù)值——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)
完畢。?
疑問待解決:
1、如何不通過經(jīng)緯度范圍,使用find函數(shù)批量提取多經(jīng)緯度對(duì)應(yīng)的行列號(hào)?
2、如何批量提取大范圍的數(shù)值?
參考鏈接:
NetCDF(.nc)文件數(shù)據(jù)的讀取_嗶哩嗶哩_bilibili
5分鐘!利用MATLAB讀取.nc(netCDF)文件并轉(zhuǎn)Excel格式_嗶哩嗶哩_bilibili文章來源:http://www.zghlxwxcb.cn/news/detail-468788.html
請(qǐng)教:matlab提取數(shù)據(jù)文件中某個(gè)經(jīng)緯度的數(shù)據(jù)?謝謝! – MATLAB中文論壇文章來源地址http://www.zghlxwxcb.cn/news/detail-468788.html
到了這里,關(guān)于利用MATLAB讀取.nc文件單像元數(shù)值并轉(zhuǎn)為Excel格式(以中國日降雨量月均數(shù)據(jù)為例)的文章就介紹完了。如果您還想了解更多內(nèi)容,請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章,希望大家以后多多支持TOY模板網(wǎng)!