回复: 大家帮忙看一看哪里出了问题呢
lz T(3,1)和T(3,2)的值并没有在实际操作中被使用 空传入而已
有两个方法
一是改变fc=[m/2+a,n/2+b,1]中旋转的中心 这样能够产生移动的效果 但是实际平移的是一个圆弧 如果你明白旋转的数学原理 应该可以理解我的意思
二是在g(ii,jj) = BilinearInterplate(f,pr);后加入平移的语句 g(ii+a,jj+b) = g(ii,jj),将你上面的10带入a,b,同时需要考虑改变pr(1)>=m||pr(1)<=1||pr(2)<=1||pr(2)>=n的门限 时旋转及平移后仍在图像可视范围的部分不要被清零 否则边角会出现“缺失现象”
PS:单独旋转是不需要建立3*3的模板的 根据数学上的知识 T=[cos(theta),sin(theta);-sin(theta),cos(theta)]即可 matlab建立这样的模板在于 同时可以设置其它参数以实现多种spatial transform 例如 lz的主函数代码文件直接改为如下 亦可以执行旋转操作
function image_transform2 %只进行旋转
f=imread('D:\lena256.bmp');
[m,n]=size(f);
imshow(f);
fc=[m/2,n/2]
theta=pi/18;
T=[cos(theta),sin(theta);-sin(theta),cos(theta)];
ii=0;
for i=1:m
ii=ii+1;
jj=0;
for j=1:n
jj=jj+1;
pf=[i,j];
pr=T*(pf-fc)'+fc';
if pr(1)>=m||pr(1)<=1||pr(2)<=1||pr(2)>=n
g(ii,jj)=0 ;
else
g(ii,jj) = BilinearInterplate(f,pr);
end
end
end
Thx for reading.
PS:若还算满意,直接点击“Thanks”,再次登陆时亦便于查看回答是否真的帮到你了。
个人观点 仅供参考 多多交流 相互学习
|