Hi,<br><br>I`m trying to read set of dicom files, extract one slice
and write it in jpeg file. But the jpeg file is very strange, it
appears in low resolution (or something like that).<br>Could anybody help me? Sorry for my english<br>
<br>My code<br><br><br>typedef itk::Image < unsigned char, 2 > WriteImageType;<br><br>typedef itk::Image < float, 3> ImageType;<br>typedef itk::Image < float, 2> ImageType2D;<br>
<br>typedef itk::ImageToVTKImageFilter <ImageType> ConnectorType;<br>typedef itk::ImageToVTKImageFilter <ImageType2D> ConnectorType2D;<br><br>typedef itk::ImageSeriesReader< ImageType > ReaderType;<br>
<br>typedef itk::ExtractImageFilter< ImageType, ImageType2D > ExtractFilterType;<br><br>typedef itk::<div>RescaleIntensityImageFilter< ImageType2D, WriteImageType > RescalerType;<br><br>
typedef itk::ImageFileWriter< WriteImageType > WriterType;<br>
<br><br>//-----------------------------------------------------------------------------<br>void ExtractSlice(ReaderType *reader, ExtractFilterType *extractFilter, int slice, int axis )<br>{<br><br> ImageType::RegionType inputRegion = reader->GetOutput()->GetLargestPossibleRegion();<br>
ImageType::SizeType size = inputRegion.GetSize();<br> <br> ImageType::IndexType start = inputRegion.GetIndex();<br> const unsigned int sliceNumber = slice; <br><br> if((axis > 2) || (axis < 0))<br> {<br>
cout << "In file imAbstractFilter.cxx at method ExtractSlice().
Variable axis must be between 2 and 0.\n" <<endl;<br>
}<br> <br> if( (sliceNumber > size[axis]) || (sliceNumber < start[axis]) )<br> {<br> cout << "In file imAbstractFilter.cxx at method ExtractSlice(). Variable slice must be between "<br> << size[axis] << " and " << start[axis] << ".\n" <<endl;<br>
}<br><br> size[axis] = 0; // Setando e armazenado o eixo de interesse.<br> start[axis] = sliceNumber; // Setando e armazenando o slice de interesse.<br> <br> ImageType::RegionType desiredRegion;<br> desiredRegion.SetSize( size );<br>
desiredRegion.SetIndex( start );<br> <br> extractFilter->SetExtractionRegion( desiredRegion );<br> extractFilter->SetInput( reader->GetOutput() );<br> <br> try<br> {<br> extractFilter->Update();<br>
}<br> catch( itk::ExceptionObject & excp )<br> {<br> cout << "In file imAbstractFilter.cxx at method ExtractSlice()\n" << endl;<br> }<br>}<br><br><br>//-----------------------------------------------------------------------------<br>
void DICOMReader_Reslice(ConnectorType *connector, bool x, bool y, bool z)<br>{<br> vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New(); <br> <br> vtkImageViewer2 * viewer = vtkImageViewer2::New();<br>
viewer->SetupInteractor( iren );<br> viewer->SetInput( connector->GetOutput() );<br> viewer->SetColorWindow( 256 );<br> viewer->SetColorLevel( 144 );<br> viewer->SetSlice(200);<br>
<br> if(x)<br> {<br> viewer->SetSliceOrientationToXY(); // orientação em X<br> }<br> else if(y)<br> {<br> viewer->SetSliceOrientationToYZ(); // orientação em Y<br> }<br> else if(z)<br>
{<br> viewer->SetSliceOrientationToXZ(); // orientação em Z<br> }<br> <br> viewer->GetRenderer()->ResetCamera();<br> viewer->Render();<br> iren->Start();<br> <br> iren->Delete();<br>
viewer->Delete();<br>}<br><br><br><br>//-----------------------------------------------------------------------------<br>int main( int argc, char **argv )<br>{<br> cout << "Lendo DICOM..." << endl;<br>
char inputDirectory[256] = {"/home/workspace2/ImageLab/data/Ricardo"};<br><br> ReaderType::Pointer reader = ReaderType::New();<br><br> typedef itk::GDCMImageIO ImageIOType;<br> ImageIOType::Pointer dicomIO = ImageIOType::New();<br>
reader->SetImageIO( dicomIO );<br><br> typedef itk::GDCMSeriesFileNames NamesGeneratorType;<br> NamesGeneratorType::Pointer nameGenerator = NamesGeneratorType::New();<br> nameGenerator->SetUseSeriesDetails( true );<br>
nameGenerator->SetDirectory( inputDirectory );<br><br> typedef std::vector< std::string > SeriesIdContainer;<br> const SeriesIdContainer & seriesUID = nameGenerator->GetSeriesUIDs();<br><br> SeriesIdContainer::const_iterator seriesItr = seriesUID.begin();<br>
SeriesIdContainer::const_iterator seriesEnd = seriesUID.end();<br> while( seriesItr != seriesEnd )<br> {<br> seriesItr++;<br> }<br><br> std::string seriesIdentifier;<br> seriesIdentifier = seriesUID.begin()->c_str();<br>
<br> typedef std::vector< std::string > FileNamesContainer;<br> FileNamesContainer fileNames;<br> fileNames = nameGenerator->GetFileNames( seriesIdentifier );<br><br> reader->SetFileNames( fileNames );<br>
reader->Update();<br><br> //######################## Filtro<br> cout << "Escrevendo JPG..." << endl;<br><br> ExtractFilterType::Pointer ExtractFilter = ExtractFilterType::New();<br> ExtractSlice(reader, ExtractFilter, 200, 0 );<br>
<br> RescalerType::Pointer rescaleFilter = RescalerType::New();<br> rescaleFilter->SetOutputMinimum(0);<br> rescaleFilter->SetOutputMaximum(255);<br> rescaleFilter->SetInput( ExtractFilter->GetOutput() );<br>
<br> char outputFileName[256] = {"/home/workspace2/ImageLab/data/Ricardo.png"};<br><br> WriterType::Pointer writer = WriterType::New();<br> writer->SetFileName( outputFileName );<br> writer->SetInput( rescaleFilter->GetOutput() );<br>
writer->Update();<br><br><br> //######################## Conector para o Pipeline ITK => VTK <br> cout << "Conectando ITK ao VTK..." <<endl;<br> ConnectorType::Pointer connector = ConnectorType::New();<br>
connector->SetInput( reader->GetOutput() );<br><br> DICOMReader_Reslice(connector, false, true, false);<br> <br> return 0;<br>}</div><br clear="all"><br>-- <br>Atenciosamente,<br><br>Eduardo Camargo<br>Analista de Sistemas<br>
HeMoLab - <a href="http://hemolab.lncc.br/equipe/desenv/eduardo.htm" target="_blank">http://hemolab.lncc.br/equipe/desenv/eduardo.htm</a><br>Laboratório Nacional de Computação Científica - LNCC<br>