[Insight-users] A question about least_square method !
zhen zonglei
amoroma at gmail.com
Tue Jan 11 05:45:00 EST 2005
Hi everyone.
I am a new user of ITK.
I want to konw whether the least_square method have been implemented
in ITK as a common Numerical method, and which class can provide the
LS method?
Thanks!
On Tue, 11 Jan 2005 01:31:20 -0500 (EST),
insight-users-request at itk.org <insight-users-request at itk.org> wrote:
> Send Insight-users mailing list submissions to
> insight-users at itk.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://www.itk.org/mailman/listinfo/insight-users
> or, via email, send a message with subject or body 'help' to
> insight-users-request at itk.org
>
> You can reach the person managing the list at
> insight-users-owner at itk.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Insight-users digest..."
>
> Today's Topics:
>
> 1. Simple image program help (RJ Stevens)
> 2. Re: Simple image program help (Jens Eickmeyer)
> 3. Re: test image for DeformableRegistration7.cxx and
> DeformableRegistration8.cxx (Jim)
> 4. Re: test image for DeformableRegistration7.cxx and
> DeformableRegistration8.cxx (Luis Ibanez)
> 5. Re: One+One_Revolutionary (Luis Ibanez)
> 6. Inversion of pixel value within a fixed intensity window
> (Amitesh Agarwal)
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 10 Jan 2005 13:16:58 -0800
> From: RJ Stevens <rj.stevens at gmail.com>
> Subject: [Insight-users] Simple image program help
> To: insight-users at itk.org
> Message-ID: <a7014a93050110131668bd6c96 at mail.gmail.com>
> Content-Type: text/plain; charset=US-ASCII
>
> Hi again, I am building a simple program to open an image based off
> the example Image2.
> I have included all the information this time, namely the program, the
> CMakeLists.txt, and the error messages.
>
> I am running SUSE linux 9.2 with Itk version 1.8.
> Thanks for any input.
>
> Richard Stevens
>
> Simple program based off of image2 example:
>
> #include "itkImage.h"
> #include "itkImageFileReader.h"
> #include <iostream>
>
> int main(int , char ** argv)
> {
>
> //define variables for image def
> typedef unsigned char PixelType;
> const unsigned int Dimension = 2;
>
> //set pixel type and dimension
> typedef itk::Image< PixelType, Dimension > ImageType;
>
> //instantiate the image reader class
> typedef itk::ImageFileReader< ImageType > ReaderType;
>
> //use reader class to create a reader object
> ReaderType::Pointer reader = ReaderType::New();
>
> //give reader oject filename of image to be loaded into memory
> const char * filename = argv[1];
> reader->SetFileName(filename);
>
> //reader object is a pipeline source object and responds to pipeline
> //update requests. Since this pipeline is of length 1, we must send the
> //update reqest manually
> reader->Update();
>
> //put the image into an image object
> ImageType::Pointer image = reader -> GetOutput();
>
> return 0;
> }
>
> The CMakeLists.txt file this program that is mostly borrowed from
> InsightApplications after my first attempt did not work:
>
> CMAKE_MINIMUM_REQUIRED(VERSION 1.6)
> PROJECT(Segment1)
> INCLUDE_REGULAR_EXPRESSION("^.*$")
>
> # Tell CMake's FIND_PACKAGE command where to look for our Find<pkg>
> # modules.
> SET(CMAKE_MODULE_PATH ${CMAKE_MODULE_PATH} ${InsightApplications_SOURCE_DIR})
>
> IF(NOT WIN32)
> OPTION(BUILD_SHARED_LIBS "Build Insight-Applications with shared
> libraries." OFF)
> ENDIF(NOT WIN32)
>
> # Don't build anything unless the version of CMake is high enough.
> # The matching ELSE/ENDIF should be the last lines in the file.
> IF("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" GREATER 1.5)
>
> FIND_PACKAGE(ITK)
> IF(ITK_FOUND)
> INCLUDE(${ITK_USE_FILE})
> ELSE(ITK_FOUND)
> MESSAGE(FATAL_ERROR
> "ITK not found. Set ITK_DIR.")
> ENDIF(ITK_FOUND)
>
> #added from CMakeLists.txt from InsightApplications
> LINK_LIBRARIES (
> ITKAlgorithms
> ITKCommon
> ITKBasicFilters
> )
>
> ADD_EXECUTABLE( Segment1 Segment1.cxx)
>
> # If the version of CMake was too old, complain and build nothing.
> # These should be the last lines in this file.
> ELSE("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" GREATER 1.5)
> MESSAGE(SEND_ERROR
> "This version of CMake is too old to build InsightApplications. "
> "Please upgrade to CMake 1.6.")
> ENDIF("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" GREATER 1.5)
>
> Errors Generated:
>
> :/VTKITK/work/picture2 # make
> Building dependencies. cmake.depends...
> Building object file Segment1.o...
> Building executable ./Segment1...
> Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE25GenerateOutputInformationEv+0x28b):
> In function `itk::ImageFileReader<itk::Image<unsigned char,
> (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> >::GenerateOutputInformation()':
> : undefined reference to `itk::ImageIOFactory::CreateImageIO(char
> const*, itk::ImageIOFactory::FileModeType)'
> Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE12GenerateDataEv+0x3ad):
> In function `itk::ImageFileReader<itk::Image<unsigned char,
> (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> >::GenerateData()':
> : undefined reference to `itk::operator<<(std::basic_ostream<char,
> std::char_traits<char> >&, itk::ImageIORegion const&)'
> Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE12GenerateDataEv+0x932):
> In function `itk::ImageFileReader<itk::Image<unsigned char,
> (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> >::GenerateData()':
> : undefined reference to `itk::ImageIOBase::GetImageSizeInBytes() const'
> Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE15DoConvertBufferEPvm+0x66b):
> In function `itk::ImageFileReader<itk::Image<unsigned char,
> (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> >::DoConvertBuffer(void*, unsigned long)':
> : undefined reference to
> `itk::ImageIOBase::GetComponentTypeAsString(itk::ImageIOBase::IOComponentType)
> const'
> Segment1.o(.gnu.linkonce.t._ZN3itk13ImageIORegionC1Ej+0x1a): In
> function `itk::ImageIORegion::ImageIORegion[in-charge](unsigned)':
> : undefined reference to `vtable for itk::ImageIORegion'
> Segment1.o(.gnu.linkonce.t._ZN3itk13ImageIORegionD1Ev+0xc): In
> function `itk::ImageIORegion::~ImageIORegion [in-charge]()':
> : undefined reference to `vtable for itk::ImageIORegion'
> Segment1.o(.gnu.linkonce.t._ZN3itk13ImageIORegionC1ERKS0_+0x1a): In
> function `itk::ImageIORegion::ImageIORegion[in-charge](itk::ImageIORegion
> const&)':
> : undefined reference to `vtable for itk::ImageIORegion'
> collect2: ld returned 1 exit status
> make[1]: *** [Segment1] Error 1
> make: *** [default_target] Error 2
>
> ------------------------------
>
> Message: 2
> Date: Mon, 10 Jan 2005 22:42:49 +0100
> From: Jens Eickmeyer <eickmeyer at gmail.com>
> Subject: Re: [Insight-users] Simple image program help
> To: RJ Stevens <rj.stevens at gmail.com>
> Cc: insight-users at itk.org
> Message-ID: <68b64933050110134241971fe2 at mail.gmail.com>
> Content-Type: text/plain; charset=US-ASCII
>
> Hi!
>
> I think that you have forgotten to link to ITKIO, too. So adding ITKIO
> to the LINK_LIBRARIES part in CMakeLists.txt should fix the problem.
>
> Jens
>
> On Mon, 10 Jan 2005 13:16:58 -0800, RJ Stevens <rj.stevens at gmail.com> wrote:
> > Hi again, I am building a simple program to open an image based off
> > the example Image2.
> > I have included all the information this time, namely the program, the
> > CMakeLists.txt, and the error messages.
> >
> > I am running SUSE linux 9.2 with Itk version 1.8.
> > Thanks for any input.
> >
> > Richard Stevens
> >
> > Simple program based off of image2 example:
> >
> > #include "itkImage.h"
> > #include "itkImageFileReader.h"
> > #include <iostream>
> >
> > int main(int , char ** argv)
> > {
> >
> > //define variables for image def
> > typedef unsigned char PixelType;
> > const unsigned int Dimension = 2;
> >
> > //set pixel type and dimension
> > typedef itk::Image< PixelType, Dimension > ImageType;
> >
> > //instantiate the image reader class
> > typedef itk::ImageFileReader< ImageType > ReaderType;
> >
> > //use reader class to create a reader object
> > ReaderType::Pointer reader = ReaderType::New();
> >
> > //give reader oject filename of image to be loaded into memory
> > const char * filename = argv[1];
> > reader->SetFileName(filename);
> >
> > //reader object is a pipeline source object and responds to pipeline
> > //update requests. Since this pipeline is of length 1, we must send the
> > //update reqest manually
> > reader->Update();
> >
> > //put the image into an image object
> > ImageType::Pointer image = reader -> GetOutput();
> >
> > return 0;
> > }
> >
> > The CMakeLists.txt file this program that is mostly borrowed from
> > InsightApplications after my first attempt did not work:
> >
> > CMAKE_MINIMUM_REQUIRED(VERSION 1.6)
> > PROJECT(Segment1)
> > INCLUDE_REGULAR_EXPRESSION("^.*$")
> >
> > # Tell CMake's FIND_PACKAGE command where to look for our Find<pkg>
> > # modules.
> > SET(CMAKE_MODULE_PATH ${CMAKE_MODULE_PATH} ${InsightApplications_SOURCE_DIR})
> >
> > IF(NOT WIN32)
> > OPTION(BUILD_SHARED_LIBS "Build Insight-Applications with shared
> > libraries." OFF)
> > ENDIF(NOT WIN32)
> >
> > # Don't build anything unless the version of CMake is high enough.
> > # The matching ELSE/ENDIF should be the last lines in the file.
> > IF("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" GREATER 1.5)
> >
> > FIND_PACKAGE(ITK)
> > IF(ITK_FOUND)
> > INCLUDE(${ITK_USE_FILE})
> > ELSE(ITK_FOUND)
> > MESSAGE(FATAL_ERROR
> > "ITK not found. Set ITK_DIR.")
> > ENDIF(ITK_FOUND)
> >
> > #added from CMakeLists.txt from InsightApplications
> > LINK_LIBRARIES (
> > ITKAlgorithms
> > ITKCommon
> > ITKBasicFilters
> > )
> >
> > ADD_EXECUTABLE( Segment1 Segment1.cxx)
> >
> > # If the version of CMake was too old, complain and build nothing.
> > # These should be the last lines in this file.
> > ELSE("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" GREATER 1.5)
> > MESSAGE(SEND_ERROR
> > "This version of CMake is too old to build InsightApplications. "
> > "Please upgrade to CMake 1.6.")
> > ENDIF("${CMAKE_MAJOR_VERSION}.${CMAKE_MINOR_VERSION}" GREATER 1.5)
> >
> > Errors Generated:
> >
> > :/VTKITK/work/picture2 # make
> > Building dependencies. cmake.depends...
> > Building object file Segment1.o...
> > Building executable ./Segment1...
> > Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE25GenerateOutputInformationEv+0x28b):
> > In function `itk::ImageFileReader<itk::Image<unsigned char,
> > (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> > >::GenerateOutputInformation()':
> > : undefined reference to `itk::ImageIOFactory::CreateImageIO(char
> > const*, itk::ImageIOFactory::FileModeType)'
> > Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE12GenerateDataEv+0x3ad):
> > In function `itk::ImageFileReader<itk::Image<unsigned char,
> > (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> > >::GenerateData()':
> > : undefined reference to `itk::operator<<(std::basic_ostream<char,
> > std::char_traits<char> >&, itk::ImageIORegion const&)'
> > Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE12GenerateDataEv+0x932):
> > In function `itk::ImageFileReader<itk::Image<unsigned char,
> > (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> > >::GenerateData()':
> > : undefined reference to `itk::ImageIOBase::GetImageSizeInBytes() const'
> > Segment1.o(.gnu.linkonce.t._ZN3itk15ImageFileReaderINS_5ImageIhLj2EEENS_25DefaultConvertPixelTraitsIhEEE15DoConvertBufferEPvm+0x66b):
> > In function `itk::ImageFileReader<itk::Image<unsigned char,
> > (unsigned)2>, itk::DefaultConvertPixelTraits<unsigned char>
> > >::DoConvertBuffer(void*, unsigned long)':
> > : undefined reference to
> > `itk::ImageIOBase::GetComponentTypeAsString(itk::ImageIOBase::IOComponentType)
> > const'
> > Segment1.o(.gnu.linkonce.t._ZN3itk13ImageIORegionC1Ej+0x1a): In
> > function `itk::ImageIORegion::ImageIORegion[in-charge](unsigned)':
> > : undefined reference to `vtable for itk::ImageIORegion'
> > Segment1.o(.gnu.linkonce.t._ZN3itk13ImageIORegionD1Ev+0xc): In
> > function `itk::ImageIORegion::~ImageIORegion [in-charge]()':
> > : undefined reference to `vtable for itk::ImageIORegion'
> > Segment1.o(.gnu.linkonce.t._ZN3itk13ImageIORegionC1ERKS0_+0x1a): In
> > function `itk::ImageIORegion::ImageIORegion[in-charge](itk::ImageIORegion
> > const&)':
> > : undefined reference to `vtable for itk::ImageIORegion'
> > collect2: ld returned 1 exit status
> > make[1]: *** [Segment1] Error 1
> > make: *** [default_target] Error 2
> > _______________________________________________
> > Insight-users mailing list
> > Insight-users at itk.org
> > http://www.itk.org/mailman/listinfo/insight-users
> >
>
> ------------------------------
>
> Message: 3
> Date: Mon, 10 Jan 2005 16:24:52 -0600
> From: Jim <itkvtk at gmail.com>
> Subject: Re: [Insight-users] test image for
> DeformableRegistration7.cxx and DeformableRegistration8.cxx
> To: Luis Ibanez <luis.ibanez at kitware.com>
> Cc: Insight-users at itk.org
> Message-ID: <5a4be7e05011014247e9c9602 at mail.gmail.com>
> Content-Type: text/plain; charset=US-ASCII
>
> Hi Luis,
>
> The DeformableRegistration7 and DeformableRegistration8 are meant for
> 3D right? I modified the line "typedef float PixelType;" to
> "typedef unsigned short PixelType; in DeformableRegistration7 for my
> 512x512x40 CT images. When I run it with the following command:
> DeformableRegistration7 image1.img image2.img out.1 out.2 out.3 out.4,
> I got an error message saying "This application has requested the
> Runtime to terminate it in an unusual way. Please contact the
> application's support team for more information."
>
> I think there is lots information needs to give to
> DeformableRegistration7, such as the image dimensions, pixel size,
> spacing, etc. I have look into the InsightDocumentations try to see
> any instructions for this example without success. I thought there
> must be some test images available somewhere I can test with before I
> test it on my own image data. If it runs fine, I can then either
> modify my image format to match test image or modify the example
> program to read my images.
>
> Thanks
>
> Jim
>
> On Sat, 08 Jan 2005 09:40:27 -0500, Luis Ibanez <luis.ibanez at kitware.com> wrote:
> > Hi Jim,
> >
> > For DeformableRegistration7.cxx do:
> >
> > DeformableRegistration7
> > fixedImage
> > movingImage
> > outputResampledMovingImage
> > [differenceImageAfterRegistration]
> > [differenceImageBeforeRegistration]
> > [deformationField]
> >
> > All in a single command line. Note that for the
> > deformation field you must use an image file format
> > that supports multiplecomponents, such as VTK or
> > MetaImage. The entries in-between "[]" are optional.
> >
> > For visualizing the deformation field you can use
> > ParaView:
> >
> > http://www.paraview.org
> >
> > For input images you coulde use:
> >
> > fixedImage = Insight/Examples/Data/
> > BrainProtonDensitySliceBorder20.png
> >
> > movingImage = Insight/Examples/Data/
> > BrainProtonDensitySliceBSplined10.png
> >
> > The same arguments are expected by DeformableRegistration8.
> >
> > Please let us know if you have further questions.
> >
> > Thanks
> >
> > Luis
> >
> > -------------------------------------------
> > Jim wrote:
> >
> > > Hi,
> > >
> > > Where can I find instruction and test images for
> > > DeformableRegistration7.cxx and DeformableRegistration8.cxx.
> > >
> > > Thanks.
> > >
> > > Jim
> > > _______________________________________________
> > > Insight-users mailing list
> > > Insight-users at itk.org
> > > http://www.itk.org/mailman/listinfo/insight-users
> > >
> > >
> >
> >
>
> ------------------------------
>
> Message: 4
> Date: Mon, 10 Jan 2005 17:41:16 -0500
> From: Luis Ibanez <luis.ibanez at kitware.com>
> Subject: Re: [Insight-users] test image for
> DeformableRegistration7.cxx and DeformableRegistration8.cxx
> To: Jim <itkvtk at gmail.com>
> Cc: Insight-users at itk.org
> Message-ID: <41E3048C.5000201 at kitware.com>
> Content-Type: text/plain; charset=us-ascii; format=flowed
>
> Hi Jim,
>
> Yes, these examples are intended for 3D images.
>
> The problems is that you are passing filenames that are not
> associated with any accepted fileformat:
>
> out.1
> out.2
> out.3
>
> you should use something like
>
> out1.mhd
> out2.mhd
> out3.mhd
>
> for using MetaImage fileformat
>
> or
>
> out1.hdr
> out2.hdr
> out3.hdr
>
> for using Analyze fileformats
>
> Please look at the Wiki FAQ for a full list of the
> fileformat supported in ITK.
>
> ---
>
> About your images: Are your ".img" files RAW data ?
> if so, you *MUST* create MetaImage headers for those
> files. In that way you will specify pixel type, number
> of pixels, pixel spacing... etc as part of the content
> of the MetaImage header.
>
> -----------
>
> How about trying with the BrainWeb images available
> in the FTP data directory:
>
> ftp://public.kitware.com/pub/itk/Data/
>
> Regards,
>
> Luis
>
> -----------------
> Jim wrote:
>
> > Hi Luis,
> >
> > The DeformableRegistration7 and DeformableRegistration8 are meant for
> > 3D right? I modified the line "typedef float PixelType;" to
> > "typedef unsigned short PixelType; in DeformableRegistration7 for my
> > 512x512x40 CT images. When I run it with the following command:
> > DeformableRegistration7 image1.img image2.img out.1 out.2 out.3 out.4,
> > I got an error message saying "This application has requested the
> > Runtime to terminate it in an unusual way. Please contact the
> > application's support team for more information."
> >
> > I think there is lots information needs to give to
> > DeformableRegistration7, such as the image dimensions, pixel size,
> > spacing, etc. I have look into the InsightDocumentations try to see
> > any instructions for this example without success. I thought there
> > must be some test images available somewhere I can test with before I
> > test it on my own image data. If it runs fine, I can then either
> > modify my image format to match test image or modify the example
> > program to read my images.
> >
> > Thanks
> >
> > Jim
> >
> > On Sat, 08 Jan 2005 09:40:27 -0500, Luis Ibanez <luis.ibanez at kitware.com> wrote:
> >
> >>Hi Jim,
> >>
> >>For DeformableRegistration7.cxx do:
> >>
> >> DeformableRegistration7
> >> fixedImage
> >> movingImage
> >> outputResampledMovingImage
> >> [differenceImageAfterRegistration]
> >> [differenceImageBeforeRegistration]
> >> [deformationField]
> >>
> >>All in a single command line. Note that for the
> >>deformation field you must use an image file format
> >>that supports multiplecomponents, such as VTK or
> >>MetaImage. The entries in-between "[]" are optional.
> >>
> >>For visualizing the deformation field you can use
> >>ParaView:
> >>
> >> http://www.paraview.org
> >>
> >>For input images you coulde use:
> >>
> >> fixedImage = Insight/Examples/Data/
> >> BrainProtonDensitySliceBorder20.png
> >>
> >> movingImage = Insight/Examples/Data/
> >> BrainProtonDensitySliceBSplined10.png
> >>
> >>The same arguments are expected by DeformableRegistration8.
> >>
> >>Please let us know if you have further questions.
> >>
> >> Thanks
> >>
> >> Luis
> >>
> >>-------------------------------------------
> >>Jim wrote:
> >>
> >>
> >>>Hi,
> >>>
> >>>Where can I find instruction and test images for
> >>>DeformableRegistration7.cxx and DeformableRegistration8.cxx.
> >>>
> >>>Thanks.
> >>>
> >>>Jim
> >>>_______________________________________________
> >>>Insight-users mailing list
> >>>Insight-users at itk.org
> >>>http://www.itk.org/mailman/listinfo/insight-users
> >>>
> >>>
> >>
> >>
> > _______________________________________________
> > Insight-users mailing list
> > Insight-users at itk.org
> > http://www.itk.org/mailman/listinfo/insight-users
> >
> >
>
> ------------------------------
>
> Message: 5
> Date: Mon, 10 Jan 2005 20:34:23 -0500
> From: Luis Ibanez <luis.ibanez at kitware.com>
> Subject: [Insight-users] Re: One+One_Revolutionary
> To: "Li, George (NIH/NCI)" <ligeorge at mail.nih.gov>
> Cc: insight-users at itk.org
> Message-ID: <41E32D1F.2080204 at kitware.com>
> Content-Type: text/plain; charset=us-ascii; format=flowed
>
> Hi George,
>
> These repeats on the Metric values are normal for the 1+1 optimizer.
>
> They reflect the way in which this family of Evolutionary
> optimizers work.
>
> Evolutionary optimizers bred a population of points in the
> parametric space. At every iteration, the population is
> renew by mixing among the best individuals, and removing
> some of the ones that are not well fitted.
>
> The values that the optimizer reports at every iteration is the
> one of the best fitted individual among the current population.
> It is therefore common to see the same value appear in multiple
> iterations until a better fitted individual appears in the
> population.
>
> The fact that you see the same value in several iterations is not
> a problem. The population of samples in changing internally, and
> simply has not found a better fit yet.
>
> For detail on the rational behind this optimizer you may want
> to look at the references indicated in the Doxygen documentation:
>
> http://www.itk.org/Insight/Doxygen/html/classitk_1_1OnePlusOneEvolutionaryOptimizer.html
>
> "Parametric estimate of intensity inhomogeneities applied to MRI" Martin
> Styner, G. Gerig, Christian Brechbuehler, Gabor Szekely, IEEE
> TRANSACTIONS ON MEDICAL IMAGING; 19(3), pp. 153-165, 2000,
> (http://www.cs.unc.edu/~styner/docs/tmi00.pdf)
>
> "Evaluation of 2D/3D bias correction with 1+1ES-optimization" Martin
> Styner, Prof. Dr. G. Gerig (IKT, BIWI, ETH Zuerich), TR-197
> (http://www.cs.unc.edu/~styner/docs/StynerTR97.pdf)
>
> Regards,
>
> Luis
>
> BTW: it is "Evolutionary", ... not "Revolutionary" :-)
>
> ---------------------------------
> Li, George (NIH/NCI) wrote:
>
> > Hi Luis:
> >
> > I have tested OnePlusOneEvolutionaryOptimizer, using the same set
> > of dicom series with initial shifts in rotation and translation.
> > The recorded iteration looks strange to me, however. It has many
> > repeats with the same versor and vector for the transformation,
> > as well as the same metric value. FYI, Mattes' MI metric is used.
> >
> > Could you please take a look at the text file? And please let me
> > know if it is supposed to be or not?
> >
> > Thanks for your advice.
> >
> > George
> >
>
> ------------------------------
>
> Message: 6
> Date: Tue, 11 Jan 2005 00:30:59 -0600
> From: Amitesh Agarwal <amiteshagarwal at gmail.com>
> Subject: [Insight-users] Inversion of pixel value within a fixed
> intensity window
> To: Insight-users at itk.org
> Cc: luis.ibanez at kitware.com
> Message-ID: <dd63bb7905011022306eff5e51 at mail.gmail.com>
> Content-Type: text/plain; charset=US-ASCII
>
> Hello,
>
> I wanted to know if there is any method available which can let me
> perform a mathematical operation on each pixel of a 3D .vtk image.
>
> the function would be of type
> G(I) = -1* f(I) + (Imax - Imin) .... where I is the intensity value of
> a pixel, Imax is the intensity maxima in an image and Imin is the
> intensity minima in a image.
>
> The reason for this pixel intensity transformation is that I wanted to
> create a negative of an Image or in other words I wanted to invert the
> intensity value within a maximum-minimum intensity range.
>
> Thanks
>
> Amitesh
>
> ------------------------------
>
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
>
> End of Insight-users Digest, Vol 9, Issue 26
> ********************************************
>
More information about the Insight-users
mailing list