本发明专利技术涉及一种基于加权信息熵和马尔可夫随机场的红外目标分割方法，技术特征在于：首先利用加权信息熵得到感兴趣区域，继而利用高斯马尔可夫随机场在感兴趣区域内建模来捕捉图像的上下文信息，并在迭代中采用ICM优化算法，使得解空间大大减小、噪声干扰大大降低，从而迭代步数减少、分割精度有明显的改善、分割速度大大提高。许多红外目标分割实验，都验证了本发明专利技术的有效性。
【技术实现步骤摘要】
本专利技术涉及一种，属于红 外图像处理方法，具体是红外目标的分割方法。
技术介绍
图像目标分割是图像分析以及模式识别的关键技术，也是目标识别、跟踪以及其 它计算机视觉相关应用的基础，其最终目的是为后续的处理和场景描述提供基础。红外图 像反映目标及周围背景的温度分布情况，但是由于探测器本身的固有的特性及外界成像环 境的影响，导致红外图像SNR(信噪比)较低，呈现出目标与背景对比度较差、目标的边缘模 糊、噪声较大等特点，因此红外目标分割是一个艰巨的任务。基于阈值的分割方法是在红外 目标分割中被广泛应用的，例如一维Otsu、二维Otsu、二维最大熵法等，但是由于红外图 像的成像特点所致，传统的分割方法的分割效果不能令人满意。马尔可夫随机场(Markov Random Field, MRF)可以捕获图像中像素间的相 互作用和约束，并能将上下文信息引入到图像的分割过程，并具有较强的噪声抑制能 力。文献"Unsupervised image segmentation using a simple MRF model with a newimplementation scheme, Pattern Recognition，vol. 37，pp. 2323-2335，2004. ” 在分析 传统方法固定加权缺点的基础上提出了一种可变的权值。该方法利用简化的MRF对图像建 模，并在区域成分和特征成分之间采用指数形式下降的权重，最后利用模拟退火算法来优 化目标函数，从而得到最终的分割结果。但是此方法由于是在整幅图上分割，噪声干扰大， 分割的准确性不高，而且利用模拟退火优化算法，使得分割速度大大降低。红外图像SNR低，噪声干扰大，目标边缘模糊，背景和目标对比度差；基于MRF的分 割算法存在着模型复杂、缺乏自适应性、收敛速度慢等缺点。
技术实现思路
要解决的技术问题为了避免现有技术的不足之处，本专利技术提出一种基于加权信息熵和马尔可夫随机 场的红外目标分割方法，解决红外图像SNR低，噪声干扰大，目标边缘模糊，背景和目标对 比度差，造成目标分割存在着模型复杂、缺乏自适应性、收敛速度慢的缺点。本专利技术的思想在于基于加权信息熵(weighted information entropy,WIE) 了感 兴趣区域(region of interest,R0I)提取，在ROI内利用高斯马尔可夫随机场(Gaussian MarkovRandom Field, GMRF)建模，并利用迭代条件模型(Iterated Condition Mode, I CM) 优化算法，缩小了解空间，极大的提高分割速度。技术方案一种，其特征在于步骤如 下步骤1基于加权信息熵的感兴趣区域提取步骤a:以大小为mX η的图像块来平分MX N图像，m = 2k，η = 21，k，1 = 1，2， 3，…，得到平分后的若干个子图像；步骤b 对每幅子图像F (U，ν)，其中0≤u≤M/m-1,0 ≤ ν ≤N/n_l，利用 when ps = O, let ps*log(ps) = O计算每幅子图像的加权信息熵WIE值H(u，v)，得到M/(m-l) XN/(n_l)大小的熵 图像，并计算出熵图像的均值μ和方差s ；步骤c 当最大的熵值的熵图像H(u，ν)≥Ht作为种子，其中Ht = μ+6Χ δ ；如果 H(u，ν) <HT = m/2，n = n/2，重复步骤a b；步骤d 在熵图像H(u，v)上，结合种子利用四邻域进行区域增长，得到包含目标的 矩形区域；步骤2 在矩形区域内利用马尔可夫随机场建模，然后做随机初始分割，利用ICM 进行迭代优化得到目标区域与背景区域的分割图像；所述迭代优化中的权值修正方案为 ，其中t表示迭代步数。有益效果本专利技术提出的，首先利用 加权信息熵做感兴趣区域提取，得到包含目标的矩形框，缩小了解空间、剔除了大量的噪声 干扰；然后在感兴趣区域内利用高斯马尔可夫随机场建模，可以捕获红外图像中的上下文 信息，有效克服红外噪声的影响，同时利用ICM算法迭代优化，可以在很少的分割步骤内达 到收敛。附图说明图1 本专利技术方法流程图 具体实施例方式现结合实施例、附图对本专利技术作进一步描述1.基于加权信息熵的感兴趣区域提取为了估计红外图像强度分布的复杂度，我们定义了 WIE，设s表示灰度值，s ∈，ps表示每种灰度级对应的概率，5表示红外图像的灰度平均值，则图像的加权信息熵 可以定义为 when ps = 0, let ps*log(ps) = 0根据公式可知，目标所在区域的熵值必定大。基于这点，我们可以把图像分成一定 区域大小的块，然后根据熵值的大小来判断种子块，最后在熵图像上利用区域增长来获取 ROI。设红外图像I，大小是MX N，则ROI提取算法描述如下(1)以大小为mXn的子图像块来平分图像，m= Z1Sn = Z1, k，l = 1，2，3，…，子 图像为F(u，ν)，其中0≤u≤M/m-1，0≤ν≤Ν/n-l。(2)对每幅子图像F (u，ν)，利用公式(1)计算WIE值H (u，ν)，得到M/(m_l) XN/ (n-1)大小的熵图像，并计算出熵图像的均值μ和方差δ。(3)获取种子。假设熵图像上最大的熵值为H(u，v)，如果H(u，v) ，其中Ht = μ+α* δ，这里 α 取 6，则 F(u，ν)即为种子。如果 H(u，ν) <乐，取111 = 111/2，η = η/2, M 复步骤(1)、(2)，直到熵图像的最大值满足H(u，ν)彡乐或!11 = 2或11 = 2。(4) ROI提取。在熵图像上，结合种子利用四邻域进行区域增长，得到包含目标的矩 形区域。2.利用高斯马尔可夫随机场对红外图像建模在Bayes框架下该后验概率可以表示为， 「 1 W lv 、P(YX^x)*P(X = X) 在给定观测图像的条件下P(Y = y)是常数，P(Y = y|X = x)是条件概率，它体现 了分割过程中的特征模型成分；Ρ(χ = χ)描述了分割结果标记的分布，它体现了分割过程 中的区域标记成分。对于马尔可夫场X，在Hammersley-Clifford定理的保证下，可以通过定义Gibbs 分布中势能函数来确定其分布P (X = X)，大多数基于MRF的分割模型都利用多层逻辑模型 (multi-level logistic, MLL)来定义势能函数，对于一般的分割任务来说，我们仅定义二 阶双点基团的势能而将其余二阶基团的势能都定义为0 (3)这里如果Xs = Xt，则 ；如果 Xs 乒 Xt，则 是事先 定义好的常量。然后Gibbs分布的形式的就可以表示成如下形式 (4)其中Z = 2L。exP[-+A⑴堤划分函数，T是温度参数。GMRF模型中各点特征服从独立的高斯分布，则采用伪似然的形式将条件分布表示 为， 其中μ m，σ m分别代表第m类的均值和标准差，这两个参数需要通过EM算法来估 计，第t+Ι次迭代的估计公式如(6) (7)所示 则后验能量可以表示为(7)式， ,Ef=Y^J +In(VI^xm)]表示能量中的特征成分，而用公式⑶表示能量中的区域成分，则能量可以表示为这两部分加权和的形式E = W*Ef+Er(9)由前述可知，分本文档来自技高网...
【技术保护点】
一种基于加权信息熵和马尔可夫随机场的红外目标分割方法，其特征在于步骤如下：步骤１基于加权信息熵的感兴趣区域提取：步骤ａ：以大小为ｍ×ｎ的图像块来平分Ｍ×Ｎ图像，ｍ＝２↑［ｋ］，ｎ＝２ｌ，ｋ，ｌ＝１，２，３，．．．，得到平分后的若干个子图像；步骤ｂ：对每幅子图像Ｆ（ｕ，ｖ），其中０≤ｕ≤Ｍ／ｍ－１，０≤ｖ≤Ｎ／ｎ－１，利用Ｈ（ｓ）＝－＊（ｓ－＊）↑［２］＊ｐ↓［ｓ］＊ｌｏｇ（ｐ↓［ｓ］）ｗｈｅｎｐ↓［ｓ］＝０，ｌｅｔｐ↓［ｓ］＊ｌｏｇ（ｐ↓［ｓ］）＝０计算每幅子图像的加权信息熵ＷＩＥ值Ｈ（ｕ，ｖ），得到Ｍ／（ｍ－１）×Ｎ／（ｎ－１）大小的熵图像，并计算出熵图像的均值μ和方差δ；步骤ｃ：当最大的熵值的熵图像Ｈ（ｕ，ｖ）≥Ｈ↓［Ｔ］作为种子，其中Ｈ↓［Ｔ］＝μ＋６×δ；如果Ｈ（ｕ，ｖ）＜Ｈ↓［Ｔ］，将ｍ＝ｍ／２，ｎ＝ｎ／２，重复步骤ａ～ｂ；步骤ｄ：在熵图像Ｈ（ｕ，ｖ）上，结合种子利用四邻域进行区域增长，得到包含目标的矩形区域；步骤２：在矩形区域内利用马尔可夫随机场建模，然后做随机初始分割，利用ＩＣＭ进行迭代优化得到目标区域与背景区域的分割图像；所述迭代优化中的权值修正方案为Ｗ＝８００＊０．９７′＋１，其中：ｔ表示迭代步数。...
【技术特征摘要】
【专利技术属性】
技术研发人员：李映，毛星锦，张艳宁，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：87
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。