21 import SimpleITK
as sitk
24 fixed_image =
sitk.Image(100, 100, sitk.sitkUInt8)
25 fixed_image[11:20, 11:20] = 200
28 moving_image =
sitk.Image(100, 100, sitk.sitkUInt8)
29 moving_image[51:60, 51:60] = 69
34 fixed_landmarks = [10, 10, 20, 10, 20, 20]
35 moving_landmarks = [50, 50, 60, 50, 60, 60]
40 landmark_initializer.SetFixedLandmarks(fixed_landmarks)
41 landmark_initializer.SetMovingLandmarks(moving_landmarks)
46 output_transform = landmark_initializer.Execute(transform)
48 print(output_transform)
53 moving_image, fixed_image, transform=output_transform, defaultPixelValue=150
62 out_name = sys.argv[1]
64 out_name =
"landmark_transform.tfm"