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