<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'><div style="text-align: left;">Hi Luis,<br><br>Thanks for the input.<br><br>I don't have any extra filters apart from the standard ones from the Confidence Connected example. CurvatureFlowImageFilter, ConfidenceConnectedImageFilter, and a CastImageFilter. I played around with the casting a little bit but it didn't seem to effect it (casting before and after the confidence connected filter). Maybe it does have something to do with that though?? I don't know.<br><br>&nbsp; smoothing-&gt;SetInput( reader-&gt;GetOutput() );<br>&nbsp; outputCaster-&gt;SetInput(smoothing-&gt;GetOutput()); &nbsp;&nbsp;  //had this line 3rd previously<br>&nbsp; confidenceConnected-&gt;SetInput( outputCaster-&gt;GetOutput() );<br>&nbsp; writer-&gt;SetInput( confidenceConnected-&gt;GetOutput() );<br><br><br>The type of DICOM file is CT (i presume that's what you meant by modality). The number of pixels, well I'm guessing that it's [7fe0|0000 = 524296] but that's just a wild stab in the dark. I'm presuming you would know exactly what to look for in the header so i've included it. This is the header for the first slice which I managed to extract through ITK-Snap...<br><br>0002|0000 = 190<br>0002|0001 = AAE=<br>0002|0002 = 1.2.840.10008.5.1.4.1.1.2<br>0002|0003 = 1.2.392.200036.9116.2.2.2.1762576288.1175672259.918312<br>0002|0010 = 1.2.840.10008.1.2<br>0002|0012 = 1.2.840.10008.5.1.4.1.1.7.0.1.0<br>0002|0013 = C-PACS<br>0008|0000 = 358<br>0008|0008 = ORIGINAL\PRIMARY\AXIAL<br>0008|0016 = 1.2.840.10008.5.1.4.1.1.2<br>0008|0018 = 1.2.392.200036.9116.2.2.2.1762576288.1175672259.918312<br>0008|0020 = 20070404<br>0008|0022 = 20070404<br>0008|0023 = 20070404<br>0008|0030 = 163042.000<br>0008|0032 = 163307.450<br>0008|0033 = 163309.450<br>0008|0060 = CT<br>0008|0070 = TOSHIBA <br>0008|0080 = Radiol. UK S-H Campus Luebeck <br>0008|1010 = 0001<br>0008|1040 = ID_DEPARTMENT <br>0008|1090 = Aquilion<br>0010|0000 = 46<br>0010|0010 = HUEFTE^^^^<br>0010|0020 = 555 <br>0018|0000 = 212<br>0018|0010 = CE<br>0018|0022 = HELICAL_CT<br>0018|0050 = 1.0 <br>0018|0060 = 120 <br>0018|0090 = 240.00<br>0018|1000 = B0622156<br>0018|1020 = V2.05ER001<br>0018|1100 = 240.00<br>0018|1120 = +0.0<br>0018|1130 = +132.00 <br>0018|1140 = CW<br>0018|1150 = 1000<br>0018|1151 = 230 <br>0018|1152 = 230 <br>0018|1210 = FC30<br>0018|5100 = HFS <br>0020|0000 = 370<br>0020|000d = 1.2.392.200036.9116.2.2.2.1762576288.1175671825.193444<br>0020|000e = 1.2.392.200036.9116.2.2.2.1762576288.1175672259.904331<br>0020|0010 = 3282<br>0020|0011 = 3 <br>0020|0012 = 3 <br>0020|0013 = 1 <br>0020|0020 = L\P <br>0020|0032 = -120.000000\-120.000000\-1198.500000<br>0020|0037 = 1.000000\0.000000\0.000000\0.000000\1.000000\0.000000 <br>0020|0052 = 1.2.392.200036.9116.2.2.2.1762576288.1175671854.644030<br>0020|1041 = -168.00 <br>0028|0000 = 154<br>0028|0002 = 1<br>0028|0004 = MONOCHROME2 <br>0028|0010 = 512<br>0028|0011 = 512<br>0028|0030 = 0.468\0.468 <br>0028|0100 = 16<br>0028|0101 = 16<br>0028|0102 = 15<br>0028|0103 = 1<br>0028|1050 = 300 <br>0028|1051 = 1600<br>0028|1052 = 0 <br>0028|1053 = 1 <br>0040|0000 = 116<br>0040|0002 = 20070404<br>0040|0003 = 163042.000<br>0040|0004 = 20070404<br>0040|0005 = 170042.000<br>0040|0244 = 20070404<br>0040|0245 = 163042.000<br>0040|0253 = 15553 <br>7fe0|0000 = 524296<br><br><br><br><br>Thanks heaps for your help, its greatly appreciated,<br><br>Cameron.<br></div><br><div></div><br><br><hr id="stopSpelling">&gt; Date: Sun, 6 Apr 2008 17:23:56 -0400<br>&gt; From: luis.ibanez@kitware.com<br>&gt; To: w_e_b_m_a_s_t_e_r_6_9@hotmail.com<br>&gt; CC: insight-users@itk.org<br>&gt; Subject: Re: [Insight-users] Problem with Saving a segmented object<br>&gt; <br>&gt; <br>&gt; Hi Cameron,<br>&gt; <br>&gt; It shouldn't be a problem to handle a 170Mb image, either for<br>&gt; reading, writing, and running a region growing segmentation.<br>&gt; <br>&gt;   ...You may have additional filters in the pipeline....<br>&gt; <br>&gt; Could you please describe the full pipeline ?<br>&gt; <br>&gt; and also provide the information of your image:<br>&gt; <br>&gt;    a) number of pixels<br>&gt;    b) Image Modality / pixel type<br>&gt; <br>&gt; <br>&gt; Thanks<br>&gt; <br>&gt; <br>&gt;      Luis<br>&gt; <br>&gt; <br>&gt; -------------------------<br>&gt; Cameron Burnett wrote:<br>&gt; &gt; <br>&gt; &gt; I just tried a smaller DICOM data set, and it worked (well not the way I <br>&gt; &gt; want it to, but it generated a vtk file). Is a 170MB DICOM data set <br>&gt; &gt; normally too big to handle with region growing? Is it possible to <br>&gt; &gt; perhaps buffer some to the file, clear the memory, then do the rest of <br>&gt; &gt; it?? I'm not at all sure what I should do. I'd consider doing it slice <br>&gt; &gt; by slice, but I read somewhere that all the 2D algorithms can be used <br>&gt; &gt; for 3D as well.<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; This is the output i'm getting. I was actually pleased that it did <br>&gt; &gt; SOMETHING, but its obviously not what I want from my segmentation....<br>&gt; &gt; <br>&gt; &gt; http://img371.imageshack.us/img371/6677/dicomtovtk1ju8.jpg<br>&gt; &gt; <br>&gt; &gt; You can see the dicom object and the vtkPointWidget in the background. <br>&gt; &gt; I'm using the point widget for the seed point, but for some reason it <br>&gt; &gt; wants to colour all the way to the boundary of the widget. *shrug*. I'll <br>&gt; &gt; just play around with the variables a bit I guess, unless anyone can <br>&gt; &gt; recommend that I try anything else.<br>&gt; &gt; <br>&gt; &gt; Thanks.<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt;     ------------------------------------------------------------------------<br>&gt; &gt;     From: w_e_b_m_a_s_t_e_r_6_9@hotmail.com<br>&gt; &gt;     To: insight-users@itk.org<br>&gt; &gt;     Date: Fri, 4 Apr 2008 01:11:15 +1000<br>&gt; &gt;     Subject: [Insight-users] Problem with Saving a segmented object<br>&gt; &gt; <br>&gt; &gt;     Hi,<br>&gt; &gt; <br>&gt; &gt;     I've just implemented a mix of the DicomSeriesReadImageWrite2 and<br>&gt; &gt;     ConfidenceConnected algorithm examples for region growing<br>&gt; &gt;     segmentation. So basically, I want to segment a DICOM file and then<br>&gt; &gt;     save it as vtk or a similar file format (Unless there is a better<br>&gt; &gt;     way that anyone knows of).<br>&gt; &gt; <br>&gt; &gt;     The problem is that once I run the segmentation it takes 10 minutes<br>&gt; &gt;     in order to spit an error message back.<br>&gt; &gt; <br>&gt; &gt;     So what is it doing? Well, my program gets through to the execution<br>&gt; &gt;     of the writer pipeline quite quickly (writer-&gt;Update() ). After that<br>&gt; &gt;     the cpu is maxed out at 100% and my RAM is practically maxed out<br>&gt; &gt;     also (I have 2GB, and 360MB is taken up once I load my DICOM file<br>&gt; &gt;     onto the screen via the raycast method). Then after about 10 minutes<br>&gt; &gt;     I get an error in itkProcessObject.cxx. This is in the<br>&gt; &gt;     ProcessObject::UpdateOutputData function, where it calls <br>&gt; &gt;     this-&gt;GenerateData(). When I mouse-over the exception in vc++ it<br>&gt; &gt;     says something about failing to allocate memory for the file.<br>&gt; &gt; <br>&gt; &gt;     I originally thought it had something to do with pixel type casting,<br>&gt; &gt;     because the pre-segmentation smoothing requires a float pixel type<br>&gt; &gt;     and I had just left it as float for the rest also. Then I tried<br>&gt; &gt;     casting it back to "signed short" before the segmentation stage....<br>&gt; &gt;     but that did the same. 10 minutes is a long time to wait to see if a<br>&gt; &gt;     small change works.<br>&gt; &gt; <br>&gt; &gt;     So the main question is.... should it take 10 minutes for a 2GHz<br>&gt; &gt;     dual core cpu (with 2GB ram) to do region growing segmentation (then<br>&gt; &gt;     die, lol). Btw, the DICOM files are about 170MB for the whole 342 of<br>&gt; &gt;     them (about 500KB each file).<br>&gt; &gt; <br>&gt; &gt;     Also, if anyone knows whats going on, or would like to request more<br>&gt; &gt;     information, then that would be great. I'm not expecting much though.<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt;     ------------------------------------------------------------------------<br>&gt; &gt;     Find out: SEEK Salary Centre Are you paid what you're worth?<br>&gt; &gt;     &lt;http://a.ninemsn.com.au/b.aspx?URL=http://ninemsn.seek.com.au/career-resources/salary-centre/?tracking%3Dsk:het:sc:nine:0:hot:text&amp;_t=764565661&amp;_r=OCT07_endtext_salary&amp;_m=EXT&gt;<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; Grab it. You dream job is up for grabs. <br>&gt; &gt; &lt;http://mycareer.com.au/?s_cid=596065 &gt;<br>&gt; &gt; <br>&gt; &gt; <br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; <br>&gt; &gt; _______________________________________________<br>&gt; &gt; Insight-users mailing list<br>&gt; &gt; Insight-users@itk.org<br>&gt; &gt; http://www.itk.org/mailman/listinfo/insight-users<br><br /><hr />Grab it.  <a href='http://mycareer.com.au/?s_cid=596065 ' target='_new'>You dream job is up for grabs.</a></body>
</html>