@@ -124,6 +124,11 @@ runRieszWaveletPhaseAnalysisTest( const std::string& inputImage,
124124 for ( unsigned int i = 0 ; i < forwardWavelet->GetNumberOfOutputs (); ++i )
125125 {
126126 std::cout << " Output #: " << i << " / " << numberOfOutputs - 1 << std::endl;
127+ // if( i == numberOfOutputs - 1 ) // TODO Held does not modify approx image, but it does not generate better results.
128+ // {
129+ // modifiedWavelets.push_back( analysisWavelets[i] );
130+ // continue;
131+ // }
127132 typename MonogenicSignalFrequencyFilterType::Pointer monoFilter =
128133 MonogenicSignalFrequencyFilterType::New ();
129134 typename VectorInverseFFTType::Pointer vecInverseFFT =
@@ -142,10 +147,6 @@ runRieszWaveletPhaseAnalysisTest( const std::string& inputImage,
142147
143148 phaseAnalyzer->SetInput ( vecInverseFFT->GetOutput () );
144149 phaseAnalyzer->SetApplySoftThreshold ( applySoftThreshold );
145- // if( i == numberOfOutputs - 1 ) // TODO Held does not modify approx image, but it does not generate better results.
146- // {
147- // phaseAnalyzer->SetApplySoftThreshold( false );
148- // }
149150 phaseAnalyzer->Update ();
150151
151152 fftForwardPhaseFilter->SetInput ( phaseAnalyzer->GetOutputCosPhase () );
@@ -193,7 +194,7 @@ runRieszWaveletPhaseAnalysisTest( const std::string& inputImage,
193194#endif
194195
195196 // Cast To Float for save as tiff.
196- typedef itk::Image< float , Dimension > ImageFloatType;
197+ typedef itk::Image< float , Dimension > ImageFloatType;
197198 typedef itk::CastImageFilter< ImageType, ImageFloatType> CastFloatType;
198199 typename CastFloatType::Pointer caster = CastFloatType::New ();
199200 caster->SetInput (inverseFFT->GetOutput ());
0 commit comments