- 积分
- 48
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-9-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
新手写了个程序,但是有些命令搞不懂...附的这个程序中,提示出错,但是我不会改...刚接触功底太差,错误很低级...但是百度神马的看命令不理解...望批评指正
clear; %从内存中清除变量和函数
filesample=hdfinfo('A20132742013304.L3m_MO_SST4_9.hdf') %将sss文件赋值给'filesample'
data1 = hdfread(filesample.SDS(1)); %读取filesample文件
data1=data1*0.00071718-2; %新赋值
data1(data1>44)=nan; %将data1中大于44的值赋予nan
data1=flipud(data1); %上下翻转数据
lat=-89.958333:0.8333:89.958333; %设置图像纬度范围
lon=-179.958333:0.8333:179.958333; %设置图像经度范围
[lon1,lat1]=meshgrid(lon,lat); %生成网格
lonmin=105;
lonmax=118;
latmin=4;
latmax=21;
goodlons=lon>=lonmin & lon<=lonmax;
goodlats=lat>=latmin & lat<=latmax;
subset_data1=data1(goodlats,goodlons); %取出绘图区域范围
subset_data1=double(subset_data1); %双精度
m_proj('miller','lon',[105,118],'lat',[4,21]); %绘图区域为南海范围
m_pcolor(lon1,lat1,subset_data1); %使用m_pcolor绘图
shading flat; %去掉网格
m_coast('patch','w'); %绘制海岸线,陆地用白色覆盖
m_grid('linestyle','none','box','fancy','tickdir','out'); %绘制网格
caxis([-2,1]); %控制伪彩色坐标刻度
set(colorbar,'ytick',[-2,-1,0,1],'yticklabel',{'0.01','0.1','1','10'}); %设置色标
title(['modis[sst] 201304 South China Sea']); %图形标题
print(gcf,'-dpng','e:\matlabdata\modis_[sst]_201304.png') %图形保存路径
提示错误:
|
|