now storing points and landmarks directly

Josh Klontz
1 parent 8726eb4e
Showing 15 changed files with 118 additions and 144 deletions
app/examples/face_recognition.cpp
app/openbr-gui/templateviewer.cpp
sdk/openbr_plugin.cpp
sdk/openbr_plugin.h
sdk/plugins/cascade.cpp
sdk/plugins/crop.cpp
sdk/plugins/draw.cpp
sdk/plugins/eyes.cpp
sdk/plugins/hist.cpp
sdk/plugins/keypoint.cpp
sdk/plugins/meta.cpp
sdk/plugins/random.cpp
sdk/plugins/regions.cpp
sdk/plugins/register.cpp
sdk/plugins/wavelet.cpp
@@ -30,10 +30,9 @@
 static void printTemplate(const br::Template &t)
 {
-    printf("%s eyes: (%d, %d) (%d, %d)\n",
-           qPrintable(t.file.fileName()),
-           t.file.get<int>("Affine_0_X"), t.file.get<int>("Affine_0_Y"),
-           t.file.get<int>("Affine_1_X"), t.file.get<int>("Affine_1_Y"));
+    const QPoint firstEye = t.file.get<QPoint>("Affine_0");
+    const QPoint secondEye = t.file.get<QPoint>("Affine_1");
+    printf("%s eyes: (%d, %d) (%d, %d)\n", qPrintable(t.file.fileName()), firstEye.x(), firstEye.y(), secondEye.x(), secondEye.y());
 }
 int main(int argc, char *argv[])
@@ -39,10 +39,8 @@ void TemplateViewer::setFile(const File &amp;file_)
     // Update landmarks
     landmarks.clear();
-    if (file.contains("Affine_0_X") && file.contains("Affine_0_Y"))
-        landmarks.append(QPointF(file.get<float>("Affine_0_X"), file.get<float>("Affine_0_Y")));
-    if (file.contains("Affine_1_X") && file.contains("Affine_1_Y"))
-        landmarks.append(QPointF(file.get<float>("Affine_1_X"), file.get<float>("Affine_1_Y")));
+    if (file.contains("Affine_0")) landmarks.append(file.get<QPointF>("Affine_0"));
+    if (file.contains("Affine_1")) landmarks.append(file.get<QPointF>("Affine_1"));
     while (landmarks.size() < NumLandmarks)
         landmarks.append(QPointF());
     nearestLandmark = -1;
@@ -148,84 +148,72 @@ float File::label() const
     return ok ? val : -1;
 }
