#include <ipPic/mitkIpPicTypeMultiplex.h>
#include "ipSegmentation.h"

Defines
#define	SEGSET(ofs)
#define	ADD_CONTOUR_POINT
Functions
template<typename PicType >
float *	tmGetContour4N (const mitkIpPicDescriptor seg, int startOfs, int &numPts, int &resSize, float result)
float *	ipMITKSegmentationGetContour4N (const mitkIpPicDescriptor seg, int startOfs, int &numPoints, int &sizeBuffer, float pointBuffer)
template<typename PicType >
float *	tmGetContour8N (const mitkIpPicDescriptor seg, int startOfs, int &numPts, int &resSize, float result)
float *	ipMITKSegmentationGetContour8N (const mitkIpPicDescriptor seg, int startOfs, int &numPoints, int &sizeBuffer, float pointBuffer)

Define Documentation

#define ADD_CONTOUR_POINT

Value:

result[numPts*2]   = xPos+xCorner[dir];      \
result[numPts*2+1] = yPos+yCorner[dir];      \
if (result[numPts*2]==result[0] && result[numPts*2+1]==result[1] && numPts>0) finished = true; \
numPts++;                    \
if (numPts==resSize) {              \
  resSize+=16+resSize/2;            \
  result = (float*)realloc( result, resSize*2*sizeof(float) ); \
  if (!result) finished = true;        \
}

Definition at line 33 of file ipSegmentationContourExtraction.cpp.

Referenced by tmGetContour4N(), and tmGetContour8N().

#define SEGSET ( ofs )

Value:

(              \
((ofs)>=0 && (ofs)<maxOfs) &&          \
( *(((PicType*)seg->data) + (ofs)) != 0 ) &&  \
( (ofs)%line != 0 || (pos+1)%line != 0 ) &&    \
( (ofs+1)%line != 0 || (pos)%line != 0 )      \
)

Definition at line 23 of file ipSegmentationContourExtraction.cpp.

Referenced by tmGetContour4N(), and tmGetContour8N().

Function Documentation

float* ipMITKSegmentationGetContour4N	(	const mitkIpPicDescriptor *	seg,
		int	startOfs,
		int &	numPoints,
		int &	sizeBuffer,
		float *	pointBuffer
	)

Definition at line 105 of file ipSegmentationContourExtraction.cpp.

References tmGetContour4N().

{
  float *newBuffer = NULL;
  mitkIpPicTypeMultiplexR4( tmGetContour4N, seg, newBuffer, startOfs, numPoints, sizeBuffer, pointBuffer );
  return newBuffer;
}

float* ipMITKSegmentationGetContour8N	(	const mitkIpPicDescriptor *	seg,
		int	startOfs,
		int &	numPoints,
		int &	sizeBuffer,
		float *	pointBuffer = `0`
	)

Same as above, but for the 8 neighbourhood contour.

Definition at line 194 of file ipSegmentationContourExtraction.cpp.

References tmGetContour8N().

Referenced by ipMITKSegmentationGetCutPoints(), mitk::SetRegionTool::OnMousePressed(), mitk::RegionGrowingTool::PerformRegionGrowingAndUpdateContour(), mitk::CorrectorAlgorithm::TobiasHeimannCorrectionAlgorithm(), and mitk::PaintbrushTool::UpdateContour().

{
  float *newBuffer = NULL;
  mitkIpPicTypeMultiplexR4( tmGetContour8N, seg, newBuffer, startOfs, numPoints, sizeBuffer, pointBuffer );
  return newBuffer;
}

template<typename PicType >

float* tmGetContour4N	(	const mitkIpPicDescriptor *	seg,
		int	startOfs,
		int &	numPts,
		int &	resSize,
		float *	result
	)

Definition at line 50 of file ipSegmentationContourExtraction.cpp.

References ADD_CONTOUR_POINT, and SEGSET.

