<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">I want to show saggital side of a CT from a dicom folder. I've three Dicom folder to test the project. Two of them are working ( 35KB * 488 dicom files and 35KB * 341 dicom files) but the last one ( 515KB * 250 dicom files) gives runtime error such that;<br><br><span style="font-weight: bold;">ERROR: In c:\vtk\common\vtkDataArrayTemplate.txx, line 142</span><br style="font-weight: bold;"><span style="font-weight: bold;">vtkIdTypeArray [00A8FD90]: Unable to allocate 131072000 elements of size</span><br><br>what can be the solution ?<br><br>Dicom()<br> {<br> reader = vtkDICOMImageReader::New();<br> viewer = vtkImageViewer::New();<br> iren = vtkRenderWindowInteractor::New();<br> saggital =
vtkImageDataGeometryFilter::New();<br> coronal = vtkImageDataGeometryFilter::New();<br> axial = vtkImageDataGeometryFilter::New();<br> saggitalMapper = vtkPolyDataMapper::New();<br> saggitalActor = vtkActor::New();<br> saggitalRenderer = vtkRenderer::New();<br> saggitalRenWin = vtkRenderWindow::New();<br> saggitalIren = vtkRenderWindowInteractor::New();<br> saggitalCam = vtkCamera::New();<br><br> }<br><br> // To Show saggital axis of CT<br> void ReadDICOMFolder(string DICOMFolder)<br> {<br> viewer->SetInput(reader->GetOutput());<br>
reader->SetDirectoryName(DICOMFolder.c_str());<br> reader->Update();<br><br> saggital->SetInputConnection(reader->GetOutputPort());<br><br> saggitalMapper->SetInput(saggital->GetOutput());<br> saggitalMapper->ScalarVisibilityOn();<br> saggitalMapper->SetScalarRange(1450, 2450);<br><br> saggitalActor->SetMapper(saggitalMapper);<br><br> saggitalRenderer->AddActor(saggitalActor);<br> saggitalRenWin->AddRenderer(saggitalRenderer);<br> saggitalIren->SetRenderWindow(saggitalRenWin);<br> saggitalIren->SetSize(640,
480);<br><br> saggitalCam->SetViewUp(0, 0, -1);<br> saggitalCam->SetPosition(1, 0, 0);<br> saggitalCam->SetFocalPoint(0, 0, 0);<br> saggitalCam->ComputeViewPlaneNormal();<br><br> saggitalRenderer->SetActiveCamera(saggitalCam);<br> saggitalRenderer->ResetCamera();<br> saggitalRenderer->SetBackground(0, 0, 0);<br><br> int dataExtent[6];<br><br> reader->GetDataExtent(dataExtent);<br><br> for(int i = 0; i <= dataExtent[1]; i++)<br> {<br> saggital->SetExtent(i, i, 0,
dataExtent[3], 0, dataExtent[5]);<br> saggitalRenderer->ResetCamera();<br> saggitalRenWin->Render();<br> } <br> }<br></td></tr></table><br>