数字图像处理第三版中文答案冈萨雷斯.doc
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第二章
(第二版是和* 的矩形,第三版是和圆形)
对应点的视网膜图像的直径x 可通过如下图题所示的相似三角形几何关系得到,即
d 2 x 2
0.30.017
解得x=。
根据节内容,我们知道:如果把中央凹处想象为一个有337000 个成像单元的圆形传感器阵列,它转换成一个大小327.52 成像单元的阵列。
假设成像单元之间的间距相
等,这表明在总长为 1.5 mm(直径)的一条线上有655 个成像单元和654 个成像单元间隔。
则每个成像单元和成像单元间隔的大小为s=[(1.5 mm)/1309]=× 10-6 m。
如果在中央凹处的成像点的大小是小于一个可分辨的成像单元,在我们可以认为改点对于眼
睛来说不可见。
换句话说,眼睛不能检测到以下直径的点:
x 0.06d 1.1 10 6 m ,即 d 18.3 10 6 m
当我们在白天进入一家黑暗剧场时,在能看清并找到空座时要用一段时间适应。
节描述的
视觉过程在这种情况下起什么作用
亮度适应。
虽然图中未显示,但交流电的却是电磁波谱的一部分。
美国的商用交流电频率是 77HZ。
问这一波谱分量的波长是多少
光速 c=300000km/s ,频率为 77Hz。
因此λ =c/v= * 10 8(m/s)/77(1/s) = *10 6m = 3894 Km.
根据图得:设摄像机能看到物体的长度为x (mm),则有 :500/x=35/14; 解得: x=200 ,所以相机的分辨率为: 2048/200=10; 所以能解析的线对为:10/2=5 线对 /mm.
假设中心在( x0,y0 )的平坦区域被一个强度分布为:
i (x, y) Ke [( x x 0) 2 ( y y 0) 2 ] 的光源照射。
为简单起见,假设区域的反射是恒定的,并等于,令 K=255。
如果图像用 k 比特的强度分辨率进行数字化,并且眼睛可检测相邻像
素间 8 种灰度的突变,那么 k 取什么值将导致可见的伪轮廓解:题中的图像是由:
f x, y i x, y r x, y 255e x x02 y y0 2
255e
x x0 2 y y0 2
1.0
一个截面图像见图(a)。
如果图像使用 k 比特的强度分辨率,然后我们有情况见图(b),其中G 255 1 2k。
因为眼睛可检测 4 种灰度突变,因此,G 4 256 2k ,K= 6。
也就是说,2k 小于 64 的话,会出现可见的伪轮廓。
(a)传输数据包 ( 包括起始比特和终止比特 ) 为: N=n+m=10bits 。
对于一幅 2048× 2048 大小
的图像,其总的数据量为 M 2048 2 N ,故以56K 波特的速率传输所需时间为:
T M 56000 2048 2 8 2 56000 748.98s 12.48 min
(b)以 3000K 波特的速率传输所需时间为
T M 3000000 2048 28 2 3000000 13.98s
解:图像宽高比为
=1920 像素 / 线。
16:9 ,且水平电视线的条数是1080 条,则:竖直电视线为1080 ×(16/9 )
由题意可知每场用则该系统每 1/30 的图像。
又因为
1s 的 1/60 ,则:每帧用时2× 1/60=1/30秒。
秒的时间形成一幅1920× 1080 分辨率的红、绿、蓝每个像素都有
90min 为 5400 秒,故储存90min 的电视节目所需的空间是:
8 比特
1080 1920 8 3 30 5400 8.0621012
bits
1.001 1012 bytes 解: p和 q如图所示:
(a)S1和 S2不是4邻接,因为q不在 N 4 p 集中。
(b)S1和 S2是8连接,因为q在 N8 p 集。
(c) S 和 S 是
m 连接,因为 q 在集合
N D p
中,且 N
4
p N
4
q 没有V值的像素。
1 2
提出将一个像素宽度的8 通路转换为 4 通路的一种算法。
解:找出一个像素点的所有邻接情况,将对角元素转化成相应的四邻接元素。
如下图所示:
提出将一个像素宽度的m通路转换为 4 通路的一种算法。
解:把 m 通道转换成 4 通道仅仅只需要将对角线通道转换成 4 通道,由于 m 通道是 8 通道与 4 通道的混合通道, 4 通道的转换不变,将8 通道转换成 4 通道即可。
如图所示:
(1)4 邻域关系不变
(2)8 领域关系变换如下图所示
(没答案,自己做的,看对不对)
(1)在V={0,1,2
距离为 6。
}时, p 和q 之间通路的D4距离为8(两种情况均为8), D8距离为4, D m
(2) 在 V={ 2,3,4 }时, p 和 q 之间通路的D4距离为∞, D8距离为p 和 q 之间不存在 4 邻接路径,因为不同时存在从p 到 q 像素的4, D m距离为 5。
4 毗邻像素和具备V 的
值,情况如图 (a) 所示。
p不能到达q。
解:
(a) 点 p(x , y)和点q(s , t) 两点之间最短 4 通路如下图所示,其中假设所有点沿路径V。
路径段长度分别为x s
和
y t ,由D4 距离的定义可知,通路总长度| X-S|+| Y-T| ,(这
个距离是独立于任何点之间可能存在的任何路径),显然 D4距离是等于这两点间的最短4
通路。
所以当路径的长度是x s y t
,满足这种情况。
(b)路径可能未必惟一的,取决于V 和沿途的点值。
由公式H [f(x,y)]=g(x,y),
让 H表示相邻的和操
作
, 让S1和 S2表示两个不同子图像区的小值, 并让S1 + S2表示相应的
总数S1和 S2像素,如
在
2.5.4 节里的解释. 注意到附近的大小( 即像素数字) 并没有随着这
总和的改变而改变。
H计算像素值是一个给定的区域。
然
后
,
H aS1 bS2
意味着 :
(1)在每个子区域里乘像素 ,
(2) 从 aS1到 bS2每个像素值相加(首先产生一个单独的子区域)
(3) 在单独的子图像区域里计算所有像素值的和。
让ap1和ap2表示两个任意(但相应的)像
素 aS1bS2。
然后我们可以依据Eq. - 1),表明H是一个线性算子。
(两个版本答案,一个意思)
(1)中值ζ表示,数集的一半数值比它大,另一半比它
小。
一个简单的例子能够表明 ,Eq. - 1) 的平均算子操作。
让 S1 = {1,-2,3}, S2 = {4,5,6}, a = b = 1.在这种情况下,H是平均算子。
然后有 H(S1 + S2)= 中值 { 5,3,9 } = 5,S1 + S2是S1和S2的和。
接下来 , 计算 H(S1)=中值 { 1 、-2 、3 } =1和H(S2)=中值{ 4、5、6 } = 5。
然后 , 从 H(aS1 + bS2) ≠ aH(S1)+ bH(S2),因此,子图像区域S 中值的算子是非线性的。
(2)
因为 g x , y f x, y x, y g x ,y 1 K
g i(x , y ) K i 1
E g x, y
1 K
g i( x, y)
1 K
x, y i x, y E E f i
K i 1 K i 1
1 E K 1 E K
f i x,y i x ,y f x, y
K i 1 K i 1
2 g x , y 2 1 K
i
K i g (x, y )
1
1 2 K
f i x, y
K 2 i 1
(没答案看看做的对不对)
(a)为 A的补集
(b) A B C
1 2 K
f i x, y i x, y K 2 i 1
1 K 1
2 2 K
2 i
x , y
i 1 K
A B B C A C 2A B C
A C
B A B B C
(看看翻的对不对)
答:使用三角区即三个约束点,所以我们可以解决以下的系数为 6 的线性方程组:
x c1x c2 y c3
y c4 x c5 y c6
实施空间变换。
插值强度可使用 2.4.4 节的方法。
(看看翻的对不对)
傅里叶变换核是可分的,因为:
r x, y,u,v e j 2 ux / M vy / N e j 2 ux / M e j 2 vy / N r1 x,u r2 y ,v
傅里叶变换核是对称的,因为:
e j 2 ux / M vy / N e j 2 ux / M e j 2 vy / N r1 x,u r1 y,v
(看看翻的对不对)
由可分离变换核的定义知其中:
当 x 值固定时,可看作 f(x,y) 某一行的一维变换,当 x 从 0 变换到 M-1 时计算出整个数组 T (x,v ),然后,通过替换这个数组的最后一行以前的方程我们可以得到T( x,v )按列的一维变换。
也就是说,当一个图像是内核可分的,我们可以计算图像沿行的一维变换,然后我
们计算中间的一列得到最终的二维变换T(u,v). 这和先计算列的一维变换再计算中间行得
到二维变换最终结果是相同的。
从式(),二维傅里叶变换是由:
它很容易验证,傅立叶变换核是可分离的(参见题),所以我们可以写这个方程:
是沿着 f(x,y)行的一维傅里叶变换,X= 0 , 1,, M-1。
第三章
( a)由 s T (r ) Ae Kr 2 ,Ae KL 02 A/3得: KL02 ln(1/ 3) , K 1.0986 / L20
1.0986r2
s T (r ) Ae L20
(b)、由,B(1 e KL 02 ) B/4得:
KL 20 ln( 3/ 4) , K 0.2877 / L 02
0.2877r 2
s T (r ) B(1 e L 0
2 )
(c )、
逐次查找像素值,如( x ,y )=(0,0)点的 f ( x ,y )值。
若该灰度值的 4 比特的第 0
位是 1,则该位置的灰度值全部置 1,变为 15;否则全部置 0,变为 0。
因此第 7 位平
面 [0,7] 置 0, [7,15] 置 1,第 6 位平面 [0,3]
, [4,7] 置 0, [8,11] , [12,15] 置 15。
依
次对图像的全部像素进行操作得到第
0 位平面, 若是第 i 位平面,则该位置的第
i 位值
是 0 还是 1,若是 1,则全置 1,变为 15,若是 0,则全置 0 设像素的总数为 n ,是输入图像的强度值,由, rk 对应 sk ,所以,由
和得由此得知,第
二次直方图均衡化处理的结果与第一次直方图均衡化处理的结果相同, 这里我们假设忽略不
计四舍五入的误差。
z 4w
0 w 0.5
v G ( z)
p z (w)dw , p z (w)
4 4 w
0.5 w 1
v G ( z)
令 s v 得
z 2 z 2
0 z 0.5 p z (w) dw
1 2 z 2 z
2
0 .5 z 1
v r 2
2 r
0.5
所以 z
G 1
(v)
2
2 0 r 1
2v 1 2( r 2
2 r ) 1
1 0.5 r 1
0.5
0.5
2
2
第 k 个点邻域内的局部增强直方图的值为:
P (r k )=n /n (k=0,1,2, K-1) 。
这里 n 是灰度级为 r
的像素个数, n 是邻域内像素的总个
r
k
k
k
数, k 是图像中可能的灰度级总数。
假设此邻域从左以一个像素为步长向右移动。
这样最左 面的列将被删除的同时在后面又产生一个新的列。
变化后的直方图则变成: (k=0,1,2,
K-1)
这里 n lk 是灰度级 r k 在左面的列出现的次数, n rk 则为在右面出现的次数。
上式也可以改写成: (k=0,1,2,
K-1)
同样的方法也适用于其他邻域的移动: 这里 a k 是灰度级 r k 在邻域内在移动中被删除的像素数,
b k 则是在移动中引入的像素数:
(k=0,1,2,
K-1)
上式等号右边的第一项为
0(因为 f 中的元素均为常数) 。
变量是噪声的简单抽样,它的方
2 差是。
因此 并且我们可以得到。
上述过程证明了式
g ( x, y)
1
K
2
( x, y)
的有效性。
(A )中值是
[( n 2 1) / 2] 的最大值
(B )一旦中值被找出,我们简单的删除邻域边缘的值,在合适的位置插入合适的值
旋转前坐标的拉普拉斯定义为
2
f
2
f
2
f
x 2
y 2 ,旋转后坐标的拉普拉斯定义为
2
f
2
f
2
f x ,
cos
y , sin 和 y ,
sin
,
cos
x '2
y '2 ,现在给出 x
x
y ,其中 指轴
旋转的角度,若想证明拉普拉斯变换是各向同性的,只需证明
2
f
2
f
2
f
2
f ,
x 2
y 2
x '2
y ' 2
首先, f
f x
f y f cos f
sin
x '
x x ,
y x ,
x
y
两边对 x ' 求导得,
2
f
2
f
2
( f
cos
( f
sin
2
f
sin 2
( 1)
x ,2
x 2 cos
) sin
) cos y 2
x y
y
x
同理可得,
f
f x f y f sin f cos
y ' x y ,
y y ,
x
y
两边对 y , 求导得,
2
f
2
f
2
f
f
2
f
2
( 2)
y ,2
x 2
cos
x
(
y ) sin
cos
y
(
x ) cos sin
y 2 sin
(1)和( 2)式相加得,
2
f
2
f 2
f
2 f ,所以拉普拉斯变换是各向同性的。
x 2
y 2
x ' 2
y ' 2
使用式()给出的拉普拉斯定义,证明从一幅图像中减去相应的拉普拉斯图像等同于对图 像进行非锐化模板处理。
2
f
[ f ( x 1, y) f ( x 1, y) f ( x, y 1)
f ( x, y 1)] 4 f ( x, y)
(3.6.6 )
考虑到下列公式
_
其中 f ( x, y) 是 f ( x, y) 预先确定的临域的平均数,更确切的说就是以
( x, y) 为中心并且包
括中心像素以及四个相邻像素。
把上面的等式的最后一行的常量视为均衡因子(或比例因子),我们可以写出
2
f ( x, y)
_
f (x, y)
f ( x, y) f ( x, y)
_
等式的右端就是等式
f s ( x, y) f ( x, y ) f ( x, y) 给出的非锐化掩膜处理的定义。
因此验证
了从一幅图像中间取相应的拉普拉斯图像等同于对图像做非锐化掩膜处理。
题
f mag( f )
[G x
2
G y 2 ]1 / 2
[( f
) 2
( f )2 ]1/ 2
( 3.6.11 )
x
y
f |G x | | G y |
( 3.6.12 )
(a )由
f f
cos
f
sin 和
f
f
sin
f
cos
x '
x y
y '
x y
2
f
2
f
2
f
2
f
2
f
2
f 1 / 2 (
2
f
2 2
' 2
'2或(
2
2 ) '2
x
y
x y x
y
x
因此 , 我们看到的梯度向量的模值是一种各向同性梯度算子
(b )从上面的结果得 | G x | | f |, | G y | |
f
|
x
y
| G x' | | f ' | | f
cos
f sin |, | G y ' | | f ' | | f sin
x x
y
y x
显然得到 |G x '
| |G y ' | | G x | |G y |
2
'f 2 )
1/ 2
y
f
cos |
重复例,但是用函数 f (t) 2A( W / 4 W / 4) 和 f (t ) 0 ,对于其他所有的t
值。
对你的结果和例子中的结果之间的任何不同,解释原因。
解:
F f t e j 2 t dt
W
4
W 2 Ae j 2t dt
4
A j e
W
4
j 2 t
W
4
A j W j W
e
j e 2 2
A j W j W
e 2 e 2
j
Q sin e j e j
2 j
F
2A
sin W
2
W
sin
2
AW
W
2
傅立叶变换的幅值是不变的;由于周期不同,
~ ~
F f t e j 2 t dt
f t t n T e j 2 t dt
n
证明式() f t t n T e j 2 t dt 中的F 在
~
n
f n e j 2 n Tt
n
两个方向上是无限周期的,周期为1/T
证明:
(1) 要证明两个方向上是无限周期1/ T ,只需证明
根据如下式子:
可得:
其中上式第三行,由于k, n是整数,且和的极限是关于原点对称。
(2)同样的需要证明
根据如下式子:
~ ~
F f t e j 2 t dt
f t t n T e j 2 t dt
n
f t t n T e j 2 t dt
n
f n e j 2 n Tt
n
可得:
j 2 kn
其中第三行由于k, n都为整数,所以e1。
可以证明( Brancewell[2000])1(t)和1(t) 。
使用前一个性质和表中的平移性质,
证明连续函数 f (t ) cos(2 nt ) 的傅立叶变换
是 F
1/ 2 n n ,其中是一个实数。
证明:
根据一维傅里叶变换公式:
可得:
F(u)
f (t )e
j 2 ut
dt
cos( 2
) j 2 ut
dt
nt e
1
[ e j 2
nt
e j 2
nt
]
j 2 ut
dt
2
e
1
e j 2 nt e j 2
ut dt
1
e j 2 nt e j 2 ut dt
2
2
根据傅里叶变换性质可得:
根据一个常数 f(t)=1
的傅里叶变换是一个脉冲响应可得:
所以可得如下两个等式:
(1)e j 2 nt ( n) (1)e - j 2 nt
( +n)
所以:
F(u)
1 (u n ) (u n )
2
考虑连续函数 f (t ) cos(2
nt )
(a)
f t 的周期是多少 (b) f t 的频率是多少
(a) 根据 2 nt 2 ,所以周期为 t 1/ n
(b) 频率为 n ,给定的正弦波的连续傅立叶变换如在图。
(a )(见习题),采样数据(示出
了几个期间)的变换所示的一般形式的如图( b )(虚线框是一个理想的过滤器,将允许重建 如果该正弦函数进行采样,采样定理满意)。
解:
(a) 根据正交性,将式直接代入式得
最后一步是根据问题的陈述中给出的正交条件,将式代入式应用同样的过程生成
f n 的相似
特性。
(b) 如上小题,根据正交性,将式直接代入式得
最后一步是根据问题的陈述中给出的正交条件, 将式代入式应用同样的过程生成
f x 的相
似特性。
证明式 F u
kM
F u 和式 f
x kM
f x 的正确性。
证明:
(1) 证明等式 F
u kM
F u
k 0, 1,
2...
M 1
将 u u kM 代入式 F u
f ( x)e j 2
ux / M
, u 0,1,2,K ,M 1 :
n 0
M 1
F u kM
f ( x)e j 2 ( u kM ) x/ M
n 0
M 1
e j 2 kx
f ( x)e j 2 ux/ M
n 0
F(u)
最后一步因为 k 和 x 都是整数, e j 2
kx
1。
(2) 同理可以对式周期性的证明,将
u u kM 代入式
f x
1 M 1
F (u)e j 2 ux / M
, u 0,1,2,K , M 1
M n
f x kM 1 M 1 u e j 2 ( u kM ) x/ M
M n
F
1 M 1
u e j 2 ux/ M
e j 2 kx
M
F
n 0
=f x
证明一个变量的离散卷积定理的正确性
[见式、式 和式 ] 。
证明:
证明卷积定理等价于证明
f ( x) h( x) F(u)H(u)
和
f ( x)h( x) F(u) H(u)
M 1 M 1
从式 f ( x) h( x) f (m) h( x m) 和式 F u f ( x) e j 2 ux / M , u 0,1,2,K , M 1 m 0 n 0
离散傅里叶变换的定义,得到:
M 1 M 1
m) e j 2 ux/ M
f ( x) h( x) f (m)h( x
x 0 m 0
M 1 M 1
j 2 ux/ M
f (m) h(x m)e
m 0 x 0
M 1
f (m) H(u) e j 2 um/ M
m 0
M 1
f (m) e j 2 um/ M
H(u)
m 0
H(u) F(u)
同理可以证明 f ( x) h(x) F(u) H(u)
1 F u H u M 1 M 1
F (m)H (t m) e j 2 ux/ M x 0 m 0
M 1 M 1
m)e j 2 ux/ M
F ( m) H (t
m 0 x 0
M 1
F ( m) h( x)e j 2 um/ M
m 0
M 1
h( x) f (m) e j 2 um/ M
m 0
h x f x
写出二维连续卷积的表达式
对式进行卷积运算得到:
f(t,z) h(t,z) f ( , )h(t, z)d d
证明一维连续和离散傅里叶变换都是线性操作
解:
若连续傅里叶变换
是线性的,只需证明:
代入傅立叶变换定义
其中第二步由于积分的分配率。
同样的,离散傅里叶变换:
证明连续和离散傅里叶变换都是平移和旋转不变的。
证明:
平移不变:根据二维离散傅立叶变换
可得
旋转不变:根据二维离散傅立叶反变换
证明离散函数 f x, y cos 2 u 0 x 2 v 0 y 的DFT 是
F u, v
1
u Mu 0 , v Nv 0
u Mu 0 , v Nv 0
2
证明: 根据欧拉公式
M 1 N 1
v y e j 2
ux/ M vy/ N
F u,v
cos 2 u 0 x 2 x 0 y 0
1M1N1
e j 2 u 0 x v 0 y
e j 2 u 0 x v 0 y e j 2 ux/ M vy/ N
2 x 0 y 0
1
M 1 N 1
e j 2 Mu 0 x/ M Nv 0 y / N e j 2 ux/ M vy/ N
2 x 0 y 0
1M1N1
j 2 Mu 0 x / M Nv 0 y/ N
e j 2 ux/ M vy / N
e
2 x 0 y 0
1
1 e
j 2Mu 0
x / M Nv 0
y/ N
1 1 e
j 2 Mu 0
x/ M Nv 0
y / N
2
2
1 u Mu 0 , v Nv
u Mu 0, v Nv 0
2
其中最后一步由于
1 u, v ,根据 DFT 平移性
1 e j
2 u 0x/ M v0 y / N u u0 , v v0。
找出一个等价的滤波器H u, v ,在他的频率域实现使用图 (a) 中拉普拉斯模版
执行的空间操作。
解:
滤波后的函数为
g x, y f x 1, y f x 1, y f x, y 1 f x, y 1 4 f x, y
又因为 G u, v H u,v F u,v ,其中
将滤波器变换为频率中心对称
当 u,v 2 M / 2, N / 2 (变换后滤波器中心) 时,H u, v0 。
对于远离中心的值,H u, v 降低。
重要的一点这是一个高通滤波器的特性,消除了直流分量,留下了高频分量。
解:
共轭复数只是从j 变成了 -j在逆变换中,所以右边的图像可以通过下述过程求出:
a
x y
1 f x, y
M 1 N 1
y
b F u, v
x
x, y e j 2 ux/ M vy/ N
1 f
x 0 y 0
M 1 N 1
y
c F u, v
x
x, y e j 2 ux/ M vy / N
1 f
x 0 y 0
d F 1 F u,v 1
MN
实部为
x y
f x, y
1
M 1N 1 M 1N1
x y
e j 2 ux / M vy/ N
1 f x, y e j
2 ux / M vy/ N
u 0 v 0 x 0 y 0
e 结果为 1 x y 1 x y
f x, y =f x, y
可以知道整个过程只是将 f x, y 上下左右颠倒,从而产生了右边的图像
解:
(a)以卷积的形式给出滤波表达式,来减少空间域的处理过程。
然后滤波后的图像由下式给
出:
其中 h 是空间滤波函数, f 是输入图像。
直方图处理结果为:
T表示直方图均衡化。
如果先进行直方图均衡化,再
与
总体来说, T 是由图像像素的属性决定的非线性的函数。
因此,
,并且先后顺序是有影响的。
(b)正如在第节,高通滤波严重削弱了图像的对比度。
虽然高频率的改进一些,但并不显著(见图)。
因此,如果对一个图像先直方图均衡化,均衡化中对对比度的改进会在滤波过程
中严重损失。
因此,该过程一般是先滤波再直方图均衡化。
证明:
因为,我们可以写出等式和,分别为
与
用归纳法证明开始显示两个方程对于n = 1 成立;
1
2 1 1 与 a 12 1 2
m 1
2
n 成立,那么可以我们从进行讨论的部分中知道这些结果是正确的,然后我们假定方程对于
得出方程对于 n+1 也成立。
从等式中,
将 m(n) 从上式替换得到 ,
因此,等式对所有的n 都成立。
从等式中,
将 a(n) 从上式替换得到,
则证明了等式成立。
第五章
给出与表中带阻滤波器对应的高斯和巴特沃斯带通滤波器的公式。
一个带通滤波通过从相应的带阻滤波而获得:
然后:
(a)理想带通滤波:
(b)巴特带通滤波:
(c)高斯带通滤波:
以式()的形式给出高斯、巴特沃斯和理想陷波带阻滤波器的公式。
带阻滤波器公式可以通过带通滤波器的公式得到。
两者的和为 1.
H br (u, v) 1H bp (u, v)
(a) 理想陷波带阻滤波:
D (u, v)
D 或 D ( u , v )
D
1 0,
2 0
H (u, v)
1
其他
(b )巴特沃斯带阻滤波:
1-
巴特沃斯带通
巴特带通滤波: (c )高斯带阻滤波:
1-
高斯带通滤波
高斯带通滤波:
二维连续余弦函数的傅里叶变换
F ( u , v )
f ( x , y )e j 2 ( ux vy )
dxdy
A cos( u 0 x v 0 y )e
j 2
( ux
vy )
dxdy
余弦的变换
cos
1 ( e i e
j )
2
带入得到
F (u, v)
A [e
j ( u 0
x v 0
y )
e
j ( u 0
x v 0
y )
]e
j 2
(ux
vy)
dxdy
A
2
A
j 2 ( u 0 x / 2
v 0 y / 2 )
e j 2 ( ux vy ) dxdy] [ e
j 2 (u 0 x / 2
v 0 y / 2 )
e j 2 ( ux vy) dxdy]
[ e
2
2
这些都是傅里叶变换的功能
并且
结果变换成
A u0
, v v0
)(u
u0
, v
v0
)] 即可
F (u, v) [ (u
2 2 2
2 2
从例子()
即因此
得出
当
这是一个持续的形式, 一个高斯密度方差
或者
减去的整体从无限数量的加上括号里面是1,因此
这是一个模糊的版本的原始图像
解决这一问题的关键是下面的函数
其中,是此函数的拉普拉斯( 对 r 的二次导数 )
那是 ,等于给定的函数。
然后我们知道从式得到函数f(x,y)
因此,我们简化了求高斯函数中的傅里叶变换。
从表格中, 我们从高斯
对可以得到函数的傅里叶变换, 其变换形式是
因此,退化函数的傅里叶变换是
这是一个简单的扩展问题。
它的目的是为了熟悉维纳滤波器的各种条件。
从式 5.8.3 得
其中
然后
从式 5.9.4 得
其中, P(u,v) 是拉普拉斯算子的傅氏变换。
这是至于这个问题
拉斯算子的变换的表达式通过问题中得到的。
然而,对P(u,v) 要求 , 并且不会简化表达式。
, 我们可以合理地解答。
拉普的代替 , 这只会增加滤波器的
5.24
因为这个系统是假定的线性和位置不变, 因此可以用式子 5.5.17 。
举行。
此外 , 我们可以用叠加问题 , 得到了系统响应的F(u,v) 和 N(u,v) 。
两个响应的和是完整的响应。
首先,仅用F(u,v)
然后,仅仅用N(u,v)
所以
第六章
给出用于产生图中标为“日光”的点的红光、绿光、蓝光的百分比。
从图中可知,x=, y=,由 x+y+z=1 可得 z=,这是三色值系数。
我们感兴趣的是三色值XYZ。
由他们的变换公式:x = X/(X+Y+Z)
相同的,故可得:X=, Y=, Y=
, y=Y/(X/Y/Z) ,z=Z/(X/Y/Z) ,可知他们的比例是
c1
用 c
和 c2
表示给定的颜色,并且给出它的坐标,用
(的距离分别为:
x0,y0 )表示, c 和c1 之间的距离以及
c1 占c 的百分比表示为:
c2 的百分比用p2 表示: p2=100-p1, 由上面的等式我们知道,作为例子,当c=c1 时,那么d(c,c1)=0, 并且 p1=100%,p2=0%,同样当 d(c,c1)=d(c1,c2)时,p1=0%,p2=100%,从它们简单的关系中可以容易地得出它们的值。
在中心点有R/2+ B/2+G= R+G+B /2 + G /2=midgray+G/2
使人们看起来像纯绿色。
,由于增加了灰色分量和强度
在每幅 12 比特图像中有
212 4096 种可能值。
对于灰度色彩,所有的RGB分量必须相等,所以有4096 种不同的灰度。
(a) R 图像中的所有像素值都是255。
在 G图像中,第一列全是0,第二列全是 1,最后一列全由 255 组成。
在 B 图像中,第一行全为 255,第二行全为254,直到最后一行全为 0。
(b)(令坐标轴编号同书中图(RGB彩色立方体示意图)相同。
)则:(0,0,0 ) =白色,(1,1,1 ) =黑色,( 1,0,0 )=青色,( 1,1,0 )=蓝色,( 1,0,1 ) =绿色,( 0,1,1 ) =红色,(0,0,1 ) =黄色,( 0,1,0 )=深红色。
从式()的RGB亮度映射函数推导出式()的CMY亮度映射函数。
s i ks i(i=1,2,3)
s i ks i(1 k) (i=1,2,3)
C R
由公式M1G
Y B
可知, CMY图像中的每个分量都是响应RGB图像单一分量的函
数。
C是 R 的函数, M是 G的函数, Y 是 B 的函数。
为清楚起见,我们使用素数标示CMY分量。
有公式
C R (i=1,2,3 )得,
s i ks i( i1,2,3)(对应RGB分量),并且有公式M 1G
Y B 得,对应于 r i 和 s i 的 CMY分量是(用素数表示),
r 1 r
i s i 1 s i
i
从而有,r
i 1 r
i
s 1 s 1 kr 1 k(1 r )
i i i i
因此,s i kr i (1 k )
最纯的红色是FF0000,对应元素( 6,6)
最纯的黄色是FFFF00,对应元素( 1, 6)
【没有答案,个人理解】推导产生一幅彩色图像的补色的CMY变换
一幅 RGB图像的补色变换为:
s i 1 r i (i 1,2,3) ( 对应 RGB分量 ) ,
C R
和 s i的CMY分量(用素数表示)是由 CMY空间定义公式M 1 G 可知,对应于
r i
Y B
r i 1 r i s i 1 s i
得:r
i 1 r
i
s
i
1 s 1 (1 r ) 1 (1 (1 r ))
i i i
因此s 1 r i
证明当 C=I(单位矩阵)时,式简化为式 .
1
D (z, a) z a [( z a)T
( z a )]
2
1
[( z a ) 2 ( z a )2 ( z
B a
B
)2]2
R R G G
1
D(z,a) [(z a)
T C 1(z a)]2
这是一个简单的问题,当 C 为单位矩阵时, C 的逆矩阵也是单位矩阵,所以式()就变成
1
了 D(z,a) [( z a)T (z a)] 2。
括号中的部分被认为是向量(z-a )与其自身的内积,
所以它与式()的右边部分是相等的。
【20、26 题个人翻译,大家参考】。