nc文件转tif代码

最新推荐文章于 2023-11-16 13:46:33 发布

dropout_z

最新推荐文章于 2023-11-16 13:46:33 发布

阅读量290

点赞数

CC 4.0 BY-SA版权

分类专栏：遥感文章标签：图像处理

本文链接：https://blog.youkuaiyun.com/jiekeheiguanglan/article/details/130178181

遥感专栏收录该内容

15 篇文章

订阅专栏

该代码段使用Matlab读取NetCDF文件中的经度、纬度、时间及覆盖度数据，通过imagesc和pcolor展示数据，然后将数据转换为tif格式，利用georasterref和geotiffwrite进行地理栅格化和保存。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

clc
clear 
ncFilePath = ['D:\homework\nc\c_gls_FCOVER-RT0_202006300000_GLOBE_PROBAV_V2.0.1.nc']; %设定NC路径
ncdisp(ncFilePath); %查看nc文件变量
 %% 读取变量值
    %根据ncdisp函数读取到的nc文件变量相应替换
    lon=ncread(ncFilePath,'lon');     %读取经度信息（范围、精度）
    lat=ncread(ncFilePath,'lat');      %读取维度信息
    time=ncread(ncFilePath,'time');         %读取时间序列
    DUS=ncread(ncFilePath,'FCOVER');   %获取覆盖度数据 
 %% 展示数据内部结构等信息
     %figure(i);
     imagesc(lat,lon,DUS);
     shading flat;                        %移除网格线，否则图上一片黑什么都没有
     [x,y]=meshgrid(lon,lat);             %根据经纬度信息产生格网，3600列（经度），1800列（纬度）
     %phandle=pcolor(x,y,DUS');        %显示一个矩阵，其中x,y,DUS的行列数必须一致
     shading flat;
     colorbar
%% 存为tif格式
    data=flipud(rot90(DUS,1));   %镜像反转，不反转的话最后的图像的南北朝向是错的
                                 %逆时针旋转90°*n，此时n为1，可自行设置调整
    data1=flipud(data);  %图像上下翻转
    %georasterref 绘制地理栅格数据（raster data）；double双精度浮点数，保留16位有效数字
    R = georasterref('RasterSize', size(data1),'Latlim', [double(min(lat)) double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
    %geotiffwrite geotiff栅格影像保存
    geotiffwrite(['D:\homework\nc\fcover.tif'],data1,R);
disp('Done!')