| /* |
| * 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.transform; |
| |
| import org.opengis.referencing.operation.TransformException; |
| import org.apache.sis.referencing.operation.matrix.MatrixSIS; |
| import org.apache.sis.referencing.operation.matrix.Matrices; |
| import org.apache.sis.referencing.operation.matrix.Matrix4; |
| import org.apache.sis.internal.referencing.ExtendedPrecisionMatrix; |
| import org.apache.sis.internal.util.DoubleDouble; |
| |
| import org.apache.sis.test.DependsOnMethod; |
| import org.apache.sis.test.DependsOn; |
| import org.opengis.test.Assert; |
| import org.junit.Test; |
| |
| import static org.opengis.test.Assert.*; |
| |
| |
| /** |
| * Tests the {@link ScaleTransform} class. |
| * |
| * @author Martin Desruisseaux (Geomatys) |
| * @version 0.7 |
| * @since 0.7 |
| * @module |
| */ |
| @DependsOn(AbstractMathTransformTest.class) |
| public final strictfp class ScaleTransformTest extends MathTransformTestCase { |
| /** |
| * Sets the {@link #transform} field to the {@link ScaleTransform} instance to test. |
| * |
| * @param sourceDimensions expected number of source dimensions. |
| * @param targetDimensions expected number of source dimensions. |
| * @param matrix the data to use for creating the transform. |
| */ |
| private void create(final int sourceDimensions, final int targetDimensions, final MatrixSIS matrix) { |
| final double[] elements = matrix.getElements(); |
| final ScaleTransform tr = new ScaleTransform(matrix.getNumRow(), matrix.getNumCol(), elements); |
| assertEquals("sourceDimensions", sourceDimensions, tr.getSourceDimensions()); |
| assertEquals("targetDimensions", targetDimensions, tr.getTargetDimensions()); |
| Assert.assertMatrixEquals("matrix", matrix, tr.getMatrix(), 0.0); |
| assertArrayEquals("elements", elements, tr.getExtendedElements(), 0.0); |
| transform = tr; |
| validate(); |
| } |
| |
| /** |
| * Tests a transform created from a square matrix with no error terms. |
| * In this test, no dimension are dropped. |
| * |
| * @throws TransformException should never happen. |
| */ |
| @Test |
| public void testConstantDimension() throws TransformException { |
| create(3, 3, new Matrix4( |
| 2, 0, 0, 0, |
| 0, 3, 0, 0, |
| 0, 0, 8, 0, |
| 0, 0, 0, 1)); |
| |
| verifyTransform(new double[] {1,1,1, 6, 0, 2, 2, Double.NaN, 6}, |
| new double[] {2,3,8, 12, 0, 16, 4, Double.NaN, 48}); |
| } |
| |
| /** |
| * Tests a transform with less output dimensions than input dimensions. |
| * This transform drops the last dimension. |
| * |
| * @throws TransformException should never happen. |
| */ |
| @Test |
| @DependsOnMethod("testConstantDimension") |
| public void testDimensionReduction() throws TransformException { |
| isInverseTransformSupported = false; // Because matrix is not square. |
| create(3, 2, Matrices.create(3, 4, new double[] { |
| 2, 0, 0, 0, |
| 0, 3, 0, 0, |
| 0, 0, 0, 1})); |
| |
| verifyTransform(new double[] {1,1,1, 6, 0, 2, 2, Double.NaN, 6}, |
| new double[] {2,3, 12, 0, 4, Double.NaN}); |
| } |
| |
| /** |
| * Tests a transform with more output dimensions than input dimensions. |
| * The extra dimension has values set to 0. This kind of transform happen |
| * in the inverse of <cite>"Geographic 3D to 2D conversion"</cite> (EPSG:9659). |
| * |
| * @throws TransformException should never happen. |
| */ |
| @Test |
| @DependsOnMethod("testDimensionReduction") |
| public void testDimensionAugmentation() throws TransformException { |
| transform = new ProjectiveTransform(Matrices.create(4, 3, new double[] { |
| 2, 0, 0, |
| 0, 3, 0, |
| 0, 0, 0, |
| 0, 0, 1})); |
| |
| assertInstanceOf("inverse", ScaleTransform.class, transform.inverse()); |
| verifyTransform(new double[] {1,1, 6, 0, 2, Double.NaN}, |
| new double[] {2,3,0, 12, 0, 0, 4, Double.NaN, 0}); |
| } |
| |
| /** |
| * Verifies that {@link ScaleTransform} stores the error terms when they exist. |
| */ |
| @Test |
| @DependsOnMethod("testConstantDimension") |
| public void testExtendedPrecision() { |
| final Number O = 0; |
| final Number l = 1; |
| final DoubleDouble r = DoubleDouble.createDegreesToRadians(); |
| final MatrixSIS matrix = Matrices.create(4, 4, new Number[] { |
| r, O, O, O, |
| O, r, O, O, |
| O, O, l, O, |
| O, O, O, l |
| }); |
| final double[] elements = ((ExtendedPrecisionMatrix) matrix).getExtendedElements(); |
| assertTrue (r.value > r.error); |
| assertFalse(r.error == 0); // Paranoiac checks for making sure that next assertion will test something. |
| assertArrayEquals(new double[] { // Paranoiac check for making sure that getExtendedElements() is not broken. |
| r.value, 0, 0, 0, |
| 0, r.value, 0, 0, |
| 0, 0, 1, 0, |
| 0, 0, 0, 1, |
| r.error, 0, 0, 0, |
| 0, r.error, 0, 0, |
| 0, 0, 0, 0, |
| 0, 0, 0, 0}, elements, 0); |
| |
| final ScaleTransform tr = new ScaleTransform(4, 4, elements); |
| assertEquals("sourceDimensions", 3, tr.getSourceDimensions()); |
| assertEquals("targetDimensions", 3, tr.getTargetDimensions()); |
| Assert.assertMatrixEquals("matrix", matrix, tr.getMatrix(), 0.0); |
| assertArrayEquals("elements", elements, tr.getExtendedElements(), 0.0); |
| transform = tr; |
| validate(); |
| } |
| } |