レジストレーションの前処理手順としての位相相関の使用

この例では、自動イメージ レジストレーションの準備手順として位相相関を使用する方法を説明します。このプロセスでは、関数 imregcorr を使用して位相相関を実行し、そのレジストレーション結果を関数 imregister に渡して最適化ベースのレジストレーションの初期条件として使用します。位相相関と最適化ベースのレジストレーションは相補アルゴリズムです。位相相関は全体的な配置を検出するのに適しており、イメージの位置が大幅にずれている場合でも使用できます。最適化ベースのレジストレーションは、初期条件が適切であれば正確な位置合わせの解を求めるのに適しています。

レジストレーションの参照イメージとなるイメージを読み取ります。

fixed = imread('cameraman.tif');
imshow(fixed);

回転、等方性スケーリングおよび y 方向のせん断を使用して、このイメージを意図的に歪めることにより、レジストレーションされていないイメージを作成します。

theta = 170;
rot = [cosd(theta) sind(theta) 0;... 
    sind(theta) cosd(theta) 0;... 
    0 0 1]; 
sc = 2.3;
scale = [sc 0 0; 0 sc 0; 0 0 1]; 
sh = 0.5;
shear = [1 sh 0; 0 1 0; 0 0 1]; 

tform = affine2d(shear*scale*rot);
moving = imwarp(fixed,tform); 

ノイズをイメージに追加し、結果を表示します。

moving = moving + uint8(10*rand(size(moving)));
imshow(moving)

これら 2 つのイメージを位置合わせするために必要なレジストレーションを推定します。関数 imregcorr により、変換を定義する affine2d オブジェクトが返されます。

tformEstimate = imregcorr(moving,fixed);

推定された幾何学的変換を位置がずれているイメージに適用します。'OutputView' を指定して、レジストレーションされたイメージが必ず参照イメージと同じサイズになるようにします。元のイメージとレジストレーションされたイメージを左右に並べて表示します。imregcorr によって、2 つのイメージの回転やスケーリングの相違点がうまく処理されていることがわかります。レジストレーションされたイメージ movingReg は元のイメージ fixed とほぼ同じ位置にあります。しかし、いくらかずれが残っています。imregcorr では、回転やスケーリングの歪みは適切に処理できますが、せん断の歪みは十分に処理できません。

Rfixed = imref2d(size(fixed));
movingReg = imwarp(moving,tformEstimate,'OutputView',Rfixed);
imshowpair(fixed,movingReg,'montage');

関数 imshowpair を使用して元のイメージの上に位置合わせされたイメージを重ねて表示します。この表示では、imshowpair は色を使用して、位置のずれがある領域を強調表示しています。

imshowpair(fixed,movingReg,'falsecolor');

レジストレーションを完了するには、imregcorr によって返された推定された変換を関数 imregister に初期条件として渡します。imregister は演算の開始時に 2 つのイメージの位置がほぼ合っている場合、より効果的です。この情報は imregister によって推定された変換を使用して imregcorr に提供します。この例では、同一センサーで得られた 2 つのイメージ レジストレーションに対する既定のオプティマイザーおよびメトリクスの値 ('monomodal') を使用します。

[optimizer, metric] = imregconfig('monomodal');
movingRegistered = imregister(moving, fixed,...
    'affine', optimizer, metric,'InitialTransformation',tformEstimate);

このレジストレーションの結果を表示します。imregcorr によって適切な初期条件が提供されたため、imregister では正確なレジストレーションが行われています。

imshowpair(fixed, movingRegistered,'Scaling','joint');

参考

imregcorr | imregister | imregconfig | imwarp

×

例を開く

この例の変更されたバージョンがあります。編集された方の例を開きますか?

いいえ、変更されたバージョンを上書きします。 はい