[Insight-users] questions about generate and display sparse point cloud

[Shaohua Fan]

--0-1225210933-1065115754=:61169
Content-Type: text/plain; charset=us-ascii
Content-Id: 
Content-Disposition: inline

Hello,

I am playing with some samples in ITK, and trying to
generate a set of sparse points and display it. I got
it run somehow, but have two questions:

1. When I displayed the 3D images in raw format with
VTK, it was automatically display with mesh instead of
just points. Can I turn the mesh display off? Or
should I output those points into a different format?

2. In the code(attached, line 199), if I initialize
the whole 3D volume with value -1.0, the generated
semi-sphere would show up. However, when I initialize
it with value 0, the whole image is black. Why it
makes difference?

Anyone could help and have a quick look and point out
the problem? Honestly, I am not sure whether I am
generating/outputing the sparse points in a good way.
It would be appreciated if someone point me to some
sample codes or how to do this.

Thanks a lot,

Shaohua

__________________________________
Do you Yahoo!?
The New Yahoo! Shopping - with improved product search
http://shopping . yahoo . com
--0-1225210933-1065115754=:61169
Content-Type: application/octet-stream; name="sparsePoints.cxx"
Content-Transfer-Encoding: base64
Content-Description: sparsePoints.cxx
Content-Disposition: attachment; filename="sparsePoints.cxx"

