[Insight-users] Region growing segmentation with a border?

[Chris Turner]

Hi,

I have an image and want to segment it into different parts. Region
growing should absolutely do the trick. But my only problem is that I
also have some borders which the algorithm shouldn't cross. The
Borders are specified as a plane. First I have to convert the plane
coordinates into image coordinates. Is there an easy way to do so or
do I have to do it by "hand" (....by myself)

About the segmentation: My first thoughts are basically write my own
region growing algorithm and add besides the intensity evaluation (if
minThreshold<pixel_intensity<maxThreshold) a second evaluation which
checks if the pixel is in the specific border/region.


What do you think about this, do you have any better ideas? And also
where should I start? I thought maybe just use the
itkConnectedThresholdImageFilter (3.20) and change the iteration loop
to something like this:

[begins line 296]:

    typedef FloodFilledImageFunctionConditionalIterator<OutputImageType,
FunctionType> IteratorType;
    IteratorType it ( outputImage, function, m_SeedList );
    it.GoToBegin();

    while( !it.IsAtEnd())
      {
       if(pixelPosition > border)
          continue;
      it.Set(m_ReplaceValue);
      ++it;
      progress.CompletedPixel();  // potential exception thrown here
      }

But I don't have the pixel position at this point - right?


So how would an ITK export approach this problem?

Thanks
Chris