一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法技术

本发明专利技术提供一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，包括如下步骤：对于给定的高分辨遥感图像，首先设定多个不同的分割尺度，获得不同尺度下的初始分割结果；建立分割对象邻接图RAG；利用光谱角计算每个分割对象的局部光谱差异，包括分割对象内与分割对象间光谱差异，获得每个分割对象的光谱角差异值；生成同一位置下，不同分割尺度的对应光谱差异值曲线，选取光谱差异值最高的尺度下对应分割对象作为该位置最优分割对象；综合所有位置下选取的分割对象，组合成为最终的分割优化结果。本发明专利技术直接从多个分割尺度下的初始分割结果选取最优分割对象，自适应地调节了遥感图像不同地物下的最优分割尺度，获得了理想的分割效果。

An Optimal Method for Multi-scale Segmentation of High Resolution Remote Sensing Images Based on Maximizing Spectral Differences

【技术实现步骤摘要】
一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法
本专利技术涉及到遥感图像的分割领域，尤其涉及一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法。
技术介绍
图像分割(章毓晋.图像分割.北京：科学出版社,2005.)是图像处理与计算机视觉领域中最为基础和重要的问题之一，也是对图像进行视觉分析和模式识别的基本前提。高分辨遥感图像以其丰富的光谱信息，包含着大量的地物信息特征。多尺度分割算法(＝multiresolutionsegmentation、均值漂移算法)已被广泛应用到高分辨遥感图像分割中。多尺度分割算法中常常需要预先设定一个尺度参数，用来控制分割对象的平均大小，尺度参数的选取对分割效果有着重要影响。因此，尺度参数的选取成为多尺度分割优化的一个重要问题。目前，大部分多尺度分割优化方法集中于从多个尺度的分割结果中选取一个单独的尺度作为最优化的参数。某一尺度下的分割结果的好坏通常由分割对象间的差异性来度量，常用欧式距离、光谱角等评价指标来度量这种差异性，最后通过这个评价指标来选取一个全局最优的尺度参数。然而，这种方法往往难以兼顾高分辨遥感图像中不同地物对分割尺度的需求，极大限制了分割效果。例如，大面积的海域、连接的养殖池塘、规则的农田等不同的地物环境，便不适用同一个分割尺度。此外，全局最优尺度的选取也忽略了其余尺度参数下的分割信息，没有充分利用。传统的遥感图像分割优化算法仅在多个分割尺度的结果中选取一个单一尺度作为最优分割结果，没有充分利用到多个分割结果中的尺度信息，依然不能保证对遥感图像中不同地物提供适应的分割尺度。
技术实现思路
针对现有技术的不足，本专利技术提出一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，该方法针对多尺度分割算法只能设置一个尺度参数的不足，自适应在不同尺度参数下的分割结果中为相应地物选取最优分割结果，可以更准确地描绘遥感图像中不同环境下的地物边缘。为达到以上目的，本专利技术具体技术方案如下：一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，其特征在于，该方法包括如下步骤：(1)设置分割尺度参数范围[1,2,…,L]，使用多尺度分割算法对待分割的高分辨遥感影像进行不同尺度的分割，获得L个不同尺度下的分割结果，作为待优化的初始分割结果；(2)对每一个尺度下的初始分割结果建立分割对象邻接图RAG，用以表示相邻分割对象的空间关系，所述的分割对象为一幅遥感图像分割后划分的子区域，这些子区域组成遥感图像的分割结果；(3)选择任意一个坐标点，该坐标点所在的分割对象为O，分割尺度为t，t∈[1,2,…,L]，进行如下步骤：(3.1)计算分割对象的局部光谱差异，包括分割对象间的光谱差异和分割对象内光谱差异，其中，分割对象间的光谱差异角计算如下：计算分割对象O的光谱均值向量s，并将该分割对象O的相邻分割对象的光谱均值向量组成集合计算分割对象间光谱差异角计算公式如下：其中，s∈RB，B是该遥感图像的波段数，每一个xi都是O的相邻分割对象的光谱均值向量，xi∈RB，u是相邻分割对象的个数；分割对象内光谱差异角计算如下：在小于尺度t的分割结果中搜索组成当前分割对象O的分割对象，组成分割对象O的每个分割对象的光谱均值向量的集合为Ω，根据下式计算分割对象内光谱差异角其中，p与q是小于当前尺度t下组成O的任意两个分割对象的光谱均值向量，r是集合Ω内的对象数目；若在小于t的尺度下没有找到组成分割对象O的分割对象，则(3.2)采用下式计算分割对象O的局部光谱差异值(3.3)计算L个分割结果在当前坐标点的局部光谱差异值EO，形成光谱差异值的曲线，选择光谱差异值最大的尺度对应分割对象作为该坐标下的最优分割对象；(4)对遥感图像上所有坐标位置重复步骤(3)，获得每个坐标下的最优分割对象；(5)将最优分割对象组合为最终的优化后分割结果。进一步地，所述光谱均值向量s计算公式如下：其中，Z为分割对象O内所有像素点的集合，zi∈RB，m是像素的个数。进一步地，所述的多尺度分割算法具体为multiresolutionsegmentation。进一步地，所述的多尺度分割算法具体为均值漂移算法。与现有技术相比，本专利技术的有益效果如下：本专利技术中，通过对比同一坐标下不同尺度的分割对象的局部光谱值，来选取该坐标下的最优分割对象，可以自适应对不同地物设置相应的分割尺度。相比与传统尺度优化方法，本专利技术的优化结果能综合利用不同尺度下的分割信息，提高分割的准确率。附图说明图1是本专利技术的基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法的流程框图；图2是本专利技术中分割对象光谱间和分割对象内差异值的计算示意图；图3是本专利技术的其中一个实施例中的最优分割对象的选取示意图，其中，图3(a)-(e)分别是尺度50、100、150、350、400下该养殖鱼塘的分割结果图，图3(f)是目标养殖鱼塘的参考分割结果图，图3(g)为各个尺度下该位置的局部光谱值E的曲线图；图4是本专利技术的其中一个实施例的初始对尺度分割结果与最终分割优化后结果对比图，图4(a)是该遥感图像，图4(b)-(g)分别是尺度70、140、200、250、300、400下的分割结果图，图4(h)是分割优化后的结果图。具体实施方式下面根据附图和优选实施例详细描述本专利技术，本专利技术的目的和效果将变得更加明白，以下结合附图和实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，其特征在于，该方法包括如下步骤：(1)设置分割尺度参数范围[1,2,…,L]，使用多尺度分割算法对待分割的高分辨遥感影像进行不同尺度的分割，获得L个不同尺度下的分割结果，作为待优化的初始分割结果；(2)对每一个尺度下的初始分割结果建立分割对象邻接图(region-adjacentgraph)RAG，用以表示相邻分割对象的空间关系；所述的分割对象为一幅遥感图像分割后划分的子区域，这些子区域组成遥感图像的分割结果；(3)选择任意一个坐标点，该坐标点所在的分割对象为O，分割尺度为t，t∈[1,2,…,L]，进行如下步骤：(3.1)计算分割对象的局部光谱差异，包括分割对象间的光谱差异和分割对象内光谱差异，其中，分割对象间的光谱差异角计算如下：计算分割对象O的光谱均值向量s，并将该分割对象O的相邻分割对象的光谱均值组成集合计算分割对象间光谱差异角计算公式如下：其中，s∈RB，B是该遥感图像的波段数，每一个xi都是O的相邻分割对象的光谱均值向量，xi∈RB，u是相邻分割对象的个数；如果O与周围相邻对象xi的光谱角越大，表示O与周围是同一种地物的可能性很小；反之，O越可能与周围对象属于同一种地物类别。分割对象内光谱差异角计算如下：在小于尺度t的分割结果中搜索组成当前分割对象O的分割对象，组成分割对象O的每个分割对象的光谱均值向量的集合为Ω，根据下式计算分割对象内光谱差异角其中，p与q是小于当前尺度t下组成O的任意两个分割对象的光谱均值向量，r是集合Ω内的对象数目；若在小于t的尺度下没有找到组成分割对象O的分割对象，则类似的，公式(2)考虑了分割对象O内部光谱的差异性，如果这种内部的差异性越大，就表示本文档来自技高网...

