发明名称 一种基于高斯尺度空间构建遥感影像语义的方法
摘要 本发明涉及一种基于高斯尺度空间构建遥感影像语义的方法,包括以下步骤:首先构建遥感影像的高斯尺度空间;然后分别建立包括同一尺度下的遥感影像不同空间区域和同一遥感影像空间区域下的不同尺度的二维高斯尺度树形层次结构;基于两个不同类型的高斯尺度树形层次结构分别构建非结构性语义和结构性语义,完成基于高斯尺度空间构建遥感影像语义。
申请公布号 CN102945550B 申请公布日期 2015.09.16
申请号 CN201210390177.8 申请日期 2012.10.15
申请人 武汉大学 发明人 邵振锋;周熙然
分类号 G06T7/00(2006.01)I;G06K9/46(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 严彦
主权项 一种基于高斯尺度空间构建遥感影像语义的方法,包括以下步骤:步骤1,基于遥感影像的内容,通过高斯函数建立遥感影像的高斯尺度空间;步骤2,提取遥感影像的稳定特征区域和尺度特征区域,所述提取遥感影像的稳定特征区域,包括在步骤1所得遥感影像的高斯尺度空间中,对每一个高斯尺度分别执行以下子步骤,步骤21.1,将同一高斯尺度下影像内容划分为不同的影像空间区域;步骤21.2,根据步骤21.1所划分同一高斯尺度下影像的影像空间区域,提取该高斯尺度下影像的稳定特征点;步骤21.3,基于步骤21.2所得的稳定特征点,提取该高斯尺度的稳定特征区域;所述提取遥感影像的尺度特征区域,包括在步骤1所得遥感影像的高斯尺度空间中执行以下子步骤,步骤22.1,对每个高斯尺度下影像内容划分不同的影像空间区域;步骤22.2,根据步骤22.1所划分各高斯尺度下影像的影像空间区域,基于各个不同尺度之间具有的尺度不变特征提取影像空间区域下的尺度特征点;步骤22.3,基于步骤22.2所得的尺度特征点,提取遥感影像的尺度特征区域;步骤3,根据步骤2所得稳定特征区域建立高斯尺度的树形层次结构,根据步骤2所得尺度特征区域建立高斯尺度的树形层次结构,得到一个二维高斯尺度空间的树形层次结构;步骤4,基于步骤3根据稳定特征区域建立的高斯尺度的树形层次结构,建立非结构性层次语义;基于步骤3中根据尺度特征区域建立的高斯尺度的树形层次结构,建立结构性层次语义;步骤21.1中,通过QaR树将同一高斯尺度下影像内容划分为不同的影像空间区域,基于QaR树的计算包括以下步骤,首先,计算基于Quad树的分解结果和基于R*‑树的分解结果的差值,如下所示:I=Rφ(f)‑Qφ(f)其中,I表示基于Quad树的分解结果和基于R*‑树的分解结果的影像差值空间区域,f表示待分解的某一高斯尺度的遥感影像,Rφ(.)表示基于R*‑树的分解,Qφ(.)表示基于Quad树的分解;然后,对式I=Rφ(f)‑Qφ(f)中所得的影像差值空间区域I中的每一个像素同基于Quad树分解所得影像的最大值和最小值进行差值计算,如果和最大值或者最小值的差值任意一个小于两者之间的Kullback‑Leibler距离,则将该像素归为影像粒的一部分,影像粒是基于遥感影像空间和遥感影像数据,所具有的最小且不可再分的影像单元,公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mo>|</mo><mi>max</mi><mrow><mo>(</mo><mi>Q&phi;</mi><mrow><mo>(</mo><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>|</mo><mo>&lt;</mo><mi>KL</mi><mo>&DoubleRightArrow;</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>|</mo><mi>min</mi><mrow><mo>(</mo><mi>Q&phi;</mi><mrow><mo>(</mo><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>|</mo><mo>&lt;</mo><mi>KL</mi><mo>&DoubleRightArrow;</mo><mi>&theta;</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mi>I</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>r</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000715923160000021.GIF" wi="1199" he="164" /></maths>其中,max(.)和min(.)分别代表最大值和最小值的计算,I(i,j,r)表示影像差值空间区域I,θ(i,j,r)表示归为影像粒的某个像素,其中i、j和r分别表示影像粒的某个像素在遥感影像f的x坐标、y坐标和光谱维中的波段信息,f(i,j,r)表示遥感影像f的坐标(i,j)处像素,KL表示Kullback‑Leibler距离,计算公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>KL</mi><mo>=</mo><mfrac><msup><mrow><mo>(</mo><mi>max</mi><mrow><mo>(</mo><mi>Q&phi;</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>min</mi><mrow><mo>(</mo><mi>Q&phi;</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mrow><mi>max</mi><mrow><mo>(</mo><mi>Q&phi;</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>min</mi><mrow><mo>(</mo><mi>Q&phi;</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000715923160000022.GIF" wi="729" he="157" /></maths>最后,将归为影像粒的结果同基于Quad树的分解结果进行合并,如下所示:IA=Qφ(f)+θ,0≤θ≤I其中,IA表示对遥感影像f基于R*‑树的分解结果和基于Quad树的分解结果融合后的影像,θ表示计算所得归为影像粒的某个像素θ(i,j,r)。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学