1

下の画像のアフィン変換行列の決定に問題があります。

元の画像:

ここに画像の説明を入力

アフィン変換された画像:

ここに画像の説明を入力

アフィン変換行列を解決するために画像上の 2 点を決定しましたが、得られた結果はオリジナルを目的に変換しません。

コードは以下のとおりです。最初に行列を解き、それを使用して元のバージョンを変換されたバージョンに変換します。

% Pixel values
p1_aff = [164; 470]; 
p1_nor = [1; 512];
p2_aff = [470; 164];
p2_nor = [512; 1];
%p3_aff = [131;68];
%p3_nor = [166;61];
%p4_aff = [328;90];
%p4_nor = [456;59];

%Transformation matrix
syms a11 a12 a21 a22;
A = [a11 a12; a21 a22];

%Solving matrix
[Sx, Sy, Sz, Sk] = solve(A*p1_nor==p1_aff, A*p2_nor==p2_aff);
a11_d = double(Sx);
a12_d = double(Sy);
a21_d = double(Sz);
a22_d = double(Sk);

lena = imread('lena512.png');
new_aff = uint8(zeros(size(lena)));

for i = 1:size(lena,1)
    for j = 1:size(lena,2)
        % Applying affine transformation
        new_coord = [a11_d a12_d 0; a21_d a22_d 0; 0 0 1]*[i; j; 1];

        % Nearest-Neighbor interpolation for placing new pixels
        if(round(new_coord(1)) > 0 && round(new_coord(2)) > 0)
            new_aff(round(new_coord(1)),round(new_coord(2))) = lena(i,j);
        end
    end
end

imwrite(new_aff, 'lenaAffine_new.png');

上記のコードの最後に、次の画像が表示されます。 ここに画像の説明を入力

誰かがここで何が悪いのか理解していますか? 私は夢中になります。

4

1 に答える 1

4

アフィン変換を決定するには、2 つの対応する点では不十分です。少なくとも 6 つの一致が必要です (私が間違っていなければ)。

cpselectgui を使用して対応する点を選択します。

>> [input_points, base_points] = cpselect(oim2, oim1, 'Wait', true);

次に、次を使用して変換できます。

>> T = cp2tform( input_points, base_points, 'affine' );
>> aim2 = tformarray( oim2, T, makeresampler('cubic','fill'), [2 1], [2 1], size(oim1(:,:,1)'), [], 0 );
于 2013-10-27T13:38:22.843 に答える