ncfilepath=『d:\田晴1\esa cci\20141002-20160430active\esacci-soilmoisture-l3s-ssms-active-20141002000000-fv04.4.nc』;%% 顯示結構
% ncdisp(ncfilepath);%在命令列視窗中以文字形式顯示 netcdf 資料來源 source 中的所有組、維度、變數定義，以及所有屬性。
% ncdisp(ncfilepath,『evap』);%顯示指定變數的內容.
% ncdisp(ncfilepath,』/』,『min』);%只顯示示例檔案 example.nc 的組層次結構和變數定義。
% ncdisp(ncfilepath,』/』,『full』);%全部顯示所有結構和定義資訊
%% 讀取變數值
% ncid = netcdf.open(ncfilepath,『nowrite』); %開啟nc檔案返回索引id
% [ndims,nvars,ngglobalatts,unlimdimid] = netcdf.inq(ncid);%獲取維數，變數數，全域性屬性數量，
% [varname,xtype,dimids,natts] = netcdf.inqvar(ncid,0); %根據變數索引號獲取變數的名稱
lon=ncread(ncfilepath,『lon』);%讀取經度變數
lat=ncread(ncfilepath,『lat』);%讀取緯度變數
sm_data=ncread(ncfilepath,『sm』);%讀取sm變數
sm_data1=sm_data(1174:1194,240:262);
sm_data1(sm_data10)=;
aa=sm_data1;
[h,l]=find(aa0);
aa(h,l)=;
%% 顯示資料
% pcolor(lat,lon,sm_data);pcolor(x,y,c) 指定頂點的 x 座標和 y 座標,c 的大小必須與 x-y 座標網格的大小匹配,例如，如果 x 和 y 定義乙個 m×n 網格，則 c 必須為 m×n 矩陣.
% [x,y]=meshgrid(lon,lat);%根據經緯度資訊產生格網.
% phandle=pcolor(x,y,sm_data』);%顯示乙個矩陣，其中x,y,sm_data的行列數必須一致。類似su***ce函式.
% colorbar
% imwrite(sm_data1,『d:\田晴1\esa cci\20141002-20160430active\aa.tif』,『tif』)
%% 儲存為地理柵格tif格式，以便arcgis讀取.
data=flipud(sm_data1);
r = georasterref(『rastersize』, size(data),『latlim』, [double(min(lat)) double(max(lat))], 『lonlim』, [double(min(lon)) double(max(lon))]);%地理柵格資料參考物件(類)
geotiffwrite(『d:\田晴1\esa cci\20141002-20160430active\aa.tif』,data,r);
MATLAB開啟nc檔案並讀取nc檔案資料
ncdisp e nc 輸入該命令後命令列視窗會顯示該nc檔案的基本資訊 根據上一步中nc檔案的基本資訊提取nc檔案的相應資料。time ncread e nc time 首先，將所需要的資料提取出來 start 125,89,123 count 17,14,124 tmax ncread e nc...
製作nc檔案（Matlab）
首先看乙個nc檔案中包含哪些部分，例如乙個標準的fvcom輸入檔案wind.nc netcdf wind可以看出nc檔案包含兩個主要部分 dimensions 各個變數維度大小 variables 變數 在上面檔案中包含三個維度 nele，node，time。其中nele，node兩個是固定長度的，...
Python讀取和處理nc檔案
首先進行一下更新，在anaconda prompt 開始選單欄 anaconda3檔案下 中輸入命令 conda update all 然後安裝讀nc檔案所需的模組netcdf4 conda install netcdf4 這樣就完成了hdf4模組的安裝。coding utf 8 spyder ed...