Регистрация изображения на основе интенсивности
moving_reg = imregister(moving,fixed,transformType,optimizer,metric)moving, чтобы она была зарегистрирована с эталонным изображением, fixed. transformType определяет тип выполняемого преобразования. metric определяет количественную меру подобия между изображениями для оптимизации. optimizer описывает метод оптимизации метрики. Функция возвращает зарегистрированное изображение, moving_reg.
[ преобразует пространственно привязанное изображение moving_reg,R_reg] = imregister(moving,Rmoving,fixed,Rfixed,transformType,optimizer,metric)moving чтобы он был зарегистрирован с пространственно привязанным изображением fixed. Rmoving и Rfixed пространственные ссылочные объекты, которые описывают мировые координатные пределы и разрешение moving и fixed.
___ = imregister(___, указывает дополнительные параметры с одним или несколькими аргументами пары имя-значение.Name,Value)
Прочитайте два изображения. В этом примере используются два магнитно-резонансных (МРТ) изображения колена. Фиксированное изображение представляет собой спинное эхо-изображение, а движущееся изображение представляет собой спинное эхо-изображение с инверсионным восстановлением. Два сагиттальных фрагмента были получены одновременно, но немного не выровнены.
fixed = dicomread('knee1.dcm'); moving = dicomread('knee2.dcm');
Просмотр несовпадающих изображений.
imshowpair(fixed, moving,'Scaling','joint')
Создайте оптимизатор и метрику, установив модальность как 'multimodal' поскольку изображения поступают от различных датчиков.
[optimizer, metric] = imregconfig('multimodal')optimizer =
registration.optimizer.OnePlusOneEvolutionary
Properties:
GrowthFactor: 1.050000e+00
Epsilon: 1.500000e-06
InitialRadius: 6.250000e-03
MaximumIterations: 100
metric =
registration.metric.MattesMutualInformation
Properties:
NumberOfSpatialSamples: 500
NumberOfHistogramBins: 50
UseAllPixels: 1
Настройте свойства оптимизатора, чтобы задача сходилась на глобальных максимумах и допускала больше итераций.
optimizer.InitialRadius = 0.009; optimizer.Epsilon = 1.5e-4; optimizer.GrowthFactor = 1.01; optimizer.MaximumIterations = 300;
Выполните регистрацию.
movingRegistered = imregister(moving, fixed, 'affine', optimizer, metric);Просмотр зарегистрированных изображений.
figure imshowpair(fixed, movingRegistered,'Scaling','joint')
Оба imregtform и imregister использовать один и тот же алгоритм регистрации. imregister выполняет дополнительный этап повторной выборки moving для получения зарегистрированного выходного изображения из оценки геометрического преобразования, вычисленной imregtform. Использовать imregtform когда требуется доступ к связывающему геометрическому преобразованию moving кому fixed. Использовать imregister когда требуется зарегистрированное выходное изображение.
Создание optimizer и metric с imregconfig функция перед вызовом imregister. Получение хороших результатов от регистрации изображений на основе оптимизации обычно требует изменения параметров оптимизатора или метрики для пары регистрируемых изображений. imregconfig функция предоставляет конфигурацию по умолчанию, которая должна рассматриваться только как начальная точка. Например, если увеличить число итераций в оптимизаторе, уменьшить размер шага оптимизатора или изменить количество выборок в стохастической метрике, регистрация улучшится до точки за счет производительности. См. выходные данные imregconfig для получения дополнительной информации о различных параметрах, которые можно изменить.
Если пространственное масштабирование изображений отличается более чем на 10%, измените их размер с помощью imresize перед их регистрацией.
Использовать imshowpair или imfuse визуализировать результаты регистрации.
Вы можете использовать imregister в автоматизированном рабочем процессе для регистрации нескольких изображений.
При наличии пространственных ссылок на регистрируемое изображение укажите информацию для imregister использование пространственных ссылочных объектов. Это помогает imregister более быстрое сближение для достижения лучших результатов, поскольку можно принимать во внимание различия в масштабе.
imfuse | imregconfig | imregcorr | imregtform | imshowpair | imwarp