Lyo9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQoNCiAgUHJvZ3JhbTog
ICBJbnNpZ2h0IFNlZ21lbnRhdGlvbiAmIFJlZ2lzdHJhdGlvbiBUb29sa2l0
DQogIE1vZHVsZTogICAgJFJDU2ZpbGU6IFJHQlBvaW50U2V0LmN4eCx2ICQN
CiAgTGFuZ3VhZ2U6ICBDKysNCiAgRGF0ZTogICAgICAkRGF0ZTogMjAwMy8w
MS8xOSAwNDoyOTo0MyAkDQogIFZlcnNpb246ICAgJFJldmlzaW9uOiAxLjQg
JA0KDQogIENvcHlyaWdodCAoYykgMjAwMiBJbnNpZ2h0IENvbnNvcnRpdW0u
IEFsbCByaWdodHMgcmVzZXJ2ZWQuDQogIFNlZSBJVEtDb3B5cmlnaHQudHh0
IG9yIGh0dHA6Ly93d3cuaXRrLm9yZy9IVE1ML0NvcHlyaWdodC5odG0gZm9y
IGRldGFpbHMuDQoNCiAgICAgVGhpcyBzb2Z0d2FyZSBpcyBkaXN0cmlidXRl
ZCBXSVRIT1VUIEFOWSBXQVJSQU5UWTsgd2l0aG91dCBldmVuIA0KICAgICB0
aGUgaW1wbGllZCB3YXJyYW50eSBvZiBNRVJDSEFOVEFCSUxJVFkgb3IgRklU
TkVTUyBGT1IgQSBQQVJUSUNVTEFSIA0KICAgICBQVVJQT1NFLiAgU2VlIHRo
ZSBhYm92ZSBjb3B5cmlnaHQgbm90aWNlcyBmb3IgbW9yZSBpbmZvcm1hdGlv
bi4NCg0KPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PSovDQoNCiNpbmNs
dWRlIDxmc3RyZWFtPg0KI2luY2x1ZGUgIml0a0ltYWdlUmVnaW9uSXRlcmF0
b3IuaCINCiNpbmNsdWRlICJpdGtMZXZlbFNldEZ1bmN0aW9uLmgiDQojaW5j
bHVkZSAiaXRrU3BhcnNlRmllbGRMZXZlbFNldEltYWdlRmlsdGVyLmgiDQoj
aW5jbHVkZSAiaXRrUmF3SW1hZ2VJTy5oIg0KI2luY2x1ZGUgIml0a0ltYWdl
RmlsZVdyaXRlci5oIg0KI2luY2x1ZGUgIml0a1JHQlBpeGVsLmgiDQojaW5j
bHVkZSAiaXRrUG9pbnRTZXQuaCINCi8vIFNvZnR3YXJlIEd1aWRlIDogQmVn
aW5Db2RlU25pcHBldA0KI2luY2x1ZGUgIml0a1JHQlBpeGVsLmgiDQojaW5j
bHVkZSAiaXRrUG9pbnRTZXQuaCINCg0KZXh0ZXJuICJDIiB7DQojaW5jbHVk
ZSAic3RkaW8uaCINCn0NCg0KY29uc3QgdW5zaWduZWQgaW50IFZfV0lEVEgg
ID0gNjQ7DQpjb25zdCB1bnNpZ25lZCBpbnQgVl9IRUlHSFQgPSA2NDsNCmNv
bnN0IHVuc2lnbmVkIGludCBWX0RFUFRIICA9IDY0Ow0KDQojZGVmaW5lIFRf
UkFESVVTIChWX1dJRFRILzEwLjApDQojZGVmaW5lIFhZX1JBRElVUyAoVl9X
SURUSC8zLjApDQoNCiNkZWZpbmUgVF9DRU5URVJfWCAoMC4wKQ0KI2RlZmlu
ZSBUX0NFTlRFUl9ZICgwLjApDQojZGVmaW5lIFRfQ0VOVEVSX1ogKDAuMCkN
Cg0KI2RlZmluZSBTUEhfUl9YICgwLjJmKlZfV0lEVEgpDQojZGVmaW5lIFNQ
SF9SX1kgKDAuMmYqVl9IRUlHSFQpDQojZGVmaW5lIFNQSF9SX1ogKDAuMmYq
Vl9ERVBUSCkNCg0KZmxvYXQgc3BoZXJlKHVuc2lnbmVkIHgsIHVuc2lnbmVk
IHksIHVuc2lnbmVkIHopDQp7DQogIGZsb2F0IGRpczsNCiAgZGlzID0gKHgg
LSAoZmxvYXQpVl9XSURUSC8yLjApKih4IC0gKGZsb2F0KVZfV0lEVEgvMi4w
KQ0KICAgIC8oU1BIX1JfWCpTUEhfUl9YKSArICh5IC0gKGZsb2F0KVZfSEVJ
R0hULzIuMCkqKHkgLSAoZmxvYXQpVl9IRUlHSFQvMi4wKQ0KICAgIC8oU1BI
X1JfWSpTUEhfUl9ZKSArICh6IC0gKGZsb2F0KVZfREVQVEgvMi4wKSooeiAt
IChmbG9hdClWX0RFUFRILzIuMCkNCiAgICAvKFNQSF9SX1oqU1BIX1JfWik7
DQoNCiAgaWYoeD09NTAmJnk9PTUwJiZ6PT01MCkgcHJpbnRmKCIoJWQsICVk
LCAlZCkgLS0gJWZcbiIsIHgsIHksIHosIDEuMC1kaXMpOw0KICByZXR1cm4o
MS4wZi1kaXMpOw0KfQ0KDQoNCnZvaWQgZXZhbHVhdGVfZnVuY3Rpb24oaXRr
OjpJbWFnZTxmbG9hdCwgMz4gKmltLA0KICAgICAgICAgICAgICAgICAgICAg
ICBmbG9hdCAoKmYpKHVuc2lnbmVkIGludCwgdW5zaWduZWQgaW50LCB1bnNp
Z25lZCBpbnQpICkNCiAgDQp7DQogIGl0azo6SW1hZ2U8ZmxvYXQsIDM+OjpJ
bmRleFR5cGUgaWR4Ow0KICANCiAgZm9yKHVuc2lnbmVkIGludCB6ID0gMDsg
eiA8IFZfREVQVEg7ICsreikNCiAgICB7DQogICAgICBpZHhbMl0gPSB6Ow0K
ICAgICAgZm9yICh1bnNpZ25lZCBpbnQgeSA9IDA7IHkgPCBWX0hFSUdIVDsg
Kyt5KQ0KICAgICAgICB7DQogICAgICAgICAgaWR4WzFdID0geTsNCiAgICAg
ICAgICBmb3IgKHVuc2lnbmVkIGludCB4ID0gMDsgeCA8IFZfV0lEVEg7ICsr
eCkNCiAgICAgICAgICAgIHsNCiAgICAgICAgICAgICAgaWR4WzBdID0geDsN
CiAgICAgICAgICAgICAgaW0tPlNldFBpeGVsKGlkeCwgZih4LCB5LCB6KSAp
Ow0KICAgICAgICAgICAgfQ0KICAgICAgICB9DQogICAgICANCiAgICB9DQp9
DQoNCg0KaW50IG1haW4oKQ0Kew0KDQoNCiAgdHlwZWRlZiBmbG9hdCAgICBQ
aXhlbFR5cGU7DQoNCiAgdHlwZWRlZiBpdGs6OlBvaW50U2V0PCBQaXhlbFR5
cGUsIDMgPiBQb2ludFNldFR5cGU7DQogIFBvaW50U2V0VHlwZTo6UG9pbnRl
ciAgcG9pbnRTZXQgPSBQb2ludFNldFR5cGU6Ok5ldygpOw0KDQogIFBvaW50
U2V0VHlwZTo6UGl4ZWxUeXBlICAgcGl4ZWw7DQogIFBvaW50U2V0VHlwZTo6
UG9pbnRUeXBlICAgcG9pbnQ7DQogIHVuc2lnbmVkIGludCBwb2ludElkID0g
IDA7DQogIGNvbnN0IGRvdWJsZSByYWRpdXMgPSBTUEhfUl9YOw0KDQoJLy9n
ZW5lcmF0ZSBzcGFyc2UgcG9pbnRzIG9uIHNlbWktc3BoZXJlDQogIGZvcih1
bnNpZ25lZCBpbnQgaT0wOyBpPDE4MDsgaSsrKQ0KICAgIHsNCiAgICBjb25z
dCBkb3VibGUgYW5nbGUxID0gaSAqIGF0YW4oMS4wKSAvIDQ1LjA7DQoNCgkJ
Zm9yKHVuc2lnbmVkIGludCBqPTA7IGo8OTA7IGorKyl7DQoJCQkNCgkJCWRv
dWJsZSBhbmdsZTIgPSBqKmF0YW4oMS4wKSAvIDQ1LjA7DQoNCgkJCXBvaW50
WzBdID0gU1BIX1JfWCArIHJhZGl1cyAqIGNvcyhhbmdsZTIpICogc2luKCBh
bmdsZTEgKTsNCgkJCXBvaW50WzFdID0gU1BIX1JfWSArIHJhZGl1cyAqIGNv
cyhhbmdsZTIpICogY29zKCBhbmdsZTEgKTsNCgkJCXBvaW50WzJdID0gU1BI
X1JfWiArIHJhZGl1cyAqIHNpbihhbmdsZTIpOw0KCQkJcGl4ZWwgPSAyNTU7
IA0KCQkJcG9pbnRTZXQtPlNldFBvaW50KCBwb2ludElkLCBwb2ludCApOyAg
IA0KCQkJcG9pbnRTZXQtPlNldFBvaW50RGF0YSggcG9pbnRJZCwgcGl4ZWwg
KTsgICANCgkJCXBvaW50SWQrKzsNCgkJfQ0KICAgIH0NCg0KCS8vcHJpbnQg
b3V0IHRoZSBkYXRhDQogIHR5cGVkZWYgIFBvaW50U2V0VHlwZTo6UG9pbnRz
Q29udGFpbmVyOjpDb25zdEl0ZXJhdG9yICAgICBQb2ludEl0ZXJhdG9yOw0K
ICBQb2ludEl0ZXJhdG9yIHBvaW50SXRlcmF0b3IgPSBwb2ludFNldC0+R2V0
UG9pbnRzKCktPkJlZ2luKCk7DQogIFBvaW50SXRlcmF0b3IgcG9pbnRFbmQg
ICAgICA9IHBvaW50U2V0LT5HZXRQb2ludHMoKS0+RW5kKCk7DQogIHdoaWxl
KCBwb2ludEl0ZXJhdG9yICE9IHBvaW50RW5kICkgew0KICAgIFBvaW50U2V0
VHlwZTo6UG9pbnRUeXBlIHBvaW50ID0gcG9pbnRJdGVyYXRvci5WYWx1ZSgp
Ow0KICAgIHN0ZDo6Y291dCA8PCBwb2ludCA8PCBzdGQ6OmVuZGw7ICANCiAg
ICArK3BvaW50SXRlcmF0b3I7ICAgICAgICAgICAgICAgICAgICAgIA0KICAg
IH0NCg0KICB0eXBlZGVmICBQb2ludFNldFR5cGU6OlBvaW50RGF0YUNvbnRh
aW5lcjo6Q29uc3RJdGVyYXRvciAgICAgUG9pbnREYXRhSXRlcmF0b3I7DQog
IFBvaW50RGF0YUl0ZXJhdG9yIHBpeGVsSXRlcmF0b3IgPSBwb2ludFNldC0+
R2V0UG9pbnREYXRhKCktPkJlZ2luKCk7DQogIFBvaW50RGF0YUl0ZXJhdG9y
IHBpeGVsRW5kICAgICAgPSBwb2ludFNldC0+R2V0UG9pbnREYXRhKCktPkVu
ZCgpOw0KICB3aGlsZSggcGl4ZWxJdGVyYXRvciAhPSBwaXhlbEVuZCApIHsN
CiAgICBQb2ludFNldFR5cGU6OlBpeGVsVHlwZSBwaXhlbCA9IHBpeGVsSXRl
cmF0b3IuVmFsdWUoKTsNCiAgICBzdGQ6OmNvdXQgPDwgcGl4ZWwgPDwgc3Rk
OjplbmRsOyAgDQogICAgKytwaXhlbEl0ZXJhdG9yOyAgICAgICAgICAgICAg
ICAgICAgICANCiAgICB9DQoNCg0KCS8vb3V0cHV0IHRoZSBzcGFyc2UgcG9p
bnRzIGludG8gYW4gaW1hZ2UNCiAgdHlwZWRlZiBpdGs6OkltYWdlPGZsb2F0
LDM+IEltYWdlVHlwZTsNCg0KICB0cnkgew0KICAgIGludCBuOw0KICAgIA0K
ICAgIGl0azo6UmF3SW1hZ2VJTzxmbG9hdCwgMz46OlBvaW50ZXIgb3V0cHV0
X2lvDQogICAgICA9IGl0azo6UmF3SW1hZ2VJTzxmbG9hdCwgMz46Ok5ldygp
Ow0KDQogICAgb3V0cHV0X2lvLT5TZXRCeXRlT3JkZXJUb0xpdHRsZUVuZGlh
bigpOw0KICAgIG91dHB1dF9pby0+U2V0RmlsZVR5cGVUb0JpbmFyeSgpOw0K
ICAgIG91dHB1dF9pby0+U2V0RmlsZURpbWVuc2lvbmFsaXR5KDMpOw0KICAg
IA0KICAgIGl0azo6SW1hZ2VGaWxlV3JpdGVyPEltYWdlVHlwZT46OlBvaW50
ZXIgd3JpdGVyDQogICAgICA9IGl0azo6SW1hZ2VGaWxlV3JpdGVyPEltYWdl
VHlwZT46Ok5ldygpOyANCiAgICB3cml0ZXItPlNldEltYWdlSU8ob3V0cHV0
X2lvKTsNCiAgICANCiAgICAvLyBJbml0aWFsaXplIHRoZSBpbml0aWFsIGFu
ZCB0YXJnZXQgdm9sdW1lcy4NCiAgICBJbWFnZVR5cGU6OlBvaW50ZXIgdm9s
X2luaXQgPSBJbWFnZVR5cGU6Ok5ldygpOw0KICAgIEltYWdlVHlwZTo6UG9p
bnRlciB2b2xfdGFyZ2V0ID0gSW1hZ2VUeXBlOjpOZXcoKTsNCiAgICANCiAg
ICBJbWFnZVR5cGU6OlNpemVUeXBlIHN6Ow0KICAgIHN6WzBdID0gVl9XSURU
SDsNCiAgICBzelsxXSA9IFZfSEVJR0hUOw0KICAgIHN6WzJdID0gVl9ERVBU
SDsNCiAgICANCiAgICBJbWFnZVR5cGU6OkluZGV4VHlwZSBpZHg7DQogICAg
aWR4WzBdID0gMDsNCiAgICBpZHhbMV0gPSAwOw0KICAgIGlkeFsyXSA9IDA7
DQogICAgDQogICAgSW1hZ2VUeXBlOjpSZWdpb25UeXBlIHJlZzsgIA0KICAg
IHJlZy5TZXRTaXplKHN6KTsNCiAgICByZWcuU2V0SW5kZXgoaWR4KTsNCiAg
ICANCiAgICB2b2xfaW5pdC0+U2V0UmVnaW9ucyhyZWcpOw0KICAgIHZvbF9p
bml0LT5BbGxvY2F0ZSgpOw0KIA0KCQkvL3NldCB0aGUgd2hvbGUgaW1hZ2Ug
dm9sdW1lIHRvIHplcm8gZmlyc3QNCgkJZm9yKHVuc2lnbmVkIGludCB6ID0g
MDsgeiA8IFZfREVQVEg7ICsreikNCgkJCXsNCgkJCQlpZHhbMl0gPSB6Ow0K
CQkJCWZvciAodW5zaWduZWQgaW50IHkgPSAwOyB5IDwgVl9IRUlHSFQ7ICsr
eSkNCgkJCQkJew0KCQkJCQkJaWR4WzFdID0geTsNCgkJCQkJCWZvciAodW5z
aWduZWQgaW50IHggPSAwOyB4IDwgVl9XSURUSDsgKyt4KQ0KCQkJCQkJCXsN
CgkJCQkJCQkJaWR4WzBdID0geDsNCi8vCQkJCQkJCQl2b2xfaW5pdC0+U2V0
UGl4ZWwoaWR4LCBzcGhlcmUoeCx5LHopICk7IA0KCQkJCQkJCQl2b2xfaW5p
dC0+U2V0UGl4ZWwoaWR4LCAtMS4wICk7IA0KCQkJCQkJCQlpZih4PT01MCYm
eT09NTAmJno9PTUwKSBwcmludGYoIiglZCwgJWQsICVkKTogJWZcZiIsIHgs
IHksIHosIHZvbF9pbml0LT5HZXRQaXhlbChpZHgpKTsNCgkJCQkJCQl9DQoJ
CQkJCX0NCiAgICAgIA0KCQkJfQ0KDQoNCgkJLy9hc3NpZ24gcGl4ZWwgdmFs
dWVzIGZvciB0aGUgcG9pbnQgb24gdGhlIHNwaGVyZQ0KCQl0eXBlZGVmICBQ
b2ludFNldFR5cGU6OlBvaW50c0NvbnRhaW5lcjo6Q29uc3RJdGVyYXRvciAg
ICAgUG9pbnRJdGVyYXRvcjsNCgkJUG9pbnRJdGVyYXRvciBwb2ludEl0ZXJh
dG9yID0gcG9pbnRTZXQtPkdldFBvaW50cygpLT5CZWdpbigpOw0KCQlQb2lu
dEl0ZXJhdG9yIHBvaW50RW5kICAgICAgPSBwb2ludFNldC0+R2V0UG9pbnRz
KCktPkVuZCgpOw0KCQl0eXBlZGVmICBQb2ludFNldFR5cGU6OlBvaW50RGF0
YUNvbnRhaW5lcjo6Q29uc3RJdGVyYXRvciAgICAgUG9pbnREYXRhSXRlcmF0
b3I7DQoJCVBvaW50RGF0YUl0ZXJhdG9yIHBpeGVsSXRlcmF0b3IgPSBwb2lu
dFNldC0+R2V0UG9pbnREYXRhKCktPkJlZ2luKCk7DQoJCVBvaW50RGF0YUl0
ZXJhdG9yIHBpeGVsRW5kICAgICAgPSBwb2ludFNldC0+R2V0UG9pbnREYXRh
KCktPkVuZCgpOw0KCQl3aGlsZSggcG9pbnRJdGVyYXRvciAhPSBwb2ludEVu
ZCAmJiAgcGl4ZWxJdGVyYXRvciAhPSBwaXhlbEVuZCkgey8vdGhleSBzaG91
bGQgZW5kIGF0IHRoZSBzYW1lIHRpbWUNCgkJCVBvaW50U2V0VHlwZTo6UG9p
bnRUeXBlIHBvaW50ID0gcG9pbnRJdGVyYXRvci5WYWx1ZSgpOw0KCSAgICBQ
b2ludFNldFR5cGU6OlBpeGVsVHlwZSBwaXhlbCA9IHBpeGVsSXRlcmF0b3Iu
VmFsdWUoKTsNCg0KCQkJaWR4WzJdID0gcG9pbnRbMl07IC8vdG8gZ2V0IHRo
ZSB6IHZhbHVlIG9mIHBvaW50IGNvb3JkDQoJCQlpZHhbMV0gPSBwb2ludFsx
XTsNCgkJCWlkeFswXSA9IHBvaW50WzBdOw0KDQoJCQl2b2xfaW5pdC0+U2V0
UGl4ZWwoaWR4LCAyNTUgKTsgDQoNCgkJCSsrcG9pbnRJdGVyYXRvcjsgICAg
ICAgICAgICAgICAgICAgICAgDQoJICAgICsrcGl4ZWxJdGVyYXRvcjsgICAg
ICAgICAgICAgICAgICAgICAgDQoJCX0NCg0KICAgIA0KCQkvL3dyaXRlIG91
dCB0aGUgc3BhcnNlIDNEIHBvaW50cw0KICAgIHdyaXRlci0+U2V0RmlsZU5h
bWUoInNwYXJzZV9zcGhlcmVQb2ludHNfM0QucmF3Iik7DQogICAgd3JpdGVy
LT5TZXRJbnB1dCh2b2xfaW5pdCk7DQogICAgd3JpdGVyLT5Xcml0ZSgpOw0K
CX0gY2F0Y2ggKGl0azo6RXhjZXB0aW9uT2JqZWN0ICplKQ0KICAgIHsNCiAg
ICAgIHN0ZDo6Y291dCA8PCBlIDw8IHN0ZDo6ZW5kbDsNCiAgICB9DQogICAN
Cg0KICByZXR1cm4gMDsNCg0KfQ0KDQoNCg0K

--0-1225210933-1065115754=:61169--