log_rmse.hpp

Go to the documentation of this file.

/*

PICCANTE
The hottest HDR imaging library!
http://vcg.isti.cnr.it/piccante

Copyright (C) 2014
Visual Computing Laboratory - ISTI CNR
http://vcg.isti.cnr.it
First author: Francesco Banterle

This Source Code Form is subject to the terms of the Mozilla Public
License, v. 2.0. If a copy of the MPL was not distributed with this
file, You can obtain one at http://mozilla.org/MPL/2.0/.

*/

#ifndef PIC_METRICS_LOG_RMSE_HPP
#define PIC_METRICS_LOG_RMSE_HPP

#include <math.h>

#include "../base.hpp"
#include "../image.hpp"
#include "../util/math.hpp"
#include "../metrics/base.hpp"

namespace pic {

PIC_INLINE double logRMSE(Image *ori, Image *cmp)
{
    if(ori == NULL || cmp == NULL) {
        return -2.0;
    }

    if(!ori->isValid() || !cmp->isValid()) {
        return -4.0;
    }

    if(!ori->isSimilarType(cmp)) {
        return -1.0;
    }

    int size = ori->size();
    int counter = 0;

    double acc = 0.0;
    for(int i = 0; i < size; i++) {
        if(ori->data[i] > 0.0f && cmp->data[i] > 0.0f) {
            double val = log2(ori->data[i] / cmp->data[i]);
            acc += val * val;
            counter++;
        }
    }

    if(counter > 0) {
        acc = acc / double(counter);
        return sqrt(acc);
    } else {
        return -3.0;
    }
}

} // end namespace pic

#endif /* PIC_METRICS_LOG_RMSE_HPP */