大家帮忙看一看哪里出了问题呢

[gj9951101]

谁能帮我看看这个对图像刚性变换的程序哪里出了问题，为什么只能旋转，不能平移呢？
function image_transform
f=imread('E:\photos\128128.bmp');
[m,n]=size(f);
fc=[m/2,n/2,1]
theta=pi/18;
T=[cos(theta),sin(theta),0;
-sin(theta),cos(theta),0
10,10,1];
ii=0;
for i=1:m
ii=ii+1; jj=0;
for j=1:n
jj=jj+1;
pf=[i,j,1];
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

[silas_xue]

lz 能否把BilinearInterplate( )也贴出来
大家一起讨论一下 是只进行双线性差值的处理而已么？

[gj9951101]

function Pixel = BilinearInterplate(im,position_s)

x0 = position_s(1);
y0 = position_s(2);
x = floor(x0);
y = floor(y0);
dx = (x0 - x); %相对位置
dy = (y0 - y);

pixel1=(1-dx)*im(x,y)+dx*im(x+1,y);
pixel2=(1-dx)*im(x,y+1)+dx*im(x+1,y+1);
Pixel=(1-dy)*pixel1+dy*pixel2;
这是那个插值的函数，平移的问题我已经解决了，可是现在的问题是，旋转一次，比如说旋转10度，出来的图像是对的，可是在旋转图像的基础上在旋转任意角度，出来的图像就变形了，不对了，这是为什么呢

[silas_xue]

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”，再次登陆时亦便于查看回答是否真的帮到你了。
个人观点 仅供参考 多多交流 相互学习