classify/mfx.cpp

Go to the documentation of this file.

/* =================
          Include Files and Type Defines
 ==================== */
#include "mfdefs.h"
#include "variables.h"
#include "sigmenu.h"
#include "mfoutline.h"
#include "clusttool.h"           //NEEDED
#include "const.h"
#include "intfx.h"
#include <math.h>

/* default values for tunable knobs */
/* old numbers corresponded to 10.0 degrees and 80.0 degrees */
/* PREV DEFAULT 0.176326981 approx. 10.0 degrees */
#define MIN_SLOPE               0.414213562
/* PREV DEFAULT 5.671281820 approx. 80.0 degrees */
#define MAX_SLOPE               2.414213562
/* no noise filtering */
#define NOISE_SEGMENT_LENGTH    (0.00)
/* no feature splitting */
#define MAX_FEATURE_LENGTH      (MAXFLOAT)

/* =================
          Macros
 ==================== */
/* miscellaneous macros */
#define NormalizeAngle(A)       ( (((A)<0)?((A)+2*PI):(A)) / (2*PI) )

/*----------------------------------------------------------------------------
          Private Function Prototypes
-----------------------------------------------------------------------------*/
void ComputeBulges(MFOUTLINE Start, MFOUTLINE End, MICROFEATURE MicroFeature);

FLOAT32 ComputeOrientation(MFEDGEPT *Start, MFEDGEPT *End);

MICROFEATURES ConvertToMicroFeatures(MFOUTLINE Outline,
                                     MICROFEATURES MicroFeatures);

MICROFEATURE ExtractMicroFeature(MFOUTLINE Start, MFOUTLINE End);

void SmearBulges(MICROFEATURES MicroFeatures, FLOAT32 XScale, FLOAT32 YScale);

/*
#if defined(__STDC__) || defined(__cplusplus)
# define _ARGS(s) s
#else
# define _ARGS(s) ()
#endif*/

/* /users/danj/wiseowl/src/danj/microfeatures/mfx.c
void ComputeBulges
  _ARGS((MFOUTLINE Start,
  MFOUTLINE End,
  MICROFEATURE MicroFeature));

FLOAT32 ComputeOrientation
  _ARGS((MFEDGEPT *Start,
  MFEDGEPT *End));

MICROFEATURES ConvertToMicroFeatures
  _ARGS((MFOUTLINE Outline,
  MICROFEATURES MicroFeatures));

MICROFEATURE ExtractMicroFeature
  _ARGS((MFOUTLINE Start,
  MFOUTLINE End));

void SmearBulges
  _ARGS((MICROFEATURES MicroFeatures,
  FLOAT32 XScale,
  FLOAT32 YScale));

#undef _ARGS
*/

/* =================
        Global Data Definitions and Declarations
 ==================== */
/* tuning knobs that can be adjusted without recompilation */
static FLOAT32 MinSlope;
static FLOAT32 MaxSlope;
static FLOAT32 NoiseSegmentLength;

/* =================
            Public Code
 ==================== */
/* =============================== */
void InitMicroFxVars() { 
  VALUE dummy;

  float_variable (MinSlope, "MinSlope", MIN_SLOPE);
  float_variable (MaxSlope, "MaxSlope", MAX_SLOPE);
  float_variable (NoiseSegmentLength, "NoiseSegmentLength",
    NOISE_SEGMENT_LENGTH);
}                                /* InitMicroFxVars */


/* =============================== */
CHAR_FEATURES BlobMicroFeatures(TBLOB *Blob, LINE_STATS *LineStats) { 
  MICROFEATURES MicroFeatures = NIL;
  FLOAT32 XScale, YScale;
  LIST Outlines;
  LIST RemainingOutlines;
  MFOUTLINE Outline;
  INT_FEATURE_ARRAY blfeatures;
  INT_FEATURE_ARRAY cnfeatures;
  INT_FX_RESULT_STRUCT results;

  if (Blob != NULL) {
    Outlines = ConvertBlob (Blob);
//    NormalizeOutlines(Outlines, LineStats, &XScale, &YScale); 
    ExtractIntFeat(Blob, blfeatures, cnfeatures, &results);
    XScale = 0.2f / results.Ry;
    YScale = 0.2f / results.Rx;

    RemainingOutlines = Outlines;
    iterate(RemainingOutlines) { 
      Outline = (MFOUTLINE) first (RemainingOutlines);
      CharNormalizeOutline (Outline,
        results.Xmean, results.Ymean,
        XScale, YScale);
    }

    RemainingOutlines = Outlines;
    iterate(RemainingOutlines) {
      Outline = (MFOUTLINE) first (RemainingOutlines);
      FindDirectionChanges(Outline, MinSlope, MaxSlope);
      FilterEdgeNoise(Outline, NoiseSegmentLength);
      MarkDirectionChanges(Outline);
      SmearExtremities(Outline, XScale, YScale);
      MicroFeatures = ConvertToMicroFeatures (Outline, MicroFeatures);
    }
    SmearBulges(MicroFeatures, XScale, YScale);
    FreeOutlines(Outlines);
  }
  return ((CHAR_FEATURES) MicroFeatures);
}                                /* BlobMicroFeatures */


/* =================
              Private Macros
 ==================== */
#define angle_of(x1,y1,x2,y2)                   \
((x2-x1) ?                                    \
   (atan2 (y2-y1, x2-x1)) :                     \
   ((y2<y1) ? (- PI / 2.0) : (PI / 2.0)))   \


#define scale_angle(x)                             \
(((x<0) ? (2.0 * PI + x) : (x)) * 0.5 / PI)  \

