Application: gvSIG desktop » gvSIG projections

package org.gvsig.geotools.proj.catalog.test;

import java.util.ArrayList;

import java.util.List;

import org.gvsig.geotools.proj.catalog.DefaultTransformationDefinition;

import org.gvsig.geotools.proj.catalog.GtCRSCatalogManager;

import org.gvsig.proj.catalog.CRSCatalogManager;

import org.gvsig.proj.catalog.CRSCatalogManagerIT;

import org.gvsig.proj.catalog.CRSDefinition;

import org.gvsig.proj.catalog.TransformationDefinition;

public class GtCRSCatalogManagerIT extends CRSCatalogManagerIT {

        private static final String EPSG25830_WKT = "PROJCS[\"ETRS89 / UTM zone 30N\",GEOGCS[\"ETRS89\",DATUM[\"European Terrestrial Reference System 1989\",SPHEROID[\"GRS 1980\",6378137.0, 298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0],AUTHORITY[\"EPSG\",\"6258\"]],PRIMEM[\"Greenwich\", 0.0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.017453292519943295],AXIS[\"Geodetic latitude\",NORTH],AXIS[\"Geodetic longitude\",EAST],AUTHORITY[\"EPSG\",\"4258\"]],PROJECTION[\"Transverse_Mercator\",AUTHORITY[\"EPSG\",\"9807\"]],PARAMETER[\"central_meridian\", -3.0],PARAMETER[\"latitude_of_origin\",0.0],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000.0],PARAMETER[\"false_northing\",0.0],UNIT[\"m\",1.0],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"25830\"]]";

        private static final String WKT_25830_TO_4326 = "CONCAT_MT[INVERSE_MT[PARAM_MT[\"Transverse_Mercator\", PARAMETER[\"semi_major\", 6378137.0], PARAMETER[\"semi_minor\", 6356752.314140356], PARAMETER[\"central_meridian\", -3.0], PARAMETER[\"latitude_of_origin\", 0.0], PARAMETER[\"scale_factor\", 0.9996], PARAMETER[\"false_easting\", 500000.0], PARAMETER[\"false_northing\", 0.0]]], PARAM_MT[\"Ellipsoid_To_Geocentric\", PARAMETER[\"dim\", 2], PARAMETER[\"semi_major\", 6378137.0], PARAMETER[\"semi_minor\", 6356752.314140356]], PARAM_MT[\"Geocentric_To_Ellipsoid\", PARAMETER[\"dim\", 2], PARAMETER[\"semi_major\", 6378137.0], PARAMETER[\"semi_minor\", 6356752.314245179]], PARAM_MT[\"Affine\", PARAMETER[\"num_row\", 3], PARAMETER[\"num_col\", 3], PARAMETER[\"elt_0_0\", 0.0], PARAMETER[\"elt_0_1\", 1.0], PARAMETER[\"elt_1_0\", 1.0], PARAMETER[\"elt_1_1\", 0.0]]]";

        private static final String TRANS_SOURCE_25830 = "EPSG:25830";

        private static final String TRANS_TARGET_4326 = "EPSG:4326";

        @Override

        protected CRSCatalogManager createProjectionManager() {

                return new GtCRSCatalogManager();

        }

        @Override

        protected String getAuthorityName() {

                return "EPSG";

        }

        @Override

        protected String getCode() {

                return "EPSG:23030";

        }

        @Override

        protected String getSearchString() {

                return "WGS";

        }

        @Override

        protected String getCoordinateReferenceSystemWKT() {

                return EPSG25830_WKT;

        }

        @Override

        protected String getCoordinateTransformationWKT() {

                return WKT_25830_TO_4326;

        }

        @Override

        protected String getCoordinateTransformationSource() {

                return TRANS_SOURCE_25830;

        }

        @Override

        protected String getCoordinateTransformationTarget() {

                return TRANS_TARGET_4326;

        }

