结合形态学冠层控制和分水岭的LiDAR数据单木提取方法技术

本发明专利技术的目的在于针对CHM、从LiDAR数据中提取单木时，减少对CHM质量的依赖，提高单木提取的精度。具体通过如下方法实现：基于形态学冠层控制，在CHM中确定外部标记；对所述CHM进行全局中值滤波处理，结合所述外部标记，使用局部极值的方式获得中值滤波处理后的CHM的内部标记；结合所述内部标记，对CHM进行两次分水岭操作；利用所述外部标记，分别对分水岭操作结果进行优化；根据单木分割结果，提取单木参数。本发明专利技术提供的方法，分割的结果只存在可忽略的过分割现象，使得单木分割的精度得到了提高，且对CHM质量的依赖程度得到降低，具有很强的实用价值。

【技术实现步骤摘要】

本专利技术涉及激光雷达数据处理技术，尤其涉及一种。
技术介绍
激光雷达(Light Detecting and Ranging, LiDAR)技术是近20年来遥感领域最具革命性的成就之一。根据其特点，LiDAR技术在林业上有着极为广泛的应用，尤其在构建冠层高度模型(Canopy Height Model, CHM)等方面有着突出的优势。目前从LiDAR数据中提取单木的算法有很多，其中一些算法是针对CHM进行的。 CHM是林业研究中的一个重要模型，它表示了地面以上的冠层高度，反映了植被的垂直高度变化以及植被上表面的水平分布情况。针对CHM、从LiDAR数据中提取单木是指，根据LiDAR 数据建立CHM，从CHM中提取单木的方法。目前，从CHM中提取单木的方法主要是区域生长法和标记分水岭法。但是，现有的从CHM中提取单木的方法对CHM质量要求较高，且无法保证单木提取的精度。
技术实现思路
本专利技术的目的是提供一种从，从而解决在提取单木时对CHM质量的依赖程度过高，以及提取单木精度无法保证的问题。本专利技术的目的通过以下技术方案来实现。—种，包括步骤1、基于LiDAR数据建立CHM ；步骤2、基于形态学冠层控制，在所述CHM中确定外部标记；步骤3、对所述CHM进行全局中值滤波处理；步骤4、结合所述外部标记，使用局部极值的方式获得全局中值滤波处理后的CHM 的内部标记；步骤5、结合所述内部标记，对步骤1中所述的CHM进行第一次分水岭操作；步骤6、利用所述外部标记，剔除第一次分水岭操作结果中的非冠层水槽，得到优化的第一次分水岭操作结果；步骤7、查找所述优化的第一次分水岭操作结果中没有内部标记的水槽，对所述没有内部标记的水槽进行重建，得到重建后的CHM ；步骤8、结合所述内部标记，对重建后的CHM进行第二次分水岭操作；步骤9、利用所述外部标记，剔除第二次分水岭操作结果中的非冠层水槽，得到优化的第二次分水岭操作结果；步骤10、使用预定的结构元素，对第二次分水岭操作结果进行开操作，得到单木分割结果；步骤11、根据单木分割结果，提取单木参数。本专利技术的关键在于始终将冠层控制的概念与单木分割过程结合，并且通过重建， 进行了两次的分水岭操作，从而大大的减少了过分割现象。分水岭分割之后，该专利技术还使用基于自适应结构元素形态学操作对结果进行优化，使得分割的单木冠层更贴近于真实的冠层形状。因此，本专利技术提供的方法，分割的结果只存在可忽略的过分割现象，使得单木分割的精度得到了提高，且降低了对CHM质量的依赖，具有很强的实用价值。附图说明图1本专利技术实施例提供的一种方法流程图；图2为本专利技术实施例提供的冠幅测量示意图。具体实施例方式本专利技术提供一种，其实现方式如图1所示，包括如下操作步骤1、基于LiDAR数据建立CHM ；步骤2、基于形态学冠层控制，在步骤1中建立的CHM中确定外部标记；所述的外部标记是一个二值图像，用来表示树冠的分布情况。所述的外部标记并不参与后续分水岭运算，以最大程度的保留分割的紧凑性，避免为控制过分割而带来的欠分割现象。外部标记的作用是针对分水岭运算的结果进行冠层控制，以剔除分水岭运算结果中的非冠层水槽。可以使用基于形态学的冠层控制进行外部标记的确定。作为举例而非限定，步骤 2的具体实现方式可以是步骤21、对上述的CHM进行形态学开操作，开操作的数学表示如下CHM0pen = CHMo S公式⑴其中，CHM0pen为经过开操作后的CHM，CHM为步骤1中建立的CHM，ο是形态学开操作算子，S为近圆形结构元素。步骤22、将开操作的结果与设定的冠层阈值进行比较，根据比较结果获得外部标记的初始二值图像；这里需要设定一个冠层阈值thresho 1 d0pen,开操作的结果CHMtjpen中大于等于 thresholds的像素被认定是有冠层区域(也就是树冠部分)，而小于该thresholfQpen的像素被认定为非冠层区域，即为地表。则步骤22的数学表示如下f 1，CHM0pen > threshold0penoutMark = \“pen .、公式⑵|0，CHM0pen < threshold0pen其中，outMark是外部标记的初始二值图像，outMark为1表示冠层区域，outMark 表示非冠层区域。步骤23、通过去除所述外部标记的初始二值图像中小于阈值的独立冠层区域和非冠层区域，得到外部标记二值图像，该外部标记二值图像即外部标记的表现形式。具体可以通过填补法和去除小像素块法实现。现有技术中，存在着通过形态学梯度来去除树冠的细节。但CHM上存在较为严重的冠层数据缺陷和丢失，这是本领域技术人员公知的。基于形态学梯度去除树冠细节只能单纯的反应冠层的边界，但无法实现对冠层数据的修复和还原。而本专利技术提供的方法，基于形态学的冠层控制来确定外部标记(外部标记反映了冠层数据)，然后通过填补法和去除小像素块法对外部标记进行优化，能够保留更多的冠层细节，最大程度的修复和还原冠层的真实形状。步骤3、对步骤1中建立的CHM进行全局中值滤波处理，全局中值滤波处理的数学表达如下CHMfcdian = Median (CHM)公式(3)其中，CHMfcdian是全局中值滤波后的CHM，MedianO是中值滤波函数。步骤4、结合上述的外部标记，使用局部极值的方式获得全局中值滤波处理后的 CHM的内部标记；步骤4的数学表示如下Γ1, CHMMedian(x,y) > Max(Sub_Window(x,y)),o^Ma^(x,y) = 1 _inMark = <, .,公式(4)其中，inMark为内部标记的二值图像、是内部标记的表现形式，MaxO为取最大值函数，Sub_ffindow()为取(X，y)像素的7*7邻域、不包括中心像素。内部标记表示可能的单木位置，相当于表示可能的树顶标记。应当指出的是，单木位置存在于树冠部分，因此，在计算单木位置时，关心的是外部标记取值为1的部分。现有技术中，通常假设树顶位于冠层的中间位置。但实际上，树顶往往并不处于冠层的中间位置。专利技术人在实现本专利技术的过程中，发现CHM中的单木树顶一定处于一些局部最大值的位置上，但也有可能会出现一个树冠有多个局部最大值的情况。为了确定树顶标记，本专利技术提供的方法是，首先对CHM进行中值滤波处理，来平滑掉一些同一树冠上的多个树顶的情况；然后结合外部标记和局部极值，查找可能的树顶标记。步骤5、结合上述的内部标记，对步骤1中建立的CHM进行第一次分水岭操作；在进行第一次分水岭操作之前，需要对步骤1中建立的CHM进行平滑操作，这样可以最大程度的降低分水岭的过分割现象。具体可以使用形态学的闭操作进行平滑处理，其数学表示如下CHMclose = CHM · S(5)其中，CHMclose是经过形态学闭操作后的CHM， 是形态学闭算子，S为近圆形结构元素。步骤5的具体实现方式可以是步骤51、对平滑处理后的CHM进行值反转操作；如果基于形态学的闭操作进行平滑处理，则步骤51的具体实现方式是，对CHMcl。se 进行值反转操作，得到CHMa_—…其数学表达如下CHMclosejieverse = Max (CHMclose) -CHMclose+a公式(6)其中，a为一个常数。步骤52本文档来自技高网...

【技术保护点】

【技术特征摘要】

【专利技术属性】
技术研发人员：赵旦，庞勇，李增元，刘清旺，徐光彩，
申请(专利权)人：中国林业科学研究院资源信息研究所，
类型：发明
国别省市：