专利名称：一种结合可见光图像信息的sar图像快速分割方法
技术领域：
本发明属于图像处理领域，更具体地说涉及合成孔径雷达(Synthetic ApertureRadar，SAR)图像的分割方法。
背景技术：
SAR是一种高分辨微波成像雷达，具有全天时、全天侯工作的优点，并且对植被具 有一定的穿透能力。因此，从上个世纪五十年代SAR被提出后，其发展非常迅速，到2008年 底，民用和军用的SAR卫星达18颗，而机载SAR系统更是不计其数。随着这些SAR系统的 成功运行，SAR图像的应用和理解亟待发展。因此，SAR成像处理技术成为当前SAR领域的 研究热点之一。SAR图像分割是图像处理的重要内容之一，它是图像目标检测和分类的基 础，因此也是SAR目标识别系统的重要组成部分之一。由于SAR系统固有的相干斑效应，使得SAR图像中存在着乘性斑点噪声，因此与可 见光图像处理相比，SAR图像分割要复杂得多。现有的SAR图像分割方法主要包括基于数 据驱动和基于模型驱动两大类(1)基于数据驱动的SAR图像分割方法是首先对SAR图像进行斑点噪声抑制，然后 利用传统的可见光图像分割方法对斑点抑制后的SAR图像进行分割。比较经典的方法有 基于直方图最优阈值分割，基于边缘检测分割，基于区域分块合并分割等。(2)基于模型驱动的SAR图像分割方法是建立在对SAR图像特性分析的基础上，在 图像分割的过程中抑制斑点噪声，主要包括基于马尔可夫随机场(Markev Random Field, MRF)模型分割，基于组合优化模型分割，以及基于多尺度模型分割等。由于受到斑点噪声的影响，现有的SAR图像分割方法均存在分割区域边界粗糙， 分割区域内部均勻度不好的缺点。基于数据驱动的SAR图像分割方法，受SAR图像斑点噪声 抑制效果影响很大，而后续的图像分割处理不具有乘性噪声抑制能力，因此分割效果较差。 基于模型驱动的SAR图像分割方法，比较充分地利用了 SAR图像本身特性，因此分割精度较 高，但是分割速度非常慢，不利于SAR图像的快速处理，无法满足目前大量SAR图像信息快 速应用的需求。随着传感器技术的迅速发展，越来越多种类的传感器可以实现对同一区域的观 测，获取多种遥感信息。信息源的增加有可能提供更多的、相互补充的信息。美国、欧盟等 的现有信息获取系统就能够提供同一地区的SAR图像、可见光图像、红外图像等。例如美国 的JDL系统、英国的多平台多传感器系统等。适当的融合方法，能够很好地综合不同传感器 图像的互补信息，降低图像信息的不确定性，提高对遥感信息的感知和利用能力。因此，本 世纪以来，利用多传感器数据融合提高信息利用和图像分析性能的技术倍受关注。
发明内容
本发明的技术解决问题克服现有技术的不足，提供一种结合可见光图像信息的 SAR图像快速分割方法，该方法能够实现准确的分割，且分割速度大大提高。
本发明的技术解决方案一种结合可见光图像信息的SAR图像快速分割方法，实 现步骤如下第一步，进行可见光图像的双阈值分割(1. 1)首先进行分割阈值确定，确定方法如下利用人眼提取出可见光图像上的一块目标区域、」和一块背景区域RBkd，则目标阈 值和背景阈值TBkd的计算式分别为 其中，s0pt(x, y)为可见光图像的灰度值，(x，y)表示图像坐标；(1.2)遍历整个图像，将可见光图像灰度值大于目标阈值的像素点判定为目标像 素点，小于背景阈值的像素点判定为背景像素点，位于所述两个阈值之间的像素判定为不 确定像素点，计算公式为 式中，Sopt(X，y)为坐标(x，y)处可见光图像的灰度值，C(x，y)表示坐标(x，y)处 像素点的类别，这里像素点的类别有目标像素点、背景像素点和不确定像素点三类。第二步，将第一步得到的可见光图像的分割结果作为初始分割场以及图像的先验 信息，提供给SAR图像进行图像分割，分割方法具体为(2. 1)结合可见光图像的分割结果确定SAR图像的初始分割将可见光图像的目标像素点对应的SAR图像像素点直接标记为目标；将可见光图 像的背景像素点对应的SAR图像像素点直接标记为背景；可见光图像的不确定像素点对应 的SAR图像像素点利用最大似然方法确定每个像素sSAK(x，y)的初始分割标记兴，(x， y) G REand, REand为不确定像素的集合，有 其中，p(sSAK(X，y) z(x,y))为SAR图像的观测模型，采用瑞利分布，即 式中a (z (x, y))与瑞利分布的均值P K(z(x，y))和方差cr〗(z(x,_y))的关系为 (2. 2)不确定像素的分割随机改变可见光分割结果中标定为不确定的那些像素点对应的SAR图像中像素 点的标记，计算新标记下每个不确定像素的局部能量，有ELk(x, y) = - [lnp (sSAE (x, y) | zk (x, y)) +lnp (zk (x, y)) ], (x, y) G REand其中，zk(x, y)表示第k迭代下的新标记，k表示迭代次数，初次迭代时有k= 1 ； k = 0表示初始分割标记的本地能量，p(zk(X，y))为SAR图像的先验模型，它可以通过下式计算 其中，Z是归一化常数；β是(0，1)之间的一个常数；δ ( ·)是狄拉克函数；R(x,y) 是像素(X，y)的邻域；若新标记下每个不确定像素的局部能量低于未改变旧标记下每个不确定像素的 局部能量，即Eut(^y) > Euri(^y),接受新标记；反之，选取在
中均勻分布的随机变 量入，当exp[-(ELk (x，y)-Elh(XJ)VT]彡λ，则接受新标记；否则，保持当前标记；T为退 火温度，采用指数降温形式，有 其中，T0为初始退火温度，X为温度系数；计算当前状态和新状态的全局能量差AEk，有 其中I · I表示取绝对值，In为自然对数。如果全局能量的变化量小于收敛阈值ε，S卩AEkS ε，则得到最终分割结果；反 之，更新温度系数，令k = k+Ι，如果k < K，并重复该步骤处理过程；如果k > K，强行结束 程序，分割处理失败，ε为预先设置的收敛阈值，K为预先设置的迭代次数的上限值。本发明与现有技术相比的优点在于SAR是一种常用的微波遥感手段，它能够全 天时、全天候工作，可以获取高的空间分辨率，并且具有一定的穿透能力，但是其固有的斑 点噪声使得图像识别能力大大受限。可见光传感器具有多频谱特性，并且图像的可辨识性 能较好，但是其图像的获取受光线条件影响很大，在黑夜无法工作。本发明的目的在于利用 同一地区的可见光图像和SAR图像，提出一种结合可见光图像信息的SAR图像分割方法。本 发明利用可见光图像将SAR图像像素划分为目标像素、背景像素和不确定像素三类，然后 针对不同的图像像素采用不同的图像分割策略。本发明利用可见光图像信息给出的目标像 素和背景像素能够实现非常准确的分割；充分利用了可见光图像提供的先验信息，仅采用 基于模型驱动的SAR图像分割算法对不确定区域进行分割，因此可以大大提高分割速度。
图1为本发明实现流程图；图2为本发明的双阈值分割处理流程图；图3为本发明的SAR图像分割采用基于MRF的分割方法流程图；图4为天津地区SAR图像和可见光图像，其中a为1米分辨率SAR图像，b为可见 光图像；图5采用本发明方法得到的SAR图像分割结果，其中a为可见光图像双阈值分割 结果，b为结合可见光的SAR图像分割结果；图6为经典的基于MRF的SAR图像分割结果。
具体实施例方式本发明的实现流程如图1所示，具体的实现步骤主要包括 1.可见光图像的双阈值分割。灰度大于目标阈值的像素划分为目标像素点，小于 背景阈值的像素划分为背景像素，位于两个阈值之间的像素判定为不确定像素点。这些不 确定像素的划分类别主要依赖SAR分割来确定。可见光图像分割中两个阈值的确定必须保 证目标像素和背景像素判定的准确性。2.结合可见光图像分割结果的SAR图像分割。可见光图像的分割结果作为初始分 割场以及图像的先验信息，提供给SAR图像进行图像分割。可见光图像的双阈值分割具体处理流程如图2，主要包括以下步骤(1)分割阈值确定。利用人眼提取出可见光图像上的一块目标区域Raj和一块背 景区域RBkd，则目标阈值和背景阈值TBkd的计算式分别为 其中，s0pt (χ, y)为可见光图像的灰度值，(x, y)表示图像坐标。(2)双阈值分割。根据目标阈值和背景阈值，遍历整个图像，将可见光图像灰度值 大于目标阈值的像素点判定为目标像素点，小于背景阈值的像素点判定为背景像素点，位 于所述两个阈值之间的像素判定为不确定像素点，具体计算公式为 式中，S0pt (x，y)为坐标(x，y)处可见光图像的灰度值，C(x，y)表示坐标(x，y)处 像素点的类别，这里像素点的类别有目标像素点、背景像素点和不确定像素点三类。结合可见光图像分割结果的SAR图像分割采用基于马尔可夫随机场(Markev Random Fiel d，MRF)的分割方法，具体处理如图3，主要包括以下步骤(1)确定初始分割。结合可见光图像的分割结果确定SAR图像的初始分割将可 见光图像的目标像素点对应的SAR图像像素点直接标记为目标；将可见光图像的背景像素 点对应的SAR图像像素点直接标记为背景；可见光图像的不确定像素点对应的SAR图像像 素点利用最大似然方法确定每个像素sSAK (X，y)的初始分割标记f(xj)， 其中(x，y)e REand, REand 为不确定像素的集合，ρ (sSAE(x, y) | ζ (χ, y))为 SAR 图像 的观测模型，这里采用瑞利分布，即 式中α (ζ (χ, y))与瑞利分布的均值μκ(ζ(χ，y))和方差σ〗(Ζ(χ,>0)的关系为 它们可以通过EM方法估计获得。
(2)不确定像素的分割。随机改变可见光分割结果中标定为不确定的那些像素点 的标记，计算新标记下每个不确定像素的局部能量，有ELk(x, y) = - [lnp (sSAE (x, y) | zk (x, y)) +lnp (zk (x, y)) ], (x, y) G REand其中，zk(x，y)表示第k迭代下的新标记，k表示迭代次数，初次迭代时有k= 1 ；k =0表示初始分割标记的本地能量。p(zk(x，y))为SAR图像的先验模型，它可以通过下式
计算 其中，Z是归一化常数邛是(0，1)之间的一个常数，常选取为0.5; S ( )是狄 拉克函数；R(x,y)是像素(x，y)的邻域，常选取4邻域。(x' ,y')是R(x,y)中的像素坐标。若新标记能量低于旧标记能量，即EuU，y) > E^ (x, y)，接受新标记；反之，选取 在W，l]中均勻分布的随机变量入，当能量变化与温度之比的负指数大于等于该随机量， 即expt-^uU，y)-ELk_1(x, y))/T]彡X，则接受新标记；否则，保持当前标记；T为退火温 度，这里采用指数降温形式，有T = TQ>ck其中，I；为初始退火温度，x为温度系数。本发明给出的方法取1 = 4，x =0.9。计算当前状态和新状态的全局能量差AEk，有 如果全局能量的变化量小于收敛阈值e，即AEk< ￡，则得到最终分割结果；反 之，更新温度系数，令k = k+1，如果k < K，并重复该步骤处理过程。如果k > K，强行结束 程序，分割处理失败。￡为预先设置的收敛阈值，取值越小则收敛精度越高，但速度越慢，本 发明取￡ =0.1。K为预先设置的迭代次数的上限值，常选取K= 1000。图4中4a是天津地区的1米分辨率TerraSAR图像，4b是同一区域的快鸟卫星获 取的可见光图像。采用本发明给出的方法，对图4a进行双阈值分割得到图5a;结合图5a给 出的信息，对图4a给出的SAR图像进行快速分割处理，得到图5b。为了比较，同时给出了传 统单幅SAR图像基于MRF方法的分割结果如图6。为了更好地说明分割效果，表1给出了传统方法和本发明方法分割结果精度和速 度的评估结果。可以看出，与传统方法相比，本发明给出的方法在能够获取更好分割效果 的同时，提高了分割处理的速度，本实例中的分割精度提高了约两倍多，分割速度提高了约 50%。表1图像分割效果评估 本发明未详细阐述部分属于本领域技术公知技术。
权利要求
一种结合可见光图像信息的SAR图像快速分割方法，其特征在于实现步骤如下第一步，进行可见光图像的双阈值分割遍历整个图像，将可见光图像灰度值大于目标阈值TObj的像素点判定为目标像素点，小于背景阈值TBkd的像素点判定为背景像素点，位于所述两个阈值之间的像素判定为不确定像素点，计算公式为式中，sOpt(x，y)为坐标(x，y)处可见光图像的灰度值，C(x，y)表示坐标(x，y)处像素点的类别，这里像素点的类别有目标像素点、背景像素点和不确定像素点三类；第二步，将第一步得到的可见光图像的双阈值分割结果作为初始分割场以及图像的先验信息，提供给SAR图像进行图像分割，分割方法具体为(2.1)结合可见光图像的分割结果确定SAR图像的初始分割将可见光图像的目标像素点对应的SAR图像像素点直接标记为目标；将可见光图像的背景像素点对应的SAR图像像素点直接标记为背景；可见光图像的不确定像素点对应的SAR图像像素点利用最大似然方法确定每个像素sSAR(x，y)的初始分割标记(x，y)∈RRand，RRand为不确定像素的集合，有 <mrow><mover> <mi>z</mi> <mo>^</mo></mover><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow><mo>=</mo><mi>arg</mi><munder> <mi>max</mi> <mrow><mi>z</mi><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow> </mrow></munder><mi>p</mi><mrow> <mo>(</mo> <msub><mi>s</mi><mi>SAR</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>|</mo> <mi>z</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>,</mo> </mrow>(x，y)∈RRand其中，p(sSAR(x，y)|z(x，y))为SAR图像的观测模型，采用瑞利分布，即 <mrow><mi>p</mi><mrow> <mo>(</mo> <msub><mi>s</mi><mi>SAR</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>|</mo> <mi>z</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><msub> <mi>s</mi> <mi>SAR</mi></msub><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow> </mrow> <mrow><mi>α</mi><msup> <mrow><mo>(</mo><mi>z</mi><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow></mfrac><mi>exp</mi><mo>{</mo><mo>-</mo><mfrac> <mrow><msubsup> <mi>s</mi> <mi>SAR</mi> <mn>2</mn></msubsup><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow> </mrow> <mrow><mn>2</mn><mi>α</mi><msup> <mrow><mo>(</mo><mi>z</mi><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow><mo>)</mo> </mrow> <mn>2</mn></msup> </mrow></mfrac><mo>}</mo><mo>,</mo> </mrow>(x，y)∈RRand式中α(z(x，y))与瑞利分布的均值μR(z(x，y))和方差的关系为 <mrow><msub> <mi>μ</mi> <mi>R</mi></msub><mrow> <mo>(</mo> <mi>z</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mi>α</mi><mrow> <mo>(</mo> <mi>z</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><msqrt> <mfrac><mi>π</mi><mn>2</mn> </mfrac></msqrt> </mrow> <mrow><msubsup> <mi>σ</mi> <mi>R</mi> <mn>2</mn></msubsup><mrow> <mo>(</mo> <mi>z</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mfrac> <mrow><mn>4</mn><mo>-</mo><mi>π</mi> </mrow> <mn>2</mn></mfrac><msup> <mi>α</mi> <mn>2</mn></msup><mrow> <mo>(</mo> <mi>z</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>,</mo> </mrow>(x，y)∈RRand；(2.2)不确定像素的分割随机改变可见光分割结果中标定为不确定像素点对应的SAR图像像素点的标记，计算新标记下每个不确定像素的局部能量，有ELk(x，y)＝-[lnp(sSAR(x，y)|zk(x，y))+lnp(zk(x，y))]，(x，y)∈RRand其中，zk(x，y)表示第k迭代下的新标记，k表示迭代次数，初次迭代时有k＝1；k＝0表示初始分割标记的本地能量，p(zk(x，y))为SAR图像的先验模型，通过下式计算 <mrow><mi>p</mi><mrow> <mo>(</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mfrac> <mn>1</mn> <mi>Z</mi></mfrac><msup> <mi>e</mi> <mrow><mo>-</mo><mi>H</mi><mrow> <mo>(</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow> </mrow></msup><mo>,</mo> </mrow> <mrow><mi>H</mi><mrow> <mo>(</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>β</mi><munder> <mi>Σ</mi> <mrow><mrow> <mo>(</mo> <msup><mi>x</mi><mo>′</mo> </msup> <mo>,</mo> <msup><mi>y</mi><mo>′</mo> </msup> <mo>)</mo></mrow><mo>∈</mo><msub> <mi>R</mi> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow></msub> </mrow></munder><mo>[</mo><mi>δ</mi><mrow> <mo>(</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>-</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><msup> <mi>x</mi> <mo>′</mo></msup><mo>,</mo><msup> <mi>y</mi> <mo>′</mo></msup><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>-</mo><mn>1</mn><mo>]</mo> </mrow>其中，Z是归一化常数；β是(0，1)之间的一个常数；δ(·)是狄拉克函数；R(x，y)是像素(x，y)的邻域；若新标记下每个不确定像素的局部能量低于未改变旧标记下每个不确定像素的局部能量，即ELk(x，y)＞ELk-1(x，y)，接受新标记；反之，选取在
中均匀分布的随机变量λ，当exp[-(ELk(x，y)-ELk-1(x，y))/T]≥λ，则接受新标记；否则，保持当前标记；T为退火温度，采用指数降温形式，有T＝T0χk其中，T0为初始退火温度，χ为温度系数；计算当前状态和新状态的全局能量差ΔEk，有 <mrow><msub> <mi>ΔE</mi> <mi>k</mi></msub><mo>=</mo><mo>|</mo><mo>-</mo><munder> <mi>Σ</mi> <mrow><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow><mo>∈</mo><msub> <mi>R</mi> <mi>Rand</mi></msub> </mrow></munder><mi>ln</mi><mi>p</mi><mrow> <mo>(</mo> <msub><mi>s</mi><mi>SAR</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>|</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>+</mo><mi>ln</mi><mi>p</mi><mrow> <mo>(</mo> <msub><mi>z</mi><mi>k</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow> </mrow> <mrow><mo>+</mo><munder> <mi>Σ</mi> <mrow><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo></mrow><mo>∈</mo><msub> <mi>R</mi> <mi>Rand</mi></msub> </mrow></munder><mi>ln</mi><mi>p</mi><mrow> <mo>(</mo> <msub><mi>s</mi><mi>SAR</mi> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>|</mo> <msub><mi>z</mi><mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn></mrow> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>+</mo><mi>ln</mi><mi>p</mi><mrow> <mo>(</mo> <msub><mi>z</mi><mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn></mrow> </msub> <mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>|</mo> </mrow>其中|·|表示取绝对值，ln为自然对数；如果全局能量的变化量小于收敛阈值ε，即ΔEk≤ε，则得到最终分割结果；反之，更新温度系数，令k＝k+1，如果k≤K，并重复该步骤处理过程；如果k＞K，强行结束程序，分割处理失败，ε为预先设置的收敛阈值，K为预先设置的迭代次数的上限值。FSA00000181511200011.tif,FSA00000181511200012.tif,FSA00000181511200015.tif
2.根据权利要求1所述的结合可见光图像信息的SAR图像快速分割方法，其特征在于 所述第一步中目标阈值T^和背景阈值TBkd的确定方法为利用人眼提取出可见光图像上的一块目标区域‘」和一块背景区域RBkd，则目标阈值 TobJ和背景阈值TBkd的计算式分别为 其中，s0pt(x, y)为可见光图像的灰度值，(X，y)表示图像坐标。
全文摘要
一种结合可见光图像信息的SAR图像快速分割方法，实现如下可见光图像的双阈值分割。灰度大于目标阈值的像素划分为目标像素点，小于背景阈值的像素划分为背景像素，位于两个阈值之间的像素判定为不确定像素点；结合可见光图像分割结果的SAR图像分割。可见光图像的分割结果作为初始分割场以及图像分割的先验信息，提供给SAR图像进行图像分割。本发明能够实现准确的分割，且分割速度大大提高。
文档编号G06T7/00GK101872481SQ20101020494
公开日2010年10月27日 申请日期2010年6月18日 优先权日2010年6月18日
发明者刘向华, 周荫清, 徐华平 申请人:北京航空航天大学