/usr/include/ITK-4.5/itkAnchorUtilities.hxx is in libinsighttoolkit4-dev 4.5.0-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 | /*=========================================================================
*
* Copyright Insight Software Consortium
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#ifndef __itkAnchorUtilities_hxx
#define __itkAnchorUtilities_hxx
#include "itkAnchorUtilities.h"
#include "itkImageRegionConstIteratorWithIndex.h"
#include "itkImageRegionConstIterator.h"
#include "itkNeighborhoodAlgorithm.h"
namespace itk
{
/**
* AnchorUtilities
* functionality in common for anchor openings/closings and
* erosions/dilation
*
**/
template< typename TImage, typename TBres, typename TAnchor, typename TLine >
void DoAnchorFace(const TImage *input,
TImage *output,
typename TImage::PixelType border,
TLine line,
TAnchor & AnchorLine,
const typename TBres::OffsetArray LineOffsets,
std::vector<typename TImage::PixelType> & inbuffer,
std::vector<typename TImage::PixelType> & outbuffer,
const typename TImage::RegionType AllImage,
const typename TImage::RegionType face)
{
// iterate over the face
// we can't use an iterator with a region outside the image. All we need here
// is to
// iterate over all the indexes of the face, without accessing the content of
// the image.
// I can't find any cleaner way, so we use a dumb image, not even allocated,
// to iterate
// over all the indexes inside the region.
//
// typedef ImageRegionConstIteratorWithIndex<TImage> ItType;
// ItType it(input, face);
typename TImage::Pointer dumbImg = TImage::New();
dumbImg->SetRegions(face);
TLine NormLine = line;
NormLine.Normalize();
// set a generous tolerance
float tol = 1.0 / LineOffsets.size();
for ( unsigned int it = 0; it < face.GetNumberOfPixels(); it++ )
{
typename TImage::IndexType Ind = dumbImg->ComputeIndex(it);
unsigned start, end;
if ( FillLineBuffer< TImage, TBres, TLine >(input, Ind, NormLine, tol, LineOffsets,
AllImage, inbuffer, start, end) )
{
const unsigned len = end - start + 1;
// compat
inbuffer[0] = border;
inbuffer[len + 1] = border;
AnchorLine.DoLine(outbuffer, inbuffer, len + 2); // compat
CopyLineToImage< TImage, TBres >(output, Ind, LineOffsets, outbuffer, start, end);
}
}
}
} // namespace itk
#endif
|