imregconfig
強度ベースのレジストレーションの構成
説明
[ は、強度ベースのイメージ レジストレーションを実行するために optimizer,metric] = imregconfig(modality)imregister に渡す optimizer と metric の構成を作成します。modality はイメージ キャプチャ モダリティを指定します。
例
イメージをワークスペースに読み込んで表示します。明度およびコントラストが似ているため、これらのイメージはモノモーダルになります。
fixed = imread('pout.tif'); moving = imrotate(fixed, 5, 'bilinear', 'crop'); imshowpair(fixed, moving,'Scaling','joint')
モダリティを 'monomodal' に設定して、オプティマイザーとメトリクスを作成します。
[optimizer, metric] = imregconfig('monomodal')optimizer =
registration.optimizer.RegularStepGradientDescent
Properties:
GradientMagnitudeTolerance: 1.000000e-04
MinimumStepLength: 1.000000e-05
MaximumStepLength: 6.250000e-02
MaximumIterations: 100
RelaxationFactor: 5.000000e-01
metric = registration.metric.MeanSquares This class has no properties.
オプティマイザーとメトリクスを imregister に渡し、レジストレーションを実行します。
movingRegistered = imregister(moving,fixed,'rigid',optimizer, metric);レジストレーションが行われたイメージを表示します。
figure imshowpair(fixed, movingRegistered,'Scaling','joint')
2 つのイメージを読み取ります。この例では、膝の 2 個の MRI (磁気共鳴) イメージを使用します。固定イメージはスピン エコー イメージであり、移動イメージは反転回復を行ったスピン エコー イメージです。同時に 2 個の矢状スライスが取得されますが、ややずれています。
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);レジストレーションが行われたイメージを表示します。
imshowpair(fixed,movingRegistered,"Scaling","joint")
