图像处理作业4
圖像處理作業4
1. 第二版課本習題4.21
本質沒有區別,只將圖片放置在中心,而周圍填充0的個數不變時,不會影響結果。因為本質都是進行了周期延拓,使得尾部的信息不會被丟棄掉。相當于濾波前將圖像進行了平移。需要注意的是,濾波后得到的圖像也會發生平移,裁剪的時候會產生區別。
2. 假設我們有一個[0,1]上的均勻分布隨機數發生器U(0,1), 請基于它構造指數分布的隨機數發生器,推導出隨機數生成方程。若我們有一個標準正態分布的隨機數發生器N(0,1),請推導出對數正態分布的隨機數生成方程。
(1) 解答
設指數分布的隨機變量為YYY,概率密度PDF表示為:f(y)=λe?λy;y>0,λ>0,其CDF表示為G(y)f(y) = \lambda e^{-\lambda y};y >0,\lambda >0,其CDF表示為G(y)f(y)=λe?λy;y>0,λ>0,其CDF表示為G(y)。
設均勻分隨機變量為XXX。并且隨機數生成方程Y=g(X)Y = g(X)Y=g(X)。
根據CDF的定義有:
G(y)=P{Y≤y}=P{g(X)≤y}=P{X≤g?1(y)}=g?1(y)G(y)=P\{Y \le y\} = P\{g(X) \le y\} = P\{X \le g^{-1}(y)\} = g^{-1}(y)G(y)=P{Y≤y}=P{g(X)≤y}=P{X≤g?1(y)}=g?1(y)
由此可知:g(y)=G?1(y)g(y)=G^{-1}(y)g(y)=G?1(y),即g(X)=G?1(X)g(X)=G^{-1}(X)g(X)=G?1(X)
指數分布CDF為G(y)=∫0yλe?λxdx=1?e?λyG(y)=\int_0^{y}\lambda e^{-\lambda x}dx = 1-e^{-\lambda y}G(y)=∫0y?λe?λxdx=1?e?λy
G?1(X)=?1λln(1?X)G^{-1}(X) = -\frac{1}{\lambda}ln(1-X)G?1(X)=?λ1?ln(1?X)
因此 Y=G?1(X)=?1λln(1?X);XY = G^{-1}(X)=-\frac{1}{\lambda}ln(1-X);XY=G?1(X)=?λ1?ln(1?X);X服從均勻分布。
(2) 解答
設標準正態分布的隨機變量為XXX,對數正態分布的隨機變量為YYY,即是求Y=g(X)Y=g(X)Y=g(X)。
設標準正態分布的概率密度函數f(x)=12πe?x2/2f(x)=\frac{1}{\sqrt{2\pi}}e^{-x^2/2}f(x)=2π?1?e?x2/2
設對數分布的概率密度是h(y)=12πσe?(lny?μ)22σ2h(y) = \frac{1}{\sqrt{2\pi }\sigma}e^{\frac{-(lny-\mu)^2}{2\sigma^2}}h(y)=2π?σ1?e2σ2?(lny?μ)2?
G(y)=P{Y≤y}=P{g(X)≤y}=P{X≤g?1(y)}=F(g?1(y))G(y) = P\{Y \le y\}=P\{g(X) \le y\}=P\{X \le g^{-1}(y)\}=F(g^{-1}(y))G(y)=P{Y≤y}=P{g(X)≤y}=P{X≤g?1(y)}=F(g?1(y))
∫?ooy12πσe?(lnt?μ)22σ2dt=∫?oog?1(y)12πe?t2/2dt\int_{-oo}^{y}\frac{1}{\sqrt{2\pi }\sigma}e^{\frac{-(lnt-\mu)^2}{2\sigma^2}}dt=\int_{-oo}^{g^{-1}(y)}\frac{1}{\sqrt{2\pi}}e^{-t^2/2}dt∫?ooy?2π?σ1?e2σ2?(lnt?μ)2?dt=∫?oog?1(y)?2π?1?e?t2/2dt
求導可得
e?(lny?μ)22σ2=σe?g?1(y)2/2?g?1′(y)e^{\frac{-(lny-\mu)^2}{2\sigma^2}}=\sigma e^{-{g^{-1}(y)}^2/2}*g^{-1'}(y)e2σ2?(lny?μ)2?=σe?g?1(y)2/2?g?1′(y)
繼而
g(y)=eσy+μg(y) = e^{\sigma y+\mu}g(y)=eσy+μ
g(X)=eσX+μg(X) = e^{\sigma X + \mu}g(X)=eσX+μ,XXX服從均勻分布
3. 請證明第二版課本習題4.5中提及的頻域內高通濾波器與低通濾波器的關系式子。
由于原圖像f(x,y)f(x,y)f(x,y)由兩部分構成,即高頻部分和低頻部分。
那么f(x,y)=fh(x,y)+fl(x,y)f(x,y) = f_h(x,y)+f_l(x,y)f(x,y)=fh?(x,y)+fl?(x,y)
可以認為高頻部分是由高頻濾波器對原圖像進行濾波得到:
fh(x,y)=f(x,y)?Hhp(x,y)f_h(x,y) = f(x,y)*H_{hp}(x,y)fh?(x,y)=f(x,y)?Hhp?(x,y)
fl(x,y)=f(x,y)?Hlp(x,y)f_l(x,y) = f(x,y)*H_{lp}(x,y)fl?(x,y)=f(x,y)?Hlp?(x,y)
f(x,y)=fh(x,y)+fl(x,y)=f(x,y)?Hhp(x,y)+f(x,y)?Hlp(x,y)f(x,y) = f_h(x,y)+f_l(x,y)=f(x,y)*H_{hp}(x,y)+f(x,y)*H_{lp}(x,y)f(x,y)=fh?(x,y)+fl?(x,y)=f(x,y)?Hhp?(x,y)+f(x,y)?Hlp?(x,y)
進行傅立葉變換
F(u,v)=F(u,v)Hhp(u,v)+F(u,v)Hlp(u,v)=F(u,v)(Hhp(u,v)+Hlp(u,v))F(u,v)=F(u,v)H_{hp}(u,v)+F(u,v)H_{lp}(u,v) = F(u,v)(H_{hp}(u,v)+H_{lp}(u,v))F(u,v)=F(u,v)Hhp?(u,v)+F(u,v)Hlp?(u,v)=F(u,v)(Hhp?(u,v)+Hlp?(u,v))
可知
Hhp(u,v)=1?Hlp(u,v)H_{hp}(u,v)=1-H_{lp}(u,v)Hhp?(u,v)=1?Hlp?(u,v)
4.給出的逆諧波濾波回答下列問題
(a)解釋為什么當Q是正值時濾波對去除“胡椒”噪聲有效?
Q>0Q>0Q>0時對像素具有增強作用。而對于胡椒噪聲來說,其灰度值較小,所以增強較小,所以對加權平均的結果影響較少。
(b)解釋為什么當Q是負值時濾波對去除“鹽”噪聲有效?
Q<0Q<0Q<0時對像素具有削弱作用。而對于胡椒噪聲來說,其灰度值較大,所以削弱較大,所以對加權平均的結果影響較少。
總結
- 上一篇: 图像处理作业第7次
- 下一篇: 如何给视频加水印标记怎么给视频加水印