-QList<QPointF> File::landmarks() const
+QList<QPointF> File::namedPoints() const
 {
     QList<QPointF> landmarks;
-    foreach (const QVariant &landmark, value("Landmarks").toList())
-        landmarks.append(landmark.toPointF());
-    return landmarks;
-}
-
-QList<QPointF> File::namedLandmarks() const
-{
-    QList<QPointF> landmarks;
-    QStringList keys = localMetadata().keys();
-    foreach (const QString &key, keys) {
-        if (!key.endsWith("_X"))
-            continue;
-        QString keyBaseName = key.left(key.size()-2);
-        if (!keys.contains(keyBaseName+"_Y") ||
-            keys.contains(keyBaseName+"_Width") ||
-            keys.contains(keyBaseName+"_Height") ||
-            keys.contains(keyBaseName+"_Radius"))
-            continue;
-        landmarks.append(QPointF(get<float>(keyBaseName+"_X"), get<float>(keyBaseName+"_Y")));
+    foreach (const QString &key, localMetadata().keys()) {
+        const QVariant &variant = m_metadata[key];
+        if (variant.canConvert<QPointF>())
+            landmarks.append(variant.value<QPointF>());
     }
     return landmarks;
 }
-void File::appendLandmark(const QPointF &landmark)
+QList<QPointF> File::points() const
 {
-    QList<QVariant> newLandmarks = m_metadata["Landmarks"].toList();
-    newLandmarks.append(landmark);
-    m_metadata["Landmarks"] = newLandmarks;
+    QList<QPointF> points;
+    foreach (const QVariant &point, m_metadata["Points"].toList())
+        points.append(point.toPointF());
+    return points;
 }
-void File::appendLandmarks(const QList<QPointF> &landmarks)
+void File::appendPoint(const QPointF &point)
 {
-    QList<QVariant> newLandmarks = m_metadata["Landmarks"].toList();
-    foreach (const QPointF &landmark, landmarks)
-        newLandmarks.append(landmark);
-    m_metadata["Landmarks"] = newLandmarks;
+    QList<QVariant> newPoints = m_metadata["Points"].toList();
+    newPoints.append(point);
+    m_metadata["Points"] = newPoints;
 }
-void File::setLandmarks(const QList<QPointF> &landmarks)
+void File::appendPoints(const QList<QPointF> &points)
 {
-    QList<QVariant> landmarkList; landmarkList.reserve(landmarks.size());
-    foreach (const QPointF &landmark, landmarks)
-        landmarkList.append(landmark);
-    m_metadata["Landmarks"] = landmarkList;
+    QList<QVariant> newPoints = m_metadata["Points"].toList();
+    foreach (const QPointF &point, points)
+        newPoints.append(point);
+    m_metadata["Points"] = newPoints;
 }
-QList<QRectF> File::ROIs() const
+QList<QRectF> File::namedRects() const
 {
-    QList<QRectF> ROIs;
-    foreach (const QVariant &ROI, value("ROIs").toList())
-        ROIs.append(ROI.toRect());
-    return ROIs;
+    QList<QRectF> rects;
+    foreach (const QString &key, localMetadata().keys()) {
+        const QVariant &variant = m_metadata[key];
+        if (variant.canConvert<QRectF>())
+            rects.append(variant.value<QRectF>());
+    }
+    return rects;
 }
-void File::appendROI(const QRectF &ROI)
+QList<QRectF> File::rects() const
 {
-    QList<QVariant> newROIs = m_metadata["ROIs"].toList();
-    newROIs.append(ROI);
-    m_metadata["ROIs"] = newROIs;
+    QList<QRectF> rects;
+    foreach (const QVariant &rect, m_metadata["Rects"].toList())
+        rects.append(rect.toRect());
+    return rects;
 }
-void File::appendROIs(const QList<QRectF> &ROIs)
+void File::appendRect(const QRectF &rect)
 {
-    QList<QVariant> newROIs = m_metadata["ROIs"].toList();
-    foreach (const QRectF &ROI, ROIs)
-        newROIs.append(ROI);
-    m_metadata["ROIs"] = newROIs;
+    QList<QVariant> newRects = m_metadata["Rects"].toList();
+    newRects.append(rect);
+    m_metadata["Rects"] = newRects;
 }
-void File::setROIs(const QList<QRectF> &ROIs)
+void File::appendRects(const QList<QRectF> &rects)
 {
-    QList<QVariant> ROIList; ROIList.reserve(ROIs.size());
-    foreach (const QRectF &ROI, ROIs)
-        ROIList.append(ROI);
-    m_metadata["ROIs"] = ROIList;
+    QList<QVariant> newRects = m_metadata["Rects"].toList();
+    foreach (const QRectF &rect, rects)
+        newRects.append(rect);
+    m_metadata["Rects"] = newRects;
 }
 /* File - private methods */
@@ -128,8 +128,6 @@ void reset_##NAME() { NAME = DEFAULT; }
  *
  * Key             | Value          | Description
  * ---             | ----           | -----------
- * path            | QString        | Resolve complete file paths from file names
- * enrollAll       | bool           | Enroll zero or more templates per file
  * separator       | QString        | Seperate #name into multiple files
  * Index           | int            | Index of a template in a template list
  * Label           | float          | Classification/Regression class
@@ -145,8 +143,8 @@ void reset_##NAME() { NAME = DEFAULT; }
  * Roll            | float          | Pose
  * Pitch           | float          | Pose
  * Yaw             | float          | Pose
- * Landmarks       | QList<QPointF> | Landmark list
- * ROIs            | QList<Rect>    | Region Of Interest (ROI) list
+ * Points          | QList<QPointF> | List of unnamed points
+ * Rects           | QList<Rect>    | List of unnamed rects
  * Age             | QString        | Age used for demographic filtering
  * _*              | *              | Reserved for internal use
  */
@@ -225,18 +223,19 @@ struct BR_EXPORT File
     inline void setLabel(float label) { set("Label", label); } /*!< \brief Convenience function for setting the file's \c Label. */
     inline bool failed() const { return get<bool>("FTE", false) || get<bool>("FTO", false); } /*!< \brief Returns \c true if the file failed to open or enroll, \c false otherwise. */
-    QList<QPointF> landmarks() const; /*!< \brief Returns the file's landmark list. */
-    QList<QPointF> namedLandmarks() const; /*!< \brief Returns landmarks derived from metadata keys. */
-    void appendLandmark(const QPointF &landmark); /*!< \brief Adds a landmark to the file's landmark list. */
-    void appendLandmarks(const QList<QPointF> &landmarks); /*!< \brief Adds landmarks to the file's landmark list. */
-    inline void clearLandmarks() { m_metadata["Landmarks"] = QList<QVariant>(); } /*!< \brief Clears the file's landmark list. */
-    void setLandmarks(const QList<QPointF> &landmarks); /*!< \brief Assigns the file's landmark list. */
-
-    QList<QRectF> ROIs() const; /*!< \brief Returns the file's ROI list. */
-    void appendROI(const QRectF &ROI); /*!< \brief Adds a ROI to the file's ROI list. */
-    void appendROIs(const QList<QRectF> &ROIs); /*!< \brief Adds ROIs to the file's ROI list. */
-    inline void clearROIs() { m_metadata["ROIs"] = QList<QVariant>(); } /*!< \brief Clears the file's landmark list. */
-    void setROIs(const QList<QRectF> &ROIs); /*!< \brief Assigns the file's landmark list. */
+    QList<QPointF> namedPoints() const; /*!< \brief Returns points convertible from metadata keys. */
+    QList<QPointF> points() const; /*!< \brief Returns the file's points list. */
+    void appendPoint(const QPointF &point); /*!< \brief Adds a point to the file's point list. */
+    void appendPoints(const QList<QPointF> &points); /*!< \brief Adds landmarks to the file's landmark list. */
+    inline void clearPoints() { m_metadata["Points"] = QList<QVariant>(); } /*!< \brief Clears the file's landmark list. */
+    inline void setPoints(const QList<QPointF> &points) { clearPoints(); appendPoints(points); } /*!< \brief Overwrites the file's landmark list. */
+
+    QList<QRectF> namedRects() const; /*!< \brief Returns rects convertible from metadata values. */
+    QList<QRectF> rects() const; /*!< \brief Returns the file's rects list. */
+    void appendRect(const QRectF &rect); /*!< \brief Adds a rect to the file's rect list. */
+    void appendRects(const QList<QRectF> &rects); /*!< \brief Adds rects to the file's rect list. */
+    inline void clearRects() { m_metadata["Rects"] = QList<QVariant>(); } /*!< \brief Clears the file's rect list. */
+    inline void setRects(const QList<QRectF> &rects) { clearRects(); appendRects(rects); } /*!< \brief Overwrites the file's rect list. */
 private:
     QMap<QString,QVariant> m_metadata;
@@ -83,7 +83,7 @@ class CascadeTransform : public UntrainableTransform
         foreach (const Rect &rect, rects) {
             dst += src;
-            dst.file.appendROI(OpenCVUtils::fromRect(rect));
+            dst.file.appendRect(OpenCVUtils::fromRect(rect));
         }
     }
 };
