[Insight-users] Problem with ImportImageFilter

Harri Tapio Jaalinoja harri.jaalinoja@helsinki.fi
Mon, 11 Nov 2002 12:44:26 +0200 (EET) 

  This message is in MIME format.  The first part should be readable text,
  while the remaining parts are likely unreadable without MIME-aware tools.
  Send mail to mime@docserver.cac.washington.edu for more info.

--0-461582285-1037011466=:20163
Content-Type: TEXT/PLAIN; charset=US-ASCII


Hello all!

I am using BSoft function to read in an image file. The function call
gives me a struct with a pointer to the data.

  Bimage* p = read_img(argv[1], 1, -1);

  if(p->datatype != 7)
    die("Only float data supported at the moment.");

  typedef float PixelType;

  PixelType* data = (PixelType *)p->data;
  std::cerr << data[0] << std::endl;
  std::cerr << data[9898345] << std::endl;
  std::cerr << data[16777215] << std::endl;

  unsigned long num = p->x * p->y * p->z;
  std::cerr << num << std::endl;

These values print out ok, so the reading part seems to be working.
---- output ----
-0.160718
-0.286178
-0.160718
16777216
---- output ----

Next I want to pass the data on to ITK, using either ImportImageFilter or
ImportImageContainer.

  itk::ImportImageFilter<PixelType, 3>::Pointer breader =
    itk::ImportImageFilter<PixelType, 3>::New();
  breader->SetImportPointer(data ,num, false);

This does not work. I add
  std::cerr << "Try to print a pixel from reader output..." << std::endl;
  itk::Image<PixelType, 3>::Pointer image = breader->GetOutput();
  std::cerr << image->GetPixel(index) << std::endl;
  std::cerr << "OK" << std::endl;

and get a segmentation fault. index points to pixel 9898345 (I tried
{1,1,1}, too).

---- output ----
Try to print a pixel from reader output...
Segmentation fault
---- output ----

To me this indicates that SetImportPointer did not do
what I expected it to. What am I doing wrong?

Thanks,

Harri



--0-461582285-1037011466=:20163
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="bitkShiftScale.cxx"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.OSF.4.30.0211111244260.20163@sirppi.helsinki.fi>
Content-Description: 
Content-Disposition: attachment; filename="bitkShiftScale.cxx"

