本发明专利技术涉及基于局部对比度累积的水中航行器热尾迹分割方法。本发明专利技术的目的是对水中航行器产生的热尾迹对应的红外图像中的热尾迹区域实现完整有效的分割。首先采用同态滤波和均值滤波对图像进行预处理，然后对预处理得到的图像进行基于局部对比度累积的显著性检测，得到初步分割结果，最后采用形态学处理及最大连通域提取等操作对结果进行了优化，取得了更好的分割效果。说明书附图中图1为发明专利技术的具体实现流程图。
Segmentation of Thermal Wake of Underwater Vehicle Based on Local Contrast Accumulation
【技术实现步骤摘要】
基于局部对比度累积的水中航行器热尾迹分割方法
：本专利技术属于图像处理领域，具体地说是通过分析水中航行器在水面产生的热尾迹红外图像特点，实现了对信噪比及对比度较低的红外图像中热尾迹区域的有效分割。
技术介绍
：热红外探测水下航行器是利用高精度红外探测仪探测来获取水面红外图像，获取含有热尾迹的区域进而实现对该区域的分析与识别。研究水中航行器热尾迹分割算法实现对热尾迹区域完整有效的分割对红外探测水下航行器具有重要意义。水下航行器在温度均匀的海水中由于动力系统放热、螺旋桨打水和运动产生摩擦等原因产生温度高于周围海水的温差信号，由于浮力作用温差信号会传递至海面。通过对温度均匀、正温度梯度和负温度梯度等三种典型环境条件进行的实验验证及仿真分析结果表明，潜航器在温度均匀及正温度梯度条件下会形成正温差信号(尾迹区域海水温度高于周围区域)；而在负温度梯度条件下则会形成负温差信号(尾迹区域海水温度低于周围区域)，我们将这种温差信号统称为热尾迹。探测人员利用机载或者星载高精度的红外探测仪能够探测到海面的弱小温差信号，也就是说可以利用红外探测仪探测热尾流来间接判断水下航行器的存在性，这也是热红外探测水下航行器原理，它对我国海防具有重要的价值。红外探测水下航行器技术由于其观测范围大、可在夜间工作，抗电子干扰强、隐蔽性好、被动工作、实时成像、成像分辨率高等优点不断被研究发展，红外探测属于“战略性”探测、“趋势”探测。水中航行器热尾迹分割算法可以实现对热尾迹区域进行有效分割，有助于排除背景干扰，对热尾迹区域进行更加有效的分析与研究，进而提升红外探测水下航行器系统的性能。因此，该研究具有重要意义。红外图像是由红外探测器接收到的目标区域红外辐射信号转换得到的。对于水面背景，水面自身辐射、水面反射天空辐射和太阳辐射等因素都会对红外图像产生影响，导致图像信噪比低、细节信息少且轮廓模糊。本专利技术针对由这些原因导致的背景亮度不均匀、低频背景干扰及盐噪声等三个主要干扰对图像进行了预处理。另外，考虑到热尾迹生成及浮升过程，热尾迹区域会存在难以检测的弱边缘及不明显的低温区域。本专利技术采用局部对比度累积的显著性检测方法，从空间域局部显著性的角度来对热尾迹区域进行分割，并采用形态学处理及最大连通域提取等操作对分割结果进行了优化。最终实现了对水面背景红外图像热尾迹区域完整有效的分割。本专利技术通过预处理、基于局部对比度累积的显著性检测、分割结果优化等主要操作，主要从空间域局部对比度特征的角度，实现了对水面背景红外图像热尾迹区域完整有效的分割。
技术实现思路
：本专利技术的目的是对水中航行器产生的热尾迹对应的红外图像中的热尾迹区域实现完整有效的分割。首先采用同态滤波和均值滤波对图像进行预处理，然后对预处理得到的图像进行基于局部对比度累积的显著性检测，得到初步分割结果，最后采用形态学处理及最大连通域提取等操作对结果进行了优化，取得了更好的分割效果。本专利技术采用的技术方案是：基于局部对比度累积的水中航行器热尾迹分割方法，其特征在于如下步骤：第一步：读取原始热尾迹图像I0，图像大小为M×N，M是图像I0的宽度像素，N是图像I0的高度像素；第二步：若图像I0为彩色图像，则通过式(1)对图像I0中每个像素点的R、G、B三个通道的值进行加权，得到灰度图像Ig；若图像I0为灰度图像，则令Ig＝I0；Ig＝0.29900×R+0.58700×G+0.11400×B(1)第三步：对图像Ig进行灰度反转，通过对图像Ig中每个像素点的灰度值进行式(2)所示的处理来实现，将处理后得到的图像记为f，其中f(x,y)表示图像f在坐标(x,y)处像素点的取值，Ig(x,y)表示图像Ig在坐标(x,y)处像素点的取值；f(x,y)＝255-Ig(x,y)(2)第四步：对图像f进行同态滤波，通过以下步骤来实现：(a)假设图像f可以由入射光分量图像i和反射光分量图像r的和来表示，如式(3)所示，其中f(x,y)、i(x,y)和r(x,y)分别表示图像f、图像i和图像r在坐标(x,y)处像素点的取值；f(x,y)＝i(x,y)+r(x,y)(3)(b)对式(3)的等号两侧取以e为底的对数得到式(4)；ln[f(x,y)]＝ln[i(x,y)]+ln[r(x,y)](4)(c)对式(4)的等号两侧分别进行二维傅里叶变换得到式(5)，其中F(u,v)、I(u,v)和R(u,v)分别表示对ln[f(x,y)]、ln[i(x,y)]和ln[r(x,y)]进行二维傅里叶变换得到的结果，(u,v)表示进行二维傅里叶变换后的频域坐标；F(u,v)＝I(u,v)+R(u,v)(5)(d)对F(u,v)进行高通滤波，通过式(6)来实现，其中S(u,v)为对F(u,v)进行高通滤波后得到的结果，符号·表示将H(u,v)和F(u,v)中相同坐标处的两个元素的值对应相乘。H(u,v)各元素的取值如式(7)所示，其中rh、rl、c和D0为需要设定的参数(rh应大于1，rl应小于1)，D(,u)v的定义如式(8)所示，其中P和Q的取值分别为F(u,v)中u和v最大取值的2倍；S(u,v)＝H(u,v)·F(u,v)(6)D(u,v)＝[(u-P/2)2+(v-Q/2)2]1/2(8)(e)对S(u,v)进行二维傅里叶逆变换得到s(x,y)，如式(9)所示，其中表示二维傅里叶逆变换；(f)对s(x,y)进行以e为底的指数操作，如式(10)所示，得到同态滤波后的图像g，式中g(x,y)表示图像g中坐标(x,y)处的灰度值，exp表示取以e为底的指数操作；g(x,y)＝exp[s(x,y)](10)第五步：采用窗口大小为8×8个像素的中值滤波器对图像g进行去噪处理，该操作可以通过MATLAB软件中的Imedfilt＝medfilt2(g,[8,8])语句来实现，Imedfilt是图像经过窗口大小为8×8个像素的中值滤波器滤波得到的图像；第六步：设置一个大小为9×9个像素的正方形滑动窗口，将这个滑动窗口中划分为9个3×3个像素大小的不重叠的正方形区域，将位于中间的正方形区域称为中心图像块L0，并把它周围的8个正方形区域分别称为图像块L1～L8，将这些图像块简记为Li(i＝0,1,2,3...8)；在窗口内进行如下操作：(a)计算图像块Li中像素的灰度均值和灰度中值，如式(11)和式(12)所示，其中mean_i表示图像块Li的灰度均值，med_i表示图像块Li的灰度中值，表示图像块Li中第k个像素的灰度值；med(·)表示获取中值；(b)对中心图像块L0周围的8个图像块(即i＝1,2,3...8时)，计算能够代表图像块Li的灰度值gray_i，通过式(13)实现，其中max_i表示图像块Li内像素点的最大灰度值，threshold为需要设定的阈值常数；(c)通过式(14)计算窗口内s_LCA的取值，并将其作为灰度值赋给中心图像块L0中的每一个像素(若s_LCA的值大于255，则令s_LCA的值为255)，其中mean_0表示中心图像块L0中的像素灰度值的均值，可由式(11)求得，|·|表示取模值；第七步：将第六步的窗口以3个像素为步长在图像Imedfilt上从左至右、从上至下滑动，每滑动到一个新的位置则进行一次第六步的操作，当窗口在I本文档来自技高网...
【技术保护点】
1.基于局部对比度累积的水中航行器热尾迹分割方法，其特征在于如下步骤：第一步：读取原始热尾迹图像I0，图像大小为M×N，M是图像I0的宽度像素，N是图像I0的高度像素；第二步：若图像I0为彩色图像，则通过式(1)对图像I0中每个像素点的R、G、B三个通道的值进行加权，得到灰度图像Ig；若图像I0为灰度图像，则令Ig＝I0；Ig＝0.29900×R+0.58700×G+0.11400×B (1)第三步：对图像Ig进行灰度反转，通过对图像Ig中每个像素点的灰度值进行式(2)所示的处理来实现，将处理后得到的图像记为f，其中f(x,y)表示图像f在坐标(x,y)处像素点的取值，Ig(x,y)表示图像Ig在坐标(x,y)处像素点的取值；f(x,y)＝255‑Ig(x,y) (2)第四步：对图像f进行同态滤波，通过以下步骤来实现：(a)假设图像f可以由入射光分量图像i和反射光分量图像r的和来表示，如式(3)所示，其中f(x,y)、i(x,y)和r(x,y)分别表示图像f、图像i和图像r在坐标(x,y)处像素点的取值；f(x,y)＝i(x,y)+r(x,y) (3)(b)对式(3)的等号两侧取以e为底的对数得到式(4)；ln[f(x,y)]＝ln[i(x,y)]+ln[r(x,y)] (4)(c)对式(4)的等号两侧分别进行二维傅里叶变换得到式(5)，其中F(u,v)、I(u,v)和R(u,v)分别表示对ln[f(x,y)]、ln[i(x,y)]和ln[r(x,y)]进行二维傅里叶变换得到的结果，(u,v)表示进行二维傅里叶变换后的频域坐标；F(u,v)＝I(u,v)+R(u,v) (5)(d)对F(u,v)进行高通滤波，通过式(6)来实现，其中S(u,v)为对F(u,v)进行高通滤波后得到的结果，符号·表示将H(u,v)和F(u,v)中相同坐标处的两个元素的值对应相乘。H(u,v)各元素的取值如式(7)所示，其中rh、rl、c和D0为需要设定的参数(rh应大于1，rl应小于1)，D(,u)v的定义如式(8)所示，其中P和Q的取值分别为F(u,v)中u和v最大取值的2倍；S(u,v)＝H(u,v)·F(u,v) (6)...
【技术特征摘要】
1.基于局部对比度累积的水中航行器热尾迹分割方法，其特征在于如下步骤：第一步：读取原始热尾迹图像I0，图像大小为M×N，M是图像I0的宽度像素，N是图像I0的高度像素；第二步：若图像I0为彩色图像，则通过式(1)对图像I0中每个像素点的R、G、B三个通道的值进行加权，得到灰度图像Ig；若图像I0为灰度图像，则令Ig＝I0；Ig＝0.29900×R+0.58700×G+0.11400×B(1)第三步：对图像Ig进行灰度反转，通过对图像Ig中每个像素点的灰度值进行式(2)所示的处理来实现，将处理后得到的图像记为f，其中f(x,y)表示图像f在坐标(x,y)处像素点的取值，Ig(x,y)表示图像Ig在坐标(x,y)处像素点的取值；f(x,y)＝255-Ig(x,y)(2)第四步：对图像f进行同态滤波，通过以下步骤来实现：(a)假设图像f可以由入射光分量图像i和反射光分量图像r的和来表示，如式(3)所示，其中f(x,y)、i(x,y)和r(x,y)分别表示图像f、图像i和图像r在坐标(x,y)处像素点的取值；f(x,y)＝i(x,y)+r(x,y)(3)(b)对式(3)的等号两侧取以e为底的对数得到式(4)；ln[f(x,y)]＝ln[i(x,y)]+ln[r(x,y)](4)(c)对式(4)的等号两侧分别进行二维傅里叶变换得到式(5)，其中F(u,v)、I(u,v)和R(u,v)分别表示对ln[f(x,y)]、ln[i(x,y)]和ln[r(x,y)]进行二维傅里叶变换得到的结果，(u,v)表示进行二维傅里叶变换后的频域坐标；F(u,v)＝I(u,v)+R(u,v)(5)(d)对F(u,v)进行高通滤波，通过式(6)来实现，其中S(u,v)为对F(u,v)进行高通滤波后得到的结果，符号·表示将H(u,v)和F(u,v)中相同坐标处的两个元素的值对应相乘。H(u,v)各元素的取值如式(7)所示，其中rh、rl、c和D0为需要设定的参数(rh应大于1，rl应小于1)，D(,u)v的定义如式(8)所示，其中P和Q的取值分别为F(u,v)中u和v最大取值的2倍；S(u,v)＝H(u,v)·F(u,v)(6)D(u,v)＝[(u-P/2)2+(v-Q/2)2]1/2(8)(e)对S(u,v)进行二维傅里叶逆变换得到s(x,y)，如式(9)所示，其中表示二维傅里叶逆变换；(f)对s(x,y)进行以e为底的指数操作，如式(10)所示，得到同态滤波后的图像g，式中g(x,y)表示图像g中坐标(x,y)处的灰度值，exp表示取以e为底的指数操作；g(x,y)＝exp[s(x,y)](10)第五步：采用窗口大小为8×8个像素的中值滤波器对图像g进行去噪处理，该操作可以通过MA...
【专利技术属性】
技术研发人员：王好贤，董衡，周志权，
申请(专利权)人：哈尔滨工业大学威海，
类型：发明
国别省市：山东,37
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。