[Insight-users] Image deformation filter question

[Lydia Ng]

This is a multi-part message in MIME format.

------_=_NextPart_001_01C34CB9.8CE2D04A
Content-Type: text/plain;
	charset="us-ascii"
Content-Transfer-Encoding: quoted-printable

Hi Phillip,

As per your suggestion, I've added a SetCoefficientImage() API so that
the B-spline coefficients can be set as an array images. Let me know you
run into any problems with it.

I've attached a small program that illustrates how to go from vector
field images to the B-spline coefficients images required by the
transform.
If you are interested, the end of the code shows how to create a higher
resolution field by from the B-spline coefficient images using the
ResampleImageFilter.

Hope this helps,
Lydia


> -----Original Message-----
> From: Phillip Cheng [mailto:pmcheng at u . washington . edu]
> Sent: Friday, July 11, 2003 4:41 PM
> To: insight-users at public . kitware . com
> Subject: [Insight-users] Image deformation filter question
>=20
> Hello,
>=20
> I would like to perform the following simple image deformation using
> ITK.
>=20
> 1.  Generate a 3D mesh containing a vector displacement in mm for each
> grid point.
> 2.  Apply the displacements to an image whose resolution is finer than
> the mesh, using B-splines for interpolation.
>=20
> I'm having difficulty connecting these steps efficiently in ITK.  In
> particular, it seems I need to use  itkBSplineDecompositionFilter to
> generate spline coefficients for each vector dimension in the mesh,
and
> then copy all these coefficients to a 3x3D concatenated array for use
> in itkBSplineDeformableTransform (probably inverting things somewhere
> so that I have an output->input mapping).  But the copying step wastes
> time and memory, especially if the mesh is large and I need to
generate
> many such transforms.  Is there a way to connect these pieces without
> copying the 3x3D arrays?
>=20
> I looked at itkBSplineDeformableTransformTest.cxx which wraps the
> parameter array into several itkImages, and that seems close to what I
> need, but I'm a little unclear as to how to direct the output of
> itkBSplineDecompositionFilter into the parameter images for the
> deformable transform.
>=20
> Thanks,
>=20
> Phillip Cheng
>=20
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk . org
> http://www . itk . org/mailman/listinfo/insight-users

------_=_NextPart_001_01C34CB9.8CE2D04A
Content-Type: application/octet-stream;
	name="BSplineDeformableTransformExample2.cxx"
Content-Transfer-Encoding: base64
Content-Description: BSplineDeformableTransformExample2.cxx
Content-Disposition: attachment;
	filename="BSplineDeformableTransformExample2.cxx"

