Application: gvSIG desktop » gvSIG projections

import org.gdal.osr.CoordinateTransformation;

import org.gdal.osr.SpatialReference;

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));

    }