[Insight-users] Read info extract write bug?

[Kevin H. Hobbs]

I have some images that are fairly large and I do not want to read them
in completely. So I call reader->UpdateOutputInformation() so that I can
get their region. Then I split the region and extract the sub-region
that I want.

With some sub-regions the output of the extraction is junk.

If instead of calling reader->UpdateOutputInformation() I call
reader-Update() then the output of these same sub-regions is fine.

I've attached a view from ParaView of these two cases with Update on the
left and UpdateOutputInformation on the right.

I've also attached the code that generated both images.

Am I doing something wrong or is this a bug?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ExtractBug.cxx
Type: text/x-c++src
Size: 1641 bytes
Desc: not available
URL: <http://www.itk.org/pipermail/insight-users/attachments/20091022/3e6fa1eb/attachment-0001.cxx>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: slices.png
Type: image/png
Size: 20390 bytes
Desc: not available
URL: <http://www.itk.org/pipermail/insight-users/attachments/20091022/3e6fa1eb/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 189 bytes
Desc: This is a digitally signed message part
URL: <http://www.itk.org/pipermail/insight-users/attachments/20091022/3e6fa1eb/attachment-0001.pgp>