@@ -26,7 +26,7 @@ namespace br
 /*!
  * \ingroup transforms
- * \brief Crops the regions of interest.
+ * \brief Crops the rectangular regions of interest.
  * \author Josh Klontz \cite jklontz
  */
 class ROITransform : public UntrainableTransform
@@ -35,8 +35,8 @@ class ROITransform : public UntrainableTransform
     void project(const Template &src, Template &dst) const
     {
-        foreach (const QRectF ROI, src.file.ROIs())
-            dst += src.m()(OpenCVUtils::toRect(ROI));
+        foreach (const QRectF &rect, src.file.rects())
+            dst += src.m()(OpenCVUtils::toRect(rect));
     }
 };
@@ -47,19 +47,19 @@ class DrawTransform : public UntrainableTransform
         const Scalar verboseColor(255, 255, 0);
         dst = src.m().clone();
-        QList<Point2f> landmarks = OpenCVUtils::toPoints(src.file.landmarks());
+        QList<Point2f> landmarks = OpenCVUtils::toPoints(src.file.points());
         if (unnamed) {
             foreach (const Point2f &landmark, landmarks)
                 circle(dst, landmark, 3, color, -1);
         }
         if (named) {
-            QList<Point2f> namedLandmarks = OpenCVUtils::toPoints(src.file.namedLandmarks());
+            QList<Point2f> namedLandmarks = OpenCVUtils::toPoints(src.file.namedPoints());
             foreach (const Point2f &landmark, namedLandmarks)
                 circle(dst, landmark, 3, color);
         }
         if (ROI) {
-            QList<Rect> ROIs = OpenCVUtils::toRects(src.file.ROIs());
+            QList<Rect> ROIs = OpenCVUtils::toRects(src.file.rects());
             foreach (const Rect ROI, ROIs)
                 rectangle(dst, ROI, color);
         }