Referenced by ipMITKSegmentationGetContour4N().

                   : DON'T TOUCH THE CODE (removed smiley, this is not funny but cruel for any maintainer!)
{
  numPts = 0; 
  int line   = seg->n[0];
  int maxOfs = seg->n[0] * seg->n[1];
  
  int straight[4]  = { 1, -line, -1, line };
  int right[4]  = { line, 1, -line, -1 };
  float xMod[4]    = { 1.0, 0.0, -1.0, 0.0 };
  float yMod[4]    = { 0.0, -1.0, 0.0, 1.0 };
  float xCorner[4]  = { 1.0, 1.0, 0.0, 0.0 };
  float yCorner[4]  = { 1.0, 0.0, 0.0, 1.0 };
  
  int dir = 0;
  int pos = startOfs;
  float xPos = (float)(pos % line);
  float yPos = (float)(pos / line);

  while ( dir<4 && SEGSET( pos+right[dir] ) ) dir++;
  if (dir==4) return result;  // no contour pixel
  
  bool finished = false;
  if (result==0) {
    resSize = 2048;
    result = (float*)malloc( resSize*2*sizeof(float) );
  }

  do {
    if ( SEGSET( pos+right[dir] ) ) {
      // modify direction (turn right):
      dir = (dir-1) & 3;
      // modify position:
      pos  += straight[dir];
      xPos += xMod[dir];
      yPos += yMod[dir];
    }
    else if ( SEGSET( pos+straight[dir] ) ) {
      ADD_CONTOUR_POINT
      // modify position:
      pos += straight[dir];
      xPos += xMod[dir];
      yPos += yMod[dir];
    }
    else {
      ADD_CONTOUR_POINT
      // modify direction (turn left):
      dir = (dir+1) & 3;
    }
  } while (!finished);
  return result;
}

template<typename PicType >

float* tmGetContour8N	(	const mitkIpPicDescriptor *	seg,
		int	startOfs,
		int &	numPts,
		int &	resSize,
		float *	result
	)

Definition at line 114 of file ipSegmentationContourExtraction.cpp.

References ADD_CONTOUR_POINT, and SEGSET.

Referenced by ipMITKSegmentationGetContour8N().

                   : DON'T TOUCH THE CODE ;-)
{
  numPts = 0; 
  int line   = seg->n[0];               // width of segmentation in pixels
  int maxOfs = seg->n[0] * seg->n[1];  // memory offset of pixel just beyond bottom-right-most pixel of segmentation (not a valid offset)
  
  int straight[4]  = { 1, -line, -1, line };  // right, top, left, down   (memory offsets)
  int right[4]  = { line, 1, -line, -1 };     // down, right, top, left
  float xMod[4]    = { 1.0, 0.0, -1.0, 0.0 };
  float yMod[4]    = { 0.0, -1.0, 0.0, 1.0 };
  float xCorner[4]  = { 1.0, 1.0, 0.0, 0.0 };
  float yCorner[4]  = { 1.0, 0.0, 0.0, 1.0 };
  
  int dir = 0;
  int pos = startOfs;    // initial position, this is where the contour search starts
  float xPos = (float)(pos % line); // calculate x and y from the memory offset
  float yPos = (float)(pos / line);

  while ( dir<4 && SEGSET( pos+right[dir] ) ) dir++;
  if (dir==4) {
    // check diagonal pixels:
    dir = 0;
    while ( dir<4 && SEGSET( pos+right[dir]+straight[dir] ) ) dir++;
    if (dir==4) return result;  // no contour pixel
    // chose next suitable neighbour:
    pos  += straight[dir];
    xPos += xMod[dir];
    yPos += yMod[dir];
  }
  
  bool finished = false;
  if (result==0) {
    resSize = 2048;
    result = (float*)malloc( resSize*2*sizeof(float) );
  }

  // here xPos,yPos are on some pixel next to a segmentation pixel. Where is "next to"? This is defined by the value of dir and the offsets in right[dir].

  // tries to complete the contour until a point is added that is identical to the first point
  do {
    if ( SEGSET( pos+right[dir] ) ) {      // valgrind complaint: jump dependent on uninitialized value (from my current understanding this could only be some pixel value outside the image
      // modify direction (turn right):    // this if will evaluate to true during the first iteration
      dir = (dir-1) & 3;   // ok, some weird logic selects a new direction
      // modify position:
      pos  += straight[dir];   // definitions of right and straight (plus xMod and yMod) lead to a counter-clockwise movement around the contour
      xPos += xMod[dir];
      yPos += yMod[dir];
    }
    else if ( SEGSET( pos+straight[dir] ) ) {
      ADD_CONTOUR_POINT
      // modify position:
      pos += straight[dir];
      xPos += xMod[dir];
      yPos += yMod[dir];
    }
    else if ( SEGSET( pos+right[dir]+straight[dir] ) ) {   // valgrind complaint: jump dependent on uninitialized value
      ADD_CONTOUR_POINT
      // modify position:
      pos += straight[dir];
      xPos += xMod[dir];
      yPos += yMod[dir];
      // modify direction (turn right):
      dir = (dir-1) & 3;
      // modify position second time:
      pos += straight[dir];
      xPos += xMod[dir];
      yPos += yMod[dir];
    }
    else {
      ADD_CONTOUR_POINT
      // modify direction (turn left):
      dir = (dir+1) & 3;
    }
  } while (!finished);
  return result;
}

ipSegmentationContourExtraction.cpp File Reference

Defines

Functions

Define Documentation

Function Documentation