gradient->TransformIndexToPhysicalPoint( neighbor, neighbor_P );
if( double( i_P.EuclideanDistanceTo( neighbor_P ) ) > this->m_MaxRadius )
continue;
-
+
// Normalize offset in terms of a l1 (manhattan) distance
TOffset offset = neighbor - i;
// std::cout << "Offset: " << offset << std::endl;
/*
- int max_off = 0;
- for( unsigned int o = 0; o < Self::Dimension; ++o )
+ int max_off = 0;
+ for( unsigned int o = 0; o < Self::Dimension; ++o )
max_off += std::abs( offset[ o ] );
- if( max_off == 0 )
+ if( max_off == 0 )
continue;
- bool normalized = true;
- TOffset normalized_offset;
- for( unsigned int o = 0; o < Self::Dimension && normalized; ++o )
- {
+ bool normalized = true;
+ TOffset normalized_offset;
+ for( unsigned int o = 0; o < Self::Dimension && normalized; ++o )
+ {
if( offset[ o ] % max_off == 0 )
- normalized_offset[ o ] = offset[ o ] / max_off;
+ normalized_offset[ o ] = offset[ o ] / max_off;
else
- normalized = false;
+ normalized = false;
- } // rof
- if( !normalized )
+ } // rof
+ if( !normalized )
normalized_offset = offset;
- std::cout << "offset: " << normalized_offset << std::endl;
+ std::cout << "offset: " << normalized_offset << std::endl;
- // Update profiles
- profiles[ normalized_offset ].push_back( neighbor );
+ // Update profiles
+ profiles[ normalized_offset ].push_back( neighbor );
*/
// Update queue