matlab：nc转tif

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

clc;                                     %清除命令行窗口
clear;                                   %清除工作区

filename='F:\Global_P_ET\ET\GLEAM\1982\E_1982_GLEAM_v3.2a.nc';
ncdisp (filename);                      %展示文件信息

%% 提取变量
lon = ncread(filename,'lon')';           %读取经度范围和精度
lat = ncread(filename,'lat');            %读取纬度范围和精度
time = length(ncread(filename,'time'));  %读取时间序列长度

data = nc_varget(filename,'E');          %提取变量

data = flipud(permute(data,[3 2 1]));   
%permute置换数组维度，flipud垂直翻转（往往拿到的nc文件的维度方向
    ...不是我们传统上习惯的【经度,纬度,时间】格式，所以要调整）
data = data(:,:,1);                      %为方便，本例只绘制第一层的数据

%% 绘图
R = georasterref('RasterSize', size(data),'Latlim', [double(min(lat))...
    double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]); ...
    %不需要改动，关于georasterref函数，请查看官网介绍
geotiffwrite(['F:\example','.tif'],data,R);     %保存绘制的地理栅格图像

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB matla atlab Atl Lab matlab 操作策略