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);     %保存绘制的地理栅格图像