实验2分析和总结.docxVIP

现代光电子技术实验报告

实验二 红外热像仪噪声统计特性分析

专 业： 电子与通信工程

姓

名：

何

超

学

号：

1405122310

组

号：

6

—、实验目的

了解红外热像仪的工作原理、操作步骤及相关的图像噪声处理方法；掌握红外热像仪的正确使用方法，独立完成相关的图像处理并通过拍摄于不同场景的红外图像分析红外热像仪的噪声统计特性。

二、实验步骤

拍摄不同场景的红外图像分析红外热像仪的噪声统计特性

拍摄红外场景

选取均匀背景拍摄均匀背景红外图像；

选取普通场景拍摄普通场景红外图像。

针对以上场景，计算噪声图像

去噪得到理想图像

均值滤波获取理想图像，滤波窗口取3?3或5?5；

改进均值滤波获取理想图像，滤波过程中首先去掉滤波窗口中的最大值点和最小值点，然后进行滤波，滤波窗口取3?3或5?5；

利用中值滤波法，获取理想图像，滤波窗口取3?3或5?5；

改进中值滤波法，获取理想图像，选取中值滤波法的中间三个值的平均作为输出，滤波窗口取3?3或5?5；

利用原始图像减去理想图像得到噪声图像。

噪声统计特性

要求计算噪声图像的直方图，均值，方差，丰度，偏度并拟合直方图曲线。三、基本原理

直方图

设n(r)为图中灰度为r的象素数（i?[0,k?1]），图像的大小为M行、N列，则r出

i

现的概率为：

i

n(r)

P(r)???i

i

(1)

r i MN

且

?k?1P(r)?1

r i

(2)

i?0

均值

设（X,X

1 2

,X ）是总体X容量为n的样本，则

n

称为样本均值（Samplemean）。

方差

X?1?n X

n i

i?1

(3)

设（X1,X2 ,Xn）是总体X容量为n的样本，则

S2?

1

n?1

?n(X

i

i?1

?X)2

(4)

称为样本方差（Samplevariance），

偏度

随机变量x的偏度定义为

? ?E?x?E?x?

? ?

??Var

??Var?x???32

1

(5)

它度量了分布的偏斜程度及偏向，是一个无量纲的数值。若g

1

?0，则称x的分布是正偏（或

右偏）的；若g

1

?0,则称x的分布是负偏（或左偏）的。g

1

越大，说明分布偏斜得越厉

害。

设x, ,x

是来自总体x的一个样本，则总体x的偏度可估计为

1 n

g???

?n(x

?x)3

(6)

n1 n?1??n?2?s3

n

i

i?1

其中x?1?nx

n i

i?1

是样本均值，s?

1n

1

n?1

?n

i?1

(x?x)2

i

峰度

峰度是另一个反映随机变量分布形状的量，随机变量x的峰度定义为

E??x?E?x???4

g ?

2 ??

Var?x

???2

?3 （3）

它度量了分布尾部的厚度。同偏度一样，峰度也是一个没有量纲的数值。峰度g2的取值范

围是［-2,∞］。正态分布的峰度为零。人们以正态分布为标准，若g

2

?0，则说明随机变

量X分布的尾部比正态分布的尾部粗，并且g

2

值越大，倾向认为尾部越粗；若g

2

?0，则

说明X分布的尾部比正态分布的尾部细，且g 值越大，倾向认为尾部越细。峰度g 可用

2 2

来比较已标准化了的各随机变量分布的尾部厚度。

设x, ,x

是一组样本数据或一组有限总体数据，则其峰度的计算公式为

1 n

n?n?1? ?n

?n?1?2

g? ??

?1??

?2??

?3?

(x?x)4?3?

?2??

?3?

（4）

2 n n

n s4 i n n

i?1

其中x和s的含义同前。

四、核心代码

在第6组拍摄的图像中，选取IMAGE060.IMG作为均匀背景拍摄均匀背景红外图像；

选取IMAGE074.MG作为普通场景拍摄普通场景红外图像，分析红外热像仪的噪声统计特性。用高德热像仪分别转换为1.raw和2.raw，进行处理。

代码如下：clearclc

%%读取均匀背景图像

F=fopen(E:\1.raw,r);

F=fread(F,[320,240],int16=int16);F=mat2gray(F);

%%均值滤波

A=fspecial(average,[33]);F_average=imfilter(F,A,replicate);

figure,subplot(4,3,1),imshow(F_average)