数字图像处理课程设计报告图像去雾

1、-. z数字图像处理课设题目：图像去雾学 院：信息与电气工程学院专 业：电子信息工程班 级：姓 名：学 号：指导教师：工业大学威海 年月日目 录 TOC o 1-3 h z u HYPERLINK l _Toc266026732一. 课程设计任务 PAGEREF _Toc266026732 h 1HYPERLINK l _Toc266026733二. 课程设计原理及设计方案 PAGEREF _Toc266026733 h 2HYPERLINK l _Toc266026734三. 课程设计的步骤和结果 PAGEREF _Toc266026734 h 3HYPERLINK l _Toc266026

2、735四. 课程设计总结 PAGEREF _Toc266026735 h 4HYPERLINK l _Toc266026736五. 设计体会 PAGEREF _Toc266026736 h 5HYPERLINK l _Toc266026737六. 参考文献 PAGEREF _Toc266026737 h 6课程设计任务由于大气的散射作用，雾天的大气退化图像具有比照度低、景物不清晰的特点，给交通系统及户外视觉系统的应用带来严重的影响。但由于成像系统聚焦模糊、拍摄场景存在相对运动以及雾天等不利环境，使得最终获取的图像往往无法使用。为了解决这一问题，设计图像复原处理软件。要求完成功能：采用直方图均衡

3、化方法增强雾天模糊图像，并比拟增强前后的图像和直方图；查阅文献，分析雾天图像退化因素，设计一种图像复原方法，比照该复原图像与原始图像以及直方图均衡化后的图像；设计软件界面课程设计原理及设计方案2.1 设计原理 在雾、霾等天气条件下, 大气中悬浮的大量微小水滴、气溶胶的散射作用导致捕获的图像严重降质,随着物体到成像设备的距离增大, 大气粒子的散射作用对成像的影响逐渐增加. 这种影响主要由两个散射过程造成: 1) 物体外表的反射光在到达成像设备的过程中, 由于大气粒子的散射而发生衰减;2) 自然光因大气粒子散射而进入成像设备参与成像. 它们的共同作用造成捕获的图像比照度、饱和度降低, 以及色调偏移

4、, 不仅影响图像的视觉效果, 而且影响图像分析和理解的性能.在计算机视觉领域中, 常用大气散射模型来描述雾、霾天气条件下场景的成像过程.Narasimhan 等给出雾、霾天气条件下单色大气散射模型(Monochrome atmospheric scat-tering model), 即窄波段摄像机所拍摄的图像灰度值可表示为 1式中, * 为空间坐标, A 表示天空亮度(Skylight),为场景反照率, d 为场景的景深,为大气反射系数。图像去雾的过程就是根据获得的有降质的图像，即I*来推算出的过程。 但由于该物理模型包含3 个未知参数, 从本质上讲, 这是一个病态反问题.在只有单幅图像的条件

5、下，我们可以考虑用假设以及推算的方式使其中的几个量固定，然后求解。2.2 设计方案2.2.1 白平衡 WP (White point) 算法, 也称为Ma*-RGB 算法, 利用R、G、B 颜色分量的最大值来估计光照的颜色。我们用天空亮度A来取代最大值。对于A，如果直接用图像中最亮的灰度估计的话会受到高亮噪声或白色物体的干扰，因此我们首先对图像颜色分量进展最小滤波，然后选择阀值为0.99，大于此值的认定为天空区域，然后取平均值为我们估计的天空亮度A。然前方程1变为 22.2.2 估计大气耗散函数 为了简化2，我们设则 2变为3V*称为大气耗散函数。根据3容易发现V*满足0V(*)I,如果对I进

6、展最小滤波得到W*，则V*也必定满足0V(*)W(*)。求V*的过程其实是滤波的过程，分为三步：其中p为恢复系数，去0,1中的值一般为0.90.95。2.2.3 场景反射率的恢复 通过以上步骤我们获得了A以及V*，至此方程1仅含有一个未知量，即我们要求的场景反射率2.2.4 色调映射 映射关系为其中G为求的去雾后的图像灰度最大值。2.2.5 与直方图去雾效果比照使用matlab中的histeq函数对图像的亮度和比照度进展直方图均衡化处理课程设计的步骤和结果3.1原图像与直方图global im;global orig;filename, pathname=uigetfile(*.jpg,ins

7、ert image);im=pathname filename;orig = imread(im);Q = rgb2gray(orig);figure; subplot(1, 2, 1); imshow(orig); title(原图像)subplot(1, 2, 2); imhist(Q, 64); title(原图像直方图)3.2原图像与均衡后图像global orig;global U;I = rgb2hsv(orig);H = I(:,:,1);S = I(:,:,2);V = I(:,:,3);M = histeq(S); % 对饱和度直方图均衡化N = histeq(V); % 对

8、亮度直方图均衡化U = hsv2rgb(H, M, N);figure; subplot(2, 2, 1); imshow(orig); title(原图像)subplot(2, 2, 2); imshow(U); title(均衡后图像)Q = rgb2gray(orig);W = rgb2gray(U);subplot(2, 2, 3); imhist(Q, 64); title(原图像直方图)subplot(2, 2, 4); imhist(W, 64); title(均衡变换后的直方图)3.3原图与去雾图像比照global im;im1=double(imread(im)/255.0;

9、sv=2*floor(ma*(size(im1)/25)+1;res=visibresto1(im1,sv,0.95,-1);figure;imshow(im1, res,0,1);课程设计总结首先，由于本方法是基于单幅图像输入来进展去雾处理，因此无法像多幅图像条件下那样综合同一场景的不同降质条件来估算天空亮度，因此难免会存在偏差。而且白平衡过程是对整幅图像进展的，这样对于简单场景天空与景物有明显的分界来说误差较小，但对于存在很多细节的图像如茂盛的树木，天空也景物难以区分以及场景中有很多白色景物时，会在边缘出现模糊。其次，本方法使用的是变形的中值滤波器来估算大气耗散函数，而中值滤波器虽然能较好

10、的保持边缘，但对于直角以及线这样的细节的处理存在固有缺陷，因此当场景中存在这样的景物如高楼、车辆等，处理后的图像会在景物的边缘出现明显的模糊。设计体会在本次课程设计过程中，收获颇丰。首先，这次实践让我又重新温习了一遍数字图像处理课程的相关知识，并在此过程中对所学知识进展了一次检验，加深了度数字图像处理这门课程的理解，让我深感数字图像处理这门课程的博大精深。然后，通过本次课程设计的体验，提高了综合运用数字图像处理各方面知识的能力，并从中发现了自己的许多缺乏之处，有待提高。最后，通过对matlab这个工具软件的运用，我又掌握了它的一些功能，这对我以后的学习和工作是一个莫大的帮助。参考文献1数字图像处理第二版 冈萨雷斯 电子工业2基于物理模型的单幅图像快速去雾方法 禹晶、大鹏、廖庆敏 自动化学报3Fast Visibility Restoration from a Single Color or Gray Level Image Jean-Philippe Tarel、Nicolas Hautiere IEEE课程设计成绩评定表设计上机验收成绩表*课题名称序号验收工程分值得分1设计