/* -*- mia-c++  -*-
 *
 * This file is part of MIA - a toolbox for medical image analysis 
 * Copyright (c) Leipzig, Madrid 1999-2015 Gert Wollny
 *
 * MIA is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with MIA; if not, see <http://www.gnu.org/licenses/>.
 *
 */

#include <stdexcept>
#include <climits>


#include <mia/core/minimizer.hh>
#include <mia/2d/correlation_weight.hh>
#include <mia/2d/ground_truth_evaluator.hh>

#include <mia/internal/autotest.hh>
using namespace mia; 
using namespace std; 

struct FixPluginpath {
	FixPluginpath(); 
};


struct PseudoGroundTruthFixture {
	PseudoGroundTruthFixture();

	static FixPluginpath fix; 
	
	C2DBounds size;
	size_t N;
	vector<P2DImage> input_series;
};

BOOST_FIXTURE_TEST_CASE( test_correlation, PseudoGroundTruthFixture )
{
	CCorrelationEvaluator ce(0.3);
	auto result = ce(input_series);

	BOOST_CHECK_EQUAL(result.horizontal.get_size().x, 2u);
	BOOST_CHECK_EQUAL(result.horizontal.get_size().y, 3u);
	BOOST_CHECK_EQUAL(result.vertical.get_size().x, 3u);
	BOOST_CHECK_EQUAL(result.vertical.get_size().y, 2u);


	BOOST_CHECK_CLOSE(result.horizontal(0,1), 1.0, 0.1);
	BOOST_CHECK_CLOSE(result.horizontal(0,2)+1, 1.0, 0.1);

	// negative
	BOOST_CHECK_CLOSE(result.horizontal(1,1), 0.94388, 0.1);

	// negative
	BOOST_CHECK_CLOSE(result.vertical(1,0) +1 ,  1,0.1);

	// positive but below threshold
	BOOST_CHECK_CLOSE(result.vertical(1,1)+1,   1.0,0.1);

	BOOST_CHECK_CLOSE(result.vertical(0,1),   0.88659,0.1);

};

BOOST_FIXTURE_TEST_CASE( test_pseudo_ground_truth_no_extra_gradients, PseudoGroundTruthFixture )
{
	const double alpha = 0.0;
	const double beta = 0.0;
	const double rho_th = 0.3;

	C2DGroundTruthEvaluator gte(alpha, beta, rho_th);

	vector<P2DImage> pgt;
	for (size_t i = 0; i < N; ++i)
		pgt.push_back(P2DImage(new C2DFImage(size)));

	gte(input_series, pgt);


	BOOST_CHECK_EQUAL(pgt.size(), N);
	BOOST_REQUIRE(pgt.size() == N);

	for (size_t i = 0; i < N; ++i) {
		BOOST_CHECK_EQUAL(pgt[i]->get_size(), input_series[i]->get_size());
		// what should the output look like ??
		const C2DFImage& res = dynamic_cast<const C2DFImage&>(*pgt[i]);
		const C2DFImage& inp = dynamic_cast<const C2DFImage&>(*input_series[i]);

		BOOST_CHECK_EQUAL(res.get_size(), inp.get_size());
		BOOST_REQUIRE(res.size() == inp.size());

		for(auto r = res.begin(), i = inp.begin(); r != res.end(); ++r, ++i)
			BOOST_CHECK_CLOSE(*r, *i, 0.1);
	}
}


BOOST_FIXTURE_TEST_CASE( test_run_pseudo_ground_truth_evaluation, PseudoGroundTruthFixture )
{
	const double alpha = 1.0;
	const double beta = 1.0;
	const double rho_th = 0.3;

	C2DGroundTruthEvaluator gte(alpha, beta, rho_th);

	vector<P2DImage> pgt;

	gte(input_series, pgt);


	BOOST_CHECK_EQUAL(pgt.size(), N);
	BOOST_REQUIRE(pgt.size() == N);

	for (size_t i = 0; i < N; ++i) {
		BOOST_CHECK_EQUAL(pgt[i]->get_size(), input_series[i]->get_size());
		// what should the output look like ??
		const C2DFImage& res = dynamic_cast<const C2DFImage&>(*pgt[i]);
		const C2DFImage& inp = dynamic_cast<const C2DFImage&>(*input_series[i]);

		BOOST_CHECK_EQUAL(res.get_size(), inp.get_size());
		BOOST_REQUIRE(res.size() == inp.size());

		cvmsg() << "GTimg(" << i << ")=[";
		for(auto r = res.begin(), i = inp.begin(); r != res.end(); ++r, ++i)
			cverb << *r << " ";
		cverb << "]\n";

	}
}

PseudoGroundTruthFixture::PseudoGroundTruthFixture():
	size(3,3),
	N(4)
{
	const float init[4][3*3] = {
		{ 0, 4, 0,  1, 1, 2,  3, 1, 1.1},
		{ 0, 3, 0,  2, 2, 3,  3, 1, 1.2},
		{ 0, 2, 0,  3, 3, 4,  4, 2, 1.3},
		{ 0, 1, 0,  4, 4, 4,  7, 1, 1.4}
	};
	for (size_t i = 0; i < N; ++i)
		input_series.push_back(P2DImage(new C2DFImage(size, init[i])));

}
namespace bfs=boost::filesystem; 

FixPluginpath::FixPluginpath()
{
}

FixPluginpath PseudoGroundTruthFixture::fix; 
