[Insight-users] itkSpatialObjectToImageStatisticsCalculator problem

[Benjamin King]

------------gx3jN7VRnCEODLqg66lYrZ
Content-Type: text/plain; charset=iso-8859-15; format=flowed

Hi Julien,

> Your code looks good.
> I think the problem comes from the spacing. Can you try to remove the 
> line:
>
>> ellipse->SetSpacing(imageReader->GetOutput()->GetSpacing());

indeed there was a problem with the spacing. It wasn't initialized when I 
called imageReader->GetOutput()->GetSpacing() because I didn't call 
imageReader->Update() before that.
But I'm afraid that my problem still isn't fixed. For a lot of parameters, 
I still get "-1.#IND 0" as output.
It doesn't help to remove the line either. Also I'd like to keep the 
spacing of the image and the ellipse consistent. It may not be obligatory, 
but it just makes me feel better =)

I clarified my code a bit, added some comments and fixed a misconception 
that I had about the meaning of spacing.
I also added some debugging code that writes an Analyze file (Debug.hdr + 
Debug.img) containing the sphere I intended to do statistics with. As my 
mailclient messed with the inlined code, this time I made it an attachment.

If you like, I can also send the anisotropic dataset (spacing 0.7, 0.7, 
2.0) and a list of parameters. That would be a 5MB zipfile.

Thanks for your help,
  Benjamin




-- 
Benjamin King
Institut für Medizinische Informatik
Medizinische Hochschule Hannover
Tel.: +49  511  532-2663
------------gx3jN7VRnCEODLqg66lYrZ
Content-Disposition: attachment;
	filename="volumestatistics.cxx"
Content-Type: application/octet-stream;
	name="volumestatistics.cxx"
Content-Transfer-Encoding: Base64

