[Insight-users] itkResampleImageFilter Up sample take too long?

Feng Ma mafeng at hotmail . com
Thu, 23 Oct 2003 16:58:30 -0400 

This is a multi-part message in MIME format.

------=_NextPart_000_d0d_7581_74e1
Content-Type: text/plain; format=flowed

Hi,

  I am using itkResampleImageFilter to do resampling both down sampling and 
up sampling. I found out that down sampling almost take no time, while up 
sampling takes significantly longer. Attached is my test code, with results 
as below:


       Start down sampling: 0
       End down sampling: 0.33
       Start up sampling: 0.33
       End up sampling: 107.77

  The code works as:

   First downsample a volume 512x512x251 to 64x64x50. Then up sample the 
small volume to original size.

  Any idea about this?

  Thanks a lot.

-Feng

_________________________________________________________________
See when your friends are online with MSN Messenger 6.0. Download it now 
FREE! http://msnmessenger-download . com

------=_NextPart_000_d0d_7581_74e1
Content-Type: application/octet-stream; name="test.cpp"
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="test.cpp"

I2luY2x1ZGUgPGZzdHJlYW0+CiNpbmNsdWRlIDxzdHJpbmc+CiNpbmNsdWRl
ICJpdGtDYXN0SW1hZ2VGaWx0ZXIuaCIKI2luY2x1ZGUgIml0a1Jlc2FtcGxl
SW1hZ2VGaWx0ZXIuaCIKI2luY2x1ZGUgIml0a05lYXJlc3ROZWlnaGJvcklu
dGVycG9sYXRlSW1hZ2VGdW5jdGlvbi5oIgojaW5jbHVkZSAiaXRrTGluZWFy
SW50ZXJwb2xhdGVJbWFnZUZ1bmN0aW9uLmgiCgpleHRlcm4gIkMiIHsKI2lu
Y2x1ZGUgPHN0cmluZy5oPgojaW5jbHVkZSA8c3RkaW8uaD4KI2luY2x1ZGUg
PHRpbWUuaD4KfQoKaW50IG1haW4oKQp7CiAgaW50IGksIHpvb21bM107CiAg
ZG91YmxlIGRpbXNbM10sIHNwYWNpbmdbM10sIG9yaWdpblszXSwgc21hbGxE
aW1zWzNdLCBzbWFsbFNwYWNpbmdbM10sIHNtYWxsT3JpZ2luWzNdOwogIGNs
b2NrX3Qgc3RhcnRjOwoKICBjb25zdCB1bnNpZ25lZCBzaG9ydCBWSW1hZ2VE
aW1lbnNpb24gPSAzOwogIHR5cGVkZWYgaXRrOjpJbWFnZTxpbnQsIDM+IElt
YWdlVHlwZTsKCiAgdHlwZWRlZiBpdGs6OlJlc2FtcGxlSW1hZ2VGaWx0ZXI8
SW1hZ2VUeXBlLCBJbWFnZVR5cGU+ICBSZXNhbXBsZUZpbHRlclR5cGU7CiAg
UmVzYW1wbGVGaWx0ZXJUeXBlOjpQb2ludGVyIHJlc2FtcGxlRmlsdGVyID0g
UmVzYW1wbGVGaWx0ZXJUeXBlOjpOZXcoKTsKCiAgdHlwZWRlZiAgaXRrOjpO
ZWFyZXN0TmVpZ2hib3JJbnRlcnBvbGF0ZUltYWdlRnVuY3Rpb248IAoJSW1h
Z2VUeXBlLCBkb3VibGUgPiAgTk5JbnRlcnBvbGF0b3JUeXBlOwogIE5OSW50
ZXJwb2xhdG9yVHlwZTo6UG9pbnRlciBuZWFyZXN0TmVpZ2hib3JJbnRlcnBv
bGF0b3IgPSAKCU5OSW50ZXJwb2xhdG9yVHlwZTo6TmV3KCk7CiAgCiAgdHlw
ZWRlZiBpdGs6OkxpbmVhckludGVycG9sYXRlSW1hZ2VGdW5jdGlvbjwKCUlt
YWdlVHlwZSwgZG91YmxlID4gIExpbmVhckludGVycG9sYXRvclR5cGU7CiAg
TGluZWFySW50ZXJwb2xhdG9yVHlwZTo6UG9pbnRlciBsaW5lYXJJbnRlcnBv
bGF0b3IgPSBMaW5lYXJJbnRlcnBvbGF0b3JUeXBlOjpOZXcoKTsKCiAgdHlw
ZWRlZiBpdGs6OkFmZmluZVRyYW5zZm9ybTxkb3VibGUsIFZJbWFnZURpbWVu
c2lvbj4gVHJhbnNmb3JtVHlwZTsKICBUcmFuc2Zvcm1UeXBlOjpQb2ludGVy
IHRyYW5zZm9ybSA9IFRyYW5zZm9ybVR5cGU6Ok5ldygpOwogIHRyYW5zZm9y
bS0+U2V0SWRlbnRpdHkoKTsKCiAgcmVzYW1wbGVGaWx0ZXItPlNldFRyYW5z
Zm9ybSh0cmFuc2Zvcm0pOwogIHJlc2FtcGxlRmlsdGVyLT5TZXRJbnRlcnBv
bGF0b3IobGluZWFySW50ZXJwb2xhdG9yKTsKICByZXNhbXBsZUZpbHRlci0+
U2V0RGVmYXVsdFBpeGVsVmFsdWUoMCk7CiAgCiAgSW1hZ2VUeXBlOjpQb2lu
dGVyIG91dHB1dEltYWdlID0gTlVMTDsKICBJbWFnZVR5cGU6OlBvaW50ZXIg
aW5wdXRJbWFnZSA9IEltYWdlVHlwZTo6TmV3KCk7CiAgSW1hZ2VUeXBlOjpS
ZWdpb25UeXBlIGlucHV0UmVnaW9uOwogIEltYWdlVHlwZTo6UmVnaW9uVHlw
ZSBvdXRwdXRSZWdpb247CiAgSW1hZ2VUeXBlOjpSZWdpb25UeXBlOjpJbmRl
eFR5cGUgaW5wdXRJbmRleDsKICBJbWFnZVR5cGU6OlJlZ2lvblR5cGU6Oklu
ZGV4VHlwZSBvdXRwdXRJbmRleDsKICBJbWFnZVR5cGU6OlJlZ2lvblR5cGU6
OlNpemVUeXBlIGlucHV0U2l6ZTsKICBJbWFnZVR5cGU6OlJlZ2lvblR5cGU6
OlNpemVUeXBlIG91dHB1dFNpemU7CgogIGRpbXNbMF0gPSA1MTI7IGRpbXNb
MV0gPSA1MTI7IGRpbXNbMl0gPSAyNTE7CiAgc3BhY2luZ1swXSA9IDU3NDsg
c3BhY2luZ1sxXSA9IDU3NDsgc3BhY2luZ1syXSA9IDEwMDA7CiAgb3JpZ2lu
WzBdID0gMDsgb3JpZ2luWzFdID0gMDsgb3JpZ2luWzJdID0gMDsKICB6b29t
WzBdID0gODsgem9vbVsxXSA9IDg7IHpvb21bMl0gPSA1OwoKICBmb3IoaSA9
IDA7IGkgPCBWSW1hZ2VEaW1lbnNpb247IGkgKyspewoJc21hbGxEaW1zW2ld
ID0gZGltc1tpXS96b29tW2ldOwoJc21hbGxTcGFjaW5nW2ldID0gc3BhY2lu
Z1tpXSp6b29tW2ldOwoJc21hbGxPcmlnaW5baV0gPSBvcmlnaW5baV07CiAg
fQoKICBmb3IoaSA9IDA7IGkgPCBWSW1hZ2VEaW1lbnNpb247IGkgKyspewoJ
aW5wdXRJbmRleFtpXSA9IDA7CglvdXRwdXRJbmRleFtpXSA9IDA7CglpbnB1
dFNpemVbaV0gPSAoaW50KWRpbXNbaV07CglvdXRwdXRTaXplW2ldID0gKGlu
dClzbWFsbERpbXNbaV07CiAgfQoKICBpbnB1dFJlZ2lvbi5TZXRJbmRleChp
bnB1dEluZGV4KTsKICBpbnB1dFJlZ2lvbi5TZXRTaXplKGlucHV0U2l6ZSk7
CiAgb3V0cHV0UmVnaW9uLlNldEluZGV4KG91dHB1dEluZGV4KTsKICBvdXRw
dXRSZWdpb24uU2V0U2l6ZShvdXRwdXRTaXplKTsKICBpbnB1dEltYWdlLT5T
ZXRSZWdpb25zKGlucHV0UmVnaW9uKTsKICBpbnB1dEltYWdlLT5TZXRTcGFj
aW5nKHNwYWNpbmcpOwogIGlucHV0SW1hZ2UtPlNldE9yaWdpbihvcmlnaW4p
OwogIGlucHV0SW1hZ2UtPkFsbG9jYXRlKCk7CgogIHN0YXJ0YyA9IGNsb2Nr
KCk7CiAgc3RkOjpjZXJyIDw8ICJTdGFydCBkb3duIHNhbXBsaW5nOiAiIDw8
IAoJKGZsb2F0KShjbG9jaygpLXN0YXJ0YykvQ0xPQ0tTX1BFUl9TRUMgPDwg
c3RkOjplbmRsOwoKICByZXNhbXBsZUZpbHRlci0+U2V0T3V0cHV0U3BhY2lu
ZyhzbWFsbFNwYWNpbmcpOwogIHJlc2FtcGxlRmlsdGVyLT5TZXRPdXRwdXRP
cmlnaW4oc21hbGxPcmlnaW4pOwogIHJlc2FtcGxlRmlsdGVyLT5TZXRTaXpl
KG91dHB1dFNpemUpOwogIHJlc2FtcGxlRmlsdGVyLT5TZXRJbnB1dChpbnB1
dEltYWdlKTsKICByZXNhbXBsZUZpbHRlci0+VXBkYXRlKCk7CiAgc3RkOjpj
ZXJyIDw8ICJFbmQgZG93biBzYW1wbGluZzogIiA8PCAKCShmbG9hdCkoY2xv
Y2soKS1zdGFydGMpL0NMT0NLU19QRVJfU0VDIDw8IHN0ZDo6ZW5kbDsKCiAg
b3V0cHV0SW1hZ2UgPSByZXNhbXBsZUZpbHRlci0+R2V0T3V0cHV0KCk7Cgog
IHN0ZDo6Y2VyciA8PCAiU3RhcnQgdXAgc2FtcGxpbmc6ICIgPDwgCgkoZmxv
YXQpKGNsb2NrKCktc3RhcnRjKS9DTE9DS1NfUEVSX1NFQyA8PCBzdGQ6OmVu
ZGw7CgogIHJlc2FtcGxlRmlsdGVyLT5TZXRJbnB1dChvdXRwdXRJbWFnZSk7
CiAgcmVzYW1wbGVGaWx0ZXItPlNldFNpemUoaW5wdXRTaXplKTsKICByZXNh
bXBsZUZpbHRlci0+U2V0T3V0cHV0U3BhY2luZyhzcGFjaW5nKTsKICByZXNh
bXBsZUZpbHRlci0+U2V0T3V0cHV0T3JpZ2luKG9yaWdpbik7CiAgcmVzYW1w
bGVGaWx0ZXItPlVwZGF0ZSgpOwoKICBzdGQ6OmNlcnIgPDwgIkVuZCB1cCBz
YW1wbGluZzogIiA8PCAKCShmbG9hdCkoY2xvY2soKS1zdGFydGMpL0NMT0NL
U19QRVJfU0VDIDw8IHN0ZDo6ZW5kbDsKCn0KCgkJCQkJCQkJIAo=


------=_NextPart_000_d0d_7581_74e1--