I want to compute blur kernel of image using radon transform, but i am not getting correct result. is my code correct?

1 次查看（过去 30 天）

Pornima Birajdar 2015-7-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/230907-i-want-to-compute-blur-kernel-of-image-using-radon-transform-but-i-am-not-getting-correct-result-i

评论： Image Analyst 2015-7-23

divides image in to overlapping blocks I=imread('runman.jpg'); I=rgb2gray(I);

        %function computes blur extent
fun2=@BlurEstimation;
fun3=@BlurEstimationMag;

dir=blkproc(I,[20 20],[2 2],fun2)%process overlapping block mag=blkproc(I,[20 20],[10 10],fun3)%process overlapping block

figure;quiver(dir,dir);

function [theta] = BlurEstimation(h) len=5; %length of kernel n=60; %noise p=1/len*[1 1 1 1 1]; % blurring kernel for linear motion

                   %..........for directional motion rotate p by %theta
 i1=imfilter(h,p,'conv')+n;    %....(1)
 di=diff(i1);         %differentiating (1)
 i1=log(fft2(di));          %fourier transform of (1).

theta = 0:180; [r,xp] = radon(i1,theta); %radon transform of i1

r=fft2(r); %fourier transform of radon transform

% Find the location of the peak of the radon transform image. maxR= max(r(:)); [row col]=find(r==maxR); theta=90-radtodeg(angle(maxR)); % minR=min(r(:)) % thetamin=radtodeg(angle(minR)) theta=90-col end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Image Analyst 2015-7-23

In the MATLAB editor, type control-a then control-i to properly indent your code. Then type control-c to copy it to the clipboard, then paste it here. Then highlight it and click {}Code. Then we'll be able to read it.

http://www.mathworks.com/matlabcentral/answers/13205-tutorial-how-to-format-your-question-with-markup

请先登录，再进行评论。

请先登录，再回答此问题。