发明名称 一种快速去雾方法
摘要 本发明涉及一种快速去雾方法,其包括如下步骤:对每个像素点求RGB的最小值得到灰度图像,对灰度图像采用最小值滤波得到暗原色图像,利用暗原色图像估计大气参数并均匀选取图像中一些点求其透射率,利用二元二次线性回归分析得到图像所有点的透射率,根据透射率和大气参数恢复无雾图像,本算法具有处理效果好和速度快的特点,能够很容易的到实时性要求较高的系统中。
申请公布号 CN102063706B 申请公布日期 2012.08.08
申请号 CN201010620645.7 申请日期 2010.12.23
申请人 哈尔滨工业大学(威海) 发明人 王好贤;毛兴鹏;李方
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 代理人
主权项 1.一种快速去雾方法,其特征在于,该方法包括如下步骤:第一步:读取原始有雾图像I;第二步:求原始图像I中每个像素点的R、G、B三个通道的最小值,并赋值给当前像素点,得到的灰度图像记为I<sub>g</sub>,若I为灰度图像,则I<sub>g</sub>=I;第三步:利用公式(1)对第二步得到的I<sub>g</sub>图像进行滤波,定义Ω(x,y)为以坐标(x,y)的像素点为中心的大小为N*N的模版区域,其中取N∈(9,21),N越大求得的大气参数越小,得到无雾图像的整体亮度相对越亮;<maths num="0001"><![CDATA[<math><mrow><msup><msub><mi>I</mi><mi>g</mi></msub><mi>dark</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>min</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>&Element;</mo><mi>&Omega;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow></munder><mrow><mo>(</mo><msub><mi>I</mi><mi>g</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>第四步:得到I<sub>g</sub><sup>dark</sup>图像中强度最大的像素点的位置坐标(x<sub>max</sub>,y<sub>max</sub>),若I为彩色图像,则大气参数A<sub>R</sub>=I<sub>R</sub>(x<sub>max</sub>,y<sub>max</sub>)、A<sub>G</sub>=I<sub>G</sub>(x<sub>max</sub>,y<sub>max</sub>)和A<sub>B</sub>=I<sub>B</sub>(x<sub>max</sub>,y<sub>max</sub>),其中I<sub>B</sub>、I<sub>G</sub>和I<sub>B</sub>分别为彩色图像I的R、G和B通道,若I为灰度图像,其大气参数A<sub>gray</sub>=I(x<sub>max</sub>,y<sub>max</sub>);第五步:求取图像的长度h和宽度w,然后利用公式(2)对I<sub>g</sub><sup>dark</sup>图像中满足坐标条件的像素点求其透射率的值,坐标(x,y)的取值满足以下条件中x和y的所有组合:x=1、1+S、1+2*S、…、h,y=1、1+S、1+2*S、…、w,步长S∈(40,100),S越大,第六步的处理速度越快,但是估计的回归方程的系数相对不太准确,去雾效果相对差些,若I为彩色图像,则式(2)中<img file="FSB00000821197800012.GIF" wi="416" he="122" />若I为灰度图像,则式(2)中A=A<sub>gray</sub>,ω0为调节参数,其中ω0∈(0,1),ω0越大去雾效果越好;<maths num="0002"><![CDATA[<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mi>&omega;</mi><mn>0</mn><mrow><mo>(</mo><mfrac><mrow><msup><msub><mi>I</mi><mi>g</mi></msub><mi>dark</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><mi>A</mi></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>第六步:利用二元二次线性回归分析对第五步涉及的像素点及对应的透射率进行回归分析,得到如式(3)的方程的系数b<sub>0</sub>、b<sub>1</sub>、b<sub>2</sub>、b<sub>3</sub>、b<sub>4</sub>和b<sub>5</sub>;t=b<sub>0</sub>+b<sub>1</sub>*x+b<sub>2</sub>*y+b<sub>3</sub>*xy+b<sub>4</sub>*x<sup>2</sup>+b<sub>5</sub>*y<sup>2</sup>    (3)第七步:利用公式(3)得到所有像素点的透射率t;第八步:利用公式(4)恢复无雾图像,若I为彩色图像,则分别在R、G、B三个通道利用公式(4)恢复无雾图像J(x),其中A分别取A<sub>R</sub>、A<sub>G</sub>和A<sub>B</sub>,I分别取I<sub>R</sub>、I<sub>G</sub>和I<sub>B</sub>,则分别得到J的三个通道J<sub>R</sub>、J<sub>G</sub>和J<sub>B</sub>,若I为灰度图像,则A=A<sub>gray</sub>。<maths num="0003"><![CDATA[<math><mrow><mi>J</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mi>A</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>t</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mi>t</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>
地址 264209 山东省威海市文化西路2号哈尔滨工业大学(威海)