No subject
Wed Oct 7 22:37:18 EDT 2009
However, I just created and image using a Matlab code below.
***** Star Matlab code ******
% Create single square 1340x1300x720
clc;
clear all;
close all;
matrix = zeros(1340, 1300, 720);
matrix(600:809, 600:809, :) = 1000;
% imshow(matrix);
fid = fopen('Bar_3Dsquare_I1000_1340x1300x720.bin', 'wb');
fwrite(fid, matrix, 'float');
status = fclose(fid);
******* Stop Matlab code ************
This code exactly creates an 4.7 GB image with white square extruded along
the Z-derection.
Now I use the ITK program ImageReadRegionOfInterestWrite to sequentially
extract series of 100 slices as below.
dzyubak at helium: /Square$ ./ImageReadRegionOfInterestWrite
Bar_3Dsquare_I1000_1340x1300x720_float.hdr
Bar_3Dsquare_I1000_1340x1300x720_float_s1_s100.hdr 0 0 0 1340 1300 100
and
dzyubak at helium: /Square$ ./ImageReadRegionOfInterestWrite
Bar_3Dsquare_I1000_1340x1300x720_float.hdr
Bar_3Dsquare_I1000_1340x1300x720_float_s100_s200.hdr 0 0 100 1340 1300 100
I checked the Bar_3Dsquare_I1000_1340x1300x720_float.hdr image and it is OK,
the square goes all the way through.
The first set is also OK and the "white square" shows up but the second set
is empty. Any ideas what could cause this problem?
Thanks,
Alex
******** Start *******
root at helium: /BUILD# ctest -R LargeImageWrite -V
Test command: /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests
itkLargeImageWriteReadTest
/mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImage01.mhd
30000L
Test timeout computed to be: 1500
Trying to allocate an image of size 1716 Mb
Initializing pixel values
Trying to write the image to disk
Trying to read the image back from disk
Comparing the pixel values.. :
Probe Tag Starts Stops Time (s)
Allocate 1 1 3.19481e-05
Compare 1 1 4.85127
Initializing 1 1 8.84029
Read 1 1 2.71803
Write 1 1 4.17017
Test PASSED !
-- Process completed
Passed
863/1701 Testing itkLargeImageWriteReadTest2
Test command: /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests
itkLargeImageWriteReadTest
/mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImage02.mhd
40000L
Test timeout computed to be: 1500
Trying to allocate an image of size 3051 Mb
Initializing pixel values
Trying to write the image to disk
Trying to read the image back from disk
Comparing the pixel values.. :
Probe Tag Starts Stops Time (s)
Allocate 1 1 3.00407e-05
Compare 1 1 7.89822
Initializing 1 1 12.6574
Read 1 1 4.8622
Write 1 1 12.6875
Test PASSED !
-- Process completed
Passed
864/1701 Testing itkLargeImageWriteConvertReadT
Test command: /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests
itkLargeImageWriteConvertReadTest
/mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImageConvert.nrrd
50000L
Test timeout computed to be: 1500
Initializing pixel values
Trying to write the image to disk
Trying to read the image back from disk
Probe Tag Starts Stops Time (s)
Allocate 1 1 3.69549e-05
Initializing 1 1 11.7354
Read 1 1 15.7987
Write 1 1 10.5514
Test PASSED !
865/1701 Testing itkLargeImageWriteReadTest3
Test command: /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests
itkLargeImageWriteReadTest
/mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImage03.mhd
70000L
Test timeout computed to be: 1500
Trying to allocate an image of size 9346 Mb
Initializing pixel values
Trying to write the image to disk
Trying to read the image back from disk
Comparing the pixel values.. :
Probe Tag Starts Stops Time (s)
Allocate 1 1 3.19481e-05
Compare 1 1 24.2805
Initializing 1 1 40.3377
Read 1 1 14.4607
Write 1 1 32.6605
Test PASSED !
The following tests passed:
itkLargeImageWriteReadTest1
itkLargeImageWriteReadTest2
itkLargeImageWriteConvertReadTest
itkLargeImageWriteReadTest3
100% tests passed, 0 tests failed out of 4
********* End ********
On Fri, Dec 11, 2009 at 6:50 PM, Luis Ibanez <luis.ibanez at kitware.com>wrote:
> Hi Oleksandr,
>
>
> The components of the itk::Index are of type "long",
> while the components of the itk::Size are of type
> "unsigned long".
>
> Therefore, On Linux 64 bits, you can manage
> images up to the size of your memory.
>
> With your small machine of 16 CPUs and
> only 128 Gb of RAM you should do fine.
>
> ... and yes, I'm saying that with envy :-)
>
>
> On Windows 64 bits, on the other hand
> we are working on fixing a problem with
> the integer representation used in Windows.
>
>
> Back to your report on Linux:
>
> In the directory Insight/Testing/Code/IO
> we have a test for large memory that
> will attempt to create, write and read back
> and image of large size (above 2 Gb).
> Let's try running it in your computer.
>
> Please do the following:
>
> 1) go to the binary directory where you
> built ITK and do
>
> make edit_cache
>
>
> 2) Turn BUILD_TESTING ON
>
>
> 3) Go to the Advanced mode,
> search for
>
> ITK_COMPUTER_MEMORY_SIZE
>
> and set it to 128
>
> (yes, the units are in Gb).
>
>
>
> 4) Do, "c" for configure
>
> 5) Do, "g" to generate
>
>
> 6) make -j32
>
> (that shouldn't take long) :-)
>
>
> 7) Run the LargeImage test with
> the command:
>
> ctest -R LargeImageWrite -V
>
> and see if they pass.
>
> 8) Now you can take it further,
> by looking at the command
> that is used to run the tests,
> The last numeric argument is
> the number of pixels along each
> dimension, of the image that is
> created by the test.
>
> You could force this number
> until the total size of this 2D
> image is comparable to your 3D
> image.
>
> For example, by 70,000 the
> image will be of size 1,2 Gb.
>
>
>
> Please let us know what you find
> after running these tests.
>
>
> Thanks
>
>
> Luis
>
>
> --------------------------------------------------------------
> On Fri, Dec 11, 2009 at 4:37 PM, Oleksandr Dzyubak <adzyubak at gmail.com>
> wrote:
> >
> > Dear users,
> >
> > I tried to process the 1340x1300x720 pixel image and after the slice #104
> I
> > am getting all zeros.
> > I checked the original image and it's good. It's full of info.
> > So I first decided to split that file using the
> > itkRegionOfInterestImageFilter class.
> > What I found is that all the chunks after the slice #104 are zeros
> > thus it's not the image processing filter that causes the problem.
> >
> > Here is my question.
> > Is there any limitation on an image index range?
> > Or in other words, what is the largest pixel volume that ITK still can
> > handle?
> >
> > Thanks,
> >
> > Alex
> >
> > PS.
> > I am using 64 bit Debian Linux.
> > 16 CPUs and 128 GB.
> >
> > _____________________________________
> > Powered by www.kitware.com
> >
> > Visit other Kitware open-source projects at
> > http://www.kitware.com/opensource/opensource.html
> >
> > Kitware offers ITK Training Courses, for more information visit:
> > http://www.kitware.com/products/protraining.html
> >
> > Please keep messages on-topic and check the ITK FAQ at:
> > http://www.itk.org/Wiki/ITK_FAQ
> >
> > Follow this link to subscribe/unsubscribe:
> > http://www.itk.org/mailman/listinfo/insight-users
> >
> >
>
--000e0cd154aa105d29047ab4554d
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
Hi Luis,<br><br>I followed the procedure 1) through 8) that you suggested.<=
br>All four tests have successfully passed and the results are below.<br>Ac=
cording to the results, the parameter 70,000L provides the image size <br>
"Trying to allocate an image of size 9346 Mb " which is what I ac=
tually need.<br>From these tests it looks like it should work.<br><br>Howev=
er, I just created and image using a Matlab code below.<br><br>***** Star M=
atlab code ******<br>
<br>% Create single square 1340x1300x720
<br>clc;
<br>clear all;
<br>close all;
<br>=A0<br>=A0<br>matrix =3D zeros(1340, 1300, 720);
<br>matrix(600:809, 600:809, :) =3D 1000;
<br>=A0<br>% imshow(matrix);
<br>=A0<br>fid =3D fopen('Bar_3Dsquare_I1000_1340x1300x720.bin', &#=
39;wb');
<br>fwrite(fid, matrix, 'float');
<br>status =3D fclose(fid);<br><br>******* Stop Matlab code ************<br=
><br>This code exactly creates an 4.7 GB image with white square extruded a=
long the Z-derection.<br>Now I use the ITK program ImageReadRegionOfInteres=
tWrite to sequentially extract series of 100 slices as below.<br>
<br>dzyubak at helium: /Square$ ./ImageReadRegionOfInterestWrite Bar_3Dsquare_=
I1000_1340x1300x720_float.hdr Bar_3Dsquare_I1000_1340x1300x720_float_s1_s10=
0.hdr 0 0 0 1340 1300 100<br><br>and<br><br>dzyubak at helium: /Square$ ./Imag=
eReadRegionOfInterestWrite Bar_3Dsquare_I1000_1340x1300x720_float.hdr Bar_3=
Dsquare_I1000_1340x1300x720_float_s100_s200.hdr 0 0 100 1340 1300 100<br>
<br>I checked the Bar_3Dsquare_I1000_1340x1300x720_float.hdr image and it i=
s OK, the square goes all the way through.<br>The first set is also OK and =
the "white square" shows up but the second set is empty. Any idea=
s what could cause this problem?<br>
<br>Thanks,<br><br>Alex<br><br><br><br>******** Start *******<br><br><br>ro=
ot at helium: /BUILD# ctest -R LargeImageWrite -V <br><br>Test command: /mnt/h=
elium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests itkLargeImageWriteReadTest =
/mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImage01.mhd 3=
0000L<br>
Test timeout computed to be: 1500<br>Trying to allocate an image of size 17=
16 Mb <br>Initializing pixel values <br>Trying to write the image to disk<b=
r>Trying to read the image back from disk<br>Comparing the pixel values.. :=
<br>
=A0=A0=A0=A0=A0=A0=A0=A0=A0 Probe Tag=A0=A0=A0 Starts=A0=A0=A0 Stops=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Time (s)<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0 Allocate=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0 1=A0=A0=A0=A0=A0=A0 3.19481e-05<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
Compare=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 4.85127<br>=A0=A0=A0=A0=A0=A0=A0 Initializin=
g=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0 8.84029<br>
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Read=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 2=
.71803<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Write=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0 4.17017<br><br>Test PASSED !<br>-- Process completed<br>=A0=A0 Pa=
ssed<br>863/1701 Testing itkLargeImageWriteReadTest2=A0=A0 <br>
Test command: /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests itkLarg=
eImageWriteReadTest /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporar=
y/LargeImage02.mhd 40000L<br>Test timeout computed to be: 1500<br>Trying to=
allocate an image of size 3051 Mb <br>
Initializing pixel values <br>Trying to write the image to disk<br>Trying t=
o read the image back from disk<br>Comparing the pixel values.. :<br>=A0=A0=
=A0=A0=A0=A0=A0=A0=A0 Probe Tag=A0=A0=A0 Starts=A0=A0=A0 Stops=A0=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0 Time (s)<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Allo=
cate=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=
=A0=A0=A0=A0=A0 3.00407e-05<br>
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Compare=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =
1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 7.89822<=
br>=A0=A0=A0=A0=A0=A0=A0 Initializing=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 12.6574<br>=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Read=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 4=
.8622<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Write=A0=A0=A0=A0=A0=A0=
=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0 12.6875<br>
<br>Test PASSED !<br>-- Process completed<br>=A0=A0 Passed<br>864/1701 Test=
ing itkLargeImageWriteConvertReadT<br>Test command: /mnt/helium/Root_BUILDs=
/ITK_3_16/BUILD/bin/itkIOTests itkLargeImageWriteConvertReadTest /mnt/heliu=
m/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImageConvert.nrrd 50000=
L<br>
Test timeout computed to be: 1500<br>Initializing pixel values <br>Trying t=
o write the image to disk<br>Trying to read the image back from disk<br>=A0=
=A0=A0=A0=A0=A0=A0=A0=A0 Probe Tag=A0=A0=A0 Starts=A0=A0=A0 Stops=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0 Time (s)<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 A=
llocate=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=
=A0=A0=A0=A0=A0=A0 3.69549e-05<br>
=A0=A0=A0=A0=A0=A0=A0 Initializing=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 11.7354<br>=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Read=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 =
1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 15.7987<=
br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Write=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=
0.5514<br><br>Test PASSED !<br>
<br>865/1701 Testing itkLargeImageWriteReadTest3=A0=A0 <br>Test command: /m=
nt/helium/Root_BUILDs/ITK_3_16/BUILD/bin/itkIOTests itkLargeImageWriteReadT=
est /mnt/helium/Root_BUILDs/ITK_3_16/BUILD/Testing/Temporary/LargeImage03.m=
hd 70000L<br>
Test timeout computed to be: 1500<br>Trying to allocate an image of size 93=
46 Mb <br>Initializing pixel values <br>Trying to write the image to disk<b=
r>Trying to read the image back from disk<br>Comparing the pixel values.. :=
<br>
=A0=A0=A0=A0=A0=A0=A0=A0=A0 Probe Tag=A0=A0=A0 Starts=A0=A0=A0 Stops=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Time (s)<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0 Allocate=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
=A0 1=A0=A0=A0=A0=A0=A0 3.19481e-05<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=
Compare=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 24.2805<br>=A0=A0=A0=A0=A0=A0=A0 Initializin=
g=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=
=A0=A0=A0=A0=A0=A0=A0=A0 40.3377<br>
=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Read=A0=A0=A0=A0=A0=A0=A0=A0=
=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=
4.4607<br>=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 Write=A0=A0=A0=A0=A0=
=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0=A0 1=A0=A0=A0=A0=A0=A0=A0=
=A0=A0=A0 32.6605<br><br>Test PASSED !<br><br>The following tests passed:<b=
r>=A0=A0=A0 itkLargeImageWriteReadTest1<br>
=A0=A0=A0 itkLargeImageWriteReadTest2<br>=A0=A0=A0 itkLargeImageWriteConver=
tReadTest<br>=A0=A0=A0 itkLargeImageWriteReadTest3<br><br>100% tests passed=
, 0 tests failed out of 4<br><br><br>********* End ********<br><br><br><br>=
<div class=3D"gmail_quote">
On Fri, Dec 11, 2009 at 6:50 PM, Luis Ibanez <span dir=3D"ltr"><<a href=
=3D"mailto:luis.ibanez at kitware.com">luis.ibanez at kitware.com</a>></span> =
wrote:<br><blockquote class=3D"gmail_quote" style=3D"border-left: 1px solid=
rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Oleksandr,<br>
<br>
<br>
The components of the itk::Index are of type "long",<br>
while the components of the itk::Size are of type<br>
"unsigned long".<br>
<br>
Therefore, On Linux 64 bits, you can manage<br>
images up to the size of your memory.<br>
<br>
With your small machine of 16 CPUs and<br>
only 128 Gb of RAM you should do fine.<br>
<br>
... and yes, I'm saying that with envy =A0 :-)<br>
<br>
<br>
On Windows 64 bits, on the other hand<br>
we are working on fixing a problem with<br>
the integer representation used in Windows.<br>
<br>
<br>
Back to your report on Linux:<br>
<br>
In the directory Insight/Testing/Code/IO<br>
we have a test for large memory that<br>
will attempt to create, write and read back<br>
and image of large size (above 2 Gb).<br>
Let's try running it in your computer.<br>
<br>
Please do the following:<br>
<br>
1) go to the binary directory where you<br>
=A0 =A0built ITK and do<br>
<br>
=A0 =A0 =A0 =A0make =A0 =A0edit_cache<br>
<br>
<br>
2) Turn BUILD_TESTING =A0 ON<br>
<br>
<br>
3) Go to the Advanced mode,<br>
=A0 =A0search for<br>
<br>
ITK_COMPUTER_MEMORY_SIZE<br>
<br>
=A0 =A0and set it to 128<br>
<br>
=A0(yes, the units are in Gb).<br>
<br>
<br>
<br>
4) Do, "c" for configure<br>
<br>
5) Do, "g" to generate<br>
<br>
<br>
6) =A0make -j32<br>
<br>
=A0 =A0 (that shouldn't take long) =A0 =A0:-)<br>
<br>
<br>
7) Run the LargeImage test with<br>
=A0 =A0the command:<br>
<br>
=A0 =A0 =A0ctest -R LargeImageWrite -V<br>
<br>
=A0 =A0and see if they pass.<br>
<br>
8) =A0Now you can take it further,<br>
=A0 =A0 by looking at the command<br>
=A0 =A0 that is used to run the tests,<br>
=A0 =A0 The last numeric argument is<br>
=A0 =A0 the number of pixels along each<br>
=A0 =A0 dimension, of the image that is<br>
=A0 =A0 created by the test.<br>
<br>
=A0 =A0 =A0You could force this number<br>
=A0 =A0 =A0until the total size of this 2D<br>
=A0 =A0 =A0image is comparable to your 3D<br>
=A0 =A0 =A0image.<br>
<br>
=A0 =A0 =A0For example, by 70,000 the<br>
=A0 =A0 =A0image will be of size 1,2 Gb.<br>
<div class=3D"im"><br>
<br>
<br>
Please let us know what you find<br>
</div>after running these tests.<br>
<br>
<br>
=A0 =A0Thanks<br>
<br>
<br>
=A0 =A0 =A0 =A0 =A0Luis<br>
<br>
<br>
--------------------------------------------------------------<br>
<div><div></div><div class=3D"h5">On Fri, Dec 11, 2009 at 4:37 PM, Oleksand=
r Dzyubak <<a href=3D"mailto:adzyubak at gmail.com">adzyubak at gmail.com</a>&=
gt; wrote:<br>
><br>
> Dear users,<br>
><br>
> I tried to process the 1340x1300x720 pixel image and after the slice #=
104 I<br>
> am getting all zeros.<br>
> I checked the original image and it's good. It's full of info.=
<br>
> So I first decided to split that file using the<br>
> itkRegionOfInterestImageFilter class.<br>
> What I found is that all the chunks after the slice #104 are zeros<br>
> thus it's not the image processing filter that causes the problem.=
<br>
><br>
> Here is my question.<br>
> Is there any limitation on an image index range?<br>
> Or in other words, what is the largest pixel volume that ITK still can=
<br>
> handle?<br>
><br>
> Thanks,<br>
><br>
> Alex<br>
><br>
> PS.<br>
> I am using 64 bit Debian Linux.<br>
> 16 CPUs and 128 GB.<br>
><br>
</div></div><div><div></div><div class=3D"h5">> ________________________=
_____________<br>
> Powered by <a href=3D"http://www.kitware.com" target=3D"_blank">www.ki=
tware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href=3D"http://www.kitware.com/opensource/opensource.html" target=
=3D"_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Kitware offers ITK Training Courses, for more information visit:<br>
> <a href=3D"http://www.kitware.com/products/protraining.html" target=3D=
"_blank">http://www.kitware.com/products/protraining.html</a><br>
><br>
> Please keep messages on-topic and check the ITK FAQ at:<br>
> <a href=3D"http://www.itk.org/Wiki/ITK_FAQ" target=3D"_blank">http://w=
ww.itk.org/Wiki/ITK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href=3D"http://www.itk.org/mailman/listinfo/insight-users" target=
=3D"_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
><br>
><br>
</div></div></blockquote></div><br>
--000e0cd154aa105d29047ab4554d--
More information about the Insight-users
mailing list