后缀名为nc的文件,是数控程序的后缀。一些温度、降雨等数据的文件格式后缀就是nc。本文使用matlab将其转换为tif格式的文件,以便于在Arcgis中进行分析。
代码如下:
%********************************%
%微信公众号:GIS攻略
%********************************%
clc
clear
% nc文件所在根目录
Path='E:\GOAST_L4A\GOASTL4A\';
% 输出文件根目录
tifPath='E:\GOAST_L4A\tif\';
% 获取文件夹下所有nc后缀文件
imageList=dir(strcat(Path,'*.nc'));
for fileindex=1:length(imageList)
% 完整文件名
filename=imageList(fileindex).name;
% 去掉后缀后的文件名
filename1=filename(1:find(filename=='.')-1);
% nc文件的完整路径
fileFullPath=strcat(Path,filename);
ncinf = ncinfo(fileFullPath);
ETaSets=ncread(fileFullPath,'all_flx_grd_opt');
SizeInfo=size(ETaSets);
for subsetIndex=1:SizeInfo(3)
disp(subsetIndex);
MonthData=ETaSets(:,:,subsetIndex);
MonthData(isnan(MonthData))=-8888;
tifOutputFullPath=strcat(tifPath,filename1,num2str(subsetIndex,'%02d'),'.tif');
MonthData = rot90(MonthData);
Refference=georasterref('RasterSize',size(MonthData),'Latlim',[-89.5 89.5],'Lonlim',[-179.5 179.5]);
Refference.ColumnsStartFrom = 'north';
geotiffwrite(tifOutputFullPath,MonthData,Refference);
end
end
disp('finish');
代码讲解:
此句代码中,需要根据自己的文件更改“ ”中的属性,查询自己文件的属性代码如下:
注:在无法跑通代码时,要想输出此代码,要去掉分号,以达到查看属性后修改代码的目的
注:如图所见GOAST L4A数据中属性分别是人工排放通量、生物量排放、陆地生态系统地表碳通量、海气通量。Rgn代表了区域,grid代表了1°×1°格网。
此句代码需要查看经纬度范围进行修改,方法同上。
最终结果图如下:
提取中国区域结果图如下: