Skip to content

Commit 2a40128

Browse files
jbvimortjbvimort
committed
ENH: Allowing the usage of float as input image
1 parent eb45e5d commit 2a40128

13 files changed

+49
-49
lines changed

‎include/itkCoocurrenceTextureFeaturesImageFilter.h‎

Lines changed: 7 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -103,13 +103,16 @@ class ITK_TEMPLATE_EXPORT CoocurrenceTextureFeaturesImageFilter
103103
typedef TInputImage InputImageType;
104104
typedef TOutputImage OutputImageType;
105105
typedef TMaskImage MaskImageType;
106-
typedef TInputImage DigitizedImageType;
107106

108107
typedef typename InputImageType::PixelType PixelType;
109108
typedef typename MaskImageType::PixelType MaskPixelType;
110109
typedef typename InputImageType::IndexType IndexType;
111110
typedef typename InputImageType::PointType PointType;
112111

112+
113+
typedef int HistogramIndexType;
114+
typedef itk::Image< HistogramIndexType, TInputImage::ImageDimension > DigitizedImageType;
115+
113116
typedef typename InputImageType::OffsetType OffsetType;
114117
typedef VectorContainer< unsigned char, OffsetType > OffsetVector;
115118
typedef typename OffsetVector::Pointer OffsetVectorPointer;
@@ -118,9 +121,9 @@ class ITK_TEMPLATE_EXPORT CoocurrenceTextureFeaturesImageFilter
118121
typedef typename InputImageType::RegionType InputRegionType;
119122
typedef typename OutputImageType::RegionType OutputRegionType;
120123

121-
typedef typename itk::ConstNeighborhoodIterator< InputImageType > NeighborhoodIteratorType;
122-
typedef typename NeighborhoodIteratorType::RadiusType NeighborhoodRadiusType;
123-
typedef typename NeighborhoodIteratorType::NeighborIndexType NeighborIndexType;
124+
typedef typename itk::ConstNeighborhoodIterator< DigitizedImageType > NeighborhoodIteratorType;
125+
typedef typename NeighborhoodIteratorType::RadiusType NeighborhoodRadiusType;
126+
typedef typename NeighborhoodIteratorType::NeighborIndexType NeighborIndexType;
124127

125128
typedef typename NumericTraits<PixelType>::RealType MeasurementType;
126129
typedef typename NumericTraits<PixelType>::RealType RealType;
@@ -198,8 +201,6 @@ class ITK_TEMPLATE_EXPORT CoocurrenceTextureFeaturesImageFilter
198201

199202
#ifdef ITK_USE_CONCEPT_CHECKING
200203
// Begin concept checking
201-
itkConceptMacro( InputPixelTypeCheck,
202-
( Concept::IsInteger< typename InputImageType::PixelType>) );
203204
itkConceptMacro( OutputPixelTypeCheck,
204205
( Concept::IsFloatingPoint< OutputRealType > ) );
205206
// End concept checking

‎include/itkCoocurrenceTextureFeaturesImageFilter.hxx‎

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -96,7 +96,7 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
9696

9797
DigitizerFunctorType digitalizer(m_NumberOfBinsPerAxis, m_InsidePixelValue, m_HistogramMinimum, m_HistogramMaximum);
9898

99-
typedef BinaryFunctorImageFilter< MaskImageType, InputImageType, InputImageType, DigitizerFunctorType> FilterType;
99+
typedef BinaryFunctorImageFilter< MaskImageType, InputImageType, DigitizedImageType, DigitizerFunctorType> FilterType;
100100
typename FilterType::Pointer filter = FilterType::New();
101101
if (this->GetMaskImage() != ITK_NULLPTR)
102102
{
@@ -143,10 +143,10 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
143143
typename TOutputImage::PixelType outputPixel;
144144

145145
// Separation of the non-boundary region that will be processed in a different way
146-
NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< TInputImage > boundaryFacesCalculator;
147-
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< InputImageType >::FaceListType
146+
NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType > boundaryFacesCalculator;
147+
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType >::FaceListType
148148
faceList = boundaryFacesCalculator( this->m_DigitizedInputImage, outputRegionForThread, m_NeighborhoodRadius );
149-
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< InputImageType >::FaceListType::iterator fit = faceList.begin();
149+
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType>::FaceListType::iterator fit = faceList.begin();
150150

151151
// Declaration of the variables useful to iterate over the all image region
152152
bool isInImage;
@@ -166,10 +166,10 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
166166
vnl_matrix<unsigned int> hist(m_NumberOfBinsPerAxis, m_NumberOfBinsPerAxis);
167167

168168
// Declaration of the variables useful to iterate over the all neighborhood region
169-
PixelType currentInNeighborhoodPixelIntensity;
169+
HistogramIndexType currentInNeighborhoodPixelIntensity;
170170

171171
// Declaration of the variables useful to iterate over the run
172-
PixelType pixelIntensity( NumericTraits<PixelType>::ZeroValue() );
172+
HistogramIndexType pixelIntensity( NumericTraits<HistogramIndexType>::ZeroValue() );
173173
OffsetType tempOffset;
174174

175175
/// ***** Non-boundary Region *****

‎include/itkDigitizerFunctor.h‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -76,7 +76,7 @@ class Digitizer
7676
}
7777
else
7878
{
79-
return ((inputPixel - m_Min)/((m_Max-m_Min)/ (float)m_NumberOfBinsPerAxis));
79+
return Math::Floor< TOutput >((inputPixel - m_Min)/((m_Max-m_Min)/ (float)m_NumberOfBinsPerAxis));
8080
}
8181
}
8282

‎include/itkRunLengthTextureFeaturesImageFilter.h‎

Lines changed: 17 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -116,13 +116,15 @@ class ITK_TEMPLATE_EXPORT RunLengthTextureFeaturesImageFilter
116116
typedef TInputImage InputImageType;
117117
typedef TOutputImage OutputImageType;
118118
typedef TMaskImage MaskImageType;
119-
typedef TInputImage DigitizedImageType;
120119

121120
typedef typename InputImageType::PixelType PixelType;
122121
typedef typename MaskImageType::PixelType MaskPixelType;
123122
typedef typename InputImageType::IndexType IndexType;
124123
typedef typename InputImageType::PointType PointType;
125124

125+
typedef int HistogramIndexType;
126+
typedef itk::Image< HistogramIndexType, TInputImage::ImageDimension > DigitizedImageType;
127+
126128
typedef typename InputImageType::OffsetType OffsetType;
127129
typedef VectorContainer< unsigned char, OffsetType > OffsetVector;
128130
typedef typename OffsetVector::Pointer OffsetVectorPointer;
@@ -131,9 +133,9 @@ class ITK_TEMPLATE_EXPORT RunLengthTextureFeaturesImageFilter
131133
typedef typename InputImageType::RegionType InputRegionType;
132134
typedef typename OutputImageType::RegionType OutputRegionType;
133135

134-
typedef typename itk::ConstNeighborhoodIterator< InputImageType > NeighborhoodIteratorType;
135-
typedef typename NeighborhoodIteratorType::RadiusType NeighborhoodRadiusType;
136-
typedef typename NeighborhoodIteratorType::NeighborIndexType NeighborIndexType;
136+
typedef typename itk::ConstNeighborhoodIterator< DigitizedImageType > NeighborhoodIteratorType;
137+
typedef typename NeighborhoodIteratorType::RadiusType NeighborhoodRadiusType;
138+
typedef typename NeighborhoodIteratorType::NeighborIndexType NeighborIndexType;
137139

138140
typedef typename NumericTraits<PixelType>::RealType MeasurementType;
139141
typedef typename NumericTraits<PixelType>::RealType RealType;
@@ -222,8 +224,6 @@ class ITK_TEMPLATE_EXPORT RunLengthTextureFeaturesImageFilter
222224

223225
#ifdef ITK_USE_CONCEPT_CHECKING
224226
// Begin concept checking
225-
itkConceptMacro( InputPixelTypeCheck,
226-
( Concept::IsInteger< typename InputImageType::PixelType>) );
227227
itkConceptMacro( OutputPixelTypeCheck,
228228
( Concept::IsFloatingPoint< OutputRealType > ) );
229229
// End concept checking
@@ -237,7 +237,7 @@ class ITK_TEMPLATE_EXPORT RunLengthTextureFeaturesImageFilter
237237
void NormalizeOffsetDirection(OffsetType &offset);
238238
bool IsInsideNeighborhood(const OffsetType &iteratedOffset);
239239
void IncreaseHistogram(vnl_matrix<unsigned int> &hist, unsigned int &totalNumberOfRuns,
240-
const PixelType &currentInNeighborhoodPixelIntensity,
240+
const HistogramIndexType &currentInNeighborhoodPixelIntensity,
241241
const OffsetType &offset, const unsigned int &pixelDistance);
242242
void ComputeFeatures( vnl_matrix<unsigned int> &hist, const unsigned int &totalNumberOfRuns,
243243
typename TOutputImage::PixelType &outputPixel);
@@ -251,16 +251,16 @@ class ITK_TEMPLATE_EXPORT RunLengthTextureFeaturesImageFilter
251251
virtual void GenerateOutputInformation() ITK_OVERRIDE;
252252

253253
private:
254-
typename InputImageType::Pointer m_DigitizedInputImage;
255-
NeighborhoodRadiusType m_NeighborhoodRadius;
256-
OffsetVectorPointer m_Offsets;
257-
unsigned int m_NumberOfBinsPerAxis;
258-
PixelType m_HistogramValueMinimum;
259-
PixelType m_HistogramValueMaximum;
260-
RealType m_HistogramDistanceMinimum;
261-
RealType m_HistogramDistanceMaximum;
262-
MaskPixelType m_InsidePixelValue;
263-
typename TInputImage::SpacingType m_Spacing;
254+
typename DigitizedImageType::Pointer m_DigitizedInputImage;
255+
NeighborhoodRadiusType m_NeighborhoodRadius;
256+
OffsetVectorPointer m_Offsets;
257+
unsigned int m_NumberOfBinsPerAxis;
258+
PixelType m_HistogramValueMinimum;
259+
PixelType m_HistogramValueMaximum;
260+
RealType m_HistogramDistanceMinimum;
261+
RealType m_HistogramDistanceMaximum;
262+
MaskPixelType m_InsidePixelValue;
263+
typename TInputImage::SpacingType m_Spacing;
264264
};
265265
} // end of namespace Statistics
266266
} // end of namespace itk

‎include/itkRunLengthTextureFeaturesImageFilter.hxx‎

Lines changed: 8 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -97,7 +97,7 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
9797

9898
DigitizerFunctorType digitalizer(m_NumberOfBinsPerAxis, m_InsidePixelValue, m_HistogramValueMinimum, m_HistogramValueMaximum);
9999

100-
typedef BinaryFunctorImageFilter< MaskImageType, InputImageType, InputImageType, DigitizerFunctorType> FilterType;
100+
typedef BinaryFunctorImageFilter< MaskImageType, InputImageType, DigitizedImageType, DigitizerFunctorType> FilterType;
101101
typename FilterType::Pointer filter = FilterType::New();
102102
if (this->GetMaskImage() != ITK_NULLPTR)
103103
{
@@ -168,10 +168,10 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
168168
alreadyVisitedImage->Allocate();
169169

170170
// Separation of the non-boundary region that will be processed in a different way
171-
NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< TInputImage > boundaryFacesCalculator;
172-
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< InputImageType >::FaceListType
171+
NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType > boundaryFacesCalculator;
172+
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType >::FaceListType
173173
faceList = boundaryFacesCalculator( this->m_DigitizedInputImage, outputRegionForThread, m_NeighborhoodRadius );
174-
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< InputImageType >::FaceListType::iterator fit = faceList.begin();
174+
typename NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType >::FaceListType::iterator fit = faceList.begin();
175175

176176
// Declaration of the variables useful to iterate over the all image region
177177
bool isInImage;
@@ -186,10 +186,10 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
186186

187187

188188
// Declaration of the variables useful to iterate over the all neighborhood region
189-
PixelType currentInNeighborhoodPixelIntensity;
189+
HistogramIndexType currentInNeighborhoodPixelIntensity;
190190

191191
// Declaration of the variables useful to iterate over the run
192-
PixelType pixelIntensity( NumericTraits<PixelType>::ZeroValue() );
192+
HistogramIndexType pixelIntensity( NumericTraits<HistogramIndexType>::ZeroValue() );
193193
OffsetType iteratedOffset;
194194
OffsetType tempOffset;
195195
unsigned int pixelDistance;
@@ -263,7 +263,7 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
263263
// Special attention paid to boundaries of bins.
264264
// For the last bin, it is left close and right close (following the previous
265265
// gerrit patch). For all other bins, the bin is left close and right open.
266-
if ( pixelIntensity == currentInNeighborhoodPixelIntensity )
266+
if (pixelIntensity == currentInNeighborhoodPixelIntensity )
267267
{
268268
alreadyVisitedImage->SetPixel( boolCurentInNeighborhoodIndex + iteratedOffset, true );
269269
++pixelDistance;
@@ -359,7 +359,7 @@ template<typename TInputImage, typename TOutputImage, typename TMaskImage>
359359
void
360360
RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
361361
::IncreaseHistogram(vnl_matrix<unsigned int> &histogram, unsigned int &totalNumberOfRuns,
362-
const PixelType &currentInNeighborhoodPixelIntensity,
362+
const HistogramIndexType &currentInNeighborhoodPixelIntensity,
363363
const OffsetType &offset, const unsigned int &pixelDistance)
364364
{
365365
float offsetDistance = 0;
@@ -428,7 +428,6 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
428428
shortRunHighGreyLevelEmphasis += ( frequency * i2 / j2 );
429429
longRunLowGreyLevelEmphasis += ( frequency * j2 / i2 );
430430
longRunHighGreyLevelEmphasis += ( frequency * i2 * j2 );
431-
432431
}
433432
}
434433
greyLevelNonuniformity =

‎test/CoocurrenceTextureFeaturesImageFilterTest.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -39,7 +39,7 @@ int CoocurrenceTextureFeaturesImageFilterTest( int argc, char *argv[] )
3939
const unsigned int VectorComponentDimension = 8;
4040

4141
// Declare types
42-
typedef int InputPixelType;
42+
typedef float InputPixelType;
4343
typedef float OutputPixelComponentType;
4444
typedef itk::Vector< OutputPixelComponentType, VectorComponentDimension >
4545
OutputPixelType;

‎test/CoocurrenceTextureFeaturesImageFilterTestVectorImageSeparateFeatures.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -47,7 +47,7 @@ int CoocurrenceTextureFeaturesImageFilterTestVectorImageSeparateFeatures( int ar
4747
const unsigned int VectorComponentDimension = 8;
4848

4949
// Declare types
50-
typedef int InputPixelType;
50+
typedef float InputPixelType;
5151
typedef float OutputPixelType;
5252

5353
typedef itk::Image< InputPixelType, ImageDimension > InputImageType;

‎test/CoocurrenceTextureFeaturesImageFilterTestWithoutMask.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -43,7 +43,7 @@ int CoocurrenceTextureFeaturesImageFilterTestWithoutMask( int argc, char *argv[]
4343
const unsigned int VectorComponentDimension = 8;
4444

4545
// Declare types
46-
typedef int InputPixelType;
46+
typedef float InputPixelType;
4747
typedef float OutputPixelComponentType;
4848
typedef itk::Vector< OutputPixelComponentType, VectorComponentDimension >
4949
OutputPixelType;

‎test/RunLengthTextureFeaturesImageFilterTest.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -46,7 +46,7 @@ int RunLengthTextureFeaturesImageFilterTest( int argc, char *argv[] )
4646
const unsigned int VectorComponentDimension = 10;
4747

4848
// Declare types
49-
typedef int InputPixelType;
49+
typedef float InputPixelType;
5050
typedef float OutputPixelComponentType;
5151
typedef itk::Vector< OutputPixelComponentType, VectorComponentDimension >
5252
OutputPixelType;

‎test/RunLengthTextureFeaturesImageFilterTestSeparateFeatures.cxx‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -47,7 +47,7 @@ int RunLengthTextureFeaturesImageFilterTestSeparateFeatures( int argc, char *arg
4747
const unsigned int VectorComponentDimension = 10;
4848

4949
// Declare types
50-
typedef int InputPixelType;
50+
typedef float InputPixelType;
5151
typedef float OutputPixelComponentType;
5252
typedef itk::Vector< OutputPixelComponentType, VectorComponentDimension >
5353
OutputPixelType;

0 commit comments

Comments
 (0)