topo_interp文档
的topo_interp函数篡改海拔相对海平面在任何地理点。的数据ETOPO5世界数字高程模型,这是在5分钟(或1/12度)提供网格分辨率。
内容
语法
子= topo_interp(意大利,loni)子= topo_interp(拉提、loni‘方法’,InterpMethod)
描述
loni子= topo_interp(意大利)返回海拔相对海平面的地理点意大利,loni。
子= topo_interp(拉提、loni‘方法’,InterpMethod)指定一个插值法,可以接受任何方法interp2。默认的方法“线性”。
例1:一个网格
这是一个全球网格分辨率0.25度:
(纬度、经度)= cdtgrid (0.25);
海拔高度在每个网格点了是这样的:
Z = topo_interp(纬度、经度);
这样的情节高程网格。下面我用cmocean地形colormap:
pcolor(经度,纬度,Z)阴影插值函数cb = colorbar;ylabel (cb、的高度(米))cmocean (“威尼斯平底渔船”,“主”)%设置colormap
示例2:海拔概要文件
假设你想从巴黎旅行(48.8567 n, 2.3508 e)到曼谷(13.7525 n, 100.494167 e)。这里有一个粗线组成的1000点这两个城市之间:
lat = linspace (48.8567, 13.7525, 1000);朗= linspace (2.3508, 100.494167, 1000);持有在情节(经度、纬度、的r -,“线宽”,2)
这是海拔剖面沿线,容易:
z = topo_interp(纬度、经度);图绘制(经度,z)轴紧盒子从包含“经”ylabel的高度(米)
或者描述土地和水更直观,使用异常功能和使用指定的颜色rgb:
图异常(经度,z,“topcolor”rgb (“土”),…“bottomcolor”rgb (“蓝色海洋”)轴紧盒子从包含“经”ylabel的高度(米)
示例3:加载原始数据
访问原始数据,只是类型
负载(“global_topography.mat”);
它包含的变量纬度,朗,Z。把它们像这样:
图显示亮度图像(经度,纬度,Z)轴xy图像包含“经”ylabel“纬度”cb = colorbar;ylabel (cb、的高度(米))cmocean (“威尼斯平底渔船”,“主”)
示例4:海平面上升的影响
让我们来快速浏览的地方可能会受到海平面上升的影响。因为我们已经加载全分辨率数据集在例3中,我们将只使用数据集,除了我们将开始通过设置所有海洋值NaN。所以请所有小于或等于零的海拔南:
Z (Z < = 0) =南;
所有格陵兰和南极的冰有潜力使全球海平面上升约67。当然这种情况不太可能很快,但让我们来看看这样一个事件的影响。我们将不得不做出一个原油(错误的)假设海平面会均匀分布在今天的海洋,但这是一个开始。
调查的影响增加约67米的海平面全球海洋,情节蒙面地形,并设置颜色轴范围从0到67:
图imagescn(经度,纬度,Z) caxis (67 [0]) cb = colorbar;ylabel (cb、“海平面上升(m)”)cmoceanamp%设置colormap
一阶近似,红色是危险的地方将水下的冰原的全面崩溃。深红色意味着更加脆弱。让我们放大在东部沿海地区。使用边界函数绘制状态边界。
轴([-100 -70 24 41])边界(“州”)
注意在混叠
ETOPO5数据集提供分辨率5分钟(或1/12的程度),但是topo_interp函数不插补前做任何反走样。以满足尼奎斯特要求以防止混叠,你应该在技术上interploate网格点至少两次间隔5分钟解决底层的数据集。这将是一个非常密集的网格为全球数据集,所以您可能希望接受一些潜在的混叠。或者你可以加载原始数据,使用imresize减少网格(在执行反锯齿)和插入。这完全取决于你。
作者信息
这个函数的一部分气候数据为Matlab工具箱。的功能和支持文档是乍得a·格林写的金宝app德州大学奥斯丁分校。