NC后缀文件转tif,以GOSAT卫星数据为例

后缀名为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');

代码讲解:

此句代码中,需要根据自己的文件更改“  ”中的属性,查询自己文件的属性代码如下:

注:在无法跑通代码时,要想输出此代码,要去掉分号,以达到查看属性后修改代码的目的

       NC后缀文件转tif,以GOSAT卫星数据为例_第1张图片                          NC后缀文件转tif,以GOSAT卫星数据为例_第2张图片

注:如图所见GOAST L4A数据中属性分别是人工排放通量、生物量排放、陆地生态系统地表碳通量、海气通量。Rgn代表了区域,grid代表了1°×1°格网。

此句代码需要查看经纬度范围进行修改,方法同上。

NC后缀文件转tif,以GOSAT卫星数据为例_第3张图片

最终结果图如下:

NC后缀文件转tif,以GOSAT卫星数据为例_第4张图片

提取中国区域结果图如下:

NC后缀文件转tif,以GOSAT卫星数据为例_第5张图片

你可能感兴趣的:(matlab,linq,p2p)