【技术保护点】
1.一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，其特征在于，该方法包括如下步骤：(1)设置分割尺度参数范围[1，2，...，L]，使用多尺度分割算法对待分割的高分辨遥感影像进行不同尺度的分割，获得L个不同尺度下的分割结果，作为待优化的初始分割结果；(2)对每一个尺度下的初始分割结果建立分割对象邻接图RAG，用以表示相邻分割对象的空间关系，所述的分割对象为一幅遥感图像分割后划分的子区域，这些子区域组成遥感图像的分割结果。(3)选择任意一个坐标点，该坐标点所在的分割对象为O，分割尺度为t，t∈[1，2，...，L]，进行如下步骤：(3.1)计算分割对象的局部光谱差异，包括分割对象间的光谱差异和分割对象内光谱差异，其中，分割对象间的光谱差异角计算如下：计算分割对象O的光谱均值向量s，并将该分割对象O的相邻分割对象的光谱均值向量组成集合

【技术特征摘要】
1.一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，其特征在于，该方法包括如下步骤：(1)设置分割尺度参数范围[1，2，...，L]，使用多尺度分割算法对待分割的高分辨遥感影像进行不同尺度的分割，获得L个不同尺度下的分割结果，作为待优化的初始分割结果；(2)对每一个尺度下的初始分割结果建立分割对象邻接图RAG，用以表示相邻分割对象的空间关系，所述的分割对象为一幅遥感图像分割后划分的子区域，这些子区域组成遥感图像的分割结果。(3)选择任意一个坐标点，该坐标点所在的分割对象为O，分割尺度为t，t∈[1，2，...，L]，进行如下步骤：(3.1)计算分割对象的局部光谱差异，包括分割对象间的光谱差异和分割对象内光谱差异，其中，分割对象间的光谱差异角计算如下：计算分割对象O的光谱均值向量s，并将该分割对象O的相邻分割对象的光谱均值向量组成集合计算分割对象间光谱差异角计算公式如下：其中，s∈RB，B是该遥感图像的波段数，每一个xi都是O的相邻分割对象的光谱均值向量，xi∈RB，u是相邻分割对象的个数。分割对象内光谱差异角计算如下：在小于尺度t的分割结果中搜索组成当前分割对象O的分割对象，组成分割对象O的每个分...

【专利技术属性】
技术研发人员：陈建裕，沈煜，毛志华，陈宁华，黄清波，
申请(专利权)人：自然资源部第二海洋研究所，
类型：发明
国别省市：浙江,33