[Insight-developers] FFT Patch / Tcon Today Afternoon

[Kevin H. Hobbs]

On 09/09/2011 10:48 AM, Luis Ibanez wrote:
> Let's find a time later today or early next week
> to have a technical discussion about the details
> of FFT in ITK, and in particular about this patch:
> 
> http://review.source.kitware.com/#change,2625
> 

I'm trying to make my fft->mask->ifft program work with this patch.

and rushing this e-mail before the tcon

I run the attached code and the mask filter complains that the images do
not occupy the same space yet when I print them they seem to.



[kevin at bubbles ConsistentFFT]$ ./FFTMaskIFFT ../Sarcomere/032911/032911\
p1\ unstretched\ middle\ 40x.lsm mask.mhd masked.mhd
Image (0xf643f0)
  RTTI typeinfo:   itk::Image<std::complex<float>, 2u>
  Reference Count: 1
  Modified Time: 242
  Debug: Off
  Observers:
    none
  Source: (0xf64300)
  Source output index: 0
  Release Data: Off
  Data Released: False
  Global Release Data: Off
  PipelineMTime: 221
  UpdateMTime: 243
  RealTimeStamp: 1.31559e+09 seconds
  LargestPossibleRegion:
    Dimension: 2
    Index: [0, 0]
    Size: [2048, 2048]
  BufferedRegion:
    Dimension: 2
    Index: [0, 0]
    Size: [2048, 2048]
  RequestedRegion:
    Dimension: 2
    Index: [0, 0]
    Size: [2048, 2048]
  Spacing: [1.09863e-07, 1.09863e-07]
  Origin: [0, 0]
  Direction:
1 0
0 1

  IndexToPointMatrix:
  1.09863e-07 0
0 1.09863e-07

  PointToIndexMatrix:
  9.10222e+06 0
0 9.10222e+06

  PixelContainer:
    ImportImageContainer (0xf645f0)
      RTTI typeinfo:   itk::ImportImageContainer<unsigned long,
std::complex<float> >
      Reference Count: 1
      Modified Time: 241
      Debug: Off
      Observers:
        none
      Pointer: 0x7f2335fff010
      Container manages memory: true
      Size: 4194304
      Capacity: 4194304
Image (0xf64da0)
  RTTI typeinfo:   itk::Image<unsigned char, 2u>
  Reference Count: 1
  Modified Time: 400
  Debug: Off
  Observers:
    none
  Source: (0xf64140)
  Source output index: 0
  Release Data: Off
  Data Released: False
  Global Release Data: Off
  PipelineMTime: 385
  UpdateMTime: 401
  RealTimeStamp: 1.31559e+09 seconds
  LargestPossibleRegion:
    Dimension: 2
    Index: [0, 0]
    Size: [2048, 2048]
  BufferedRegion:
    Dimension: 2
    Index: [0, 0]
    Size: [2048, 2048]
  RequestedRegion:
    Dimension: 2
    Index: [0, 0]
    Size: [2048, 2048]
  Spacing: [1.09863e-07, 1.09863e-07]
  Origin: [0, 0]
  Direction:
1 0
0 1

  IndexToPointMatrix:
  1.09863e-07 0
0 1.09863e-07

  PointToIndexMatrix:
  9.10225e+06 0
0 9.10225e+06

  PixelContainer:
    ImportImageContainer (0xf64f50)
      RTTI typeinfo:   itk::ImportImageContainer<unsigned long, unsigned
char>
      Reference Count: 1
      Modified Time: 399
      Debug: Off
      Observers:
        none
      Pointer: 0x14d9850
      Container manages memory: true
      Size: 4194304
      Capacity: 4194304
Exception caught !

itk::ExceptionObject (0xf761a0)
Location: "void itk::ImageToImageFilter<TInputImage,
TOutputImage>::VerifyInputInformation() [with TInputImage =
itk::Image<std::complex<float>, 2u>, TOutputImage =
itk::Image<std::complex<float>, 2u>]"
File:
/home/kevin/kitware/ITK_Dev/ITK/Modules/Core/Common/include/itkImageToImageFilter.hxx
Line: 253
Description: itk::ERROR: MaskImageFilter(0xf72ce0): Inputs do not occupy
the same physical space!
InputImage Spacing: [1.09863e-07, 1.09863e-07], InputImage1 Spacing:
[1.09863e-07, 1.09863e-07]
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FFTMaskIFFT.cxx
Type: text/x-c++src
Size: 2887 bytes
Desc: not available
URL: <http://www.itk.org/mailman/private/insight-developers/attachments/20110909/0d8a51ca/attachment.cxx>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 253 bytes
Desc: OpenPGP digital signature
URL: <http://www.itk.org/mailman/private/insight-developers/attachments/20110909/0d8a51ca/attachment.pgp>