今天，在我的系列在计算美国大陆的鸿沟时，我们将看一种叫做最小值实施．之前我们看了数据导入和显示，分水岭变换和标签矩阵，不同种类的最小值，以及操纵二值图像形成海洋掩膜。

在第三部分中，我们发现我们的美国大陆DEM数据集包含了超过50万个区域最小值，每个区域最小值都通过分水岭转换变成了一个集水盆地。我们只需要找到两个集水区。

最小值实施是一种修改表面高度的技术，将一些高度值向上推，使结果表面只在指定的位置有区域最小值。Minima拼版执行这种修改，同时保持相邻像素之间的局部高度排序关系。

这里有一个一维的例子来说明。

A = [1 2 3 4 3 4 2 1]

A = 1 2 3 4 3 3 2 1

向量一个有三个不同的区域最小位置:

imregionalmin (a)

Ans = 1 0 0 0 1 0 0 1

让我们做一个蒙版，表示两个我们想要区域最小值的位置:

掩码= [1 0 0 0 0 0 1]

掩码= 1 0 0 0 0 0 1

这个模板说的是我们只希望在向量的每一端都有区域极小值，所以我们希望中间的区域极小值以某种方式消失。我们使用图像处理工具箱函数来实现这一点imimposemin．

imimposemin(面具)

ans = -Inf 2.0030 3.0030 4.0030 4.0030 4.0030 3.0030 2.0030 -Inf

新的结果只在期望的位置有区域最小值。向量中间的局部最小值被推到包围元素的级别。(当我写这篇文章的时候，我注意到这个小的突起imimposemin给予所有中间元素似乎是不必要的。这是我将在另一个时间调查的事情。)

在我之前的文章中，我创建了一个代表海洋面具的二进制图像:

s =负载(“continental_divide”）;民主党= s.dem_cropped;ocean_mask = s.ocean_mask;imshow (ocean_mask“InitialMagnification”，“健康”）

我可以使用imimposemin，连同ocean_mask图像，以修改DEM，使两个海洋是只有区域最小值。

Dem_modified = imimposemin(dem, ocean_mask);

它看起来没什么不同:

imshow (dem_modified (-500 3000),“InitialMagnification”，“健康”）

但我们可以证实现在只有两个区域极小值:

regional_minima_labeled = bwlabel (imregionalmin (dem_modified));max (regional_minima_labeled (:))

ans = 2

imshow (label2rgb (regional_minima_labeled),“InitialMagnification”，“健康”）

我将把修改后的DEM添加到我在continental_divide.mat中收集的数据中:

保存continental_dividedem_modified附加

下一次，我将执行计算大陆分水岭的最后一步，并展示一种使用透明度将其可视化的方法。

关于这个系列

本系列探讨了计算美国大陆分水岭位置的问题。这道分水岭将大西洋和太平洋流域划分为北美大陆。

作为一个算法开发问题，计算分水岭让我们可以探索数据导入和可视化、操作二值图像掩模、标签矩阵、区域最小值和分水岭变换等方面。

数据来源:全球任务团队和其他人(Hastings, David A, Paula K. Dunbar, Gerald M. Elphingstone, Mark Bootz, Hiroshi Murakami, Hiroshi Maruyama, Hiroshi Masaharu, Peter Holland, John Payne, Nevin A. Bryant, Thomas L. Logan, j.p。穆勒，冈特·施莱尔，约翰·s·麦克唐纳)编。, 1999年。全球陆地一公里基高程(GLOBE)数字高程模型，版本1.0。美国国家地球物理数据中心，科罗拉多州博尔德市百老汇325号，80305-3328http://www.ngdc.noaa.gov/mgg/topo/globe.html)和cd - rom。

获取MATLAB代码

MATLAB®7.8发布