同潮图等值线在岸上的问题

飞翔的小熊

画同潮图时插值之后画出来的图形一部分等值线都在岸上了。因为原振幅与迟角数据皆为一维的数据，所以用了griddata函数进行插值，但是画出来的图形如下图所示。求助给位大神有什么解决办法吗？最后，附上我的程序，请大家指正。


lon=ncread('sbqt_0001.nc','lon');
lat=ncread('sbqt_0001.nc','lat');
h=ncread('sbqt_0001.nc','zeta');
lon1=double(linspace(min(lon),max(lon),102));
lat1=double(linspace(min(lat),max(lat),102));
lon = double(lon);lat = double(lat);
m_proj('miller','lat',[29 34],'lon',[119 124]);
m_grid('linewidh',1,'xaxislocation','bottom');
m_gshhs_f('patch',[0 0 0],'edgecolor','k');
hold on;
theta=[];
for i=1:77326;
hh=h(i,:);
tide_names='M2';
[NAME,FREQ,TIDECON,XOUT]=t_tide(hh,'latitude',33,'start time',[2015,11,01,0,0,0],'rayleigh',tide_names);
theta(i,1)=TIDECON(:,3);
end
[Lon1,Lat1,theta1]=griddata(lon,lat,theta,lon1,lat1');
cc=0:30:360;
[cs,h]=m_contour(Lon1,Lat1,theta1,cc,'b');
clabel(cs,h,'manual');
xlabel('Longitude');
ylabel('Latitude');

二爷名声在外

先画等值线，再画海岸

飞翔的小熊

二爷名声在外 发表于 2017-11-18 18:51
先画等值线，再画海岸

我试过了，不能直接覆盖掉

飞翔的小熊

求助各位matlab大神~

二爷名声在外

那你就继续提问吧……