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');

代码讲解:

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

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

                                 

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

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

最终结果图如下:

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

评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值