[Insight-users] Why MultiResolutionPyramidImageFilter does not produce right image?
Zhiyong Xie
zyxie98 at yahoo.com
Thu, 5 Feb 2004 08:51:26 -0800 (PST)
--0-2105917344-1075999886=:45247
Content-Type: text/plain; charset=us-ascii
Lydia:
In MultiResolutionPDEDeformableRegistration, it seems we cannot set the schedule
of image pyramid (which is provide by SetSchedule method in MultiResolutionPyramidImageFilter ). It will be convenient if we can set different shrink
factor for each dimensions and for each image (fixed and moving) through
MultiResolutionPDEDeformableRegistration. I am not sure whether
we can do it in current implmentation. otherwise, it will be nice if the code can be
modified by implementing method SetFixedImagePyramid and SetMovingImagePyramid,
as well as a little change in GenerateData, where moving and fixed images can be
retrived from image pyramid instead of creating pyramid form image, if the pyramid
have already been set.
Thanks
Sean
Lydia Ng <lng at insightful.com> wrote:
Hi Fucang,
I think the filter requires the output pixel type to be floating point. Could
you try it out as an experiment?
- Lydia
> -----Original Message-----
> From: Fucang Jia [mailto:jiafucang at hotmail.com]
> Sent: Wednesday, November 19, 2003 1:17 AM
> To: insight-users at itk.org
> Subject: [Insight-users] Why MultiResolutionPyramidImageFilter does not
> produce right image?
>
> Hi, everyone,
>
> I am using ITK 1.4 and want to do some experiments based on
> MultiResolutionPyramid. But the following program does not give correct
> images. It does produce 8 images, but the first level image is all zeros.
> So
> the error spreads to the next level images. Is there some error in this
> program?
>
> Thanks a lot!
>
> Fucang
>
>
> =============================
> int main(int argc, char *argv[] )
> {
>
>
> if(argc<2)
> {
> std::cerr<> return -1;
> }
>
> typedef unsigned short PixelType;
> const unsigned short Dimension = 2;
>
> typedef itk::Image
InputImageType;
> typedef itk::Image
OutputImageType;
> typedef itk::ImageFileReader [input] ReaderType;
> typedef itk::ImageFileWriter WriterType;
> typedef
> itk::MultiResolutionPyramidImageFilter [input]
> PyramidType;
> ReaderType::Pointer reader = ReaderType::New();
> WriterType::Pointer writer = WriterType::New();
> PyramidType::Pointer pyramid = PyramidType::New();
>
>
> reader->SetFileName(argv[1]);
> reader->Update();
>
> pyramid->SetInput( reader->GetOutput() );
>
> unsigned int numLevels;
> itk::Vector factors;
>
> numLevels = 8;
> factors[0] = 256; factors[1] = 256;
>
> pyramid->SetNumberOfLevels( numLevels );
> pyramid->SetStartingShrinkFactors( 256 );
> pyramid->Print( std::cout );
>
> pyramid->Update();
>
> for (int i=0;i> {
> char str[4096];
> sprintf(str,"test_%d.png",i);
> writer->SetFileName(str);
> writer->SetInput(pyramid->GetOutput(i));
> writer->Update();
> }
>
> return 0;
> }
>
> _________________________________________________________________
> The new MSN 8: advanced junk mail protection and 2 months FREE*
> http://join.msn.com/?pagežatures/junkmail
>
> _______________________________________________
> 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
---------------------------------
Do you Yahoo!?
Yahoo! Finance: Get your refund fast by filing online
--0-2105917344-1075999886=:45247
Content-Type: text/html; charset=us-ascii
<DIV>Lydia:</DIV>
<DIV> </DIV>
<DIV>In MultiResolutionPDEDeformableRegistration, it seems we cannot set the schedule</DIV>
<DIV>of image pyramid (which is provide by SetSchedule method in MultiResolutionPyramidImageFilter ). It will be convenient if we can set different shrink</DIV>
<DIV>factor for each dimensions and for each image (fixed and moving) through</DIV>
<DIV>MultiResolutionPDEDeformableRegistration. I am not sure whether</DIV>
<DIV>we can do it in current implmentation. otherwise, it will be nice if the code can be</DIV>
<DIV>modified by implementing method SetFixedImagePyramid and SetMovingImagePyramid,</DIV>
<DIV>as well as a little change in GenerateData, where moving and fixed images can be </DIV>
<DIV>retrived from image pyramid instead of creating pyramid form image, if the pyramid </DIV>
<DIV>have already been set.</DIV>
<DIV> </DIV>
<DIV>Thanks</DIV>
<DIV> </DIV>
<DIV>Sean</DIV>
<DIV> </DIV>
<DIV><BR><B><I>Lydia Ng <lng at insightful.com></I></B> wrote:</DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">Hi Fucang,<BR><BR>I think the filter requires the output pixel type to be floating point. Could<BR>you try it out as an experiment?<BR><BR>- Lydia<BR><BR>> -----Original Message-----<BR>> From: Fucang Jia [mailto:jiafucang at hotmail.com]<BR>> Sent: Wednesday, November 19, 2003 1:17 AM<BR>> To: insight-users at itk.org<BR>> Subject: [Insight-users] Why MultiResolutionPyramidImageFilter does not<BR>> produce right image?<BR>> <BR>> Hi, everyone,<BR>> <BR>> I am using ITK 1.4 and want to do some experiments based on<BR>> MultiResolutionPyramid. But the following program does not give correct<BR>> images. It does produce 8 images, but the first level image is all zeros.<BR>> So<BR>> the error spreads to the next level images. Is there some error in this<BR>> program?<BR>> <BR>> Thanks a lot!<BR>> <BR>> Fucang<BR>> <BR>> <BR>>
=============================<BR>> int main(int argc, char *argv[] )<BR>> {<BR>> <BR>> <BR>> if(argc<2)<BR>> {<BR>> std::cerr<<ARGV[0]<<" Input2DImage?<<std::endl;<br>> return -1;<BR>> }<BR>> <BR>> typedef unsigned short PixelType;<BR>> const unsigned short Dimension = 2;<BR>> <BR>> typedef itk::Image<PIXELTYPE,DIMENSION> InputImageType;<BR>> typedef itk::Image<PIXELTYPE,DIMENSION> OutputImageType;<BR>> typedef itk::ImageFileReader<INPUTIMAGETYPE> ReaderType;<BR>> typedef itk::ImageFileWriter<OUTPUTIMAGETYPE> WriterType;<BR>> typedef<BR>> itk::MultiResolutionPyramidImageFilter<INPUTIMAGETYPE,OUTPUTIMAGETYPE><BR>> PyramidType;<BR>> ReaderType::Pointer reader = ReaderType::New();<BR>> WriterType::Pointer writer = WriterType::New();<BR>> PyramidType::Pointer pyramid = PyramidType::New();<BR>> <BR>> <BR>> reader->SetFileName(argv[1]);<BR>> reader->Update();<BR>> <BR>> pyramid->SetInput(
reader->GetOutput() );<BR>> <BR>> unsigned int numLevels;<BR>> itk::Vector<UNSIGNED int,Dimension> factors;<BR>> <BR>> numLevels = 8;<BR>> factors[0] = 256; factors[1] = 256;<BR>> <BR>> pyramid->SetNumberOfLevels( numLevels );<BR>> pyramid->SetStartingShrinkFactors( 256 );<BR>> pyramid->Print( std::cout );<BR>> <BR>> pyramid->Update();<BR>> <BR>> for (int i=0;i<NUMLEVELS;I++)<BR>> {<BR>> char str[4096];<BR>> sprintf(str,"test_%d.png",i);<BR>> writer->SetFileName(str);<BR>> writer->SetInput(pyramid->GetOutput(i));<BR>> writer->Update();<BR>> }<BR>> <BR>> return 0;<BR>> }<BR>> <BR>> _________________________________________________________________<BR>> The new MSN 8: advanced junk mail protection and 2 months FREE*<BR>> http://join.msn.com/?pagežatures/junkmail<BR>> <BR>> _______________________________________________<BR>> Insight-users mailing list<BR>>
Insight-users at itk.org<BR>> http://www.itk.org/mailman/listinfo/insight-users<BR>_______________________________________________<BR>Insight-users mailing list<BR>Insight-users at itk.org<BR>http://www.itk.org/mailman/listinfo/insight-users</BLOCKQUOTE><p><hr SIZE=1>
Do you Yahoo!?<br>
Yahoo! Finance: <a href="http://us.rd.yahoo.com/evt=22055/*http://taxes.yahoo.com/filing.html">Get your refund fast by filing online</a>
--0-2105917344-1075999886=:45247--