@@ -154,11 +154,11 @@ class EditTransform : public UntrainableTransform
     {
         (void) event;
         if (flags) {
-            QList<QRectF> ROIs = currentTemplate.file.ROIs();
-            for (int i=ROIs.size()-1; i>=0; i--)
-                if (ROIs[i].contains(x,y))
-                    ROIs.removeAt(i);
-            currentTemplate.file.setROIs(ROIs);
+            QList<QRectF> rects = currentTemplate.file.rects();
+            for (int i=rects.size()-1; i>=0; i--)
+                if (rects[i].contains(x,y))
+                    rects.removeAt(i);
+            currentTemplate.file.setRects(rects);
         }
         Template temp;
@@ -147,7 +147,7 @@ public:
 private:
     void project(const Template &src, Template &dst) const
     {
-        Rect roi = OpenCVUtils::toRect(src.file.ROIs().first());
+        Rect roi = OpenCVUtils::toRect(src.file.rects().first());
         Mat gray;
         OpenCVUtils::cvtGray(src.m()(roi), gray);
@@ -183,12 +183,10 @@ private:
         float second_eye_y = (right_rect.y + maxLoc.y)*gray.rows/height+roi.y;
         dst = src;
-        dst.file.appendLandmark(QPointF(first_eye_x, first_eye_y));
-        dst.file.appendLandmark(QPointF(second_eye_x, second_eye_y));
-        dst.file.set("ASEF_Right_Eye_X", first_eye_x);
-        dst.file.set("ASEF_Right_Eye_Y", first_eye_y);
-        dst.file.set("ASEF_Left_Eye_X", second_eye_x);
-        dst.file.set("ASEF_Left_Eye_Y", second_eye_y);
+        dst.file.appendPoint(QPointF(first_eye_x, first_eye_y));
+        dst.file.appendPoint(QPointF(second_eye_x, second_eye_y));
+        dst.file.set("ASEF_Right_Eye", QPointF(first_eye_x, first_eye_y));
+        dst.file.set("ASEF_Left_Eye", QPointF(second_eye_x, second_eye_y));
     }
 };
@@ -174,7 +174,7 @@ class VarianceChangeDetectorTransform : public UntrainableTransform
         int *buffer = new int[bins];
         float bestRatio = -std::numeric_limits<float>::max();
-        QRect bestROI;
+        QRectF bestRect;
         const int rows = m.rows;
         const int cols = m.cols/bins;
@@ -202,7 +202,7 @@ class VarianceChangeDetectorTransform : public UntrainableTransform
                     if (ratio > bestRatio) {
                         bestRatio = ratio;
-                        bestROI = QRect(j*radius, i*radius, scale*radius, scale*radius);
+                        bestRect = QRect(j*radius, i*radius, scale*radius, scale*radius);
                     }
                 }
             }
@@ -210,7 +210,7 @@ class VarianceChangeDetectorTransform : public UntrainableTransform
         }
         delete[] buffer;
-        dst.file.appendROI(bestROI);
+        dst.file.appendRect(bestRect);
         dst.file.setLabel(bestRatio);
     }
 };
@@ -57,10 +57,10 @@ class KeyPointDetectorTransform : public UntrainableTransform
             dst.file.set("FTE", true);
         }
-        QList<Rect> ROIs;
+        QList<Rect> rects;
         foreach (const KeyPoint &keyPoint, keyPoints)
-            ROIs.append(Rect(keyPoint.pt.x, keyPoint.pt.y, keyPoint.size, keyPoint.size));
-        dst.file.setROIs(OpenCVUtils::fromRects(ROIs));
+            rects.append(Rect(keyPoint.pt.x, keyPoint.pt.y, keyPoint.size, keyPoint.size));
+        dst.file.setRects(OpenCVUtils::fromRects(rects));
     }
 };
