Application: gvSIG desktop » gvSIG projections

/* gvSIG. Sistema de Informaci?n Geogr?fica de la Generalitat Valenciana

 *

 * Copyright (C) 2006 Instituto de Desarrollo Regional and Generalitat Valenciana.

 *

 * This program 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 2

 * 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 this program; if not, write to the Free Software

 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307,USA.

 *

 * For more information, contact:

 *

 *  Generalitat Valenciana

 *   Conselleria d'Infraestructures i Transport

 *   Av. Blasco Ib??ez, 50

 *   46010 VALENCIA

 *   SPAIN

 *

 *      +34 963862235

 *   gvsig@gva.es

 *      www.gvsig.gva.es

 *

 *    or

 *

 *   Instituto de Desarrollo Regional (Universidad de Castilla La-Mancha)

 *   Campus Universitario s/n

 *   02071 Alabacete

 *   Spain

 *

 *   +34 967 599 200

 */

package org.gvsig.crs.proj;

import org.gdal.osr.CoordinateTransformation;

import org.gdal.osr.SpatialReference;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

import org.gvsig.crs.CrsException;

public class CrsProj {

    private static Logger logger = LoggerFactory.getLogger(CrsProj.class);

    private SpatialReference spatialReference;

    private String _strCrs;

//    private String abrev = null;

//    private Bounds bounds = null;

//    public class Bounds {

//        double minX;

//        double maxX;

//        double minY;

//        double maxY;

//        public Bounds(double minX, double maxX, double minY, double maxY){

//            this.minX = minX;

//            this.maxX = maxX;

//            this.minY = minY;

//            this.maxY = maxY;

//        }

//

//

//        public double getMinX() {

//            return minX;

//        }

//

//

//        public double getMaxX() {

//            return maxX;

//        }

//

//

//        public double getMinY() {

//            return minY;

//        }

//

//

//        public double getMaxY() {

//            return maxY;

//        }

//

//    }

        public CrsProj(String proj4String) {

        try {

            if( proj4String == null ) {

                logger.error("Can't create CRS from 'null' value.");

                return;

            }

            createCrs(proj4String);

        } catch(CrsProjException e) {

            logger.error("Can't create CRS from '"+proj4String+"'.",e);

        }

        // FIXME: borralo

        showInfo();

    }

//FIXME: Esto era un parche para intentar generar el envelope de un sistema de referencia.

//TODO: A eliminar.

//    public CrsProj(String proj4String, String abrev) {

//        //FIXME: Esto es un parche para ver la diferencia entre las cadenas proj4

//        // que nos da gvSIG y las que nos da gdal. A Eliminar

////        if(StringUtils.startsWith(abrev, "EPSG:")){

////            SpatialReference sr = new SpatialReference();

////            sr.ImportFromEPSG(Integer.parseInt(StringUtils.removeStart(abrev, "EPSG:")));

////            proj4String = sr.ExportToProj4();

////        }

//

//         this(proj4String);

////        try {

////            if (proj4String == null) {

////                logger.error("Can't create CRS from 'null' value.");

////                return;

////            }

////            createCrs(proj4String);

////        } catch (CrsProjException e) {

////            logger.error("Can't create CRS from '" + proj4String + "'.", e);

////        }

//        this.abrev = abrev;

//        showInfo();

//    }

//

//    public Bounds getBounds() {

//        if (bounds == null) {

//            String unit = spatialReference.GetAttrValue("UNIT");

//            String unitConv = spatialReference.GetAttrValue("UNIT", 1);

//            if (StringUtils.equalsIgnoreCase(unit, osr.SRS_UA_DEGREE)) {

//                bounds = new Bounds(-180, 180, -90, 90);

//            } else if (StringUtils.equalsIgnoreCase(unit, osr.SRS_UA_RADIAN)) {

//                bounds = new Bounds(-Math.PI, Math.PI, -Math.PI / 2, Math.PI / 2);

//            } else {

//                bounds =

//                    new Bounds(-spatialReference.GetSemiMajor() * 2* Math.PI / Double.parseDouble(unitConv),

//                        +spatialReference.GetSemiMajor() * 2 * Math.PI / Double.parseDouble(unitConv),

//                        -spatialReference.GetSemiMinor() * 2 * Math.PI / Double.parseDouble(unitConv),

//                        +spatialReference.GetSemiMinor() * 2 * Math.PI / Double.parseDouble(unitConv));

//            }

//        }

//

//        return bounds;

//    }

