将图像的代数中心转移到几何中心上

function  k2  = CenPtran( k,m,n)

%UNTITLED Summary of this function goes here

%   Detailed explanation goes here

centerGPos = cenP(k);

centerGPos1 = [floor(m/2)+1,floor(n/2)+1];

if mod(abs(centerGPos(1)-centerGPos1(1)),2) == 0 && mod(abs(centerGPos(2)-centerGPos1(2)),2) ==0

    k1 = zeros(m+abs(centerGPos(1)-centerGPos1(1)),n+abs(centerGPos(2)-centerGPos1(2)));

elseif  mod((centerGPos(1)-centerGPos1(1)),2) == 0&&mod(abs(centerGPos(2)-centerGPos1(2)),2) ~=0

    k1 = zeros(m+abs(centerGPos(1)-centerGPos1(1)),n+abs(centerGPos(2)-centerGPos1(2))+1);

elseif  mod((centerGPos(2)-centerGPos1(2)),2) ==0&&mod((centerGPos(1)-centerGPos1(1)),2) ~= 0

    k1 = zeros(m+abs(centerGPos(1)-centerGPos1(1))+1,n+abs(centerGPos(2)-centerGPos1(2)));

else

    k1 = zeros(m+abs(centerGPos(1)-centerGPos1(1))+1,n+abs(centerGPos(2)-centerGPos1(2))+1);

end

[k1m k1n] = size(k1);

if centerGPos(1)-centerGPos1(1)>= 0 && centerGPos(2)-centerGPos1(2)>=0

    for ki =1:m

        for kj =1:n

            k1(ki,kj) = k(ki,kj);              

        end

    end

elseif centerGPos(1)-centerGPos1(1)<0 && centerGPos(2)-centerGPos1(2)==0

    for ki = (abs(centerGPos(1)-centerGPos1(1))+1):(abs(centerGPos(1)-centerGPos1(1))+m)

        for kj =1:k1n

            k1(ki,kj) = k(ki-abs(centerGPos(1)-centerGPos1(1)),kj);               

        end

    end

elseif centerGPos(2)-centerGPos1(2)<0&&centerGPos(1)-centerGPos1(1)== 0

    for ki = 1:k1m

        for kj =( (abs(centerGPos(2)-centerGPos1(2)))+1):((abs(centerGPos(2)-centerGPos1(2)))+n)

            k1(ki,kj) = k(ki,kj-abs(centerGPos(2)-centerGPos1(2)));               

        end

    end

elseif  centerGPos(1)-centerGPos1(1)<0&& centerGPos(2)-centerGPos1(2)>0 

     for ki = (abs(centerGPos(1)-centerGPos1(1))+1):(abs(centerGPos(1)-centerGPos1(1))+m)

        for kj =1:n

            k1(ki,kj) = k(ki-abs(centerGPos(1)-centerGPos1(1)),kj);               

        end

     end

elseif centerGPos(1)-centerGPos1(1)>0&& centerGPos(2)-centerGPos1(2)<0 

        for ki =1:m

            for kj =( (abs(centerGPos(2)-centerGPos1(2)))+1):((abs(centerGPos(2)-centerGPos1(2)))+n)

                k1(ki,kj) = k(ki,kj-abs(centerGPos(2)-centerGPos1(2)));               

            end

        end

else

    for ki = (abs(centerGPos(1)-centerGPos1(1))+1):(abs(centerGPos(1)-centerGPos1(1))+m)

        for kj =( (abs(centerGPos(2)-centerGPos1(2)))+1):((abs(centerGPos(2)-centerGPos1(2)))+n)

            k1(ki,kj) = k(ki-abs(centerGPos(1)-centerGPos1(1)),kj-abs(centerGPos(2)-centerGPos1(2)));               

        end

    end

end

centerGPos2 = cenP(k1);

x1 = centerGPos2(1)-floor(m/2);

x2 = x1+2*floor(m/2);

y1 =  centerGPos2(2)-floor(n/2);

y2 = y1+2*floor(n/2);

k2 = k1(x1:x2,y1:y2);