tx = "Affine"
pnum = 0
tx = "SyN"
pnum = 1
metrics = ["MeanSquares", "GC", "Mattes"]
asr = 1
for m in range(0, len(metrics)):
mytx1 = ants.registration(
fixed=fi,
moving=mi,
type_of_transform=tx,
aff_metric=metrics[m],
aff_random_sampling_rate=asr,
)
p1 = ants.read_transform(mytx1["fwdtransforms"][pnum])
p1 = ants.get_ants_transform_parameters(p1)
for k in range(0, 2):
mytx2 = ants.registration(
fixed=fi,
moving=mi,
type_of_transform=tx,
aff_metric=metrics[m],
aff_random_sampling_rate=asr,
)
p2 = ants.read_transform(mytx2["fwdtransforms"][pnum])
p2 = ants.get_ants_transform_parameters(p2)
metval = ants.image_mutual_information(fi, mytx2["warpedmovout"])
chg = abs(p2 - p1).sum()
print(metrics[m] + " Change: try " + str(k) + " Is " + str(chg) +
" Similarity= " + str(metval))
def test_example(self):
fi = ants.image_read(ants.get_ants_data('r16')).clone('float')
mi = ants.image_read(ants.get_ants_data('r64')).clone('float')
mival = ants.image_mutual_information(fi, mi) # -0.1796141
def test_example(self):
fi = ants.image_read(ants.get_ants_data("r16")).clone("float")
mi = ants.image_read(ants.get_ants_data("r64")).clone("float")
mival = ants.image_mutual_information(fi, mi) # -0.1796141
# -*- coding:utf-8 -*-
import ants
img1 = ants.image_read(
'../mmwhs/CT_test/205_5/ct-image_crop_man_reg_resize/ct_train_1013_image.nii.gz'
)
img2 = ants.image_read(
'../mmwhs/CT_test/205_5/ct-image_crop_man_reg_resize/ct_train_1014_image.nii.gz'
)
print(img2)
print(img1)
g1 = ants.iMath_grad(img1)
g2 = ants.iMath_grad(img2)
reg1 = ants.registration(img1, img2, 'SyNOnly')
demonsMetric = ['demons', g1, g2, 1, 1]
ccMetric = ['CC', img1, img2, 2, 4]
metrics = list()
metrics.append(demonsMetric)
print("start to registration...........")
reg2 = ants.registration(img1, img2, 'SyNOnly', multivariate_extras=metrics)
metrics.append(ccMetric)
reg3 = ants.registration(img1, img2, 'SyNOnly', multivariate_extras=metrics)
print(ants.image_mutual_information(img1, img2))
print(ants.image_mutual_information(img1, reg1['warpedmovout']))
print(ants.image_mutual_information(img1, reg2['warpedmovout']))
print(ants.image_mutual_information(img1, reg3['warpedmovout']))