Add SetCoefficientImages for B-spline transform about simpleitk HOT 8 CLOSED

import SimpleITK as sitk
import numpy as np

# Create the transformation (when working with images it is easier to use the BSplineTransformInitializer function
# or its object oriented counterpart BSplineTransformInitializerFilter).
dimension = 3
spline_order = 3
direction_matrix_row_major = [1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0] # identity, mesh is axis aligned
origin = [-1.0,-1.0,-1.0]  
domain_physical_dimensions = [2,2,2]

bspline = sitk.BSplineTransform(dimension, spline_order)
bspline.SetTransformDomainOrigin(origin)
bspline.SetTransformDomainDirection(direction_matrix_row_major)
bspline.SetTransformDomainPhysicalDimensions(domain_physical_dimensions)
bspline.SetTransformDomainMeshSize((3,4,5))

# Random displacement of the control points.
original_parameters = np.random.random(len(bspline.GetParameters()))
bspline.SetParameters(original_parameters)

image_0, image_1, image_2 = bspline.GetCoefficientImages()
retrieved_parameters = list(sitk.GetArrayViewFromImage(image_0).flatten()) + list(sitk.GetArrayViewFromImage(image_1).flatten()) + list(sitk.GetArrayViewFromImage(image_2).flatten())

print(np.max(original_parameters - retrieved_parameters))

image_0, image_1, image_2 = bspline_tx.GetCoefficientImages()
image_0_arr = sitk.GetArrayFromImage(image_0)
image_1_arr = sitk.GetArrayFromImage(image_1)
image_2_arr = sitk.GetArrayFromImage(image_2)
aggregated_params = (
            list(image_0_arr.flatten())
            + list(image_1_arr.flatten())
            + list(image_2_arr.flatten())
    )
bspline_tx.SetParameters(aggregated_params)

image_0, image_1, image_2 = bspline_tx.GetCoefficientImages()
image_0_arr = np.swapaxes(sitk.GetArrayFromImage(image_0), 0, 2)
image_1_arr = np.swapaxes(sitk.GetArrayFromImage(image_1), 0, 2)
image_2_arr = np.swapaxes(sitk.GetArrayFromImage(image_2), 0, 2)
aggregated_params = (
            list(image_0_arr.flatten())
            + list(image_1_arr.flatten())
            + list(image_2_arr.flatten())
    )
bspline_tx.SetParameters(aggregated_params)