Укажите два изображения Используя пространственную ссылку, чтобы улучшить отображение

В этом примере показано, как использовать пространственные объекты привязки, чтобы изучить пространственные отношения между двумя изображениями в регистрации изображений и отобразить их эффективно. Этот пример приносит одно из изображений, названных moving отобразите, в выравнивание с другим изображением, названным fixed изображение.

Считайте два изображения той же сцены, которые немного неправильно выравниваются.

fixed = imread('westconcordorthophoto.png');
moving = imread('westconcordaerial.png');

Отобразите движущееся (незарегистрированное) изображение.

iptsetpref('ImshowAxesVisible','on')
imshow(moving)
text(size(moving,2),size(moving,1)+30, ...
    'Image courtesy of mPower3/Emerge', ...
    'FontSize',7,'HorizontalAlignment','right');

Figure contains an axes object. The axes object contains 2 objects of type image, text.

Загрузите MAT-файл, который содержит предварительно выобранные контрольные точки для fixed и moving изображения и создают подгонку геометрического преобразования к контрольным точкам, с помощью fitgeotrans .

load westconcordpoints
tform = fitgeotrans(movingPoints, fixedPoints, 'projective');

Выполните преобразование, необходимое, чтобы указать moving отобразите с fixed отобразите, с помощью imwarp . Этот пример использует дополнительный 'FillValues' параметр, чтобы задать (белое) значение заливки, который поможет при отображении fixed отобразите по преобразованному moving отобразите, чтобы проверять регистрацию. Заметьте, что полное содержимое геометрически преобразованного moving изображение присутствует, теперь названное registered . Также обратите внимание, что нет никаких пустых строк или столбцов.

registered = imwarp(moving, tform,'FillValues', 255);
figure, imshow(registered);

Figure contains an axes object. The axes object contains an object of type image.

Наложите преобразованное изображение, registered , по fixed отобразите, с помощью imshowpair . Заметьте, как два изображения появляются misregistered. Это происходит потому что imshowpair принимает, что изображения находятся оба во внутренней системе координат по умолчанию. Следующие шаги обеспечивают два способа исправить эту проблему отображения.

figure, imshowpair(fixed,registered,'blend');

Figure contains an axes object. The axes object contains an object of type image.

Ограничьте преобразованное изображение, registered , к одинаковому числу строк и столбцам и тем же пространственным пределам как fixed изображение. Это гарантирует, что зарегистрированное изображение кажется зарегистрированным в фиксированном изображении, но области зарегистрированного изображения, которое экстраполировало бы вне степени фиксированного изображения, отбрасываются. Для этого создайте пространственный объект привязки по умолчанию, который задает размер и местоположение fixed отобразите и используйте imwarp's 'OutputView' параметр, чтобы создать ограниченное передискретизируемое изображение registered1. Отобразите зарегистрированное изображение по фиксированному изображению. В этом представлении изображения, кажется, были указаны, но не все незарегистрированное изображение отображается.

Rfixed = imref2d(size(fixed));
registered1 = imwarp(moving,tform,'FillValues', 255,'OutputView',Rfixed);
figure, imshowpair(fixed,registered1,'blend');

Figure contains an axes object. The axes object contains an object of type image.

Как альтернатива, используйте дополнительный imwarp синтаксис, который возвращает выходной параметр пространственный объект привязки, который указывает на положение полного преобразованного изображения в той же внутренней системе координат по умолчанию как fixed изображение. Отобразите зарегистрированное изображение по фиксированному изображению и обратите внимание что теперь полный registered изображение отображается.

[registered2, Rregistered] = imwarp(moving, tform,'FillValues', 255);
figure, imshowpair(fixed,Rfixed,registered2,Rregistered,'blend');

Figure contains an axes object. The axes object contains an object of type image.

Очистка.

iptsetpref('ImshowAxesVisible','off')