Application: gvSIG desktop

/**

 * gvSIG. Desktop Geographic Information System.

 *

 * Copyright (C) 2007-2013 gvSIG Association.

 *

 * 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 3

 * 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., 51 Franklin Street, Fifth Floor, Boston,

 * MA  02110-1301, USA.

 *

 * For any additional information, do not hesitate to contact us

 * at info AT gvsig.com, or visit our website www.gvsig.com.

 */

package org.gvsig.fmap.geom.jts.primitive;

import org.cresques.cts.CoordTransRuntimeException;

import org.cresques.cts.ICoordTrans;

import org.slf4j.Logger;

import org.slf4j.LoggerFactory;

import org.gvsig.fmap.geom.jts.primitive.point.Point2D;

import org.gvsig.fmap.geom.primitive.Envelope;

import org.gvsig.fmap.geom.primitive.EnvelopeNotInitializedException;

import org.gvsig.fmap.geom.primitive.Point;

import org.gvsig.tools.ToolsLocator;

import org.gvsig.tools.dynobject.DynStruct;

import org.gvsig.tools.lang.Cloneable;

import org.gvsig.tools.persistence.PersistenceManager;

/**

 * @author <a href="mailto:jpiera@gvsig.org">Jorge Piera</a>

 */

public class Envelope2D extends DefaultEnvelope implements Cloneable{

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

        public static final String PERSISTENCE_DEFINITION_NAME = "Envelope2Dimensions";

        public Envelope2D() {

                super();

        }

        public Envelope2D(Point min, Point max) {

                super(min, max);

        }

        public Envelope2D(double minX,double minY,double maxX, double maxY){

                this(new Point2D(minX, minY), new Point2D(maxX, maxY));

        }

        /**

     * @param bounds

     */

    public Envelope2D(com.vividsolutions.jts.geom.Geometry bounds) {

        this();

        com.vividsolutions.jts.geom.Envelope envelope = bounds.getEnvelope().getEnvelopeInternal();

        setLowerCorner(new Point2D(envelope.getMinX(), envelope.getMinY()));

        setUpperCorner(new Point2D(envelope.getMaxX(), envelope.getMaxY()));

    }

    /**

     * @param envelope

     */

    public Envelope2D(com.vividsolutions.jts.geom.Envelope envelope) {

        this(new Point2D(envelope.getMinX(),envelope.getMinY()), new Point2D(envelope.getMaxX(),envelope.getMaxY()));

    }

    /* (non-Javadoc)

         * @see org.gvsig.fmap.geom.primitive.Envelope#getDimension()

         */

        public int getDimension() {

                return 2;

        }

        /*

         * (non-Javadoc)

         * @see org.gvsig.fmap.geom.primitive.Envelope#convert(org.cresques.cts.ICoordTrans)

         */

        public Envelope convert(ICoordTrans trans) {

            if (isEmpty){

            throw new EnvelopeNotInitializedException();

        }

            if (trans == null) {

                // clone

                return new Envelope2D(this.getLowerCorner(), this.getUpperCorner());

            }

            if (this.getDimension() > 2) {

                        return null;

                }

            // We'll reproject by taking samples like this:

            //

            //  *---*---*---*---*

        //  |               |

        //  *   *   *   *   *

        //  |               |

        //  *   *   *   *   *

        //  |               |

        //  *   *   *   *   *

        //  |               |

        //  *---*---*---*---*

            //

            // This is because:

            //

            // - In some CRS (for example EPSG:4326) the north/south pole is a "line"

            //   while in other CRS the north/south pole is a point, so if you

            //   reproject the bounding box of the world, the result can be absurd.

            // - Sometimes the top/bottom/right/bottom of one envelope

            //   corresponds to a strange point in the the other envelope

            //   (not even a point in the perimeter)

            // - More generally, reprojecting usually implies a rotation (the result

            //   is a rotated envelope) so it's better to use a few

            //   samples along the perimeter.

            double xmin = getMinimum(0);

            double ymin = getMinimum(1);

            double step_w = 0.0625 * (getMaximum(0) - xmin);

            double step_h = 0.125 * (getMaximum(1) - ymin);

            java.awt.geom.Point2D sample = null;

            java.awt.geom.Point2D sample_trans = null;

            // Init with worst values

        java.awt.geom.Point2D res_min = new java.awt.geom.Point2D.Double(

            Double.POSITIVE_INFINITY, Double.POSITIVE_INFINITY);

        java.awt.geom.Point2D res_max = new java.awt.geom.Point2D.Double(

            Double.NEGATIVE_INFINITY, Double.NEGATIVE_INFINITY);

        int added = 0;

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

                for (int j=0; j<=8; j++) {

                        sample = new java.awt.geom.Point2D.Double(

                            xmin + i * step_w,

                            ymin + j * step_h);

                        sample_trans = new java.awt.geom.Point2D.Double(0,0);

                        try {

                            sample_trans = trans.convert(sample, sample_trans);

                        } catch (Exception exc) {

                            // Unable to convert this one: ignore

                            continue;

                        }

                        // Update max/min found

                        // X

                        if (sample_trans.getX() > res_max.getX()) {

                            res_max.setLocation(sample_trans.getX(), res_max.getY());

                        }

                    if (sample_trans.getX() < res_min.getX()) {

                        res_min.setLocation(sample_trans.getX(), res_min.getY());

                    }

                        // Y

                    if (sample_trans.getY() > res_max.getY()) {

                        res_max.setLocation(res_max.getX(), sample_trans.getY());

                    }

                    if (sample_trans.getY() < res_min.getY()) {

                        res_min.setLocation(res_min.getX(), sample_trans.getY());

                    }

                    added++;

                }

            }

            if (added == 0) {

                logger.error("Unable to reproject envelope with transf: " + trans.toString());

                return null;

            }

                return new Envelope2D(

                    res_min.getX(),

                    res_min.getY(),

                    res_max.getX(),

                    res_max.getY());

        }

        public static void registerPersistent() {

                PersistenceManager manager = ToolsLocator.getPersistenceManager();

                if( manager.getDefinition(PERSISTENCE_DEFINITION_NAME)==null ) {

                        DynStruct definition = manager.addDefinition(

                                        Envelope2D.class,

                                        PERSISTENCE_DEFINITION_NAME,

                                        "Envelope2D persistence definition",

                                        null,

                                        null

                        );

                        definition.extend(manager.getDefinition(DefaultEnvelope.PERSISTENCE_DEFINITION_NAME));

                }

        }

    public Object clone() throws CloneNotSupportedException {

        return super.clone();

    }

    private void createPoints() {

        this.min = new Point2D(0, 0);

        this.max = new Point2D(0, 0);

    }

    public void add(Envelope envelope) {

        int i;

        if( envelope==null || envelope.isEmpty() ) {

            return;

        }

        if (isEmpty){

            createPoints();

            for (i=0;i<getDimension();i++){

                this.min.setCoordinateAt(i, envelope.getMinimum(i));

                this.max.setCoordinateAt(i, envelope.getMaximum(i));

            }

            isEmpty = false;

        } else {

            for (i=0;i<getDimension();i++){

                this.min.setCoordinateAt(i,

                    Math.min(this.min.getCoordinateAt(i), envelope.getMinimum(i)));

                this.max.setCoordinateAt(i,

                    Math.max(this.max.getCoordinateAt(i), envelope.getMaximum(i)));

            }

        }

    }

}