- 积分
- 49
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-5-20
- 最后登录
- 1970-1-1
|
发表于 2019-4-9 21:38:05
|
显示全部楼层
本帖最后由 yinwangyi01 于 2019-4-10 10:31 编辑
老师,出现了这种问题,我查了所有评论,也发现了回答这个问题的部分,但是我还是不知道要怎么改。
'longitudes' 的值无效。它必须满足函数: @(x)isvector(x)&&length(x)==2。
我的程序是
startLoc = [1000 300 1]; % 起始位置
count = [283 161 1]; % 读取的数据数
stride = [1 1 1];
lon=ncread('C:\Users\YLL\Documents\Tencent Files\2763361736\FileRecv\GLDAS_NOAH025_v2.0_2002\GLDAS_NOAH025_M.A200201.020.nc4','lon',1000,283);%读取lon
lat=ncread('C:\Users\YLL\Documents\Tencent Files\2763361736\FileRecv\GLDAS_NOAH025_v2.0_2002\GLDAS_NOAH025_M.A200201.020.nc4','lat',300,161);%读取lat
%开始读取每年数据
%2006年地表温度空间平均
AvgSurfT_inst200601=ncread('C:\Users\YLL\Documents\Tencent Files\2763361736\FileRecv\GLDAS_NOAH025_v2.0_2006\GLDAS_NOAH025_M.A200601.020.nc4','AvgSurfT_inst',startLoc,count,stride);
然后直接添加了
%绘制空间图像
figure
hold on;
[X,Y]=meshgrid(lat,lon);
[cs,h]=contourf(Y,X,AvgSurfT_inst200601);
colorbar;
maskmap('C:\Users\YLL\Downloads\china\china_basic_map\bou2_4p.shp',true,'lon',lon,'lat',lat,'linewidth',0.5, 'edgecolor','b')
|
|