depth_image.zip_ETOPO1_depth image_site:www.pudn.com_水深_水深线

clear;clc; % ncdisp('ETOPO1_Bed_g_gmt4.grd'); lon=ncread('ETOPO1_Bed_g_gmt4.grd','x'); lat=ncread('ETOPO1_Bed_g_gmt4.grd','y'); d=ncread('ETOPO1_Bed_g_gmt4.grd','z'); de=double(d); %% 研究区域 % m=find(lat<=-66.2&lat>=-66.6); % n=find(lon<=110.6&lon>=110); m=find(lat<=-65.7&lat>=-66.3); n=find(lon<=110&lon>=107); dep=de(n,m); LON=lon(n); LAT=lat(m); [PT,PD]=meshgrid(LAT,LON); %% 文森湾区域 m=find(lat<=-64&lat>=-66.9); n=find(lon<=114&lon>=104); dep=de(n,m); LON=lon(n); LAT=lat(m); [PT,PD]=meshgrid(LAT,LON); dep(find(dep(:)>-10))=nan; %% 南极大陆 coast_antarctic=load('Antarctic.dat'); figure bndry_lon=[104,104,114,114,104];bndry_lat=[-64,-67.2,-67.2,-64,-64]; % m_proj('lambert','long',[109 111],'lat',[-67 -66]) m_proj('stereographic','lat',-90,'radius',30,'rotangle',0); % m_contourf(PD,PT,dep) % shading interp % hold on % m_gshhs_f('patch',[0.7 0.7 0.7],'edgecolor','k'); m_patch(coast_antarctic(:,1),coast_antarctic(:,2),[.7 .7 .7],'edgecolor','none'); m_line(bndry_lon,bndry_lat,'linewi',1,'color','r'); % Area outline ... % m_gshhs_f('speckle','color','k'); % with speckle added m_grid('xtick',12,'ytick',5,'linewidh',1,'xaxislocation','top','backcolor',[1 1 1]); % set(gca,'ylim',[30 90],'ytick',(30:30:90)) set(get(gca,'xLabel'),'Fontsize',2,'fontweight','bold') % colormap(jet) % colorbar set(gcf,'PaperUnits','inches','PaperPosition',[0 0 5 5]) % print(gcf, '-dpng','-r600',['F:\文生湾海豹CTD项目\data\AUSTRALIA\figures\论文用图2\Antarctic_area.png']) % print(gcf, '-depsc2','-r600',['F:\文生湾海豹CTD项目\data\AUSTRALIA\figures\论文用图2\Antarctic_area.eps']) %% 提取展位数据 load('E:\文生湾海豹CTD项目\data\AUSTRALIA\afterdata\data8\Lon.mat'); load('E:\文生湾海豹CTD项目\data\AUSTRALIA\afterdata\data8\Lat.mat'); load('E:\文生湾海豹CTD项目\data\AUSTRALIA\afterdata\data8\DATE.mat'); X=[Lon,Lat,DATE];X(DATE(:)<3010000|DATE(:)>5010000,:)=[]; X=unique(X,'rows');%去除重复项 X=sortrows(X,3); S=X(:,3)+2012e+07; S=datevec(num2str(S),'yyyymddHHMM'); S=datenum(S); %% coast=load('coast_shelf.dat'); figure % m_proj('miller','long',[109 111],'lat',[-67 -66]);%研究区域范围 m_proj('miller','lon',[104 114],'lat',[-67.2 -64]);%文森湾范围 % [c,h]=m_contourf(PD,PT,dep,[0,-250,-500,-1000,-1500,-2000,-2500,-3000,-3500],'ShowText','on'); % % % [c,h]=m_contourf(PD,PT,dep,[-100,-200,-300,-400,-500,-1000],'ShowText','on'); % clabel(c,h,'FontSize',10,'Color','black','labelSpacing',300); % set(h,'linestyle','-'); % m_pcolor(PD,PT,dep); % shading interp %m_contourf(PD,PT,dep,20,'linestyle','none'); % a=[195 226 214;157 211 207;119 199 214;79 189 211;53 161 194;46 122 173;35 67 114; 31 42 74; 15 19 37]/255;%colormap % colormap(flipud(a)); colormap(bone); freezeColors; % colormap([bone(64);jet(64)]) [c,h]=m_contourf(PD,PT,dep,[-250,-500,-1000,-1500,-2000,-2500,-3000,-3500],'k','ShowText','on'); % h1 = get(handle) % [c,h]=m_contourf(PD,PT,dep,[-100,-200,-300,-400,-500,-1000],'ShowText','on'); clabel(c,h,'FontSize',10,'Color','black','labelSpacing',300); % set(h,'linewidth',1,'linestyle','-'); % co=colorbar; % set(get(co,'ylabel'),'String','bathymetry/m','fontsize',12,'fontweight','bold') % % set(co,'ylim',[-3500,-250],'yTick',[-250:500:-3500],'yticklabel',{'-250','-750','-1250'},'fontweight','bold'); % set(gca,'clim',[-3500,-250]); %设定colorbar范围 % cbfreeze(co); % b=colormap(bone);freezeColors(b); hold on colormap(jet); freezeColors; h2=m_scatter(X(:,1),X(:,2),5,S,'filled'); caxis([min(S),max(S)]); x=colorbar('south'); get(x,'position'); set(x,'position',[0.68,0.28,0.18,0.015]); set(gca,'clim',[min(S) max(S)]); set(x,'ylim',[min(S),max(S)],'yTick',[min(S):(max(S)-min(S))/2:max(S)],'yticklabel',{'3/1' '3/31' '4/30'},'fontweight','bold'); cbfreeze(x); % m_gshhs_f('patch',[.6 .6 .6]) m_gshhs_f('speckle','color','k'); % with speckle added m_patch(coast(:,1),coast(:,2),[.7 .7 .7]); m_grid('linewidh',1,'box','fancy','xaxislocation','bottom'); xlabel('Longitude','fontsize', 15); %控制坐标轴名称的字体大小 ylabel('Latitude','fontsize', 15); %控制坐标轴名称的字体大小 hold on axes('position',[0.12,0.25,0.2,0.2]); coast_antarctic=load('Antarctic.dat'); % figure bndry_lon=[104,104,114,114,104];bndry_lat=[-64,-67.2,-67.2,-64,-64]; % m_proj('lambert','long',[109 111],'lat',[-67 -66]) m_proj('stereographic','lat',-90,'radius',30,'rotangle',0); % m_contourf(PD,PT,dep) % shading interp % hold on % m_gshhs_f('patch',[0.7 0.7 0.7],'edgecolor','k'); m_patch(coast_antarctic(:,1),coast_antarctic(:,2),[.7 .7 .7],'edgecolor','none'); m_line(bndry_lon,bndry_lat,'linewi',1,'color','r'); % Area outline ... % m_gshhs_f('speckle','color','k'); % with speckle added m_grid('xtick',12,'ytick',4,'linewidth',1,'xaxislocation','top','backcolor',[1 1 1],'xticklable',{ },'yticklable',{ }); % set(gca,'ylim',[30 90],'ytick',(30:30:90)) % colormap(jet) % colorbar % set(get(x,'ylabel'),'String','date','Fontsize',14,'fontweight','bold') % set(get(x,'ylabel'),'position','Southeast','String','Date'); % cbfreeze(colorbar); % set(gca,'clim',[min(S) max(S)]); % set(x,'ylim',[min(S),max(S)],'yTick',[min(S):(max(S)-min(S))/9:max(S)],'yticklabel',{'3/1' '3/5'... % '3/10' '3/15' '3/20' '3/25' '3/30' '4/5' '4/10' '4/15'}); set(gcf,'PaperUnits','inches','PaperPosition',[0 0 9 6]) % print(gcf, '-dpng','-r600',['F:\文生湾海豹CTD项目\data\AUSTRALIA\figures\论文用图2\文森湾2.png']) % print(gcf, '-depsc2','-r600',['F:\文生湾海豹CTD项目\data\AUSTRALIA\figures\论文用图2\文森湾.eps'])