m_map绘制多波束数据-程序员宅基地

之前，分享了m_map如何绘制本地高程数据：m_map导入本地地形数据。

其实，m_map还可以绘制多波束测深数据。

1. 多波束网站

多波束测深数据可以从https://www.ncei.noaa.gov/maps/grid-extract/或者https://www.ncei.noaa.gov/maps/bathymetry/下载。
多波束测深数据的分辨率为3秒，约为90 m 。

在这里插入图片描述

2. 绘图

以纽约以东的一片海域为例。

这是原始的多波束数据，里面无有效数据处，被填充为0值，显示如下：

在这里插入图片描述

将0值替换为NaN值，即显示为空白，显示如下：

在这里插入图片描述

将m_map自带的1分的分辨率地形图，融合到多波束数据中，显示如下：

在这里插入图片描述

显示带有海山的区域：

注：绘图色标参考自：https://blog.csdn.net/slandarer/article/details/127719784

3. 代码

mat_name='newyork';
fignum=1

load([mat_name,'_elevation.mat'])
elev=flipud(elev);

% 将0值替换为NaN
% elev(find(elev==0))=nan;

lonlim=[extent(1:2)]
latlim=[extent(3:4)]
[nlat,nlon]=size(elev);
Lon=[linspace(lonlim(1),lonlim(2),nlon)];
Lat=[linspace(latlim(1),latlim(2),nlat)]';


%% 高程数据融合
% 在多波束数据无数据（标志为0处，替换为自带的1分高程数据）
% 首先提取1分的高程数据

[ELEV,LON,LAT]=m_etopo2([lonlim(1)-1 lonlim(2)+1 latlim(1)-1,latlim(2)+1]);
[Lon,Lat]=meshgrid(Lon,Lat);
% 将1分数据插值为和多波束数据同样的大小
ELEV_p=interp2(LON,LAT,ELEV,Lon,Lat);

for i=1:nlat
    for j=1:nlon
        if elev(i,j)==0
            elev(i,j)=ELEV_p(i,j);
        end
    end
end

%% 绘图

m_proj('mercator','long',[-67 lonlim(2) ],'lat',[38.5,latlim(2)]);

caxis([-5000 0])
colormap(slanCM('rainbow')) 
hc=colorbar;
set(get(hc,'title'),'string','Elevation(m)')
set(hc,'tickdir','out')

m_shadedrelief(Lon(1,:),Lat(:,1),elev,'lightangle',45,'gradient',10) 

m_gshhs('ic','color','k')
m_grid('box','on','tickdir','out','gridlines','no')

set(gcf,'position',[10 10 400 400])

figname=[mat_name,num2str(fignum)];
print('-dpng','-r500',[figname,'.png'])  % 导出png图片

本文链接：https://blog.csdn.net/weixin_39509073/article/details/134909517

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

readelf_windows readelf-程序员宅基地

文章浏览阅读120次。系统里的目标文件是按照特定的目标文件格式来组织的，各个系统的目标文件格式都不相同。从贝尔实验室诞生的第一个Unix系统使用的是a.out格式（直到今天，可执行文件仍然称为a.out文件）。Windows使用可移植可执行（PortableExecutable，PE)格式。Mac OS-X使用Mach-O格式。现代x86-64Linux和Unix系统使用可执行可链接格式（Executable and Linkable Format，ELF)。_windows readelf

详解冬奥冠军背后的AI黑科技-程序员宅基地

文章浏览阅读3.6k次。用人工智能普惠体育发展。

form表单提交的几种方式_提交表单-程序员宅基地

文章浏览阅读10w+次，点赞92次，收藏495次。表单提交方式一：直接利用form表单提交html页面代码：<!DOCTYPE html><html><head><meta charset="UTF-8" /><title>Insert title here</title></head><body><form action="h..._提交表单

Unity Spine SkeletonGraphic 动画重复播放过度残影透明渐变Bug 解决方案_unity skeletongraphic-程序员宅基地

文章浏览阅读5.1k次。Unity Spine SkeletonGraphic 重复播放过度残影Bug 解决方案不推荐使用SetToSetupPose和Setup Pose相关，代码直接贴上/// <summary>/// Spine播放设置/// </summary>/// <param name="trackIndex">填写0</param>/// <param name="animationName">动画名</param>/// &l_unity skeletongraphic

高斯分布3——边缘概率与条件概率_高斯分布的条件概率-程序员宅基地

文章浏览阅读3.5k次。一、推导过程:二、结果：边缘分布x1,x2 各自依然服从 μi，写反差矩阵 Σii 的多元高斯分布；条件概率分布给定 xj 求 xi 的分布：μi|j=μi+ΣijΣ−1jj(xj−μj)Σi|j=Σjj−ΣTijΣ−1iiΣij..._高斯分布的条件概率

Ratelimitcache: Python缓存库，支持速率限制-程序员宅基地

文章浏览阅读339次，点赞8次，收藏8次。Ratelimitcache: Python缓存库，支持速率限制项目链接: https://gitcode.com/simonw/ratelimitcache?utm_source=artical_gitcode如果你正在寻找一个Python缓存库，并且希望对缓存操作进行速率限制，那么Ratelimitcache可能是你的理想选择。什么是Ratelimitcache？Ratelimitca..._python ratelimit基于什么

随便推点

基于OFDM+64QAM系统的载波同步matlab仿真,输出误码率,星座图,鉴相器,锁相环频率响应以及NCO等-程序员宅基地

文章浏览阅读454次。正交频分复用（OFDM）是一种在现代通信系统中广泛使用的调制技术，它具有高效的频谱利用和抗多径衰落等特点。64QAM（64-ary Quadrature Amplitude Modulation）是一种调制方式，可以在每个符号中传输更多的位信息。在OFDM系统中，保持载波同步对确保数据传输的可靠性至关重要。_基于ofdm+64qam系统的载波同步matlab仿真,

Springboot毕设项目超市商品销售管理系统37x2w（java+VUE+Mybatis+Maven+Mysql）_vue+springboot+mybatis商品管理系统-程序员宅基地

文章浏览阅读67次。Jdk1.8 + Tomcat8.5 + Mysql + HBuilderX（Webstorm也行）+ Eclispe（IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持）。若包含，则为maven项目，否则为非maven项目。Springboot毕设项目超市商品销售管理系统37x2w（java+VUE+Mybatis+Maven+Mysql）Springboot + mybatis + Maven + Vue 等等组成，B/S模式 + Maven管理等等。其他版本理论上也可以。_vue+springboot+mybatis商品管理系统