classify/picofeat.cpp

Go to the documentation of this file.

/* =================
          Include Files and Type Defines
 ==================== */
#include "picofeat.h"
#include "mfoutline.h"
#include "variables.h"
#include "sigmenu.h"
#include "hideedge.h"
#include "fpoint.h"

#include <math.h>

#include "ocrfeatures.h"         //Debug
#include <stdio.h>               //Debug
#include "efio.h"                //Debug
//#include "christydbg.h"

#define PICO_FEATURE_LENGTH 0.05

/*---------------------------------------------------------------------------
          Private Function Prototypes
----------------------------------------------------------------------------*/
void ConvertSegmentToPicoFeat(FPOINT *Start,
                              FPOINT *End,
                              FEATURE_SET FeatureSet);

void ConvertToPicoFeatures2(MFOUTLINE Outline, FEATURE_SET FeatureSet); 

void NormalizePicoX(FEATURE_SET FeatureSet); 

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

/* /users/danj/wiseowl/src/danj/microfeatures/picofeat.c
void ConvertSegmentToPicoFeat
  _ARGS((FPOINT *Start,
  FPOINT *End,
  FEATURE_SET FeatureSet));

void ConvertToPicoFeatures2
  _ARGS((MFOUTLINE Outline,
  FEATURE_SET FeatureSet));

void NormalizePicoX
  _ARGS((FEATURE_SET FeatureSet));

#undef _ARGS
*/

/* =================
 Global Data Definitions and Declarations
 ==================== */

/* =================
     Public Code
 ==================== */
FEATURE_SET ExtractPicoFeatures(TBLOB *Blob, LINE_STATS *LineStats) { 
  LIST Outlines;
  LIST RemainingOutlines;
  MFOUTLINE Outline;
  FEATURE_SET FeatureSet;
  FLOAT32 XScale, YScale;

  FeatureSet = NewFeatureSet (MAX_PICO_FEATURES);

  Outlines = ConvertBlob (Blob);

  NormalizeOutlines(Outlines, LineStats, &XScale, &YScale); 
  RemainingOutlines = Outlines;
  iterate(RemainingOutlines) { 
    Outline = (MFOUTLINE) first (RemainingOutlines);
    /*---------Debug--------------------------------------------------*
    OFile = fopen ("f:/ims/debug/pfOutline.logCPP", "r");
    if (OFile == NULL)
    {
      OFile = Efopen ("f:/ims/debug/pfOutline.logCPP", "w");
      WriteOutline(OFile, Outline);
    }
    else
    {
      fclose (OFile);
      OFile = Efopen ("f:/ims/debug/pfOutline.logCPP", "a");
    }
    WriteOutline(OFile, Outline);
    fclose (OFile);
    *--------------------------------------------------------------------*/
    ConvertToPicoFeatures2(Outline, FeatureSet); 
  }
  if (NormMethod == baseline)
    NormalizePicoX(FeatureSet); 
  /*---------Debug--------------------------------------------------*
  File = fopen ("f:/ims/debug/pfFeatSet.logCPP", "r");
  if (File == NULL)
  {
    File = Efopen ("f:/ims/debug/pfFeatSet.logCPP", "w");
    WriteFeatureSet(File, FeatureSet);
  }
  else
  {
    fclose (File);
    File = Efopen ("f:/ims/debug/pfFeatSet.logCPP", "a");
  }
  WriteFeatureSet(File, FeatureSet);
  fclose (File);
  *--------------------------------------------------------------------*/
  FreeOutlines(Outlines); 
  return (FeatureSet);

}                                /* ExtractPicoFeatures */


/* =============================== */
void InitPicoFXVars() { 

   VALUE dummy;
   
   float_variable (PicoFeatureLength, "PicoFeatureLength",
      PICO_FEATURE_LENGTH);
   
}  /* InitPicoFXVars */

/* =================
    Private Code
 ==================== */
void ConvertSegmentToPicoFeat(FPOINT *Start,
                              FPOINT *End,
                              FEATURE_SET FeatureSet) {
  FEATURE Feature;
  FLOAT32 Angle;
  FLOAT32 Length;
  int NumFeatures;
  FPOINT Center;
  FPOINT Delta;
  int i;

  Angle = NormalizedAngleFrom (Start, End, 1.0);
  Length = DistanceBetween (*Start, *End);
  NumFeatures = (int) floor (Length / PicoFeatureLength + 0.5);
  if (NumFeatures < 1)
    NumFeatures = 1;

  /* compute vector for one pico feature */
  Xof (Delta) = XDelta (*Start, *End) / NumFeatures;
  Yof (Delta) = YDelta (*Start, *End) / NumFeatures;

  /* compute position of first pico feature */
  Xof (Center) = Xof (*Start) + Xof (Delta) / 2.0;
  Yof (Center) = Yof (*Start) + Yof (Delta) / 2.0;

  /* compute each pico feature in segment and add to feature set */
  for (i = 0; i < NumFeatures; i++) {
    Feature = NewFeature (&PicoFeatDesc);
    ParamOf (Feature, PicoFeatDir) = Angle;
    ParamOf (Feature, PicoFeatX) = Xof (Center);
    ParamOf (Feature, PicoFeatY) = Yof (Center);
    AddFeature(FeatureSet, Feature); 

    Xof (Center) += Xof (Delta);
    Yof (Center) += Yof (Delta);
  }
}                                /* ConvertSegmentToPicoFeat */


/* =============================== */
void ConvertToPicoFeatures2(MFOUTLINE Outline, FEATURE_SET FeatureSet) { 
  MFOUTLINE Next;
  MFOUTLINE First;
  MFOUTLINE Current;

  if (DegenerateOutline (Outline))
    return;

  First = Outline;
  Current = First;
  Next = NextPointAfter (Current);
  do {
    if (IsVisible (PointAt (Next)))
      ConvertSegmentToPicoFeat (&(PositionOf (PointAt (Current))),
        &(PositionOf (PointAt (Next))), FeatureSet);

    Current = Next;
    Next = NextPointAfter (Current);
  }
  while (Current != First);

}                                /* ConvertToPicoFeatures2 */


/* =============================== */
void NormalizePicoX(FEATURE_SET FeatureSet) { 
  int i;
  FEATURE Feature;
  FLOAT32 Origin = 0.0;

  for (i = 0; i < NumFeaturesIn (FeatureSet); i++) {
    Feature = FeatureIn (FeatureSet, i);
    Origin += ParamOf (Feature, PicoFeatX);
  }
  Origin /= NumFeaturesIn (FeatureSet);

  for (i = 0; i < NumFeaturesIn (FeatureSet); i++) {
    Feature = FeatureIn (FeatureSet, i);
    ParamOf (Feature, PicoFeatX) -= Origin;
  }
}                                /* NormalizePicoX */

---