clc;
clear;
close all
% 读取nc数据
filename=['/Users/mayuetong/Documents/eee/resilience2/dataset-derived-near-surface-meteorological-variables-fedb5834-831b-45f6-b5ec-5fc7bfd7c407' ...
'/Wind_WFDE5_CRU_201809_v2.1.nc'];
ncdisp(filename);
%展示文件信息%
% 提取变量
lon = ncread(filename,'lon')';
lat = ncread(filename,'lat');%读取纬度范围和精度
wind = ncread(filename,'Wind'); %读取风速范围和精度
time=length(ncread(filename,'time'));%读取时间序列长度
%data=tmp(:,:,1336);%为方便,本例只绘制2020年8月份温度数据
data=flipud(rot90(wind,1)); %镜像反转,不反转的话最后的图像的南北朝向是错的,逆时针旋转90°%
% 绘图
R = georasterref('RasterSize', size(data),'Latlim', [double(min(lat)) ...
double(max(lat))], 'Lonlim', [double(min(lon)) double(max(lon))]);
% 不需要改动,关于georasterref函数,请查看官网介绍% 保存绘制的地理栅格图像
geotiffwrite(['/Users/mayuetong/Documents/eee/resilience2/数据处理/shanzhuanalysis/test2_wind_201809','.tif'],data,R);