@@ -92,10 +92,10 @@ class KeyPointDescriptorTransform : public UntrainableTransform
     {
         std::vector<KeyPoint> keyPoints;
         if (size == -1) {
-            foreach (const QRectF &ROI, src.file.ROIs())
+            foreach (const QRectF &ROI, src.file.rects())
                 keyPoints.push_back(KeyPoint(ROI.x(), ROI.y(), (ROI.width() + ROI.height())/2));
         } else {
-            foreach (const QPointF &landmark, src.file.landmarks())
+            foreach (const QPointF &landmark, src.file.points())
                 keyPoints.push_back(KeyPoint(landmark.x(), landmark.y(), size));
         }
         descriptorExtractor->compute(src, keyPoints, dst);
@@ -166,7 +166,7 @@ class SIFTDescriptorTransform : public UntrainableTransform
     void project(const Template &src, Template &dst) const
     {
         std::vector<KeyPoint> keyPoints;
-        foreach (const QPointF &val, src.file.landmarks())
+        foreach (const QPointF &val, src.file.points())
             keyPoints.push_back(KeyPoint(val.x(), val.y(), size));
         Mat m;
@@ -200,7 +200,7 @@ class GridTransform : public UntrainableTransform
             for (float j=column_step/2; j<src.m().cols; j+=column_step)
                 landmarks.append(QPointF(i,j));
         dst = src;
-        dst.file.setLandmarks(landmarks);
+        dst.file.setPoints(landmarks);
     }
 };
@@ -37,18 +37,18 @@ static TemplateList Expanded(const TemplateList &amp;templates)
         }
         const bool fte = t.file.get<bool>("FTE", false);
-        QList<QPointF> landmarks = t.file.landmarks();
-        QList<QRectF> ROIs = t.file.ROIs();
-        if (landmarks.size() % t.size() != 0) qFatal("Uneven landmark count.");
-        if (ROIs.size() % t.size() != 0) qFatal("Uneven ROI count.");
-        const int landmarkStep = landmarks.size() / t.size();
-        const int ROIStep = ROIs.size() / t.size();
+        QList<QPointF> points = t.file.points();
+        QList<QRectF> rects = t.file.rects();
+        if (points.size() % t.size() != 0) qFatal("Uneven point count.");
+        if (rects.size() % t.size() != 0) qFatal("Uneven rect count.");
+        const int pointStep = points.size() / t.size();
+        const int rectStep = rects.size() / t.size();
         for (int i=0; i<t.size(); i++) {
             if (!fte || !t.file.get<bool>("enrollAll", false)) {
                 expanded.append(Template(t.file, t[i]));
-                expanded.last().file.setROIs(ROIs.mid(i*ROIStep, ROIStep));
-                expanded.last().file.setLandmarks(landmarks.mid(i*landmarkStep, landmarkStep));
+                expanded.last().file.setRects(rects.mid(i*rectStep, rectStep));
+                expanded.last().file.setPoints(points.mid(i*pointStep, pointStep));
             }
         }
     }
@@ -155,7 +155,7 @@ class RndPointTransform : public Transform
     void project(const Template &src, Template &dst) const
     {
         dst = src;
-        dst.file.appendLandmark(QPointF(src.m().cols * x, src.m().rows * y));
+        dst.file.appendPoint(QPointF(src.m().cols * x, src.m().rows * y));
     }
 };
