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图片
智能推荐
qpushbutton添加menu时不显示三角下表_如何看待奶粉中添加的香精香料？不同人群这样说...-程序员宅基地
文章浏览阅读115次。婴幼儿配方奶粉是作为宝宝在无法母乳喂养时的无奈选择，也是除母乳以外最好的选择，宝妈在选择奶粉时也会格外认真，当看到奶粉中含有香精香料成分时，部分宝妈便不会购买，甚至认为此类成分是对宝宝身体有所危害的，所以坚决抵制购买，那么添加到奶粉中的香精香料到底怎么样，长期吃有没有坏处？我们来看看吧~不同人群，对香精香料有不同看法"宝妈们这样说：前段时间网上一位宝妈的评论吸引了笔者的注意：“我这奶粉买回家一开罐...
【备忘录】 iPad air3 pro 10.5 原装键盘 失灵 维修_ipad10.5键盘套没反应-程序员宅基地
文章浏览阅读2.5k次。并不是苹果说的触点问题，而是软布导线被腐蚀断线问题。用铜箔胶带粘一下就好了。不过二手直接买一个咸鱼上也就三百多，也不贵。这个明显是苹果设计缺陷。直接丢了挺可惜的。拆的时候小心点。我的QQ: 13008312096，有空的话可以帮忙有偿代处理。参考链接如下：nullhttps://www.youtube.com/watch?v=buNYHzMZJdk修复好的样子，毫无违和感：..._ipad10.5键盘套没反应
gif图用photoshop快速抽帧,压缩,去底,修改时间_gif抽帧-程序员宅基地
文章浏览阅读1.3w次，点赞3次，收藏13次。网上下载的gif图很多都不透底,而且尺寸时间什么的都不合适,用ps简单修改一下就好了1.抽帧很多gif图帧数很多,抽帧可以大幅压缩gif大小1.导入gif到ps2.勾选动作和时间轴面板这时我们发现帧和图层的隐藏显示是相对应的这样的话我们操作图层就会打乱帧,我们要让所有的图层都显示,但又不干扰帧的显示3.这时我们就要把它转换为视频时间轴4.然后再转换帧>>转换为..._gif抽帧
c++ 优先队列_面试必知必会|堆和优先队列-程序员宅基地
文章浏览阅读235次。通过本文将了解到以下内容：优先队列的概念优先队列的实现优先队列的应用1.优先队列的概念优先队列是计算机科学中的一类抽象数据类型。优先队列中的每个元素都有各自的优先级，优先级最高的元素最先得到服务；优先级相同的元素按照其在优先队列中的顺序得到服务。优先队列至少需要支持下述操作：a.插入带优先级的元素b.取出具有最高优先级的元素c.查看最高优先级的元素。综合考虑插入和删除的性能 优先队列一般采用堆来实..._优先队列 c++面试
windows系统VMware安装Linux虚拟机、Linux配置JDK环境、Linux安装tomcat、Linux安装mysql、Linux通过SSH连接Navicat 数据库可视化工具_vm中linux下载安装mysql jdk-程序员宅基地
文章浏览阅读281次。windows系统VMware安装Linux虚拟机、配置JDK环境、安装tomcat、安装mysql、安装redis第一部分VMware安装Linux—————————————————————————————————xshell、VMware工具链接：解压之后安装VMware工具xshell不需要安装解压可直接使用..._vm中linux下载安装mysql jdk
Python 加密解密_咕咕加密v4-程序员宅基地
文章浏览阅读1k次。# -*- coding: utf-8 -*-import hashlib'''加密解密'''#MD5是最常见的摘要算法，速度很快，生成结果是固定的128 bit字节，# 通常用一个32位的16进制字符串表示。md5 = hashlib.md5()update = md5.update('hhhhhaaa')print(md5.hexdigest())#SHA1的结果是160..._咕咕加密v4
随便推点
Ubuntu18.04安装教程——超详细的图文教程_ubuntu系统18.04-程序员宅基地
文章浏览阅读10w+次，点赞112次，收藏974次。Ubuntu18.04镜像_ubuntu系统18.04
二、Json对象、Json数组和Json字符串_json字符串数组-程序员宅基地
文章浏览阅读6.7k次，点赞3次，收藏17次。一、Json字符串和Json对象定义：1、Json字符串：所谓字符串：单引号或者双引号引起来，是一个String类型的字符串：如下：var person='{"name":"shily","sex":"女","age":"23"}';//json字符串console.log(person)console.log(person.name)console.log(typeof person) 2、Json对象：最显著的特征：对象的值可以用 “对象.属性” 进行访问，_json字符串数组
Linux系统100条命令：关于Ubuntu和 CentOS 7 相同功能的不同的终端操作命令_ubuntu 命令跟centos-程序员宅基地
文章浏览阅读718次。CentOS 7：ip link set interface_name up 或 ip link set interface_name down。Ubuntu：ifconfig interface_name up 或 ifconfig interface_name down。CentOS 7：编辑 /etc/sysconfig/network-scripts/ifcfg-eth0 文件。Ubuntu：编辑 /etc/network/interfaces 文件。_ubuntu 命令跟centos
windows10下VS2019编译jpegsrc.v9e.tar.gz为lib静态库（已验证）_jpeg library error vs2019-程序员宅基地
文章浏览阅读652次。jpegsr9e windows vs2019生成方法，以及库下载_jpeg library error vs2019
重磅？华为 Mate60/Pro 系列网速实测结果公布，最高 1205.57 Mbps_华为mate60pro+核实网络-程序员宅基地
文章浏览阅读647次。总的来说，华为Mate 60/Pro系列手机的高速网速表现引起了广泛的关注，这也是消费者对该系列手机购买热情高涨的一个重要因素。可以看出，华为Mate 60/Pro系列手机的网速表现非常出色，这也是消费者购买该系列手机的一个重要原因。此前，华为Mate 60 Pro的供应量已经增至1500万至1700万台，而最新消息称，华为Mate 60 Pro和Mate 60 Pro+的出货量甚至已上调至2000万台。目前，在中国市场上，手机竞争愈发激烈，不仅华为Mate 60系列，其他品牌的手机也都受到了高温的迎接。_华为mate60pro+核实网络
access查找出生日期年份_access怎样利用出生日期计算年龄呀！-程序员宅基地
文章浏览阅读7.1k次。公告： 为响应国家净网行动，部分内容已经删除，感谢读者理解。话题：access怎样利用出生日期计算年龄呀！回答：lt;%set rs = server.createobject(＂adodb.recordset＂) curid=request(＂id＂) sql = ＂UPDATE pany SET a_num=a_num+1,day_count=day_count+1 WHERE day_lda..._access出生年份表达式