【matlab程序】海图坐标轴单位的唯一

【matlab程序】海图坐标轴单位的唯一

本文写作来源，从实际出发，用于实际：

热带海洋学报，投稿须知：

其中一条关于海图制作规范中：

经度标识（E, W）置于最右端数字右侧。经纬度标识一个即可，不需要在每一个经纬度刻度后面均标识N或E。

我们平时画出的图片，每个刻度下方都有单位：

如图：

图片

按照要求改正为，每个刻度下方，只需要最后一个刻度后面出现标识N或者E；

如图：

图片

具体做法，我已做出，再算温度锋面过程中，顺带讲出几个实际问题和例子！

以下内容为，matlab代码和所需要的数据函数：包括colorbar函数，命名为：colormore_66.m;

地图所需的数据为WOA2023月平均数据2月份数据：woa23_decav91C0_t02_04.nc

地图填充颜色：color.mat

文件如下图：

图片

收费不贵，安心使用，文件可联系作者，付费后，获取，；

谢谢理解，坚持分享！感谢各位厚爱！

%% 本文目的写一个程序，寻找温度锋面：
% 温度锋面的定义或者说计算方法：Sobel算子
% 具体计算方法见：[1]王兴宇,纪棋严,彭腾腾等.东海黑潮温度锋的三维结构特征分析[J].海洋预报,2022,39(01):67-79.
% [1]王建辉. 北欧海海洋锋的时空变化研究[D].浙江海洋大学,2022.DOI:10.27747/d.cnki.gzjhy.2022.000413.
% 本文作者：海洋与大气科学
% Burger_miao
% 时间：2023年04月11日；
% 使用数据为WOA2023的年平均、季节平均和月平均；
clear;clc;close all;
% 01：读取温度数据
ncdisp ‘K:\03_data\WOA2023\月平均\woa23_decav91C0_t02_04.nc’;
file=‘K:\03_data\WOA2023\月平均\woa23_decav91C0_t02_04.nc’;%二月份温度文件
lon=double(ncread(file,‘lon’));
lat=double(ncread(file,‘lat’));
depth=double(ncread(file,‘depth’));
time=double(ncread(file,‘time’));
temp=double(ncread(file,‘t_an’));%‘sea_water_temperature’;
% 研究范围：为了方便对比本次计算结果的准确性，
% 研究范围和以往研究范围一致，120-131°E，22-33°E；
% 寻找范围内的温度数据
ln=find(lon>=119&lon<=132);
la=find(lat>=21.5&lat<=33.5);
lonx=lon(ln);laty=lat(la);
[x,y]=meshgrid(lonx,laty);x=x’;y=y’;
sst=squeeze(temp(ln,la,1,1));
%% 画个平面温度分布，看下数据分布情况；
close all
figure
m_proj(‘miller’,‘lon’,[120 131],‘lat’,[22 33])
m_contourf(x,y,sst,20,‘linestyle’,‘none’);
hold on
[c,h]=m_contour(x,y,sst,6:2:30,‘color’,‘k’,‘showtext’,‘on’);
clabel(c,h,‘fontsize’,10,‘fontname’,‘time news roman’,‘fontweight’,‘bold’);
col=load(‘color.mat’);
m_gshhs_f(‘patch’,col.color_list_temp,‘edgecolor’,‘k’);
m_grid;
map=colormore_66(100);
colormap(flip(map));
col_h=colorbar;
set(col_h,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’);
set(get(col_h,‘Title’),‘string’,‘°C’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’);
caxis([5 30])
export_fig(‘未改变的坐标轴.png’,‘-r300’)
%%
close all
figure
m_proj(‘miller’,‘lon’,[120 131],‘lat’,[22 33])
m_contourf(x,y,sst,20,‘linestyle’,‘none’);
hold on
[c,h]=m_contour(x,y,sst,6:2:30,‘color’,‘k’,‘showtext’,‘on’);
clabel(c,h,‘fontsize’,10,‘fontname’,‘time news roman’,‘fontweight’,‘bold’);
col=load(‘color.mat’);
m_gshhs_f(‘patch’,col.color_list_temp,‘edgecolor’,‘k’);
m_grid(‘xtick’,120:2:130,‘xticklabel’,[],‘ytick’,22:2:32,‘yticklabel’,[]…
,‘tickdir’,‘out’,‘Layer’,‘top’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘box’,‘on’,‘fontweight’,‘bold’);%(‘xtick’,120:5:130,‘xticklabel’,[{‘120’},{‘125’},{‘130’}]);
% xtick
m_text(119.5,21.5,‘120°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(121.5,21.5,‘122°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(123.5,21.5,‘124°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(125.5,21.5,‘126°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(127.5,21.5,‘128°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(129.5,21.5,‘130°E’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
% ytick
m_text(119,21.9,‘22°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(119,23.9,‘24°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(119,25.9,‘26°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(119,27.9,‘28°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(119,29.9,‘30°’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
m_text(119,31.9,‘32°N’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’)
%
map=colormore_66(100);
colormap(flip(map));
col_h=colorbar;
set(col_h,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’);
set(get(col_h,‘Title’),‘string’,‘°C’,‘FontName’,‘Times new roman’,‘FontSize’,10,‘fontweight’,‘bold’);
caxis([5 30])
export_fig(‘改变的坐标轴.png’,‘-r300’)