使用deconvblind估計PSF

2021-09-17 20:56:16 字數 1066 閱讀 9753

f=checkerboard(8);%產生乙個一邊為8個正方形的測試板

psf=fspecial('gaussian',7,10);%用fspecial處理影象的運動模糊,型別為高斯型,濾波器大小為7*10,psf為空間濾波器

subplot(2,2,1),imshow(pixeldup(psf,73),);title('原始psf')

sd=0.01;%sd為標準差

g=imnoise(imfilter(f,psf),'gaussian',0,sd^2);%用psf模糊影象f,在其上新增均值為0,標準差為0.01的高斯雜訊

initpsf=ones(size(psf));%initpsf是點擴散函式的初始值,令initpsf是元素全為1,大小與psf相同的矩陣

numit=5;%迭代次數為5

dampar=10*sd;%dampar指定了結果影象與原影象之間的偏離閾值

lim=ceil(size(psf,1)/2);

weight=zeros(size(g));%產生權重矩陣

weight(lim+1:end-lim,lim+1:end-lim)=1;%不良畫素賦以權重為零,其他畫素都是1

[fr,psfe]=deconvblind(g,initpsf,numit,dampar,weight);%使用deconvblind執行盲目卷積

subplot(2,2,2),imshow(pixeldup(psfe,73),);title('5次迭代估計的psf')

numit=10;%迭代次數為10

[fr,psfe]=deconvblind(g,initpsf,numit,dampar,weight);

subplot(2,2,3),imshow(pixeldup(psfe,73),);title('10次迭代估計的psf')

numit=20;%迭代次數為20

[fr,psfe]=deconvblind(g,initpsf,numit,dampar,weight);

subplot(2,2,4),imshow(pixeldup(psfe,73),);title('20次迭代估計的psf')

使用openmp加速法線估計

在pcl中使用pcl normalestimation來對法線進行估計。但在對數量較大的點雲進行法線計算時會耗時很長,實時性很差。pcl庫提供了利用openmp來進行多核 多執行緒開發,以加快計算速度。但在使用時有可能出現以下問題 num threads 子句的引數必須是正值。這種結果通常出現在32...

引數估計 點估計和區間估計

引數估計就是根據樣本統計量的數值對總體引數進行估計的過程。根據引數估計的性質不同,可以分成兩種型別 點估計和區間估計。點估計就是用樣本統計量的某一具體數值直接推斷未知的總體引數。例如,在進行有關小學生身高的研究中,隨機抽取1000名小學生並計算出他們的平均身高為1.46公尺。如果直接用這個1.46公...

引數估計之點估計和區間估計

作者 cda資料分析師 引數估計 parameter estimation 是根據從總體中抽取的樣本估計總體分布中包含的未知引數的方法。人們常常需要根據手中的資料,分析或推斷資料反映的本質規律。即根據樣本資料如何選擇統計量去推斷總體的分布或數字特徵等。統計推斷是數理統計研究的核心問題。所謂統計推斷是...