/*---------------------------------------------------------------------------
            Private Code
---------------------------------------------------------------------------*/
/* =============================== */
void ComputeBulges(MFOUTLINE Start, MFOUTLINE End, MICROFEATURE MicroFeature) { 
  MATRIX_2D Matrix;
  MFEDGEPT *Origin;
  MFOUTLINE SegmentStart, SegmentEnd;
  FPOINT CurrentPoint, LastPoint;
  FLOAT32 BulgePosition;

  /* check for simple case */
  if (End == NextPointAfter (Start))
    FirstBulgeOf (MicroFeature) = SecondBulgeOf (MicroFeature) = 0;
  else {
    Origin = PointAt (Start);

    InitMatrix(&Matrix);
    RotateMatrix (&Matrix, OrientationOf (MicroFeature) * -2.0 * PI);
    TranslateMatrix (&Matrix, -XPositionOf (Origin), -YPositionOf (Origin));

    SegmentEnd = Start;
    FillPoint (CurrentPoint, 0, 0);
    BulgePosition = LengthOf (MicroFeature) / 3;
    CopyPoint(CurrentPoint, LastPoint);
    while (Xof (CurrentPoint) < BulgePosition) {
      SegmentStart = SegmentEnd;
      SegmentEnd = NextPointAfter (SegmentStart);
      CopyPoint(CurrentPoint, LastPoint);
      MapPoint (&Matrix, PositionOf (PointAt (SegmentEnd)), CurrentPoint);
    }
    FirstBulgeOf (MicroFeature) =
      XIntersectionOf(LastPoint, CurrentPoint, BulgePosition);

    BulgePosition *= 2;

    // Prevents from copying the points before computing the bulge if
    // CurrentPoint will not change. (Which would cause to output nan
    // for the SecondBulge.)
    if (Xof (CurrentPoint) < BulgePosition)
      CopyPoint(CurrentPoint, LastPoint);
    while (Xof (CurrentPoint) < BulgePosition) {
      SegmentStart = SegmentEnd;
      SegmentEnd = NextPointAfter (SegmentStart);
      CopyPoint(CurrentPoint, LastPoint);
      MapPoint (&Matrix, PositionOf (PointAt (SegmentEnd)), CurrentPoint);
    }
    SecondBulgeOf (MicroFeature) =
      XIntersectionOf(LastPoint, CurrentPoint, BulgePosition);

    FirstBulgeOf (MicroFeature) /= BULGENORMALIZER *
      LengthOf(MicroFeature);
    SecondBulgeOf (MicroFeature) /= BULGENORMALIZER *
      LengthOf(MicroFeature);
  }
}                                /* ComputeBulges */


/* =============================== */
FLOAT32 ComputeOrientation(MFEDGEPT *Start, MFEDGEPT *End) { 
  FLOAT32 Orientation;

  Orientation = NormalizeAngle (AngleFrom (PositionOf (Start),
    PositionOf (End)));

  /* ensure that round-off errors do not put circular param out of range */
  if ((Orientation < 0) || (Orientation >= 1))
    Orientation = 0;
  return (Orientation);
}                                /* ComputeOrientation */


/* =============================== */
MICROFEATURES ConvertToMicroFeatures(MFOUTLINE Outline,
                                     MICROFEATURES MicroFeatures) {
  MFOUTLINE Current;
  MFOUTLINE Last;
  MFOUTLINE First;
  MICROFEATURE NewFeature;

  if (DegenerateOutline (Outline))
    return (MicroFeatures);

  First = NextExtremity (Outline);
  Last = First;
  do {
    Current = NextExtremity (Last);
    NewFeature = ExtractMicroFeature (Last, Current);
    if (NewFeature != NULL)
      MicroFeatures = push (MicroFeatures, NewFeature);
    Last = Current;
  }
  while (Last != First);

  return (MicroFeatures);
}                                /* ConvertToMicroFeatures */


/* =============================== */
MICROFEATURE ExtractMicroFeature(MFOUTLINE Start, MFOUTLINE End) { 
  MICROFEATURE NewFeature;
  MFEDGEPT *P1, *P2;

  P1 = PointAt (Start);
  P2 = PointAt (End);

  NewFeature = NewMicroFeature ();
  CenterX (NewFeature) = AverageOf (XPositionOf (P1), XPositionOf (P2));
  CenterY (NewFeature) = AverageOf (YPositionOf (P1), YPositionOf (P2));
  LengthOf (NewFeature) = DistanceBetween (PositionOf (P1), PositionOf (P2));
  OrientationOf (NewFeature) =
    NormalizedAngleFrom (&(PositionOf (P1)), &(PositionOf (P2)), 1.0);
  ComputeBulges(Start, End, NewFeature);
  return (NewFeature);
}                                /* ExtractMicroFeature */


/* =============================== */
void SmearBulges(MICROFEATURES MicroFeatures, FLOAT32 XScale, FLOAT32 YScale) { 
  MICROFEATURE MicroFeature;
  FLOAT32 MinSmear;
  FLOAT32 MaxSmear;
  FLOAT32 Cos, Sin;
  FLOAT32 Scale;

  iterate(MicroFeatures) {
    MicroFeature = NextFeatureOf (MicroFeatures);

    Cos = fabs (cos (2.0 * PI * OrientationOf (MicroFeature)));
    Sin = fabs (sin (2.0 * PI * OrientationOf (MicroFeature)));
    Scale = YScale * Cos + XScale * Sin;

    MinSmear = -0.5 * Scale / (BULGENORMALIZER * LengthOf (MicroFeature));
    MaxSmear = 0.5 * Scale / (BULGENORMALIZER * LengthOf (MicroFeature));

    FirstBulgeOf (MicroFeature) += UniformRandomNumber (MinSmear, MaxSmear);
    SecondBulgeOf (MicroFeature) += UniformRandomNumber (MinSmear, MaxSmear);
  }
}                                /* SmearBulges */

---