{finalbug}

3. 空间滤波
(a)噪点生成器。1. 找到（或开发）一段程序能在一个图像中加入高斯噪点,并且必须能够指出噪点的均值和方差。2.找到（或开发）一段程序能在一个图像中加入黑白相间的（脉冲式）噪点,并且必须能指出这两种噪点中每一种的出
现概率。
(b)写一段程序实现空间均值滤波，把你的滤波器应用到在(a)中含噪点的图像上。你可以把你的空间遮罩的大小固定设成3X3，但是系数要是可变的并作为你程序的输入。
(c)修改你在(b)中生成的程序以实现一个3X3 的中值滤波。比较用均值滤波和中值滤波过滤后的图像的不同之处。

结果图：

先把高斯噪声的伪随机数生成方法记下来。其他的以后有空再写。

/**
* 用雅可比变换生成一个呈高斯分布的伪随机数。
*
* @param mu 平均值μ
* @param sigma 标准差σ^2
* @return 伪随机数。
*/
private function getGaussianRandom(mu:Number = 0, sigma:Number = 1):Number
{
	var r1:Number = Math.random();
	var r2:Number = Math.random();
	var rs:Number = Math.sqrt(-2 * Math.log(r1)) * Math.cos(2 * Math.PI * r2 ) * sigma + mu;
	return rs;
}

还有整整一个月的时间就要交《数字图像处理》的作业了，一直有任何的头绪。本来是打算有了现成的程序之后自己就写点论文就好了，但是看来大家都没有太大的兴致来写。于是横下一条心，干脆自己来做吧！虽然不专业，对位图处理还是有点经验的。开始向用freemat或者opencv之类的东西来做，但是毕竟不是日常使用的东西，不熟悉，看了看作业的要求，似乎只要是能转换成位图矩阵就能使用，所以还是用Flex做吧，发布成AIR，但愿到时候老师知道我的是什么东西。。

这东西只能断断续续的做，还是将自己做的东西记录下来，免得忘掉了。

首要的问题，在AS中似乎没有直接获取灰度值的方法，因此我用的是将RGB值分解之后取平均值。顺便记录一下分解RGB的简单算法：

var red:uint = color >> 16 & 0xFF;
 var green:uint = color >> 8 & 0xFF;
 var blue:uint = color & 0xFF;

如果不将RGB的值转换成一个灰度值，可以分别对三个值进行计算，得到三个直方图，让后用同样的算法对三个值进行计算，之后再组成一个RGB值，可以作为整个图形的计算。

第一题，直方图均衡化

• 写一段程序计算一个图像的直方图
• 实现直方图均衡化
• 你的程序必须大体上能允许所有的灰阶图像作为输入

至少，在你的报告中需要包括原始图像和它的直方图，一张直方图均衡化的图，转换函数，改进后的图像和他的直方图。

实现效果：

按照我自己的理解，直方图的算法如下。未必正确，仅供参考。

// 以下代码仅描述算法，不能直接运行。
// 第一步，获取该位图从0-255的每个灰度值的点的个数。
for(i = 1 ; i < 255 ; i++)
{
	for(j = 1 ; j < 255 ; j++)
	{
		// grays数组的index为灰度值，值为该灰度值的点的个数
		grays[getGrayAt(i, j)]++;
	}
}

// 第二步，获得归一化之后的直方图
for(pi = 0 ; pi < grays.length ; pi++)
{
	p[pi] = grayList[pi] / (imgWidth * imgHeight); // p数组中保存的是每个灰度的百分比
}
// 第三步，累计的归一化直方图
for(ci = 0 ; ci < p.length ; ci++)
{
	c[ci] = 0;
	for(cj = 0 ; cj < ci ; cj++)
	{
		c[ci] += p[cj]; // 数组c中保存的是归一化之后的直方图
	}
}
//
// 最后，转换成新的位图
for(ti = 1 ; ti < 255 ; ti++)
{
	for(tj = 1 ; tj < 255 ; tj++)
	{
		t[ti][tj] = c[getGray(originaImage.getGrayAt(ti, tj))] * (maxGray - minGray) + minGray;
	}
}
// maxGray是原图中最高的灰度值，minGray是最低的。
// 数组t即为新的图的像素矩阵。

补充几个AS的方法：

public static function splitColor(color:uint):Array
{
	var red:uint = color >> 16 & 0xFF;
	var green:uint = color >> 8 & 0xFF;
	var blue:uint = color & 0xFF;
	return [red, green, blue];
}

public static function getGray(color:uint):uint
{
	var RGB:Array = splitColor(color);
	return uint((RGB[0] + RGB[1] + RGB[2]) / 3);
}

public static function getColor(gray:uint):uint
{
	return gray * 65536 + gray * 256 + gray;
}