|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我编凑了这么一个程序:
>>A=dicomread('g:\Transv01_DS0681.dcm'); %处理对象A是医学图片,dicom格式,
2维数组,大小128*128
>>PSF=fspecial('gaussian',5,8); %成像系统的空间分辨率是5微米,
点扩展函数(PSF)近似选为高斯,
不知高斯的参数这样选是否合理?
成像系统空间分辨率与图像像素有什么关系?
>>INITPSF=ones(size(PSF));
>>WT=zeros(size(A));
>>WT(5:end-4,5:end-4) = 1; %WT的取法也心存疑虑,请问怎么取才合理?
>>[J P]=deconvblind(A,INITPSF,20,[],WT); %NUMIT随便取20,DAMPAR取默认,合理吗?
>>max(max(A))
ans =
169
>>max(max(J))
ans =
176
>>subplot(221),imshow(A,[0,169]),title('blurred image');
Warning: The initial magnification of the image is set to 'fit' in a docked figure.
> In imshow at 181
>>subplot(222),imshow(PSF,[]),title('PSF');
Warning: The initial magnification of the image is set to 'fit' in a docked figure.
> In imshow at 181
>>subplot(223),imshow(J,[0,176]),title('deblurred image');
Warning: The initial magnification of the image is set to 'fit' in a docked figure.
> In imshow at 181
>>subplot(224),imshow(P,[]),title('restored PSF');
Warning: The initial magnification of the image is set to 'fit' in a docked figure.
> In imshow at 181
%出现这样的警告怎么处理?
我的图像是从一幅多帧的图像拿出来的,图像的各种属性选择都组合过了,都是一样的结果。
以上,输入图像A和结果图像J是int16,其它为double。
请大虾帮帮我看看是什么问题
[ 本帖最后由 eight 于 2008-5-4 21:18 编辑 ] |
|