[Insight-users] FW: Dicom Series reader and Image Registration
raghu venkatram
raghu_420 at hotmail.com
Sat, 06 Mar 2004 20:07:54 +0000
<html><div style='background-color:'><DIV class=RTE>
<P>Hi Luis,</P>
<P>Thank you for replying and the useful links regarding image registration.</P>
<P>Medical Image Registration edited by Joseph Hajnal, David Hawkes and Derek Hill </P>
<P>is available online at <A href="http://www-ipg.umds.ac.uk/d.hill/hhh/">http://www-ipg.umds.ac.uk/d.hill/hhh/</A></P>
<P>thought this mighy help others also.</P>
<P>Warm regards,</P>
<P>Raghu<BR><BR></P></DIV>
<DIV></DIV>>From: Luis Ibanez <LUIS.IBANEZ at KITWARE.COM>
<DIV></DIV>>To: raghu venkatram <RAGHU_420 at HOTMAIL.COM>
<DIV></DIV>>CC: insight-users at itk.org
<DIV></DIV>>Subject: Re: [Insight-users] FW: Dicom Series reader and Image Registration
<DIV></DIV>>Date: Wed, 03 Mar 2004 23:16:20 -0500
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>Hi Raghu,
<DIV></DIV>>
<DIV></DIV>>There are two implementation of DICOM readers in ITK 1.4.
<DIV></DIV>>They are:
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>> Insight/Code/IO/
<DIV></DIV>> DICOMImageIO
<DIV></DIV>> DICOMImageIO2
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>DICOMImageIO2 is more recent than DICOMImageIO and use
<DIV></DIV>>a different parser mechanism based on the library in
<DIV></DIV>>
<DIV></DIV>> Insight/Utilities/DICOMParser
<DIV></DIV>>
<DIV></DIV>>They have different levels of verfication of compliance
<DIV></DIV>>with the DICOM format. It is likely that in one machine
<DIV></DIV>>you are using DICOMImageIO and in the other DICOMImageIO2.
<DIV></DIV>>Only one of them complains about magic numbers.
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>You may want try reading your DICOM files with MRIConvert
<DIV></DIV>>http://lcni.uoregon.edu/~jolinda/MRIConvert
<DIV></DIV>>
<DIV></DIV>>or with VolView
<DIV></DIV>>http://www.kitware.com/products/volview.html
<DIV></DIV>>
<DIV></DIV>>--
<DIV></DIV>>
<DIV></DIV>>About methods for evaluating registration result you may
<DIV></DIV>>want to look at the publications by Fitzpatrick
<DIV></DIV>>http://www.vuse.vanderbilt.edu/~jmf/jmf.html
<DIV></DIV>>
<DIV></DIV>>and the retrospective image registration evaluation project
<DIV></DIV>>http://www.vuse.vanderbilt.edu/~image/registration/
<DIV></DIV>>
<DIV></DIV>>you may also want to look at the book:
<DIV></DIV>>"Medical Image Registration" by Hajnal, Hill and Hawkes
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>---
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>In the deformable case you can speed up convergence
<DIV></DIV>>(if using Demons) if you precompute a deformation field
<DIV></DIV>>and provide it as initialization.
<DIV></DIV>>
<DIV></DIV>>In the FEM case, fiducial will not directly help to speed
<DIV></DIV>>up convergence, but they will allow you to define anatomical
<DIV></DIV>>features that should match between the two images, and by
<DIV></DIV>>doing so you can restrict the evolution of the deformation.
<DIV></DIV>>
<DIV></DIV>>All coordinates in registration should be interpreted in the
<DIV></DIV>>physical coordinate system of the fixed image. You may have
<DIV></DIV>>to look at the origin and spacing of your image in order to
<DIV></DIV>>get an idea of this coordinate frame.
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>Regards,
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>> Luis
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>-----------------------
<DIV></DIV>>raghu venkatram wrote:
<DIV></DIV>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>Hi Luis,
<DIV></DIV>>>Thank you for replying, it looks like dicom series reader is rather
<DIV></DIV>>>strange. I checked with one of my co workers, and it works fine on
<DIV></DIV>>>her machine, the same computer, OS, using ITK 1.4, same data. I do
<DIV></DIV>>>see the "no magic number warning", but it reads in all the files,
<DIV></DIV>>>unlike on my computer, where it generates all the file names and
<DIV></DIV>>>gives an abnormal program termination error, the same happens with
<DIV></DIV>>>yet another coworker of mine.
<DIV></DIV>>>
<DIV></DIV>>>the registration aspect, i am using analyze files generated by
<DIV></DIV>>>MRICro for the timebeing.
<DIV></DIV>>>ct data: 512x 512x58
<DIV></DIV>>>mr data: 256x 256 x 72
<DIV></DIV>>>
<DIV></DIV>>>the ct pixel sizes are roughly half the size of mr pixels, so
<DIV></DIV>>>spatially they occupy approx. the same volume.
<DIV></DIV>>>
<DIV></DIV>>>Mattesmutual information metric with 50,000 samples, time: 6-8
<DIV></DIV>>>mins. I was wondering if you could suggest methods to evaluate
<DIV></DIV>>>registration results. are there any papers, guidlines as far as the
<DIV></DIV>>>number of samples goes?
<DIV></DIV>>>
<DIV></DIV>>>also in the deformable case would using a pre-computed mesh be
<DIV></DIV>>>possible, to reduce time.
<DIV></DIV>>>
<DIV></DIV>>>would using landmarks/fiducials help in improving time complexity?
<DIV></DIV>>>i have seen the format of the landmarks file, but i was wondering
<DIV></DIV>>>as to what origin/coordinate system was used, is the origin in the
<DIV></DIV>>>center of the volume, or the top left hand corner of the volume?
<DIV></DIV>>>
<DIV></DIV>>>I apologize if my questions are too lame or have been previously
<DIV></DIV>>>stated.
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>Thanks in advance.
<DIV></DIV>>>Raghu
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>>From: Luis Ibanez <LUIS.IBANEZ at KITWARE.COM>
<DIV></DIV>>>>To: raghu venkatram <RAGHU_420 at HOTMAIL.COM>
<DIV></DIV>>>>CC: insight-users at itk.org
<DIV></DIV>>>>Subject: Re: [Insight-users] FW:
<DIV></DIV>>>>Date: Wed, 03 Mar 2004 17:04:46 -0500
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>Hi Raghu,
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>Good news: Your code works for me.
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>Bad news: Your data seems to be the problem.
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>The fact that you are getting the warning about
<DIV></DIV>>>>the "no magic number" is a bad sign about the
<DIV></DIV>>>>compliance of your files with the DICOM standard.
<DIV></DIV>>>>
<DIV></DIV>>>>You may want to get back to your data source
<DIV></DIV>>>>and complaint about the format used for your files.
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>Regards,
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>> Luis
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>----------------------
<DIV></DIV>>>>raghu venkatram wrote:
<DIV></DIV>>>>
<DIV></DIV>>>>>Hi Luis,
<DIV></DIV>>>>>
<DIV></DIV>>>>>I wish that was the problem, but i understand that the 0/1 flag
<DIV></DIV>>>>>is controls the ordering of the images. Also I have commented out
<DIV></DIV>>>>>the writing out part. For the time being i am just trying to read
<DIV></DIV>>>>>in the series from a directory which contains only one series and
<DIV></DIV>>>>>no other files. I see the filenames in the series, also it gives
<DIV></DIV>>>>>me the 'no magic number warning', for both dicom and
<DIV></DIV>>>>>dicom2iofactory. and after it spits out all the file names, it
<DIV></DIV>>>>>gives an abnormal program termination error.
<DIV></DIV>>>>>
<DIV></DIV>>>>>Thanks
<DIV></DIV>>>>>
<DIV></DIV>>>>>Raghu
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>>
<DIV></DIV>>>>> >From: Luis Ibanez
<DIV></DIV>>>>> >To: raghu venkatram
<DIV></DIV>>>>> >CC: insight-users at itk.org
<DIV></DIV>>>>> >Subject: Re: [Insight-users] FW:
<DIV></DIV>>>>> >Date: Tue, 02 Mar 2004 17:53:09 -0500
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >Hi Raghu,
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >The problem seems to be with the arguments that you
<DIV></DIV>>>>> >are passing to the program. Right now you are using
<DIV></DIV>>>>> >the example in:
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > Insight/Examples/IO/
<DIV></DIV>>>>> > DicomSeriesReadImageWriter.cxx
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >This program expects as argument:
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > 1) The directory where the DICOM files are
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > 2) The filename of an output image file where the
<DIV></DIV>>>>> > set of image slices will be saved as a 3D dataset.
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > 3) A *string identifying a dicom series*
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >You are passing:
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > F:\itk\data\MR\temp1\ 0
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >You are probably expecting the second argument "0"
<DIV></DIV>>>>> >to refer to the first series. That's not the case,
<DIV></DIV>>>>> >the program is taking "0" as the filename that you
<DIV></DIV>>>>> >want to use for saving the image data.
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >You sould be passing something like
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > F:\itk\data\MR\temp1 myOutputImageFile.mhd
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >The reason why the program is printing out all the
<DIV></DIV>>>>> >series Ids is for making it easier for you to take
<DIV></DIV>>>>> >one of those strings and use it as second argument
<DIV></DIV>>>>> >to the program. So, in a second run you can do
<DIV></DIV>>>>> >something more specific like:
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > DicomSeriesReaderImageWriter.exe
<DIV></DIV>>>>> > F:\itk\data\MR\temp1
<DIV></DIV>>>>> > myOutputImageFile.mhd
<DIV></DIV>>>>> > 1.3.12.2.1107.5.1.4.24515.4.0.6770456214062487
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >(all in a single line...)
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >Where the last string "1.3.12.2...." is the identifier
<DIV></DIV>>>>> >of one of the series.
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >Once you save your DICOM data in a fileformat
<DIV></DIV>>>>> >such as Analyze of MetaImage, you will be ready
<DIV></DIV>>>>> >for start testing the ImageRegistration examples.
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >Deformable registration can easily take 40 minutes
<DIV></DIV>>>>> >if you are dealing with 3D datasets.
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >Did you use the Multi-Resolution framework ?
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > Regards,
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> > Luis
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >----------------------
<DIV></DIV>>>>> >raghu venkatram wrote:
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >>Hello,
<DIV></DIV>>>>> >>I am using the dicom imageseriesreader, to read in a series
<DIV></DIV>>>>>of
<DIV></DIV>>>>> >>dicom images.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>My program arguments are:
<DIV></DIV>>>>> >>F:\itk\data\MR\temp1\ 0
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>I only have one series in this directory, and nothing other
<DIV></DIV>>>>>than
<DIV></DIV>>>>> >>the dicom files
<DIV></DIV>>>>> >>are in this directory.
<DIV></DIV>>>>> >>Its reads through, gets the filenames in the series and gives
<DIV></DIV>>>>>an
<DIV></DIV>>>>> >>"abnormal
<DIV></DIV>>>>> >>program termination error". Debugging shows that the it
<DIV></DIV>>>>>crashes on
<DIV></DIV>>>>> >>the second
<DIV></DIV>>>>> >>file in the series. I would appreciate if anyone could help
<DIV></DIV>>>>>me out.
<DIV></DIV>>>>> >>also i see
<DIV></DIV>>>>> >>Chunyan had postes a similar problem.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>If i am able to read in a series of dicom images, i want to
<DIV></DIV>>>>>use one
<DIV></DIV>>>>> >>of the
<DIV></DIV>>>>> >>registration algorithms and pass the two series, I am not
<DIV></DIV>>>>>sure if
<DIV></DIV>>>>> >>this would
<DIV></DIV>>>>> >>give me 3d registration. I tried using the analyze format,
<DIV></DIV>>>>>with
<DIV></DIV>>>>> >>deformable
<DIV></DIV>>>>> >>registration, it takes close to 40 mins, which makes me think
<DIV></DIV>>>>>i am
<DIV></DIV>>>>> >>doing
<DIV></DIV>>>>> >>something wrong. I would be glad if somebody could help me
<DIV></DIV>>>>>sort
<DIV></DIV>>>>> >>things out.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>All my dicom images are viewable using the QT viewr, MRICro.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>Thanks in advance
<DIV></DIV>>>>> >>Raghu Venkatram
<DIV></DIV>>>>> >>
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>/*=========================================================================
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> Program: Insight Segmentation & Registration Toolkit
<DIV></DIV>>>>> >> Module: $RCSfile: DicomSeriesReadImageWrite.cxx,v $
<DIV></DIV>>>>> >> Language: C++
<DIV></DIV>>>>> >> Date: $Date: 2003/12/30 21:13:04 $
<DIV></DIV>>>>> >> Version: $Revision: 1.2 $
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> Copyright (c) Insight Software Consortium. All rights
<DIV></DIV>>>>>reserved.
<DIV></DIV>>>>> >> See ITKCopyright.txt or
<DIV></DIV>>>>>http://www.itk.org/HTML/Copyright.htm for
<DIV></DIV>>>>> >>details.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> This software is distributed WITHOUT ANY WARRANTY;
<DIV></DIV>>>>>without
<DIV></DIV>>>>> >>even
<DIV></DIV>>>>> >> the implied warranty of MERCHANTABILITY or FITNESS FOR A
<DIV></DIV>>>>> >>PARTICULAR
<DIV></DIV>>>>> >> PURPOSE. See the above copyright notices for more
<DIV></DIV>>>>> >>information.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>=========================================================================*/
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>// Software Guide : BeginLatex
<DIV></DIV>>>>> >>//
<DIV></DIV>>>>> >>// This example illustrates how to read a DICOM series into
<DIV></DIV>>>>>a
<DIV></DIV>>>>> >>volume
<DIV></DIV>>>>> >>// and then save this volume in another file format.
<DIV></DIV>>>>> >>//
<DIV></DIV>>>>> >>// Software Guide : EndLatex
<DIV></DIV>>>>> >>
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>/*=========================================================================
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> Program: Insight Segmentation & Registration Toolkit
<DIV></DIV>>>>> >> Module: $RCSfile: itkDICOMImageSeriesTest.cxx,v $
<DIV></DIV>>>>> >> Language: C++
<DIV></DIV>>>>> >> Date: $Date: 2003/09/10 14:30:10 $xgoto-l
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> Version: $Revision: 1.5 $
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> Copyright (c) Insight Software Consortium. All rights
<DIV></DIV>>>>>reserved.
<DIV></DIV>>>>> >> See ITKCopyright.txt or
<DIV></DIV>>>>>http://www.itk.org/HTML/Copyright.htm for
<DIV></DIV>>>>> >>details.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> This software is distributed WITHOUT ANY WARRANTY;
<DIV></DIV>>>>>without
<DIV></DIV>>>>> >>even
<DIV></DIV>>>>> >> the implied warranty of MERCHANTABILITY or FITNESS FOR A
<DIV></DIV>>>>> >>PARTICULAR
<DIV></DIV>>>>> >> PURPOSE. See the above copyright notices for more
<DIV></DIV>>>>> >>information.
<DIV></DIV>>>>> >>
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>=========================================================================*/
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>#if defined(_MSC_VER)
<DIV></DIV>>>>> >>#pragma warning ( disable : 4786 )
<DIV></DIV>>>>> >>#endif
<DIV></DIV>>>>> >>#include "itkDICOMImageIOFactory.h"
<DIV></DIV>>>>> >>#include "itkDICOMImageIO2.h"
<DIV></DIV>>>>> >>#include "itkImageSeriesReader.h"
<DIV></DIV>>>>> >>#include "itkDICOMSeriesFileNames.h"
<DIV></DIV>>>>> >>#include "itkAnalyzeImageIO.h"
<DIV></DIV>>>>> >>#include "itkImageFileWriter.h"
<DIV></DIV>>>>> >>#include
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>"F:/itk/InsightToolkit-1.4.0/Testing/Code/BasicFilters/itkFilterWatcher.h"
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>int main(int ac, char* av[])
<DIV></DIV>>>>> >>{
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> typedef itk::AnalyzeImageIO ImageIOType;
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> if(ac < 3)
<DIV></DIV>>>>> >> {
<DIV></DIV>>>>> >> std::cerr << "Usage: " << av[0] << " DicomDirectory
<DIV></DIV>>>>> >>ReverseOrder(0/1)\n";
<DIV></DIV>>>>> >> return EXIT_FAILURE;
<DIV></DIV>>>>> >> }
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> typedef itk::Image ImageNDType;
<DIV></DIV>>>>> >> typedef itk::Image Image3DType;
<DIV></DIV>>>>> >> typedef itk::ImageSeriesReader ReaderType;
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> itk::DICOMImageIO2::Pointer io = itk::DICOMImageIO2::New();
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> // Get the DICOM filenames from the directory
<DIV></DIV>>>>> >> itk::DICOMSeriesFileNames::Pointer names =
<DIV></DIV>>>>> >>itk::DICOMSeriesFileNames::New();
<DIV></DIV>>>>> >> names->SetDirectory(av[1]);
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> ReaderType::Pointer reader = ReaderType::New();
<DIV></DIV>>>>> >> reader->SetFileNames(names->GetFileNames());
<DIV></DIV>>>>> >> reader->SetImageIO(io);
<DIV></DIV>>>>> >> std::cout << names;
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> FilterWatcher watcher(reader);
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> try
<DIV></DIV>>>>> >> {
<DIV></DIV>>>>> >> if (atoi(av[2]))
<DIV></DIV>>>>> >> {
<DIV></DIV>>>>> >> reader->ReverseOrderOn();
<DIV></DIV>>>>> >> }
<DIV></DIV>>>>> >> reader->Update();
<DIV></DIV>>>>> >> reader->GetOutput()->Print(std::cout);
<DIV></DIV>>>>> >> }
<DIV></DIV>>>>> >> catch (itk::ExceptionObject &ex)
<DIV></DIV>>>>> >> {
<DIV></DIV>>>>> >> std::cout << ex;
<DIV></DIV>>>>> >> return EXIT_FAILURE;
<DIV></DIV>>>>> >> }
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>/*
<DIV></DIV>>>>> >> typedef itk::ImageFileWriter< Image3DType > WriterType;
<DIV></DIV>>>>> >> WriterType::Pointer writer = WriterType::New();
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> std::cout << "Writing the image as " << std::endl <<
<DIV></DIV>>>>> >>std::endl;
<DIV></DIV>>>>> >> std::cout << av[2] << std::endl << std::endl;
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> writer->SetFileName( "out.hdr" );
<DIV></DIV>>>>> >> ImageIOType::Pointer analyzeIO = ImageIOType::New();
<DIV></DIV>>>>> >> writer->SetImageIO(analyzeIO);
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> writer->SetInput( reader->GetOutput() );
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> try
<DIV></DIV>>>>> >> {
<DIV></DIV>>>>> >> writer->Update();
<DIV></DIV>>>>> >> }
<DIV></DIV>>>>> >> catch (itk::ExceptionObject &ex)
<DIV></DIV>>>>> >> {
<DIV></DIV>>>>> >> std::cout << ex;
<DIV></DIV>>>>> >> return EXIT_FAILURE;
<DIV></DIV>>>>> >> }
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>*/
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >> return EXIT_SUCCESS;
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>}
<DIV></DIV>>>>> >>
<DIV></DIV>>>>>
<DIV></DIV>>>>> >>_________________________________________________________________
<DIV></DIV>>>>> >>Fast. Reliable. Get MSN 9 Dial-up - 1 month FREE!
<DIV></DIV>>>>>(Limited-time
<DIV></DIV>>>>> >>Offer)
<DIV></DIV>>>>>http://click.atdmt.com/AVE/go/onm00200361ave/direct/01/
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >>_______________________________________________
<DIV></DIV>>>>> >>Insight-users mailing list
<DIV></DIV>>>>> >>Insight-users at itk.org
<DIV></DIV>>>>> >>http://www.itk.org/mailman/listinfo/insight-users
<DIV></DIV>>>>> >>
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>> >
<DIV></DIV>>>>>
<DIV></DIV>>>>>------------------------------------------------------------------------
<DIV></DIV>>>>>Store more e-mails with MSN Hotmail Extra Storage – 4 plans to
<DIV></DIV>>>>>choose from! <HTTP: 2737??PS="" 8HMAENUS g.msn.com>
<DIV></DIV>>>>>_______________________________________________ Insight-users
<DIV></DIV>>>>>mailing list Insight-users at itk.org
<DIV></DIV>>>>>http://www.itk.org/mailman/listinfo/insight-users
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>
<DIV></DIV>>>>_______________________________________________
<DIV></DIV>>>>Insight-users mailing list
<DIV></DIV>>>>Insight-users at itk.org
<DIV></DIV>>>>http://www.itk.org/mailman/listinfo/insight-users
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>>_________________________________________________________________
<DIV></DIV>>>Frustrated with dial-up? Lightning-fast Internet access for as low
<DIV></DIV>>>as $29.95/month.
<DIV></DIV>>>http://click.atdmt.com/AVE/go/onm00200360ave/direct/01/
<DIV></DIV>>>
<DIV></DIV>>>
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>
<DIV></DIV>>_______________________________________________
<DIV></DIV>>Insight-users mailing list
<DIV></DIV>>Insight-users at itk.org
<DIV></DIV>>http://www.itk.org/mailman/listinfo/insight-users
<DIV></DIV></div><br clear=all><hr> <a href="http://g.msn.com/8HMBENUS/2752??PS=">Find things fast with the new MSN Toolbar – includes FREE pop-up blocking!</a> </html>