clc

clear

ncFilePath='E:\批處理讀取NC文件-MATLAB語(yǔ)言\1NC文件\VIC_Output_China_evap_1952.nc';

%% 顯示結(jié)構(gòu)

% ncdisp(ncFilePath);%顯示nc文件的所有結(jié)構(gòu)，以便大概了解里面的內(nèi)容

% ncdisp(ncFilePath,'evap');%顯示指定變量的內(nèi)容，注意一定要是變量variables才可以

% ncdisp(ncFilePath,'/','min');%簡(jiǎn)單顯示結(jié)構(gòu)以及定義

% ncdisp(ncFilePath,'/','full');%全部顯示所有結(jié)構(gòu)和定義信息

%% 讀取變量值

% ncid = netcdf.open(ncFilePath,'NOWRITE'); %打開(kāi)nc文件返回索引ID

% [ndims,nvars,ngglobalatts,unlimdimid] = netcdf.inq(ncid);%獲取維數(shù)，變量數(shù)，全局屬性數(shù)量，

% [varname,xtype,dimids,natts] = netcdf.inqVar(ncid,0); %根據(jù)變量索引號(hào)獲取變量的名稱

lon=ncread(ncFilePath,'lon');%讀取經(jīng)度變量

lat=ncread(ncFilePath,'lat');%讀取緯度變量

time=ncread(ncFilePath,'time');%讀取時(shí)間變量

evap=ncread(ncFilePath,'evap');%獲取蒸散發(fā)變量數(shù)據(jù)

evap1=evap(:,:,200);

evap1(evap1==0)=[];

aa=evap1;

[h,l]=find(aa==0);

aa(h,l)=[];

%% 顯示數(shù)據(jù)

% pcolor(lat,lon,evap1);

% [x,y]=meshgrid(lon,lat);%根據(jù)經(jīng)緯度信息產(chǎn)生格網(wǎng)，284列(經(jīng)度)，164列(緯度)

% phandle=pcolor(x,y,evap1');%顯示一個(gè)矩陣，其中x,y,evap1的行列數(shù)必須一致。類(lèi)似surface函數(shù)

% colorbar

% imwrite(evap1','E:\批處理讀取NC文件-MATLAB語(yǔ)言\aa.tif','tif')

%% 保存為地理柵格tif格式，以便arcgis讀取

data=flipud(evap1');

R = georasterref('RasterSize', size(data),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);%地理柵格數(shù)據(jù)參考對(duì)象(類(lèi))

geotiffwrite('E:\批處理讀取NC文件-MATLAB語(yǔ)言\aa4.tif',data,R);

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

以上是生活随笔為你收集整理的matlab中nc文件,MATLAB读取.nc文件的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。