    /**

     * Test method

     * {@link org.gvsig.proj.catalog.CRSCatalogManager#getTransformationDefinition(String, CRSDefinition, CRSDefinition)}

     * for a well-known ntv2 transformation

     */

    public void testNTv2Grid() throws Exception {

            String transCode = "EPSG:15932";

        CRSDefinition sourceCrs =

            manager.getCRSDefinition("EPSG:4230");

        CRSDefinition targetCrs =

                manager.getCRSDefinition("EPSG:4258");

        TransformationDefinition transDef = manager.getTransformationDefinition(transCode);

        assertNotNull(transDef);

        assertNotNull(transDef.getName());

        assertEquals(transCode, transDef.getIdentifier());

        assertEquals(sourceCrs, transDef.getSourceDefinition());

        assertEquals(targetCrs, transDef.getTargetDefinition());

        assertEquals(getAuthorityName(), transDef.getAuthorityName());

    }

    /**

     * This test should be moved to TransformationDefinitionTests class!!!

     * 

     * Tests the instantiation of TransformationDefinition and ensures some basic information is available

     * for a set of coordinate operations using heterogeneous methods

     */

    public void testTransformationDefinition() throws Exception {

            ArrayList<String> codesToFind = new ArrayList<String>();

            codesToFind.add("EPSG:1632"); // 7-params

            codesToFind.add("EPSG:15932"); // NTv2 grid

        CRSDefinition sourceCrs =

            manager.getCRSDefinition("EPSG:23030");

        CRSDefinition targetCrs =

                manager.getCRSDefinition("EPSG:4258");

        testTransformationDefinition(sourceCrs, targetCrs, codesToFind);

            codesToFind = new ArrayList<String>();

            codesToFind.add("EPSG:1633"); // 7-params

            codesToFind.add("EPSG:15933"); // NTv2 grid

        sourceCrs =

            manager.getCRSDefinition("EPSG:23030");

        targetCrs =

                manager.getCRSDefinition("EPSG:4326");

        testTransformationDefinition(sourceCrs, targetCrs, codesToFind);

    }

    protected void testTransformationDefinition(CRSDefinition sourceCrs, CRSDefinition targetCrs, ArrayList<String> codesToFind) throws Exception {

        List<TransformationDefinition> transDefs = manager.getCoordinateTransformations(sourceCrs, targetCrs);

        for (TransformationDefinition transDef: transDefs) {

            assertNotNull(transDef);

            assertNotNull(transDef.getName());

            System.out.println(transDef.getName());

            if (transDef.getIdentifier()!=null && codesToFind.contains(transDef.getIdentifier())) {

                    codesToFind.remove(transDef.getIdentifier());

            }

            else {

                    for (String code: new ArrayList<String>(codesToFind)) {

                            if (transDef.getName().contains(code)) {

                                    codesToFind.remove(code);

                            }

                    }

            }

            assertEquals(sourceCrs, transDef.getSourceDefinition());

            assertEquals(targetCrs, transDef.getTargetDefinition());

        }

        assertEquals(codesToFind.size(), 0);

    }

    /**

     * Test method

     * {@link org.gvsig.proj.catalog.CRSCatalogManager#getTransformationDefinition(String, CRSDefinition, CRSDefinition)}

     * for a well-known transformation

     */

    public void testTransformation() throws Exception {

            String transCode = "EPSG:1632";

        CRSDefinition sourceCrs =

            manager.getCRSDefinition("EPSG:4230");

        CRSDefinition targetCrs =

                manager.getCRSDefinition("EPSG:4258");

        TransformationDefinition transDef = manager.getTransformationDefinition(transCode);

        assertNotNull(transDef);

        assertEquals(transCode, transDef.getIdentifier());

        assertEquals(sourceCrs, transDef.getSourceDefinition());

        assertEquals(targetCrs, transDef.getTargetDefinition());

        assertEquals(getAuthorityName(), transDef.getAuthorityName());

    }        

}