专利名称：一种基于空间相关性的亚像元空间定位方法
技术领域：
本发明涉及一种亚像元内部各组分的空间定位方法，属于地球空间信息技术领域。
背景技术：
利用数字遥感影像制作专题地图(影像分类)是遥感技术中的热点研究领域。然而， 由于地表复杂性，光谱传输过程中各种环境的影响及传感器探测元件本身的物理特性等 因素的作用，遥感影像中的像元很少是由单一地物组成的， 一般都包含两种或两种以上 地物。这种由多种地物组成的像元被称为混合像元，其光谱特征是所对应的不同土地覆 盖类型光谱响应的综合特征。混合像元在栅格编码影像中是普遍存在的，所占比重也很 大，尤其在不同类型地物的交界地带最为普遍，它的存在使同一类型或属性地物的遥感 影像灰度值具有一定的变化范围。解决混合像元的问题能够进一步细化影像分类结果， 使遥感技术应用进入到亚像元尺度，更加详细地表示地物细节信息，而不仅仅是停留在 像元尺度上，就成为一个必要而且有意义的研究内容。
"软分类(soft classification)"(相对于"硬分类(hard classification)"而 言)方法在一定程度上克服了混合像元的限制，例如光谱混合模型，监督模糊均值分 类，人工神经网络，支持向量机，基于知识和基于规则的推理技术等，能够提供了更为 理想的土地覆盖类型表达，但是并没有给出混合像元中每个地物类别的空间分布。像元 内部地物类别的空间分布信息能够提高一幅粗糙影像的空间分辨率，同时为不确定性的 研究提供更为丰富的信息。Atkinson在1997年首先提出亚像元成像(subpixel mapping) 的概念，即通过假设像元内部和像元之间的空间依赖性(spatial d印endence)，对亚像 元边界进行映射，确定每个像元内部出现的一定比例的地物类别的位置，以使空间依赖 性达到最大化。其后，许多学者相继对此提出了一些解决方法。例如，Aplin等在1999 年提出"per-field"方法在亚像元尺度上进行地域边界定位。Aplin和Atkinson也利 用矢量数据进行像元分割，然后获得亚像元定位。Schneider提出利用遥感影像边界特 征的先验知识实现亚像元的自动定位。Foody利用具有更高空间分辨率的锐化影像，釆 用基于简单回归的方法，对更粗分辨率影像的软分类输出结果进行锐化，得到亚像元地 物类型图。Atkinson在2001年又提出遥感影像中的超分辨率目标成像算法解决此问题。 Mertens et al.将遗传算法与空间依赖性假设进行结合，根据邻域值参数对像元内部所有可能的亚像元配置进行评价，从而为每一个亚像元分配位置。Tatem等将Hopfield神 经网形式化为能量最小化工具以检测地物目标，利用软分类的输出结果对Hopfield神 经网进行约束，对亚像元定位。Tatem在2003年利用真实的Landsat TM农作物影像对 所提出的方法进行了测试。但是，算法在PII-350计算机上的运行时间分别为约210分 钟和510分钟，是令人无法接受的。此外凌峰等还提出基于元胞自动机模型进行亚像元 定位。
上述方法在亚像元成像过程中，不仅依赖于原始影像，还应用了其它相关信息，例 如矢量数据、先验知识等，而且运算上需要多次迭代，速度较慢。鉴于此，本发明提出 一个混合像元内部地物空间定位的新方法，从而在亚像元的级别上提供了地物类型的边 界信息，与传统方法相比，此算法计算简单，解释性较强和算法效率高，更为重要的是， 除像元内部各端元组分面积比信息之外，此算法不需要任何其他辅助数据，例如矢量数 据、具有更高分辨率的影像、地物分布特征的先验知识等。
发明内容
本发明解决的技术问题提供一种基于空间相关性的亚像元空间定位方法，该方法 利用混合像元自身及其邻域像元的组分面积矢量，计算像元内部基本组分的亚像元空间 分布，进而模拟出较常规遥感分类方法更细的地物边界，能够提高遥感影像的分辨率， 更重要的是能够提高分类的精度。
本发明的技术方案 一种基于空间相关性的亚像元空间定位方法包括如下步骤 步骤1、预处理遥感数据，采集样本点，进行软分类获取被处理影像上每个像元的软分 类信息，即每个像元内部的各组分也称作端元的面积矢量或称丰度；
步骤2、设置亚像元放大系数N， N=4m, m = 1, 2, 3…，为自然数； 步骤3、根据设置的亚像元放大系数，利用混合像元自身及其邻域8个像元的组分面积 矢量，计算邻域像元及混合像元自身对混合像元内部各组分的贡献，即计算每个邻域像元 在其对应的边界线段上对中心像元各组分空间分布的贡献，这种贡献表现为在对应的线 段上的长度和位置，因此在这一步中首先计算邻域像元在对应线段上的长度，然后计算 在对应线段上的位置，进一步计算混合像元自身得组分面积矢量在混合像元内部对应的 长度和位置；
步骤4、根据邻域像元对中心混合像元的贡献和自身组分对中心混合像元的线段的位置， 确定中心混合像元内部各组分的空间分布，即得到原始影像中每个像元内部各组分的空间 分布。
所述步骤3中利用混合像元自身及其邻域像元的组分面积矢量，计算每个像元内部各基本组分的亚像元空间分布，假如混合像元内部有三类地物'A， B'邻'C，首先要计算邻域8个
像元对中心混合像元的贡献，即计算每个邻域像元在其对应的边界线段上对中心混合像
元A组分空间分布的贡献，这种贡献表现为在对应的线段上的长度和位置；其次，计算 中心像元内部线段位置。在确定了全部邻域像元和中心混合像元在其边界和中心线段上 的取值长度和位置之后，将由这些线段确定的边界点连接成多边形，得到中心混合像元 内组分A的模拟分布；同样的过程重复两次，就可以模拟出背景组分B和C在中心混合
像元内的空间分布。
本发明与现有技术相比的优点在于本发明利用空间相关性提出一个混合像元内部地 物空间定位方法，从而在亚像元的级别上提供了地物类型的边界信息，与传统方法相比，此 算法计算简单，解释性较强和算法效率高，更为重要的是，除像元内部各端元组分面积比信 息之外，此算法不需要任何其他辅助数据，例如矢量数据、具有更高分辨率的影像、地物分 布特征的先验知识等。
图1为本发明的主流程图2为中心像元及其邻域像元集合；
图3为中心像元尸。邻域像元尸,在像元^内部和边界上的对应关系。
具体实施例方式
如图1所示，本发明的具体实施步骤如下
步骤l，预处理遥感数据，采集样本点，进行软分类获取被处理影像上每个像元的软分 类信息，即每个像元内部的端元组分面积矢量或称丰度。
假定一个混合像元包含两种组分类型A、 B，当取一阶邻域时就可以得到图2中所
示的中心像元4和其邻域像元f的空间分布示意图。中心像元尸c.内部两种组分类型所占
面积比记为(ac,6^，各邻域像元《中两种组分类型面积比记为"，6,}。像元内部各组分 所占面积比满足下列条件"C+6C=1;《+6,=1, !:0,1,…，V-1 。把这些软分类信息
作为输入数据。
步骤2，设置亚像元放大系数N， N=4m, m = 1， 2, 3…，其目的在于将混合像元内部 细化的程度，如111=1， N-4，这表明每个混合像元被细分为4 x 4个亚像元，如果m^2， 则N-8，这表明每个混合像元被细分为8 x 8个亚像元。
步骤3，计算邻域像元及混合像元自身对混合像元内部各组分的贡献 (1)计算邻域像元对中心混合像元的贡献即计算每个邻域像元在其对应的边界线段上对中心像元各组分空间分布的贡献，定 义为F。这种贡献r表现为在对应的线段上的长度和位置。例如，在图3中，P。与线段
^5对应，其对中心像元的贡献p^的长度直接可以通过像元尸。内部组分计算。为后面计
算方便，再次将线段^58等分，等分点分别标以7;、 r2、石、r、 r4、 r5、 r6，如图3 中爿万线段上所示。这里组分a在A边界线段^s上的取值记为
其中，I l表示长度；[x]表示实数;c的圆整数(rounded rmmber)，即最接近于x的 整数。
当0<^<1时，需要进一步确定F^在线段^5上的取值位置。F^位置的确定过程 分为两步
(1) 初始位置rj^的确定从线段^4S的中点r开始，分别向两侧依次取等分点，直 至长度与已8相等。若^^表示的长度关于r不对称，则^s的位置更靠近两个相邻像元 中组分面积a,更大的一侧。图2中所示，^￡线段上，774标有0.25，也就是说当.a。 =0.25
时，f^)=|774|;同理，。。=0.50时，f^)=|r2r5|; 0。 = 0.75时，-|773 ;|等。
(ii)位置偏移量a^的确定初始位置f:确定之后，如果a^a，则需要将f^向 a、 ^中较大的像元一侧进行偏移，偏移量由下式确定-
<formula>formula see original document page 6</formula>
其中，省略了偏移量A^的单位。容易看出，A朋的单位为，，即线段乂S上的一
8
个8等分线段的长度。其余邻域像元对中心像元的贡献也可以通过类似的方法来确定。
(2) 计算中心混合像元内部线段位置
中心混合像元尸c对应十字交叉线段MCW、尸Og。 MCW、尸C^为过像元尸c的中心 点O,且互相垂直、长度相等的两条线段，长度都等于像元^边长的1/2，且O户-OQ， (9M-CW。将线段尸2、層八等分，等分点分别为^、 F2、 K、 (9、 K、 F5、 r6， lt3、 t/2、 t/,、 O、 f/4、 "、 f/6，如图3所示。像元尸c与线段MO/v、 POg对应，也就是说， 在水平方向上，组分a在像元4中所占面积&的值，以及与尸「相邻的两个邻域像元户6、A中A、 ^之间的大小关系，确定了组分A在线段MCW上的取值长度和位'置。同理， 在垂直方向上，组分A在像元A中所占面积+的值，以及与尸,、相邻的两个邻域像元户。、
i^中a。、 ^之间的大小关系，共同确定组分A在线段PC^上的取值长度和位置。组分A 在像元尸c中所占面积&的值确定了组分A的空间分布在线段MCW、 POg上的取值长 度。取值规则与前面所述的规则相同，只是对于中心像元^.自身来说，取值规则应该先 后应用在两条线段MCW、尸OG上，即在水平方向上，由^值的大小以及^、 ^之间的 大小关系共同确定在线段MOiV上的取值长度和位置；而在垂直方向上，则由&值的大 小以及a。、 "4之间的大小关系共同确定在线段户02上的取值长度和位置。
步骤4，根据邻域像元对中心混合像元的贡献和自身组分对中心混合像元的线段的位置， 确定中心混合像元尸。.内部各组分的空间分布，进而得到整个原始影像中每个像元内部各组 分的空间分布。
在确定了全部邻域像元S和中心像元《,在^的边界和中心线段上的取值长度和位 置之后，将由这些线段确定的边界点连接成多边形，就得到了像元g.内组分A的模拟分 布。同样的过程重复一次，可以由6,、 ^的值模拟出背景组分B在像元4内的空间分布。
依此类推，多类别组分也可以通过此方法获得各组分在亚像元的空间分布。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术 人员来说，本发听可以有各种更改和变化。凡在本发明的精神和原财之内，所作的任何 修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。
权利要求
1、一种基于空间相关性的亚像元空间定位方法，其特征在于步骤如下步骤1、预处理遥感数据，采集样本点，进行软分类获取被处理影像上每个像元的软分类信息，即每个像元内部的各组分面积矢量或称丰度；步骤2、设置亚像元放大系数N，N＝4m，m＝1，2，3…；步骤3、根据设置的亚像元放大系数，利用混合像元自身及其邻域8个像元的组分面积矢量，计算邻域像元及混合像元自身对混合像元内部各组分的贡献，即计算每个邻域像元在其对应的边界线段上对中心像元各组分空间分布的贡献，其计算方法为首先计算邻域像元在对应线段上的长度，然后计算在对应线段上的位置，再计算混合像元自身得组分面积矢量在混合像元内部对应的长度和位置；步骤4、根据邻域像元对中心混合像元的贡献和自身组分对中心混合像元的线段的位置，确定中心混合像元内部各组分的空间分布，即得到原始影像中每个像元内部各组分的空间分布。
2、 根据权利要求1所述的基于空间相关性的亚像元空间定位方法，其特征在于所述第二步中设置亚像元放大系数的m = 1， N=4;或m = 2 ， N=8。
3、 根据权利要求1所述的基于空间相关性的亚像元空间定位方法，其特征在于所述 步骤3中利用混合像元自身及其邻域像元的组分面积矢量，计算每个像元内部各基本组分的 亚像元空间分布假如混合像元内部有三类地物A， B和C，首先要计算邻域8个像元对中 心混合像元的贡献，即计算每个邻域像元在其对应的边界线段上对中心混合像元A组分 空间分布的贡献，这种贡献表现为在对应的线段上的长度和位置；其次，计算中心像元 内部线段位置。在确定了全部邻域像元和中心混合像元在其边界和中心线段上的取值长 度和位置之后，将由这些线段确定的边界点连接成多边形，得到中心混合像元内组分A 的模拟分布；同样的过程重复两次，就可以模拟出背景组分B和C在中心混合像元内的 空间分布。
全文摘要
本发明是一种基于空间相关性的亚像元空间定位方法。包括利用混合像元自身及其邻域像元的组分(端元)面积矢量或称丰度，计算像元内部基本组分的亚像元空间分布，进而模拟出较常规遥感分类方法更细的地物边界。本发明具有实用性广、模拟精度高等优点，适用于遥感数据进行目标提取、土地利用变化监测、环境污染评价、自然灾害分析等地学信息提取和地学数据挖掘工作。
文档编号G01S17/00GK101561882SQ20091008547
公开日2009年10月21日 申请日期2009年5月22日 优先权日2009年5月22日
发明者李三平, 咏 葛 申请人:中国科学院地理科学与资源研究所