/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
 * Copyright 2012 The MITRE Corporation                                      *
 *                                                                           *
 * Licensed under the Apache License, Version 2.0 (the "License");           *
 * you may not use this file except in compliance with the License.          *
 * You may obtain a copy of the License at                                   *
 *                                                                           *
 *     http://www.apache.org/licenses/LICENSE-2.0                            *
 *                                                                           *
 * Unless required by applicable law or agreed to in writing, software       *
 * distributed under the License is distributed on an "AS IS" BASIS,         *
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.  *
 * See the License for the specific language governing permissions and       *
 * limitations under the License.                                            *
 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */

#include <opencv2/imgproc/imgproc.hpp>
#include <openbr_plugin.h>

#include "core/common.h"
#include "core/opencvutils.h"

using namespace cv;

namespace br
{

/*!
 * \ingroup transforms
 * \brief Selects a random transform.
 * \author Josh Klontz \cite jklontz
 */
class RndTransformTransform : public Transform
{
    Q_OBJECT
    Q_PROPERTY(QList<br::Transform*> transforms READ get_transforms WRITE set_transforms RESET reset_transforms STORED false)
    BR_PROPERTY(QList<br::Transform*>, transforms, QList<br::Transform*>())

    int selectedIndex;
    Transform *selectedTransform;

    void train(const TemplateList &data)
    {
        selectedIndex = theRNG().uniform(0, transforms.size());
        selectedTransform = transforms[selectedIndex]->clone();
        selectedTransform->train(data);
    }

    void project(const Template &src, Template &dst) const
    {
        selectedTransform->project(src, dst);
    }

    void store(QDataStream &stream) const
    {
        stream << selectedIndex << *selectedTransform;
    }

    void load(QDataStream &stream)
    {
        stream >> selectedIndex;
        selectedTransform = transforms[selectedIndex]->clone();
        stream >> *selectedTransform;
    }
};

BR_REGISTER(Transform, RndTransformTransform)

/*!
 * \ingroup transforms
 * \brief Generates a random subspace.
 * \author Josh Klontz \cite jklontz
 */
class RndSubspaceTransform : public Transform
{
    Q_OBJECT
    Q_PROPERTY(float fraction READ get_fraction WRITE set_fraction RESET reset_fraction STORED false)
    Q_PROPERTY(bool weighted READ get_weighted WRITE set_weighted RESET reset_weighted STORED false)
    BR_PROPERTY(float, fraction, 0.5)
    BR_PROPERTY(bool, weighted, false)

    Mat map;

    void train(const TemplateList &data)
    {
        int cols = data.first().m().cols;
        int size = data.first().m().rows * cols;
        QList<float> weights; weights.reserve(size);
        if (weighted) {
            Mat flatData = OpenCVUtils::toMat(data.data());
            for (int i=0; i<size; i++) {
                Scalar mean, stddev;
                cv::meanStdDev(flatData.col(i), mean, stddev);
                weights.append(pow(stddev[0],2.0));
            }
        } else {
            for (int i=0; i<size; i++)
                weights.append(1);
        }
        const int dimsOut = weights.size()*fraction;

        QList<int> sample = Common::RandSample(dimsOut, weights);
        Mat xMap(1, dimsOut, CV_16SC1);
        Mat yMap(1, dimsOut, CV_16SC1);
        for (int j=0; j<dimsOut; j++) {
            int index = sample[j];
            xMap.at<short>(0,j) = index % cols;
            yMap.at<short>(0,j) = index / cols;
        }
        std::vector<Mat> mv;
        mv.push_back(xMap);
        mv.push_back(yMap);

        merge(mv, map);
    }

    void project(const Template &src, Template &dst) const
    {
        remap(src, dst, map, Mat(), INTER_NEAREST);
    }

    void store(QDataStream &stream) const
    {
        stream << fraction << weighted << map;
    }

    void load(QDataStream &stream)
    {
        stream >> fraction >> weighted >> map;
    }
};

BR_REGISTER(Transform, RndSubspaceTransform)

/*!
 * \ingroup transforms
 * \brief Selects a random region.
 * \author Josh Klontz \cite jklontz
 */
class RndRegionTransform : public Transform
{
    Q_OBJECT
    Q_PROPERTY(float x READ get_x WRITE set_x RESET reset_x)
    Q_PROPERTY(float y READ get_y WRITE set_y RESET reset_y)
    Q_PROPERTY(float width READ get_width WRITE set_width RESET reset_width)
    Q_PROPERTY(float height READ get_height WRITE set_height RESET reset_height)
    BR_PROPERTY(float, x, -1)
    BR_PROPERTY(float, y, -1)
    BR_PROPERTY(float, width, -1)
    BR_PROPERTY(float, height, -1)

    void train(const TemplateList &data)
    {
        (void) data;

        RNG &rng = theRNG();
        width = rng.uniform(0.f, 1.f);
        height = rng.uniform(0.f, 1.f);
        x = rng.uniform(0.f, 1.f-width);
        y = rng.uniform(0.f, 1.f-height);
    }

    void project(const Template &src, Template &dst) const
    {
        dst = src.m()(Rect(src.m().cols * x,
                           src.m().rows * y,
                           src.m().cols * width,
                           src.m().rows * height));
    }
};

BR_REGISTER(Transform, RndRegionTransform)

/*!
 * \ingroup transforms
 * \brief Generates a random landmark.
 * \author Josh Klontz \cite jklontz
 */
class RndPointTransform : public Transform
{
    Q_OBJECT
    Q_PROPERTY(float x READ get_x WRITE set_x RESET reset_x)
    Q_PROPERTY(float y READ get_y WRITE set_y RESET reset_y)
    BR_PROPERTY(float, x, -1)
    BR_PROPERTY(float, y, -1)

    void train(const TemplateList &data)
    {
        (void) data;

        RNG &rng = theRNG();
        x = rng.uniform(0.f, 1.f);
        y = rng.uniform(0.f, 1.f);
    }

    void project(const Template &src, Template &dst) const
    {
        dst = src;
        dst.file.appendLandmark(QPointF(src.m().cols * x, src.m().rows * y));
    }
};

BR_REGISTER(Transform, RndPointTransform)

} // namespace br

#include "random.moc"