    private void showInfo() {

                logger.info("=========================================================================");

                logger.info("_strCrs" + " = " + _strCrs);

                logger.info("IsProjected" + " = " + spatialReference.IsProjected());

                logger.info("IsGeographic" + " = " + spatialReference.IsGeographic());

                logger.info("IsGeocentric" + " = " + spatialReference.IsGeocentric());

                logger.info("IsVertical" + " = " + spatialReference.IsVertical());

                logger.info("SemiMajor" + " = " + spatialReference.GetSemiMajor());

                logger.info("SemiMinor" + " = " + spatialReference.GetSemiMinor());

                logger.info("InvFlattening" + " = " + spatialReference.GetInvFlattening());

                logger.info("PROJCS" + " = " + spatialReference.GetAttrValue("PROJCS"));

                logger.info("GEOGCS" + " = " + spatialReference.GetAttrValue("GEOGCS"));

                logger.info("DATUM" + " = " + spatialReference.GetAttrValue("DATUM"));

                logger.info("SPHEROID" + " = " + spatialReference.GetAttrValue("SPHEROID"));

                logger.info("PROJECTION" + " = " + spatialReference.GetAttrValue("PROJECTION"));

                logger.info("UNIT" + " = " + spatialReference.GetAttrValue("UNIT"));

                logger.info("conversion a metros" + " = " + spatialReference.GetAttrValue("UNIT",1));

                logger.info("Linear unit name" + " = " + spatialReference.GetLinearUnitsName());

                logger.info("linear units" + " = " + spatialReference.GetLinearUnits());

                logger.info("Angular units" + " = " + spatialReference.GetAngularUnits());

                logger.info("TOWGS84" + " = " + spatialReference.GetAttrValue("TOWGS84"));

                logger.info("TOWGS84" + " = " + spatialReference.GetAttrValue("TOWGS84",1));

                logger.info("TOWGS84" + " = " + spatialReference.GetAttrValue("TOWGS84",2));

                logger.info("SRS_PP_LATITUDE_OF_ORIGIN" + " = " + spatialReference.GetAttrValue("SRS_PP_LATITUDE_OF_ORIGIN"));

                logger.info("=========================================================================");

    }

    protected void createCrs(String strCrs) throws CrsProjException {

        spatialReference = new SpatialReference();

        int errNo = 0;

        try {

            errNo = spatialReference.ImportFromProj4(strCrs);

        } catch (Exception e) {

            throw new CrsProjException("Error creating CRS", e);

        }

        if (errNo!=0) {

            throw new CrsProjException("Error creating CRS. Erro code = "+errNo);

        }

        try {

            errNo = spatialReference.Fixup();

        } catch (Exception e) {

            throw new CrsProjException("Error fixing up CRS", e);

        }

        if (errNo!=0) {

            throw new CrsProjException("Error creating CRS. Erro code = "+errNo);

        }

        _strCrs=strCrs;

    }

    public void reloadCrs() throws CrsException {

        // do nothing

    }

    public void changeStrCrs(String code) {

        _strCrs += code;

    }

    public boolean isLatlong() {

        return spatialReference.IsGeographic()!=0;

    }

    public boolean isGeocent() {

        return spatialReference.IsGeocentric()!=0;

    }

    public String getUnits() {

        return spatialReference.GetAttrValue("UNIT");

    }

    public double getConversionFactorToMeters() {

        return Double.parseDouble(spatialReference.GetAttrValue("UNIT",1));

    }

    public String getStr() {

        return spatialReference.ExportToProj4();

    }

    public long getPtr() {

        return SpatialReference.getCPtr(spatialReference);

    }

    public SpatialReference getSpatialReference(){

        return this.spatialReference;

    }

    public static int operate(double[] firstCoord, double[] secondCoord, double[] thirdCoord, CrsProj srcCrs,

        CrsProj destCrs) throws OperationCrsException {

        try {

            logger.warn("CrsProj.operate: No deber?a pasar por aqu?. (Excepcion provocada para ver de d?nde viene).");

            throw new Exception();

        } catch (Exception e) {

            e.printStackTrace();

        }

        long srcCrsPtr = srcCrs.getPtr();

        long destCrsPtr = destCrs.getPtr();

        if (srcCrsPtr == 0 || destCrsPtr == 0) {

            throw new OperationCrsException(srcCrs, destCrs, "");

        }

        int points =  firstCoord.length;

        if (secondCoord.length > points){

            points = secondCoord.length;

        }

        if (thirdCoord.length > points){

            points = thirdCoord.length;

        }

        CoordinateTransformation ct = null;

        try {

            ct = CoordinateTransformation.CreateCoordinateTransformation(srcCrs.spatialReference, destCrs.spatialReference);

        } catch (Exception e) {

            throw new OperationCrsException(srcCrs, destCrs, "Can't create transformation");

        }

        try {

            if (ct!=null) {

                for (int i = 0; i < points; i++) {

                    double[] resultCoords = ct.TransformPoint(firstCoord[i], secondCoord[i], thirdCoord[i]);

                    firstCoord[i] = resultCoords[0];

                    secondCoord[i] = resultCoords[1];

                    thirdCoord[i] = resultCoords[2];

                }

            }

        } catch (Exception e) {

            throw new OperationCrsException(srcCrs, destCrs, "Can't transform points");

        }

        return 0;

    }

    public static int compareDatums(CrsProj crs1, CrsProj crs2) {

        long crs1Ptr = crs1.getPtr();

        long crs2Ptr = crs2.getPtr();

        if (crs1Ptr == 0 || crs2Ptr == 0) {

            throw new IllegalArgumentException("Coordinate compare error: " + crs1.getStr() + " to " + crs2.getStr());

        }

        return Long.valueOf(crs1Ptr).compareTo(Long.valueOf(crs2Ptr));

    }

}