matlab对等深线加单位

matlab对等深线加单位
我试了好久，能做，但不够成熟！

第一种情况：标注某一条直线：自己编写代码，更改了0.4m（波高数据）

第二种：都加上单位,使用函数

虽然地图也能行，但是呢，不能够旋转和直线一个方向，这就是我认为不够成熟的地方。

代码放进来，一起学习吧！

clear;clc;close all
%% 微信公众号：海洋与大气科学
%% 时间：2022年11月28号
%% 1：读取数据：波高数据；
datahs=load('hs.mat');
hs=datahs.hs_new;
clear datahs
hs=hs-0.6;
[x,y]=meshgrid(linspace(117,127,81),linspace(33,41,65));
%% 2 颜色设定
% colormore_2
% 调用格式三种：
% colormore_2(1)
% colormore_2([1,2])
% colormore_2([1:5])
cmap=flip(colormore_2(1:11));%% 3:画图 %% 等值线分布情况
%% 一条等深线加上单位的情况
figure
set(gcf,'position',[50 50 1000 800],'color','w')
[c,h]=contour(hs,'color','k','linewidth',1.5);
clabel(c,h,[-0.5:0.1:0.3 0.5]);
hs_fix=unique(round(hs,1));
loc1=find(c(1,:)==0.4);
loc2=find(c(1,:)==0.5);
x=c(1,loc1+1:loc2-1);
y=c(2,loc1+1:loc2-1);
hold on
plot(x,y,'-r','linewidth',2)
hold on
plot(x(1:4),y(1:4),'-w')
hold on
plot(x(90:95),y(90:95),'-w')
h1=text(x(4),y(4),'0.4m','fontsize',12,'fontweight','bold','fontname','time news roman');
k=(y(5)-y(4))/(x(5)-x(4));
h2=text(x(90),y(90),'0.4m','fontsize',12,'fontweight','bold','fontname','time news roman');
k1=(y(91)-y(90))/(x(91)-x(90));
sita1=atan(k)*180/pi;
sita2=atan(k1)*180/pi;
set(h1,'rotation',sita1)
set(h2,'rotation',sita2)
hold on;
title('中国黄渤海区域2019年平均波高图');
export_fig('4.png','-r300')
%%
close all
figure
set(gcf,'position',[50 50 1000 800],'color','w')
[c,h]=contour(hs,[0.1:0.1:0.5],'linewidth',1.5);
hh=clabel_mb(c,'Units','m')% 
colormap(cmap)
export_fig('5.png','-r300')
%% 加上地图
%% 
figure
set(gcf,'position',[50 50 1000 800],'color','w')
m_proj('lambert','lat',[33 41],'long',[117 127]);
c=m_contour(x,y,hs,[0.1:0.1:0.5],'color','k');
hh=clabel_mb(c,'Units','m')% 
title('中国黄渤海区域平均波高图','fontsize',12,'fontname','宋体');
% m_gshhs_f('patch',[.8 .8 .8]);
m_coast('patch',[.8 .8 .8]);
m_grid('box','fancy');
colormap(cmap)
export_fig('6.png','-r300')