Lyo9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQoNCiAgUHJvZ3JhbTog
ICBJbnNpZ2h0IFNlZ21lbnRhdGlvbiAmIFJlZ2lzdHJhdGlvbiBUb29sa2l0
DQogIE1vZHVsZTogICAgJFJDU2ZpbGU6IGl0a1dhdGVyc2hlZFNlZ21lbnRh
dGlvbkV4YW1wbGUuY3h4LHYgJA0KICBMYW5ndWFnZTogIEMrKw0KICBEYXRl
OiAgICAgICREYXRlOiAyMDAyLzA4LzI1IDE1OjUyOjExICQNCiAgVmVyc2lv
bjogICAkUmV2aXNpb246IDEuOSAkDQoNCiAgQ29weXJpZ2h0IChjKSAyMDAy
IEluc2lnaHQgQ29uc29ydGl1bS4gQWxsIHJpZ2h0cyByZXNlcnZlZC4NCiAg
U2VlIElUS0NvcHlyaWdodC50eHQgb3IgaHR0cDovL3d3dy5pdGsub3JnL0hU
TUwvQ29weXJpZ2h0Lmh0bSBmb3IgZGV0YWlscy4NCg0KICAgICBUaGlzIHNv
ZnR3YXJlIGlzIGRpc3RyaWJ1dGVkIFdJVEhPVVQgQU5ZIFdBUlJBTlRZOyB3
aXRob3V0IGV2ZW4gDQogICAgIHRoZSBpbXBsaWVkIHdhcnJhbnR5IG9mIE1F
UkNIQU5UQUJJTElUWSBvciBGSVRORVNTIEZPUiBBIFBBUlRJQ1VMQVIgDQog
ICAgIFBVUlBPU0UuICBTZWUgdGhlIGFib3ZlIGNvcHlyaWdodCBub3RpY2Vz
IGZvciBtb3JlIGluZm9ybWF0aW9uLg0KDQo9PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09Ki8NCiNpZiBkZWZpbmVkKF9NU0NfVkVSKQ0KI3ByYWdtYSB3
YXJuaW5nICggZGlzYWJsZSA6IDQ3ODYgKQ0KI2VuZGlmDQojaW5jbHVkZSA8
ZnN0cmVhbT4NCiNpbmNsdWRlICJpdGtJbXBvcnRJbWFnZUNvbnRhaW5lci5o
Ig0KI2luY2x1ZGUgIml0a0ltcG9ydEltYWdlRmlsdGVyLmgiDQojaW5jbHVk
ZSAiaXRrU2hpZnRTY2FsZUltYWdlRmlsdGVyLmgiDQojaW5jbHVkZSAiaXRr
UmF3SW1hZ2VJTy5oIg0KI2luY2x1ZGUgIml0a0ltYWdlRmlsZVJlYWRlci5o
Ig0KI2luY2x1ZGUgIml0a0ltYWdlRmlsZVdyaXRlci5oIg0KI2luY2x1ZGUg
PHN0cmluZz4NCiNpbmNsdWRlICJpdGtDb21tYW5kLmgiDQoNCiNpbmNsdWRl
IDxyd01SQy5oPg0KI2luY2x1ZGUgPHN0cmluZy5oPg0KI2luY2x1ZGUgPHN0
ZGlvLmg+DQoNCg0Kc3RhdGljIGlubGluZSB2b2lkIGRpZShzdGQ6OnN0cmlu
ZyBtZXNzYWdlKSB7IA0KICBzdGQ6OmNlcnIgPDwgImJXYXRlcnNoZWQgdjAu
MCIgPDwgc3RkOjplbmRsOw0KICBzdGQ6OmNlcnIgPDwgbWVzc2FnZSA8PCBz
dGQ6OmVuZGw7IA0KICBzdGQ6OmNlcnIgPDwgIkV4aXRpbmcuLiIgPDwgc3Rk
OjplbmRsOw0KICBleGl0KDEpOyANCn0NCg0Kc3RhdGljIGlubGluZSB2b2lk
IGRpZShzdGQ6OnN0cmluZyBtZXNzYWdlLCBpdGs6OkV4Y2VwdGlvbk9iamVj
dCBlKSB7IA0KICAgICAgbWVzc2FnZSA9ICJcbiI7DQogICAgICBtZXNzYWdl
ICs9IGUuR2V0TG9jYXRpb24oKTsNCiAgICAgIG1lc3NhZ2UgKz0gIlxuIjsN
CiAgICAgIG1lc3NhZ2UgKz0gZS5HZXREZXNjcmlwdGlvbigpOw0KICAgICAg
ZGllKG1lc3NhZ2UpOw0KfQ0KDQppbnQgbWFpbihpbnQgYXJnYywgY2hhciAq
YXJndltdKQ0Kew0KICBmbG9hdCBzaGlmdDsNCiAgZmxvYXQgc2NhbGU7DQoN
CiAgc3RkOjpzdHJpbmcgb3V0bmFtZTsNCiAgc3RkOjpzdHJpbmcgbWVzc2Fn
ZTsNCg0KICBpZiAoYXJnYyAhPSA0KQ0KICAgIGRpZSAoIlVzZTogYml0a1No
aWZ0U2NhbGUgaW5wdXRfZmlsZSBzaGlmdCBzY2FsZSIpOw0KICBzc2NhbmYo
YXJndlsyXSwgIiVmIiwgJnNoaWZ0KTsNCiAgc3NjYW5mKGFyZ3ZbM10sICIl
ZCIsICZzY2FsZSk7DQoNCiAgc3RkOjpjZXJyIDw8ICJyZWFkaW5nLi4uIiA8
PCBzdGQ6OmVuZGw7DQogIEJpbWFnZSogcCA9IHJlYWRfaW1nKGFyZ3ZbMV0s
IDEsIC0xKTsNCi8vICAgc3RkOjpjZXJyIDw8ICJkb25lISIgPDwgc3RkOjpl
bmRsOw0KLy8gICBzdGQ6OmNlcnIgPDwgcC0+eCA8PCBzdGQ6OmVuZGw7DQov
LyAgIHN0ZDo6Y2VyciA8PCBwLT55IDw8IHN0ZDo6ZW5kbDsNCi8vICAgc3Rk
OjpjZXJyIDw8IHAtPnogPDwgc3RkOjplbmRsOw0KLy8gICBzdGQ6OmNlcnIg
PDwgcC0+YyA8PCBzdGQ6OmVuZGw7DQovLyAgIHN0ZDo6Y2VyciA8PCBwLT5t
aW4gPDwgc3RkOjplbmRsOw0KLy8gICBzdGQ6OmNlcnIgPDwgcC0+bWF4IDw8
IHN0ZDo6ZW5kbDsNCi8vICAgc3RkOjpjZXJyIDw8IHAtPmF2ZyA8PCBzdGQ6
OmVuZGw7DQovLyAgIHN0ZDo6Y2VyciA8PCBwLT5zdGQgPDwgc3RkOjplbmRs
Ow0KDQogIGlmKHAtPmRhdGF0eXBlICE9IDcpIA0KICAgIGRpZSgiT25seSBm
bG9hdCBkYXRhIHN1cHBvcnRlZCBhdCB0aGUgbW9tZW50LiIpOw0KDQogIHR5
cGVkZWYgZmxvYXQgUGl4ZWxUeXBlOw0KDQogIFBpeGVsVHlwZSogZGF0YSA9
IChQaXhlbFR5cGUgKilwLT5kYXRhOw0KICBzdGQ6OmNlcnIgPDwgZGF0YVsw
XSA8PCBzdGQ6OmVuZGw7DQogIHN0ZDo6Y2VyciA8PCBkYXRhWzk4OTgzNDVd
IDw8IHN0ZDo6ZW5kbDsNCiAgc3RkOjpjZXJyIDw8IGRhdGFbMTY3NzcyMTVd
IDw8IHN0ZDo6ZW5kbDsNCg0KICB1bnNpZ25lZCBsb25nIG51bSA9IHAtPngg
KiBwLT55ICogcC0+ejsNCiAgc3RkOjpjZXJyIDw8IG51bSA8PCBzdGQ6OmVu
ZGw7DQoNCiAgLy8gaW5kZXggdGhhdCBjb3JyZXNwb25kcyB0byA5ODk4NDM1
ID0gMTUxKjI1Nl4yICsgOSoyNTYgKyAxMDUNCiAgaXRrOjpJbmRleDwzPiBp
bmRleCA9IHsxMDUsIDksIDE1MX07DQoNCiAgdHlwZWRlZiBpdGs6OkltYWdl
PFBpeGVsVHlwZSwgMz4gSW1hZ2VUeXBlOw0KICB0eXBlZGVmIGl0azo6SW1h
Z2VGaWxlV3JpdGVyPCBJbWFnZVR5cGUgPiBGaWxlVHlwZTsNCiAgdHlwZWRl
ZiBpdGs6OlJhd0ltYWdlSU88UGl4ZWxUeXBlLCAzPiBJbWFnZUlPVHlwZTsN
Cg0KLy8gICBpdGs6OkltcG9ydEltYWdlQ29udGFpbmVyPGxvbmcgdW5zaWdu
ZWQgaW50LCBQaXhlbFR5cGU+OjpQb2ludGVyIGltcG9ydA0KLy8gICAgID0g
aXRrOjpJbXBvcnRJbWFnZUNvbnRhaW5lcjxsb25nIHVuc2lnbmVkIGludCxQ
aXhlbFR5cGUgPjo6TmV3KCk7DQovLyAgIGltcG9ydC0+SW5pdGlhbGl6ZSgp
Ow0KLy8gICBpbXBvcnQtPlJlc2VydmUobnVtKTsNCi8vICAgaW1wb3J0LT5T
ZXRJbXBvcnRQb2ludGVyKGRhdGEsbnVtLGZhbHNlKTsNCiAgDQovLyAgIGl0
azo6SW1hZ2U8UGl4ZWxUeXBlLCAzPjo6UG9pbnRlciBpbWcgPSBpdGs6Oklt
YWdlPFBpeGVsVHlwZSwgMz46Ok5ldygpOw0KLy8gICBpbWctPlNldFBpeGVs
Q29udGFpbmVyKGltcG9ydCk7DQoNCi8vICAgc3RkOjpjZXJyIDw8IGltZy0+
R2V0UGl4ZWwoaW5kZXgpIDw8IHN0ZDo6ZW5kbDsNCg0KICBpdGs6OkltcG9y
dEltYWdlRmlsdGVyPFBpeGVsVHlwZSwgMz46OlBvaW50ZXIgYnJlYWRlciA9
IA0KICAgIGl0azo6SW1wb3J0SW1hZ2VGaWx0ZXI8UGl4ZWxUeXBlLCAzPjo6
TmV3KCk7DQogIGJyZWFkZXItPlNldEltcG9ydFBvaW50ZXIoZGF0YSAsbnVt
LCBmYWxzZSk7IA0KICANCi8vIC8vICAgLy8gU2V0IHVwIHRoZSBmaWx0ZXIN
Ci8vIC8vICAgaXRrOjpTaGlmdFNjYWxlSW1hZ2VGaWx0ZXI8SW1hZ2VUeXBl
LCBJbWFnZVR5cGU+OjpQb2ludGVyIHNzZmlsdGVyDQovLyAvLyAgICAgPSBp
dGs6OlNoaWZ0U2NhbGVJbWFnZUZpbHRlcjxJbWFnZVR5cGUsIEltYWdlVHlw
ZT46Ok5ldygpOw0KLy8gLy8gICBzc2ZpbHRlci0+U2V0U2hpZnQoc2hpZnQp
Ow0KLy8gLy8gICBzc2ZpbHRlci0+U2V0U2NhbGUoc2NhbGUpOw0KLy8gLy8g
ICBzc2ZpbHRlci0+U2V0SW5wdXQoYnJlYWRlci0+R2V0T3V0cHV0KCkpOw0K
Ly8gLy8gLy8gICBzc2ZpbHRlci0+U2V0SW5wdXQoaW1nKTsNCg0KICBzdGQ6
OmNlcnIgPDwgIlRyeSB0byBwcmludCBhIHBpeGVsIGZyb20gcmVhZGVyIG91
dHB1dC4uLiIgPDwgc3RkOjplbmRsOw0KICBpdGs6OkltYWdlPFBpeGVsVHlw
ZSwgMz46OlBvaW50ZXIgaW1hZ2UgPSBicmVhZGVyLT5HZXRPdXRwdXQoKTsN
CiAgc3RkOjpjZXJyIDw8IGltYWdlLT5HZXRQaXhlbChpbmRleCkgPDwgc3Rk
OjplbmRsOw0KICBzdGQ6OmNlcnIgPDwgIk9LIiA8PCBzdGQ6OmVuZGw7DQoN
CiAgLy8gU2V0IHVwIHRoZSBmaWxlIHdyaXRlcg0KICBJbWFnZUlPVHlwZTo6
UG9pbnRlciBvdXRwdXRfaW8gPSBJbWFnZUlPVHlwZTo6TmV3KCk7DQogIG91
dHB1dF9pby0+U2V0Qnl0ZU9yZGVyVG9MaXR0bGVFbmRpYW4oKTsNCiAgb3V0
cHV0X2lvLT5TZXRGaWxlRGltZW5zaW9uYWxpdHkoMyk7DQogIG91dHB1dF9p
by0+U2V0RmlsZVR5cGVUb0JpbmFyeSgpOw0KDQogIEZpbGVUeXBlOjpQb2lu
dGVyIHdyaXRlciA9IEZpbGVUeXBlOjpOZXcoKTsNCiAgd3JpdGVyLT5TZXRJ
bWFnZUlPKG91dHB1dF9pbyk7DQogIG91dG5hbWUgPSAic3NfIiArIHN0ZDo6
c3RyaW5nKGFyZ3ZbMl0pICsgIi5yYXciOw0KICB3cml0ZXItPlNldEZpbGVO
YW1lKG91dG5hbWUuY19zdHIoKSk7DQogIHdyaXRlci0+U2V0SW5wdXQoYnJl
YWRlci0+R2V0T3V0cHV0KCkpOw0KLy8gLy8gICB3cml0ZXItPlNldElucHV0
KHNzZmlsdGVyLT5HZXRPdXRwdXQoKSk7DQoNCiAgdHJ5IA0KICAgIHsNCiAg
ICAgIHdyaXRlci0+V3JpdGUoKTsNCiAgICB9ICANCiAgY2F0Y2goIGl0azo6
RXhjZXB0aW9uT2JqZWN0ICYgZSApDQogICAgew0KICAgICAgbWVzc2FnZSA9
ICJFeGNlcHRpb24gY2F1Z2h0IGR1cmluZyBwcm9jZXNzaW5nLiBcbiI7DQog
ICAgICBtZXNzYWdlICs9IGUuR2V0TG9jYXRpb24oKTsNCiAgICAgIG1lc3Nh
Z2UgKz0gIlxuIjsNCiAgICAgIG1lc3NhZ2UgKz0gZS5HZXREZXNjcmlwdGlv
bigpOw0KICAgICAgZGllKG1lc3NhZ2UpOw0KICAgIH0NCn0NCg==
--0-461582285-1037011466=:20163--