DQojaW5jbHVkZSAiaXRrQlNwbGluZURlZm9ybWFibGVUcmFuc2Zvcm0uaCINCiNpbmNsdWRlICJp
dGtCU3BsaW5lRGVjb21wb3NpdGlvbkltYWdlRmlsdGVyLmgiDQojaW5jbHVkZSAiaXRrUmVzYW1w
bGVJbWFnZUZpbHRlci5oIg0KI2luY2x1ZGUgIml0a0JTcGxpbmVSZXNhbXBsZUltYWdlRnVuY3Rp
b24uaCINCiNpbmNsdWRlICJpdGtJbXBvcnRJbWFnZUZpbHRlci5oIg0KDQojaW5jbHVkZSAiaXRr
SW1hZ2VGaWxlV3JpdGVyLmgiDQojaW5jbHVkZSA8c3RyaW5nPg0KDQojaW5jbHVkZSAiaXRrVGV4
dE91dHB1dC5oIg0KDQppbnQgbWFpbigpDQp7DQogIGl0azo6T3V0cHV0V2luZG93OjpTZXRJbnN0
YW5jZShpdGs6OlRleHRPdXRwdXQ6Ok5ldygpKTsNCg0KICB1bnNpZ25lZCBpbnQgajsNCg0KICAv
KioNCiAgICogRGVmaW5lIGEgdmVjdG9yIGZpZWxkIGFzIERpbWVuc2lvbiBudW1iZXIgb2YgaW1h
Z2VzDQogICAqLw0KICBjb25zdCB1bnNpZ25lZCBpbnQgRGltZW5zaW9uID0gMjsNCiAgdHlwZWRl
ZiBkb3VibGUgUGl4ZWxUeXBlOw0KICB0eXBlZGVmIGl0azo6SW1hZ2U8UGl4ZWxUeXBlLERpbWVu
c2lvbj4gIEltYWdlVHlwZTsNCg0KICAvLyBTZXQgdXAgZmllbGQgc3BhY2luZywgb3JpZ2luLCBy
ZWdpb24NCiAgZG91YmxlIHNwYWNpbmdbRGltZW5zaW9uXTsNCiAgZG91YmxlIG9yaWdpbltEaW1l
bnNpb25dOw0KICBJbWFnZVR5cGU6OlNpemVUeXBlIHNpemU7DQogIEltYWdlVHlwZTo6UmVnaW9u
VHlwZSByZWdpb247DQoNCiAgZm9yICggaiA9IDA7IGogPCBEaW1lbnNpb247IGorKyApDQogICAg
ew0KICAgIHNwYWNpbmdbal0gPSAxMC4wOw0KICAgIG9yaWdpbltqXSAgPSAtMTAuMDsNCiAgICB9
DQoNCiAgc2l6ZVswXSA9IDU7DQogIHNpemVbMV0gPSA3Ow0KDQogIHJlZ2lvbi5TZXRTaXplKCBz
aXplICk7DQoNCiAgSW1hZ2VUeXBlOjpQb2ludGVyIGZpZWxkW0RpbWVuc2lvbl07DQogIGZvciAo
IGogPSAwOyBqIDwgRGltZW5zaW9uOyBqKysgKQ0KICAgIHsNCiAgICBmaWVsZFtqXSA9IEltYWdl
VHlwZTo6TmV3KCk7DQogICAgZmllbGRbal0tPlNldFNwYWNpbmcoIHNwYWNpbmcgKTsNCiAgICBm
aWVsZFtqXS0+U2V0T3JpZ2luKCBvcmlnaW4gKTsNCiAgICBmaWVsZFtqXS0+U2V0UmVnaW9ucygg
cmVnaW9uICk7DQogICAgZmllbGRbal0tPkFsbG9jYXRlKCk7DQogICAgfQ0KDQogIC8vIGZpbGwg
dGhlIGZpZWxkIHdpdGggYSBjb25zdGFudCBkaXNwbGFjbWVudA0KICBpdGs6OlZlY3Rvcjxkb3Vi
bGUsRGltZW5zaW9uPiB2Ow0KICB2WzBdID0gNTsNCiAgdlsxXSA9IDc7DQoNCiAgZm9yICggaiA9
IDA7IGogPCBEaW1lbnNpb247IGorKyApDQogICAgew0KICAgIGZpZWxkW2pdLT5GaWxsQnVmZmVy
KCB2W2pdICk7DQogICAgfQ0KDQoNCiAgLyoqDQogICAqIFRoaXMgc2VjdGlvbiBvZiBjb2RlIGls
bHVzdHJhdGUgaG93IHRvDQogICAqIGNvbnZlcnQgdGhlIHZlY3RvciBmaWVsZCBpbWFnZXMgaW50
byBCLXNwbGluZSBjb2VmZmljaWVudHMNCiAgICogZm9yIHVzZSBpbiB0aGUgQlNwbGluZURlZm9y
bWFibGVUcmFuc2Zvcm0NCiAgICoNCiAgICogVGhlIG91dHB1dCBvZiB0aGUgQlNwbGluZURlY29t
cG9zaXRpb25JbWFnZUZpbHRlciBpcyB0aGVuDQogICAqIGNvcGllZCBpbnRvIGEgcGFyYW1ldGVy
cyBhcnJheSB0byB1c2Ugd2l0aCB0aGUgQlNwbGluZURlZm9ybWFibGVUcmFuc2Zvcm0NCiAgICoN
CiAgICovDQoNCiAgLy8gU2V0IHVwIHRoZSB0cmFuc2Zvcm0NCiAgY29uc3QgdW5zaWduZWQgaW50
IFNwbGluZU9yZGVyID0gMzsNCiAgdHlwZWRlZiBkb3VibGUgQ29vcmRSZXA7DQogIHR5cGVkZWYg
aXRrOjpCU3BsaW5lRGVmb3JtYWJsZVRyYW5zZm9ybTxDb29yZFJlcCxEaW1lbnNpb24sU3BsaW5l
T3JkZXI+IFRyYW5zZm9ybVR5cGU7DQogIA0KICBUcmFuc2Zvcm1UeXBlOjpQb2ludGVyIHRyYW5z
Zm9ybSA9IFRyYW5zZm9ybVR5cGU6Ok5ldygpOw0KDQogIFRyYW5zZm9ybVR5cGU6OkltYWdlUG9p
bnRlciBjb2VmZltEaW1lbnNpb25dOw0KDQogIHR5cGVkZWYgaXRrOjpCU3BsaW5lRGVjb21wb3Np
dGlvbkltYWdlRmlsdGVyPEltYWdlVHlwZSxJbWFnZVR5cGU+IERlY29tcG9zaXRpb25UeXBlOw0K
DQogIC8vIFNldCB1cCB0aGUgZGVjb21wb3NpdGlvbiBmaWx0ZXINCiAgRGVjb21wb3NpdGlvblR5
cGU6OlBvaW50ZXIgZGVjb21wb3NpdGlvbiA9IERlY29tcG9zaXRpb25UeXBlOjpOZXcoKTsNCg0K
ICBkZWNvbXBvc2l0aW9uLT5TZXRTcGxpbmVPcmRlciggU3BsaW5lT3JkZXIgKTsNCg0KICBmb3Ig
KCBqID0gMDsgaiA8IERpbWVuc2lvbjsgaisrICkNCiAgICB7DQogICAgZGVjb21wb3NpdGlvbi0+
U2V0SW5wdXQoIGZpZWxkW2pdICk7DQogICAgZGVjb21wb3NpdGlvbi0+VXBkYXRlKCk7DQogICAg
Y29lZmZbal0gPSBkZWNvbXBvc2l0aW9uLT5HZXRPdXRwdXQoKTsNCiAgICBjb2VmZltqXS0+RGlz
Y29ubmVjdFBpcGVsaW5lKCk7DQogICAgfQ0KDQogIHRyYW5zZm9ybS0+U2V0Q29lZmZpY2llbnRJ
bWFnZSggY29lZmYgKTsNCiAgdHJhbnNmb3JtLT5QcmludCggc3RkOjpjb3V0ICk7DQoNCiAgLyoq
DQogICAqIFRyYW5zZm9ybSBzb21lIHBvaW50cw0KICAgKi8NCiAgVHJhbnNmb3JtVHlwZTo6SW5w
dXRQb2ludFR5cGUgaW5wdXRQb2ludDsNCiAgVHJhbnNmb3JtVHlwZTo6T3V0cHV0UG9pbnRUeXBl
IG91dHB1dFBvaW50Ow0KDQogIHRyeQ0KICAgIHsNCg0KICAgIC8vIHRyeSBhIHBvaW50IGluc2lk
ZSB0aGUgdmFsaWRlIHJlZ2lvbg0KICAgIGlucHV0UG9pbnQuRmlsbCggMTAuMCApOw0KICAgIG91
dHB1dFBvaW50ID0gdHJhbnNmb3JtLT5UcmFuc2Zvcm1Qb2ludCggaW5wdXRQb2ludCApOw0KICAg
IHN0ZDo6Y291dCA8PCAiIElucHV0UG9pbnQ6ICIgPDwgaW5wdXRQb2ludDsNCiAgICBzdGQ6OmNv
dXQgPDwgIiBPdXRwdXRQb2ludDogIiA8PCBvdXRwdXRQb2ludDsNCiAgICBzdGQ6OmNvdXQgPDwg
c3RkOjplbmRsOw0KDQogICAgLy8gdHJ5IGEgcG9pbnQgb24gdGhlIHZhbGlkIHJlZ2lvbiBib3Vu
ZGFyeQ0KICAgIGlucHV0UG9pbnQuRmlsbCggMC4wICk7DQogICAgb3V0cHV0UG9pbnQgPSB0cmFu
c2Zvcm0tPlRyYW5zZm9ybVBvaW50KCBpbnB1dFBvaW50ICk7DQogICAgc3RkOjpjb3V0IDw8ICIg
SW5wdXRQb2ludDogIiA8PCBpbnB1dFBvaW50Ow0KICAgIHN0ZDo6Y291dCA8PCAiIE91dHB1dFBv
aW50OiAiIDw8IG91dHB1dFBvaW50Ow0KICAgIHN0ZDo6Y291dCA8PCBzdGQ6OmVuZGw7DQoNCiAg
ICAvLyB0cnkgYSBwb2ludCBvbiB0aGUgdmFsaWQgcmVnaW9uIGJvdW5kYXJ5DQogICAgaW5wdXRQ
b2ludFswXSA9IDE5Ljk7DQogICAgaW5wdXRQb2ludFsxXSA9IDMwLjA7DQogICAgb3V0cHV0UG9p
bnQgPSB0cmFuc2Zvcm0tPlRyYW5zZm9ybVBvaW50KCBpbnB1dFBvaW50ICk7DQogICAgc3RkOjpj
b3V0IDw8ICIgSW5wdXRQb2ludDogIiA8PCBpbnB1dFBvaW50Ow0KICAgIHN0ZDo6Y291dCA8PCAi
IE91dHB1dFBvaW50OiAiIDw8IG91dHB1dFBvaW50Ow0KICAgIHN0ZDo6Y291dCA8PCBzdGQ6OmVu
ZGw7DQoNCiAgICAvLyB0cnkgYSBwb2ludCBvdXRzaWRlIHRoZSB2YWxpZCByZWdpb24NCiAgICBp
bnB1dFBvaW50WzBdID0gMjAuMDsNCiAgICBpbnB1dFBvaW50WzFdID0gMzAuMDsNCiAgICBvdXRw
dXRQb2ludCA9IHRyYW5zZm9ybS0+VHJhbnNmb3JtUG9pbnQoIGlucHV0UG9pbnQgKTsNCiAgICBz
dGQ6OmNvdXQgPDwgIiBJbnB1dFBvaW50OiAiIDw8IGlucHV0UG9pbnQ7DQogICAgc3RkOjpjb3V0
IDw8ICIgT3V0cHV0UG9pbnQ6ICIgPDwgb3V0cHV0UG9pbnQ7DQogICAgc3RkOjpjb3V0IDw8IHN0
ZDo6ZW5kbDsNCg0KICAgIH0NCiAgY2F0Y2goIGl0azo6RXhjZXB0aW9uT2JqZWN0JiBlcnIgKQ0K
ICAgIHsNCiAgICBzdGQ6OmNvdXQgPDwgZXJyIDw8IHN0ZDo6ZW5kbDsNCiAgICByZXR1cm4gRVhJ
VF9GQUlMVVJFOw0KICAgIH0NCg0KICAvKioNCiAgICogVGhpcyBzZWN0aW9uIG9mIGNvZGUgaWxs
dXN0cmF0ZSBob3cgdG8NCiAgICogY29udmVydCB0aGUgY29lZmZpY2llbnRzIG9mIGEgQlNwbGlu
ZURlZm9ybWFibGVUcmFuc2Zvcm0NCiAgICogYmFjayBpbnRvIGEgdmVjdG9yIGZpZWxkIG9mIGEg
aGlnaGVyIHJlc29sdXRpb24uDQogICAqDQogICAqLw0KICB0eXBlZGVmIGl0azo6SW1wb3J0SW1h
Z2VGaWx0ZXI8ZG91YmxlLERpbWVuc2lvbj4gSW1wb3J0ZXJUeXBlOw0KICBJbXBvcnRlclR5cGU6
OlBvaW50ZXIgaW1wb3J0ZXIgPSBJbXBvcnRlclR5cGU6Ok5ldygpOw0KDQogIGltcG9ydGVyLT5T
ZXRSZWdpb24oIHJlZ2lvbiApOw0KICBpbXBvcnRlci0+U2V0U3BhY2luZyggc3BhY2luZyApOw0K
ICBpbXBvcnRlci0+U2V0T3JpZ2luKCBvcmlnaW4gKTsNCiAgDQogIC8vIFNldCB1cCBuZXcgb3V0
cHV0IHNwYWNpbmcsIG9yaWdpbiBhbmQgcmVnaW9uDQogIGZvciAoIGogPSAwOyBqIDwgRGltZW5z
aW9uOyBqKysgKQ0KICAgIHsNCiAgICBzcGFjaW5nW2pdID0gNS4wOw0KICAgIG9yaWdpbltqXSAg
PSAtNS4wOw0KICAgIH0NCg0KICBzaXplWzBdID0gODsNCiAgc2l6ZVsxXSA9IDEyOw0KDQogIA0K
ICAvLyBTZXQgdXAgdGhlIHJlc2FtcGxlcg0KICB0eXBlZGVmIGl0azo6UmVzYW1wbGVJbWFnZUZp
bHRlcjxJbWFnZVR5cGUsSW1hZ2VUeXBlPiBSZXNhbXBsZVR5cGU7DQogIFJlc2FtcGxlVHlwZTo6
UG9pbnRlciByZXNhbXBsZSA9IFJlc2FtcGxlVHlwZTo6TmV3KCk7DQoNCiAgdHlwZWRlZiBpdGs6
OkJTcGxpbmVSZXNhbXBsZUltYWdlRnVuY3Rpb248SW1hZ2VUeXBlLGRvdWJsZT4gRnVuY3Rpb25U
eXBlOw0KICBGdW5jdGlvblR5cGU6OlBvaW50ZXIgZnVuY3Rpb24gPSBGdW5jdGlvblR5cGU6Ok5l
dygpOw0KDQogIHJlc2FtcGxlLT5TZXRJbnRlcnBvbGF0b3IoIGZ1bmN0aW9uICk7DQogIHJlc2Ft
cGxlLT5TZXRTaXplKCBzaXplICk7DQogIHJlc2FtcGxlLT5TZXRPdXRwdXRPcmlnaW4oIG9yaWdp
biApOw0KICByZXNhbXBsZS0+U2V0T3V0cHV0U3BhY2luZyggc3BhY2luZyApOw0KICByZXNhbXBs
ZS0+U2V0RGVmYXVsdFBpeGVsVmFsdWUoIDAuMCApOw0KDQogIC8vIFVwc2FtcGxlIHRoZSBmaWVs
ZCBpbWFnZXMgb25lIGJ5IG9uZQ0KICBJbWFnZVR5cGU6OlBvaW50ZXIgdXBzYW1wbGVkRmllbGRb
RGltZW5zaW9uXTsNCg0KICBmb3IgKCBqID0gMDsgaiA8IERpbWVuc2lvbjsgaisrICkNCiAgICB7
DQoNCiAgICByZXNhbXBsZS0+U2V0SW5wdXQoIHRyYW5zZm9ybS0+R2V0Q29lZmZpY2llbnRJbWFn
ZSgpW2pdICk7DQogICAgcmVzYW1wbGUtPlVwZGF0ZSgpOw0KDQogICAgdXBzYW1wbGVkRmllbGRb
al0gPSByZXNhbXBsZS0+R2V0T3V0cHV0KCk7DQogICAgdXBzYW1wbGVkRmllbGRbal0tPkRpc2Nv
bm5lY3RQaXBlbGluZSgpOw0KDQogICAgfQ0KDQovLyBVbmNvbW1lbnQgdG8gd3JpdGUgaW1hZ2Vz
IHRvIGZpbGUNCi8qDQogIGZvciAoIGogPSAwOyBqIDwgRGltZW5zaW9uOyBqKysgKQ0KICAgIHsN
Cg0KICAgIHR5cGVkZWYgaXRrOjpJbWFnZUZpbGVXcml0ZXI8SW1hZ2VUeXBlPiBXcml0ZXJUeXBl
Ow0KICAgIFdyaXRlclR5cGU6OlBvaW50ZXIgd3JpdGVyID0gV3JpdGVyVHlwZTo6TmV3KCk7DQoN
CiAgICBjaGFyIGJ1ZmZlcls1XTsNCiAgICBzcHJpbnRmKCBidWZmZXIsICIlZCIsIGogKTsNCg0K
ICAgIHN0ZDo6c3RyaW5nIGZpbGVOYW1lID0gImZpZWxkIjsNCiAgICBmaWxlTmFtZSArPSBidWZm
ZXI7DQogICAgZmlsZU5hbWUgKz0gIi5taGQiOw0KDQogICAgc3RkOjpjb3V0IDw8ICJXcml0aW5n
IGZpbGU6ICIgPDwgZmlsZU5hbWUgPDwgc3RkOjplbmRsOw0KICAgIHdyaXRlci0+U2V0RmlsZU5h
bWUoIGZpbGVOYW1lLmNfc3RyKCkgKTsNCiAgICB3cml0ZXItPlNldElucHV0KCB1cHNhbXBsZWRG
aWVsZFtqXSApOw0KICAgIHdyaXRlci0+V3JpdGUoKTsNCg0KICAgIH0NCiovDQoNCiAgcmV0dXJu
IEVYSVRfU1VDQ0VTUzsNCn0=

------_=_NextPart_001_01C34CB9.8CE2D04A--