/*
 * Licensed to the Apache Software Foundation (ASF) under one or more
 * contributor license agreements.  See the NOTICE file distributed with
 * this work for additional information regarding copyright ownership.
 * The ASF licenses this file to You 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.
 */
package org.apache.sis.referencing.operation.builder;

import java.util.Map;
import java.util.HashMap;
import java.util.Random;
import java.util.Collections;
import java.awt.geom.AffineTransform;
import org.opengis.util.FactoryException;
import org.opengis.geometry.DirectPosition;
import org.opengis.referencing.operation.Matrix;
import org.apache.sis.referencing.operation.matrix.Matrix3;
import org.apache.sis.geometry.DirectPosition1D;
import org.apache.sis.geometry.DirectPosition2D;
import org.apache.sis.test.DependsOnMethod;
import org.apache.sis.test.TestUtilities;
import org.apache.sis.test.TestCase;
import org.junit.Test;

import static org.apache.sis.test.Assert.*;


/**
 * Tests {@link LinearTransformBuilder}.
 *
 * @author  Martin Desruisseaux (Geomatys)
 * @version 1.0
 * @since   0.5
 * @module
 */
public final strictfp class LinearTransformBuilderTest extends TestCase {
    /**
     * Tests a very simple case where an exact answer is expected.
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    public void testMinimalist1D() throws FactoryException {
        final LinearTransformBuilder builder = new LinearTransformBuilder();
        final Map<DirectPosition1D,DirectPosition1D> pos = new HashMap<>(4);
        assertNull(pos.put(new DirectPosition1D(1), new DirectPosition1D(1)));
        assertNull(pos.put(new DirectPosition1D(2), new DirectPosition1D(3)));
        builder.setControlPoints(pos);

        assertArrayEquals(new double[] {1}, builder.getControlPoint(new int[] {1}), STRICT);
        assertArrayEquals(new double[] {3}, builder.getControlPoint(new int[] {2}), STRICT);
        assertNull(                         builder.getControlPoint(new int[] {3}));

        final Matrix m = builder.create(null).getMatrix();
        assertEquals("m₀₀",  2, m.getElement(0, 0), STRICT);
        assertEquals("m₀₁", -1, m.getElement(0, 1), STRICT);
        assertArrayEquals("correlation", new double[] {1}, builder.correlation(), STRICT);
    }

    /**
     * Tests a very simple case where an exact answer is expected.
     * Tolerance threshold is set to zero because the math transform has been built from exactly 3 points,
     * in which case we expect an exact solution without rounding errors at the scale of the {@code double}
     * type. This is possible because SIS implementation uses double-double arithmetic.
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    public void testMinimalist2D() throws FactoryException {
        final Map<DirectPosition2D,DirectPosition2D> pos = new HashMap<>(8);
        assertNull(pos.put(new DirectPosition2D(1, 1), new DirectPosition2D(3, 2)));
        assertNull(pos.put(new DirectPosition2D(1, 2), new DirectPosition2D(3, 5)));
        assertNull(pos.put(new DirectPosition2D(2, 2), new DirectPosition2D(5, 5)));
        final LinearTransformBuilder builder = new LinearTransformBuilder();
        builder.setControlPoints(pos);

        assertArrayEquals(new double[] {3, 2}, builder.getControlPoint(new int[] {1, 1}), STRICT);
        assertArrayEquals(new double[] {3, 5}, builder.getControlPoint(new int[] {1, 2}), STRICT);
        assertArrayEquals(new double[] {5, 5}, builder.getControlPoint(new int[] {2, 2}), STRICT);
        assertNull(                            builder.getControlPoint(new int[] {2, 1}));

        final Matrix m = builder.create(null).getMatrix();

        // First row (x)
        assertEquals("m₀₀",  2, m.getElement(0, 0), STRICT);
        assertEquals("m₀₁",  0, m.getElement(0, 1), STRICT);
        assertEquals("m₀₂",  1, m.getElement(0, 2), STRICT);

        // Second row (y)
        assertEquals("m₁₀",  0, m.getElement(1, 0), STRICT);
        assertEquals("m₁₁",  3, m.getElement(1, 1), STRICT);
        assertEquals("m₁₂", -1, m.getElement(1, 2), STRICT);

        assertArrayEquals("correlation", new double[] {1, 1}, builder.correlation(), STRICT);
    }

    /**
     * Tests a two-dimensional case where sources coordinates are explicitly given.
     *
     * @throws FactoryException if the transform can not be created.
     *
     * @since 0.8
     */
    @Test
    public void testExplicitSource2D() throws FactoryException {
        testSetAllPoints(new LinearTransformBuilder());
        testSetEachPoint(new LinearTransformBuilder());
    }

    /**
     * Same test than {@link #testExplicitSource2D()}, but using the
     * {@link LinearTransformBuilder#LinearTransformBuilder(int...)} constructor.
     *
     * @throws FactoryException if the transform can not be created.
     *
     * @since 0.8
     */
    @Test
    @DependsOnMethod("testExplicitSource2D")
    public void testImplicitSource2D() throws FactoryException {
        testSetAllPoints(new LinearTransformBuilder(2, 3));
        testSetEachPoint(new LinearTransformBuilder(2, 3));
    }

    /**
     * Execution of {@link #testExplicitSource2D()} and {@link #testImplicitSource2D()}
     * where all control points are specified by a map.
     */
    private void testSetAllPoints(final LinearTransformBuilder builder) throws FactoryException {
        final Map<DirectPosition2D,DirectPosition2D> pos = new HashMap<>(8);
        assertNull(pos.put(new DirectPosition2D(0, 0), new DirectPosition2D(3, 9)));
        assertNull(pos.put(new DirectPosition2D(0, 1), new DirectPosition2D(4, 7)));
        assertNull(pos.put(new DirectPosition2D(0, 2), new DirectPosition2D(6, 6)));
        assertNull(pos.put(new DirectPosition2D(1, 0), new DirectPosition2D(4, 8)));
        assertNull(pos.put(new DirectPosition2D(1, 1), new DirectPosition2D(5, 4)));
        assertNull(pos.put(new DirectPosition2D(1, 2), new DirectPosition2D(8, 2)));
        builder.setControlPoints(pos);
        verify(builder);
    }

    /**
     * Execution of {@link #testExplicitSource2D()} and {@link #testImplicitSource2D()}
     * where all control points are specified one-by-one.
     */
    private void testSetEachPoint(final LinearTransformBuilder builder) throws FactoryException {
        builder.setControlPoint(new int[] {0, 0}, new double[] {3, 9});
        builder.setControlPoint(new int[] {0, 1}, new double[] {4, 7});
        builder.setControlPoint(new int[] {0, 2}, new double[] {6, 6});
        builder.setControlPoint(new int[] {1, 0}, new double[] {4, 8});
        builder.setControlPoint(new int[] {1, 1}, new double[] {5, 4});
        builder.setControlPoint(new int[] {1, 2}, new double[] {8, 2});
        verify(builder);
    }

    /**
     * Verifies the transform created by {@link #testExplicitSource2D()} and {@link #testImplicitSource2D()}.
     */
    private void verify(final LinearTransformBuilder builder) throws FactoryException {
        final Matrix m = builder.create(null).getMatrix();

        assertArrayEquals(new double[] {3, 9}, builder.getControlPoint(new int[] {0, 0}), STRICT);
        assertArrayEquals(new double[] {4, 7}, builder.getControlPoint(new int[] {0, 1}), STRICT);
        assertArrayEquals(new double[] {6, 6}, builder.getControlPoint(new int[] {0, 2}), STRICT);
        assertArrayEquals(new double[] {4, 8}, builder.getControlPoint(new int[] {1, 0}), STRICT);
        assertArrayEquals(new double[] {5, 4}, builder.getControlPoint(new int[] {1, 1}), STRICT);
        assertArrayEquals(new double[] {8, 2}, builder.getControlPoint(new int[] {1, 2}), STRICT);
        /*
         * Expect STRICT results because Apache SIS uses double-double arithmetic.
         */
        assertEquals("m₀₀",  16 / 12d, m.getElement(0, 0), STRICT);        // First row (x)
        assertEquals("m₀₁",  21 / 12d, m.getElement(0, 1), STRICT);
        assertEquals("m₀₂",  31 / 12d, m.getElement(0, 2), STRICT);
        assertEquals("m₁₀", -32 / 12d, m.getElement(1, 0), STRICT);        // Second row (y)
        assertEquals("m₁₁", -27 / 12d, m.getElement(1, 1), STRICT);
        assertEquals("m₁₂", 115 / 12d, m.getElement(1, 2), STRICT);

        assertArrayEquals("correlation", new double[] {0.9656, 0.9536}, builder.correlation(), 0.0001);
    }

    /**
     * Tests with a random number of points with an exact solution expected.
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    @DependsOnMethod("testMinimalist1D")
    public void testExact1D() throws FactoryException {
        final Random rd = TestUtilities.createRandomNumberGenerator(-6080923837183751016L);
        for (int i=0; i<10; i++) {
            test1D(rd, rd.nextInt(900) + 100, false, 1E-14, 1E-12);
        }
    }

    /**
     * Tests with a random number of points with an exact solution expected.
     *
     * <p><b>Note:</b> this test can pass with a random seed most of the time. But we fix the seed anyway
     * because there is always a small probability that truly random points are all colinear, in which case
     * the test would fail. Even if the probability is low, we do not take the risk of random build failures.</p>
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    @DependsOnMethod("testMinimalist2D")
    public void testExact2D() throws FactoryException {
        final Random rd = TestUtilities.createRandomNumberGenerator(41632405806929L);
        for (int i=0; i<10; i++) {
            test2D(rd, rd.nextInt(900) + 100, false, 1E-14, 1E-12);
        }
    }

    /**
     * Tests with a random number of points and a random errors in target points.
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    @DependsOnMethod("testExact1D")
    public void testNonExact1D() throws FactoryException {
        final Random rd = TestUtilities.createRandomNumberGenerator(8819436190826166876L);
        for (int i=0; i<4; i++) {
            test1D(rd, rd.nextInt(900) + 100, true, 0.02, 0.5);
        }
    }

    /**
     * Tests with a random number of points and a random errors in target points.
     *
     * <p><b>Note:</b> this test can pass with a random seed most of the time. But we fix the seed anyway
     * because there is always a small probability that truly random points are all colinear, or that lot
     * of errors are in the same directions (thus introducing a larger bias than expected), in which case
     * the test would fail. Even if the probability is low, we do not take the risk of such random build
     * failures.</p>
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    @DependsOnMethod("testExact2D")
    public void testNonExact2D() throws FactoryException {
        final Random rd = TestUtilities.createRandomNumberGenerator(270575025643864L);
        for (int i=0; i<4; i++) {
            test2D(rd, rd.nextInt(900) + 100, true, 0.02, 0.5);
        }
    }

    /**
     * Implementation of {@link #testExact1D()} and {@link #testNonExact1D()}.
     *
     * @param  rd              the random number generator to use.
     * @param  numPts          the number of points to generate.
     * @param  addErrors       {@code true} for adding a random error in the target points.
     * @param  scaleTolerance  tolerance threshold for floating point comparisons.
     */
    private static void test1D(final Random rd, final int numPts, final boolean addErrors,
            final double scaleTolerance, final double translationTolerance) throws FactoryException
    {
        final double scale  = rd.nextDouble() * 30 - 12;
        final double offset = rd.nextDouble() * 10 - 4;
        final Map<DirectPosition1D,DirectPosition1D> pos = new HashMap<>(numPts);
        for (int i=0; i<numPts; i++) {
            final DirectPosition1D src = new DirectPosition1D(rd.nextDouble() * 100 - 50);
            final DirectPosition1D tgt = new DirectPosition1D(src.coordinate * scale + offset);
            if (addErrors) {
                tgt.coordinate += rd.nextDouble() * 10 - 5;
            }
            assertNull(pos.put(src, tgt));
        }
        /*
         * Create the fitted transform to test.
         */
        final LinearTransformBuilder builder = new LinearTransformBuilder();
        builder.setControlPoints(pos);
        final Matrix m = builder.create(null).getMatrix();
        assertEquals("m₀₀", scale,  m.getElement(0, 0), scaleTolerance);
        assertEquals("m₀₁", offset, m.getElement(0, 1), translationTolerance);
        assertEquals("correlation", 1, StrictMath.abs(builder.correlation()[0]), scaleTolerance);
    }

    /**
     * Implementation of {@link #testExact2D()} and {@link #testNonExact2D()}.
     *
     * @param  rd              the random number generator to use.
     * @param  numPts          the number of points to generate.
     * @param  addErrors       {@code true} for adding a random error in the target points.
     * @param  scaleTolerance  tolerance threshold for floating point comparisons.
     */
    private static void test2D(final Random rd, final int numPts, final boolean addErrors,
            final double scaleTolerance, final double translationTolerance) throws FactoryException
    {
        /*
         * Create an AffineTransform to use as the reference implementation.
         */
        final AffineTransform ref = AffineTransform.getRotateInstance(
                rd.nextDouble() * (2 * StrictMath.PI),  // Rotation angle
                rd.nextDouble() * 30 - 12,              // Center X
                rd.nextDouble() * 10 - 8);              // Center Y
        final Map<DirectPosition2D,DirectPosition2D> pos = new HashMap<>(numPts);
        for (int i=0; i<numPts; i++) {
            final DirectPosition2D src = new DirectPosition2D(rd.nextDouble() * 100 - 50, rd.nextDouble() * 200 - 75);
            final DirectPosition2D tgt = new DirectPosition2D();
            assertSame(tgt, ref.transform(src, tgt));
            if (addErrors) {
                tgt.x += rd.nextDouble() * 10 - 5;
                tgt.y += rd.nextDouble() * 10 - 5;
            }
            assertNull(pos.put(src, tgt));
        }
        /*
         * Create the fitted transform to test.
         */
        final LinearTransformBuilder builder = new LinearTransformBuilder();
        builder.setControlPoints(pos);
        final Matrix m = builder.create(null).getMatrix();
        /*
         * Compare the coefficients with the reference implementation.
         */
        assertEquals("m₀₀", ref.getScaleX(),     m.getElement(0, 0), scaleTolerance);
        assertEquals("m₀₁", ref.getShearX(),     m.getElement(0, 1), scaleTolerance);
        assertEquals("m₀₂", ref.getTranslateX(), m.getElement(0, 2), translationTolerance);
        assertEquals("m₁₀", ref.getShearY(),     m.getElement(1, 0), scaleTolerance);
        assertEquals("m₁₁", ref.getScaleY(),     m.getElement(1, 1), scaleTolerance);
        assertEquals("m₁₂", ref.getTranslateY(), m.getElement(1, 2), translationTolerance);
        assertArrayEquals("correlation", new double[] {1, 1}, builder.correlation(), scaleTolerance);
    }

    /**
     * Tests {@link LinearTransformBuilder#getControlPoints()} with gridded source points.
     */
    @Test
    public void testGetControlPoints() {
        testGetControlPoints(new LinearTransformBuilder(3, 4));
    }

    /**
     * Tests {@link LinearTransformBuilder#getControlPoints()} with non-gridded source points.
     */
    @Test
    public void testGetUngriddedControlPoints() {
        testGetControlPoints(new LinearTransformBuilder());
    }

    /**
     * Tests {@link LinearTransformBuilder#getControlPoints()} with the given builder.
     * If the builder is backed by a grid, then the grid size shall be at least 3×4.
     */
    private static void testGetControlPoints(final LinearTransformBuilder builder) {
        final DirectPosition2D s12, s23, s00;
        final DirectPosition2D t12, t23, t00;
        s12 = new DirectPosition2D(1, 2);   t12 = new DirectPosition2D(3, 2);
        s23 = new DirectPosition2D(2, 3);   t23 = new DirectPosition2D(4, 1);
        s00 = new DirectPosition2D(0, 0);   t00 = new DirectPosition2D(7, 3);

        final Map<DirectPosition2D,DirectPosition2D> expected = new HashMap<>();
        final Map<DirectPosition,DirectPosition> actual = builder.getControlPoints();
        assertEquals(0, actual.size());
        assertTrue(actual.isEmpty());
        assertFalse(actual.containsKey  (s12));
        assertFalse(actual.containsKey  (s23));
        assertFalse(actual.containsKey  (s00));
        assertFalse(actual.containsValue(t12));
        assertFalse(actual.containsValue(t23));
        assertFalse(actual.containsValue(t00));
        assertMapEquals(expected, actual);

        builder.setControlPoint(new int[] {1, 2}, t12.getCoordinate());
        assertNull(expected.put(s12, t12));
        assertEquals(1, actual.size());
        assertFalse(actual.isEmpty());
        assertTrue (actual.containsKey  (s12));
        assertFalse(actual.containsKey  (s23));
        assertFalse(actual.containsKey  (s00));
        assertTrue (actual.containsValue(t12));
        assertFalse(actual.containsValue(t23));
        assertFalse(actual.containsValue(t00));
        assertMapEquals(expected, actual);

        builder.setControlPoint(new int[] {2, 3}, t23.getCoordinate());
        assertNull(expected.put(s23, t23));
        assertEquals(2, actual.size());
        assertFalse(actual.isEmpty());
        assertTrue (actual.containsKey  (s12));
        assertTrue (actual.containsKey  (s23));
        assertFalse(actual.containsKey  (s00));
        assertTrue (actual.containsValue(t12));
        assertTrue (actual.containsValue(t23));
        assertFalse(actual.containsValue(t00));
        assertMapEquals(expected, actual);

        builder.setControlPoint(new int[] {0, 0}, t00.getCoordinate());
        assertNull(expected.put(s00, t00));
        assertEquals(3, actual.size());
        assertFalse(actual.isEmpty());
        assertTrue (actual.containsKey  (s12));
        assertTrue (actual.containsKey  (s23));
        assertTrue (actual.containsKey  (s00));
        assertTrue (actual.containsValue(t12));
        assertTrue (actual.containsValue(t23));
        assertTrue (actual.containsValue(t00));
        assertMapEquals(expected, actual);
    }

    /**
     * Tests the effect of {@link LinearTransformBuilder#addLinearizers(Map, int...)}.
     *
     * @throws FactoryException if the transform can not be created.
     */
    @Test
    public void testLinearizers() throws FactoryException {
        final int width  = 3;
        final int height = 4;
        final LinearTransformBuilder builder = new LinearTransformBuilder(width, height);
        for (int y=0; y<height; y++) {
            final int[]    source = new int[2];
            final double[] target = new double[2];
            for (int x=0; x<width; x++) {
                source[0] = x;
                source[1] = y;
                target[0] = StrictMath.cbrt(3 + x*2);
                target[1] = StrictMath.sqrt(1 + y);
                builder.setControlPoint(source, target);
            }
        }
        final NonLinearTransform tr = new NonLinearTransform();
        builder.addLinearizers(Collections.singletonMap("x² y³", tr));
        builder.addLinearizers(Collections.singletonMap("x³ y²", tr), 1, 0);
        final Matrix m = builder.create(null).getMatrix();
        assertEquals("linearizer", "x³ y²", builder.linearizerID());
        assertNotSame("linearizer", tr, builder.linearizer().get());    // Not same because axes should have been swapped.
        assertMatrixEquals("linear",
                new Matrix3(2, 0, 3,
                            0, 1, 1,
                            0, 0, 1), m, 1E-15);
    }
}
