求助大神:海表面温度绘等温线图问题

wjh0123

有两个问题不解:使用的月均OISST.V2数据:sst.mnmean.nc+Ismask.nc
1.等温线图中陆地上也有SST值,如何消除?(使用了Ismask.nc,但是在绘图区域121.5-136.5E,22-34N内显示的均是海洋)
2.如何将等温线的凸点连接起来,并得到凸点的经纬度坐标?
绘制的图(其中一张):

代码如下:
clc;clear;clf;
ncdisp('sst.mnmean.nc');
sst=ncread('sst.mnmean.nc','sst');
lon=ncread('sst.mnmean.nc','lon');lon1=lon(122:137);
lat=ncread('sst.mnmean.nc','lat');lat1=lat(56:69);
t=ncread('sst.mnmean.nc','time');t1=t(2:433);
mask=ncread('lsmask.nc','mask');
[I,J]=find(mask==0);
for i=1:length(I)
        sst(I(i),J(i),:)=nan;
end;
sst1=sst(122:137,56:69,2:433);
for i=1:12
    sstm(:,:,i)=mean(sst1(:,:,i:12:end),3);
end
for i=1:12
figure(j+i)
m_proj('Equidistant Cylindrical','lon',[121.5,136.5],'lat',[21.5,34.5]);
m_coast('patch',[.01 .01 .01],'edgecolor','none');
m_grid('xlabeldir','end','fontsize',10);
hold on
[cs,h]=m_contourf(lon1,lat1,sstm(:,:,i)','levelstep',0.5,'linestyle','--','linewidth',0.1);
colorbar;
clabel(cs,h,'labelspacing',300);
hold off
saveas(gcf,num2str(i),'bmp');
end

二爷名声在外

1.陆地数据可以使用掩膜进行操作，具体示例本版里面有
2.关于凸点的寻找，需要使用相关算法自行设计。具体思路无非就是寻找极值，楼主可以思考一下

wjh0123

rabin_xu 发表于 2018-3-13 11:34
你可以看看。不会是纬度方向返了吗？

是由于数据分辨率较低导致的,用高精度数据会好很多

rabin_xu

我没用过OISST的月值场，但是从海岛周围海温分布这么连续来看，总感觉数据有问题，海岛周围应该不会这么连续，你在仔细看看数据

11zttt

                                  同问

wjh0123

rabin_xu 发表于 2018-3-13 09:40
我没用过OISST的月值场，但是从海岛周围海温分布这么连续来看，总感觉数据有问题，海岛周围应该不会这么连 ...

我也感觉数据存在问题 我再下个日均数据看看会不会好一点

rabin_xu

wjh0123 发表于 2018-3-13 09:58
我也感觉数据存在问题 我再下个日均数据看看会不会好一点

你可以看看。不会是纬度方向返了吗？

wjh0123

rabin_xu 发表于 2018-3-13 11:34
你可以看看。不会是纬度方向返了吗？

已解决,是数据精度不够导致的

rabin_xu

wjh0123 发表于 2018-3-18 15:48
已解决,是数据精度不够导致的

看看现在的图呗？

wjh0123

rabin_xu 发表于 2018-3-19 09:03
看看现在的图呗？

基本没什么变化,就是抠出了缺测值