I2luY2x1ZGUgIml0a0ltYWdlLmgiDQojaW5jbHVkZSAiaXRrSW1hZ2VGaWxl
UmVhZGVyLmgiDQojaW5jbHVkZSAiaXRrU3BhdGlhbE9iamVjdC5oIg0KI2lu
Y2x1ZGUgIml0a0VsbGlwc2VTcGF0aWFsT2JqZWN0LmgiDQojaW5jbHVkZSAi
aXRrU3BhdGlhbE9iamVjdFRvSW1hZ2VTdGF0aXN0aWNzQ2FsY3VsYXRvci5o
Ig0KI2luY2x1ZGUgIml0a0FmZmluZVRyYW5zZm9ybS5oIg0KI2luY2x1ZGUg
Iml0a0ltYWdlRmlsZVdyaXRlci5oIg0KI2luY2x1ZGUgIml0a1NwYXRpYWxP
YmplY3RUb0ltYWdlRmlsdGVyLmgiDQojaW5jbHVkZSAiaXRrQmluYXJ5RnVu
Y3RvckltYWdlRmlsdGVyLmgiDQojaW5jbHVkZSAiaXRrTXVsdGlwbHlJbWFn
ZUZpbHRlci5oIg0KI2luY2x1ZGUgPGlvc3RyZWFtPg0KDQp1c2luZyBzdGQ6
OmNvdXQ7DQp1c2luZyBzdGQ6OmVuZGw7DQoNCmludCBtYWluKGludCBhcmdj
LCBjaGFyKiogYXJndikNCnsNCiNkZWZpbmUgRElNRU5TSU9OIDMNCiAgaWYg
KGFyZ2MgIT0gNikgew0KICAgIGNvdXQgPDwgIlxuVm9sdW1lIFN0YXRpc3Rp
Y3NcblxuIg0KICAgICAgIlRoaXMgcHJvZ3JhbSBjb21wdXRlcyB0aGUgbWVh
biBhbmQgc3RhbmRhcmQgZGV2aWF0aW9uIG9mIFxuIg0KICAgICAgImEgc3Bo
ZXJpY2FsIHJlZ2lvbiBpbiBhICIgPDwgRElNRU5TSU9OIDw8ICIgZGltZW5z
aW9uYWwgaW1hZ2UuXG5cbiINCiAgICAgICJVU0FHRTpcbiINCiAgICAgICIg
IHZvbHVtZXN0YXRpc3RpY3MgPGlucHV0IGZpbGU+IDx4PiA8eT4gPHo+IDxy
YWRpdXM+XG4iDQogICAgICAiICAgIDxpbnB1dCBmaWxlPjogIiA8PCBESU1F
TlNJT048PCAiZCAxNmJpdCBBbmFseXplL0RpY29tIGZpbGVcbiINCiAgICAg
ICIgICAgPHg+OiAgICAgICAgICB4IGNvb3JkaW5hdGUgb2YgdGhlIHNwaGVy
ZSdzIGNlbnRlclxuIg0KICAgICAgIiAgICA8eT46ICAgICAgICAgIHkgY29v
cmRpbmF0ZSBvZiB0aGUgc3BoZXJlJ3MgY2VudGVyXG4iDQogICAgICAiICAg
IDx6PjogICAgICAgICAgeiBjb29yZGluYXRlIG9mIHRoZSBzcGhlcmUncyBj
ZW50ZXJcbiINCiAgICAgICIgICAgPHJhZGl1cz46ICAgICByYWRpdXMgb2Yg
dGhlIHNwaGVyZVxuXG4iDQogICAgICAiSElOVDogRGVmaW5lIHRoZSBjZW50
ZXIgY29vcmRpbmF0ZSBhbmQgdGhlIHJhZGl1cyB3aXRoIHJlc3BlY3QgdG8g
dGhlXG4iDQogICAgICAiICBzcGFjaW5nIGFuZCBvcmlnaW4gb2YgeW91ciBp
bnB1dCBpbWFnZS5cblxuIg0KICAgICAgIk9VVFBVVDpcbiINCiAgICAgICIg
IEEgc2luZ2xlIGxpbmUgd2l0aCB0d28gZmxvYXRpbmcgcG9pbnQgbnVtYmVy
cyA8bWVhbj4gYW5kIDxzdGFuZGFyZCBkZXZpYXRpb24+XG4iDQogICAgICAi
ICBzZXBhcmF0ZWQgYnkgYSBzcGFjZSINCiAgICAgIDw8IGVuZGw7DQoNCiAg
ICByZXR1cm4gMTsNCiAgfQ0KICANCiAgdHlwZWRlZiBpdGs6OkltYWdlPCB1
bnNpZ25lZCBzaG9ydCwgRElNRU5TSU9OID4gSW1hZ2VUeXBlOw0KICB0eXBl
ZGVmIGl0azo6SW1hZ2VGaWxlUmVhZGVyPCBJbWFnZVR5cGUgPiBSZWFkZXJU
eXBlOw0KICBSZWFkZXJUeXBlOjpQb2ludGVyIGltYWdlUmVhZGVyID0gUmVh
ZGVyVHlwZTo6TmV3KCk7DQogIGltYWdlUmVhZGVyLT5TZXRGaWxlTmFtZShh
cmd2WzFdKTsNCg0KLy8gREVCVUc6IFByaW50IHNwYWNpbmcNCi8vICB7DQov
LyAgICBpbWFnZVJlYWRlci0+VXBkYXRlKCk7DQovLyAgICBjb25zdCBkb3Vi
bGUgKnNwYWNpbmcgPSBpbWFnZVJlYWRlci0+R2V0T3V0cHV0KCktPkdldFNw
YWNpbmcoKTsNCi8vICAgIGNvdXQgPDwgIlNwYWNpbmc6ICI7DQovLyAgICBm
b3IgKGludCBpID0gMDsgaSA8IERJTUVOU0lPTjsgKytpKQ0KLy8gICAgICBj
b3V0IDw8IHNwYWNpbmdbaV0gPDwgIiAiOw0KLy8gICAgY291dCA8PCBlbmRs
Ow0KLy8gIH0NCg0KICAvLyBnZW5lcmF0ZSBhIHNwaGVyZSBjZW50ZXJlZCBh
dCAoMCwgMCwgMCkNCiAgdHlwZWRlZiBpdGs6OkVsbGlwc2VTcGF0aWFsT2Jq
ZWN0PCBESU1FTlNJT04+IEVsbGlwc2VUeXBlOw0KICBFbGxpcHNlVHlwZTo6
UG9pbnRlciBlbGxpcHNlID0gRWxsaXBzZVR5cGU6Ok5ldygpOw0KICAvLyBU
byBzZXQgdGhlIHNwYWNpbmcgb2YgdGhlIGVsbGlwc2UsIGltYWdlUmVhZGVy
IGhhcyB0byBiZSB1cGRhdGVkLg0KICBpbWFnZVJlYWRlci0+VXBkYXRlKCk7
DQogIGNvbnN0IGRvdWJsZSAqc3BhY2luZyA9IGltYWdlUmVhZGVyLT5HZXRP
dXRwdXQoKS0+R2V0U3BhY2luZygpOw0KICBlbGxpcHNlLT5TZXRTcGFjaW5n
KHNwYWNpbmcpOw0KDQogIC8vIFNldCB0aGUgcmFkaXVzDQogIGNvbnN0IGRv
dWJsZSByYWQgPSBhdG9mKGFyZ3ZbNV0pOw0KICBFbGxpcHNlVHlwZTo6QXJy
YXlUeXBlIHJhZGl1czsNCiAgZm9yIChpbnQgaSA9IDA7IGkgPCBESU1FTlNJ
T047ICsraSkNCiAgICByYWRpdXNbaV0gPSByYWQgLyAoc3BhY2luZ1tpXSAq
IHNwYWNpbmdbaV0pOw0KICBlbGxpcHNlLT5TZXRSYWRpdXMocmFkaXVzKTsN
Cg0KICAvLyBTaGlmdCBpdCB0byB0aGUgcmlnaHQgcGxhY2UNCiAgRWxsaXBz
ZVR5cGU6OlZlY3RvclR5cGUgb2Zmc2V0Ow0KICBvZmZzZXRbMF0gPSBhdG9m
KGFyZ3ZbMl0pIC8gc3BhY2luZ1swXTsNCiAgb2Zmc2V0WzFdID0gYXRvZihh
cmd2WzNdKSAvIHNwYWNpbmdbMV07DQogIG9mZnNldFsyXSA9IGF0b2YoYXJn
dls0XSkgLyBzcGFjaW5nWzJdOw0KICBlbGxpcHNlLT5HZXRJbmRleFRvT2Jq
ZWN0VHJhbnNmb3JtKCktPlNldE9mZnNldChvZmZzZXQpOw0KICBlbGxpcHNl
LT5Db21wdXRlT2JqZWN0VG9QYXJlbnRUcmFuc2Zvcm0oKTsNCg0KICAvLyBj
b21wdXRlIHRoZSBzdGF0aXN0aWNzDQogIHR5cGVkZWYgaXRrOjpTcGF0aWFs
T2JqZWN0VG9JbWFnZVN0YXRpc3RpY3NDYWxjdWxhdG9yPCBJbWFnZVR5cGUs
IEVsbGlwc2VUeXBlID4gU3RhdGlzdGljc1R5cGU7DQogIFN0YXRpc3RpY3NU
eXBlOjpQb2ludGVyIHN0YXRpc3RpYyA9IFN0YXRpc3RpY3NUeXBlOjpOZXco
KTsNCiAgc3RhdGlzdGljLT5TZXRJbWFnZShpbWFnZVJlYWRlci0+R2V0T3V0
cHV0KCkpOw0KICBzdGF0aXN0aWMtPlNldFNwYXRpYWxPYmplY3QoZWxsaXBz
ZSk7DQogIHN0YXRpc3RpYy0+VXBkYXRlKCk7DQoNCiAgY291dCA8PCBzdGF0
aXN0aWMtPkdldE1lYW4oKVswXSA8PCAiICIgPDwgc3FydChzdGF0aXN0aWMt
PkdldENvdmFyaWFuY2VNYXRyaXgoKVswXVswXSkgPDwgZW5kbDsNCg0KLy8g
REVCVUc6IEdlbmVyYXRlIGFuZCB3cml0ZSBhbiBpbWFnZSBvZiB0aGUgdm9s
dW1lIGVuY2xvc2VkIGJ5IHRoZSBlbGxpcHNlDQogIHR5cGVkZWYgaXRrOjpT
cGF0aWFsT2JqZWN0VG9JbWFnZUZpbHRlcjwgRWxsaXBzZVR5cGUsIEltYWdl
VHlwZSA+IE9iamVjdFRvSW1hZ2VUeXBlOw0KICBPYmplY3RUb0ltYWdlVHlw
ZTo6UG9pbnRlciBvYmplY3RUb0ltYWdlID0gT2JqZWN0VG9JbWFnZVR5cGU6
Ok5ldygpOw0KICBvYmplY3RUb0ltYWdlLT5TZXRJbnB1dChlbGxpcHNlKTsN
CiAgb2JqZWN0VG9JbWFnZS0+U2V0SW5zaWRlVmFsdWUoMSk7DQogIG9iamVj
dFRvSW1hZ2UtPlNldE91dHNpZGVWYWx1ZSgwKTsNCiAgb2JqZWN0VG9JbWFn
ZS0+U2V0U2l6ZShpbWFnZVJlYWRlci0+R2V0T3V0cHV0KCktPkdldExhcmdl
c3RQb3NzaWJsZVJlZ2lvbigpLkdldFNpemUoKSk7DQogIG9iamVjdFRvSW1h
Z2UtPlVwZGF0ZSgpOw0KDQogIHR5cGVkZWYgaXRrOjpNdWx0aXBseUltYWdl
RmlsdGVyPCBJbWFnZVR5cGUsIEltYWdlVHlwZSwgSW1hZ2VUeXBlID4gTXVs
dGlwbHlUeXBlOw0KICBNdWx0aXBseVR5cGU6OlBvaW50ZXIgbXVsdGlwbHlG
aWx0ZXIgPSBNdWx0aXBseVR5cGU6Ok5ldygpOw0KICBtdWx0aXBseUZpbHRl
ci0+U2V0SW5wdXQxKGltYWdlUmVhZGVyLT5HZXRPdXRwdXQoKSk7DQogIG11
bHRpcGx5RmlsdGVyLT5TZXRJbnB1dDIob2JqZWN0VG9JbWFnZS0+R2V0T3V0
cHV0KCkpOw0KICBtdWx0aXBseUZpbHRlci0+VXBkYXRlKCk7DQoNCiAgdHlw
ZWRlZiBpdGs6OkltYWdlRmlsZVdyaXRlcjxJbWFnZVR5cGU+IFdyaXRlclR5
cGU7DQogIFdyaXRlclR5cGU6OlBvaW50ZXIgaW1hZ2VXcml0ZXIgPSBXcml0
ZXJUeXBlOjpOZXcoKTsNCiAgaW1hZ2VXcml0ZXItPlNldElucHV0KG11bHRp
cGx5RmlsdGVyLT5HZXRPdXRwdXQoKSk7DQogIGltYWdlV3JpdGVyLT5TZXRG
aWxlTmFtZSgiRGVidWcuaGRyIik7DQogIGltYWdlV3JpdGVyLT5VcGRhdGUo
KTsNCg0KICByZXR1cm4gMDsNCn0NCg==
------------gx3jN7VRnCEODLqg66lYrZ
Content-Disposition: attachment;
	filename="CMakeLists.txt"
Content-Type: text/plain;
	name="CMakeLists.txt"

# This project is intended to be built outside the Insight source tree
PROJECT(volumestatistics)

# Find ITK.
FIND_PACKAGE(ITK)
IF(ITK_FOUND)
  INCLUDE(${ITK_USE_FILE})
ELSE(ITK_FOUND)
  MESSAGE(FATAL_ERROR
          "Cannot build without ITK.  Please set ITK_DIR.")
ENDIF(ITK_FOUND)

ADD_EXECUTABLE(volumestatistics volumestatistics.cxx )

TARGET_LINK_LIBRARIES(volumestatistics ITKCommon ITKIO ITKBasicFilters)


------------gx3jN7VRnCEODLqg66lYrZ--