@@ -148,7 +148,7 @@ class RectFromLandmarksTransform : public UntrainableTransform
     void project(const Template &src, Template &dst) const
     {
-        if (src.file.landmarks().isEmpty()) {
+        if (src.file.points().isEmpty()) {
             qWarning("No landmarks");
             dst = src;
             return;
@@ -160,12 +160,12 @@ class RectFromLandmarksTransform : public UntrainableTransform
         maxX = maxY = -std::numeric_limits<int>::max();
         foreach(int index, indices) {
-            if (src.file.landmarks().size() > index+1) {
-                if (src.file.landmarks()[index].x() < minX) minX = src.file.landmarks()[index].x();
-                if (src.file.landmarks()[index].x() > maxX) maxX = src.file.landmarks()[index].x();
-                if (src.file.landmarks()[index].y() < minY) minY = src.file.landmarks()[index].y();
-                if (src.file.landmarks()[index].y() > maxY) maxY = src.file.landmarks()[index].y();
-                dst.file.appendLandmark(src.file.landmarks()[index]);
+            if (src.file.points().size() > index+1) {
+                if (src.file.points()[index].x() < minX) minX = src.file.points()[index].x();
+                if (src.file.points()[index].x() > maxX) maxX = src.file.points()[index].x();
+                if (src.file.points()[index].y() < minY) minY = src.file.points()[index].y();
+                if (src.file.points()[index].y() > maxY) maxY = src.file.points()[index].y();
+                dst.file.appendPoint(src.file.points()[index]);
             }
         }
@@ -67,17 +67,14 @@ class AffineTransform : public UntrainableTransform
         else           dstPoints[2] = Point2f(x3*width, y3*height);
         Point2f srcPoints[3];
-        if (src.file.contains("Affine_0_X") &&
-            src.file.contains("Affine_0_Y") &&
-            src.file.contains("Affine_1_X") &&
-            src.file.contains("Affine_1_Y") &&
-            (src.file.contains("Affine_2_X") || twoPoints) &&
-            (src.file.contains("Affine_2_Y") || twoPoints)) {
-            srcPoints[0] = Point2f(src.file.get<float>("Affine_0_X"), src.file.get<float>("Affine_0_Y"));
-            srcPoints[1] = Point2f(src.file.get<float>("Affine_1_X"), src.file.get<float>("Affine_1_Y"));
-            if (!twoPoints) srcPoints[2] = Point2f(src.file.get<float>("Affine_2_X"), src.file.get<float>("Affine_2_Y"));
+        if (src.file.contains("Affine_0") &&
+            src.file.contains("Affine_1") &&
+            (src.file.contains("Affine_2") || twoPoints)) {
+            srcPoints[0] = OpenCVUtils::toPoint(src.file.get<QPointF>("Affine_0"));
+            srcPoints[1] = OpenCVUtils::toPoint(src.file.get<QPointF>("Affine_1"));
+            if (!twoPoints) srcPoints[2] = OpenCVUtils::toPoint(src.file.get<QPointF>("Affine_2"));
         } else {
-            const QList<Point2f> landmarks = OpenCVUtils::toPoints(src.file.landmarks());
+            const QList<Point2f> landmarks = OpenCVUtils::toPoints(src.file.points());
             if ((landmarks.size() < 2) || (!twoPoints && (landmarks.size() < 3))) {
                 resize(src, dst, Size(width, height));
@@ -87,14 +84,9 @@ class AffineTransform : public UntrainableTransform
                 srcPoints[1] = landmarks[1];
                 if (!twoPoints) srcPoints[2] = landmarks[2];
-                dst.file.set("Affine_0_X", landmarks[0].x);
-                dst.file.set("Affine_0_Y", landmarks[0].y);
-                dst.file.set("Affine_1_X", landmarks[1].x);
-                dst.file.set("Affine_1_Y", landmarks[1].y);
-                if (!twoPoints) {
-                    dst.file.set("Affine_2_X", landmarks[2].x);
-                    dst.file.set("Affine_2_Y", landmarks[2].y);
-                }
+                dst.file.set("Affine_0", OpenCVUtils::fromPoint(landmarks[0]));
+                dst.file.set("Affine_1", OpenCVUtils::fromPoint(landmarks[1]));
+                if (!twoPoints) dst.file.set("Affine_2", OpenCVUtils::fromPoint(landmarks[2]));
             }
         }
         if (twoPoints) srcPoints[2] = getThirdAffinePoint(srcPoints[0], srcPoints[1]);
@@ -189,11 +189,11 @@ class GaborJetTransform : public UntrainableTransform
     void project(const Template &src, Template &dst) const
     {
-        const QList<QPointF> landmarks = src.file.landmarks();
-        dst = Mat(landmarks.size(), kReals.size(), CV_32FC1);
-        for (int i=0; i<landmarks.size(); i++)
+        const QList<QPointF> points = src.file.points();
+        dst = Mat(points.size(), kReals.size(), CV_32FC1);
+        for (int i=0; i<points.size(); i++)
             for (int j=0; j<kReals.size(); j++)
-                    dst.m().at<float>(i,j) = response(src, landmarks[i], kReals[j], kImaginaries[j], component);
+                    dst.m().at<float>(i,j) = response(src, points[i], kReals[j], kImaginaries[j], component);
     }
 };