| gvSIG desktop 1

    public static final int EQUALS = 0;

    public static final int DISJOINT = 1;

    public static final int INTERSECTS = 2;

    public static final int TOUCHES = 3;

    public static final int CROSSES = 4;

    public static final int WITHIN = 5;

    public static final int CONTAINS = 6;

    public static final int OVERLAPS = 7;

        // TODO: OJO!!!!. Est? implementado como un rectangulo.

        // Lo correcto deber?a ser calculando las distancias reales

        // es decir, con un c?rculo.

        Rectangle2D recPoint = new Rectangle2D.Double(p.getX() - (tolerance / 2),

                p.getY() - (tolerance / 2), tolerance, tolerance);

		return queryByRect(recPoint);

import java.util.List;

import com.iver.cit.gvsig.fmap.layers.FBitSet;

import com.iver.cit.gvsig.fmap.layers.FLyrVect;

import com.vividsolutions.jts.geom.Coordinate;

import com.vividsolutions.jts.geom.Envelope;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.geom.IntersectionMatrix;

    /* (non-Javadoc)

     * @see com.iver.cit.gvsig.fmap.operations.strategies.Strategy#queryByShape(com.iver.cit.gvsig.fmap.core.IGeometry, int)

     */

    public FBitSet queryByShape(IGeometry g, int relationship)

    throws DriverException, VisitException {

        // Si hay un ?ndice espacial, lo usamos para hacer el query.

        FLyrVect lyr = (FLyrVect) capa;

        if (lyr.getSpatialIndex() == null)

            return super.queryByShape(g, relationship);

        VectorialAdapter va = lyr.getSource();

        ICoordTrans ct = lyr.getCoordTrans();

        Rectangle2D bounds = g.getBounds2D();

        Coordinate c1 = new Coordinate(bounds.getMinX(), bounds.getMinY());

        Coordinate c2 = new Coordinate(bounds.getMaxX(), bounds.getMaxY());

        Envelope env = new Envelope(c1, c2);

        List lstRecs = lyr.getSpatialIndex().query(env);

        Integer idRec;

        FBitSet bitset = new FBitSet();

        Geometry jtsShape = g.toJTSGeometry();

        IntersectionMatrix m;

        int index;

        try {

            va.start();

            for (int i=0; i < lstRecs.size(); i++)

            {

                idRec = (Integer) lstRecs.get(i);

                index = idRec.intValue();

                IGeometry geom = va.getShape(index);

                if (ct != null) {

                    geom.reProject(ct);

                }                

                Geometry jtsGeom = geom.toJTSGeometry();

                switch (relationship) {

                case CONTAINS:

                    m = jtsShape.relate(jtsGeom);

                    if (m.isContains()) {

                        bitset.set(index, true);

                    }

                    break;

                case CROSSES:

                    m = jtsShape.relate(jtsGeom);

                    if (m.isCrosses(jtsGeom.getDimension(), jtsShape.getDimension())) {

                        bitset.set(index, true);

                    }

                    break;

                case DISJOINT:

                    m = jtsShape.relate(jtsGeom);            

                    if (m.isDisjoint()) {

                        bitset.set(index, true);

                    }

                    break;

                case EQUALS:

                    m = jtsShape.relate(jtsGeom); 

                    if (m.isEquals(jtsGeom.getDimension(), jtsShape.getDimension())) {

                        bitset.set(index, true);

                    }

                    break;

                case INTERSECTS:

                    m = jtsShape.relate(jtsGeom);

                    if (m.isIntersects()) {

                        bitset.set(index, true);

                    }

                    break;

                case OVERLAPS:

                    m = jtsShape.relate(jtsGeom);

                    if (m.isOverlaps(jtsGeom.getDimension(), jtsShape.getDimension()))

                    {

                        bitset.set(index, true);

                    }

                    break;

                case TOUCHES:

                    m = jtsShape.relate(jtsGeom);

                    if (m.isTouches(jtsGeom.getDimension(), jtsShape.getDimension()))

                    {

                        bitset.set(index, true);

                    }

                    break;

                case WITHIN:

                    m = jtsShape.relate(jtsGeom);

                    if (m.isWithin()) {

                        bitset.set(index, true);

                    }

                    break;

                }

            }

            va.stop();

        } catch (DriverIOException e) {

            // TODO Auto-generated catch block

            e.printStackTrace();

        }

        return bitset;

    }

    public FBitSet queryByRect(Rectangle2D rect) throws DriverException {

        // Si hay un ?ndice espacial, lo usamos para hacer el query.

        FLyrVect lyr = (FLyrVect) capa;

        if (lyr.getSpatialIndex() == null)

            return super.queryByRect(rect);

        VectorialAdapter va = lyr.getSource();

        ICoordTrans ct = lyr.getCoordTrans();

        Rectangle2D bounds = rect;

        Coordinate c1 = new Coordinate(bounds.getMinX(), bounds.getMinY());

        Coordinate c2 = new Coordinate(bounds.getMaxX(), bounds.getMaxY());

        Envelope env = new Envelope(c1, c2);

        List lstRecs = lyr.getSpatialIndex().query(env);

        Integer idRec;

        FBitSet bitset = new FBitSet();

        int index;

        try {

            va.start();

            for (int i=0; i < lstRecs.size(); i++)

            {

                idRec = (Integer) lstRecs.get(i);

                index = idRec.intValue();

                IGeometry geom = va.getShape(index);

                if (ct != null) {

                    geom.reProject(ct);

                }

                if (geom.intersects(rect))

                    bitset.set(index, true);

            }

            va.stop();

        } catch (DriverIOException e) {

            // TODO Auto-generated catch block

            e.printStackTrace();

        }

        return bitset;

    }

    /* (non-Javadoc)

     * @see com.iver.cit.gvsig.fmap.operations.strategies.Strategy#queryByPoint(java.awt.geom.Point2D, double)

     */

    public FBitSet queryByPoint(Point2D p, double tolerance)

    throws DriverException {

        // TODO: OJO!!!!. Est? implementado como un rectangulo.

        // Lo correcto deber?a ser calculando las distancias reales

        // es decir, con un c?rculo.

        Rectangle2D recPoint = new Rectangle2D.Double(p.getX() - (tolerance / 2),

                p.getY() - (tolerance / 2), tolerance, tolerance);

        return queryByRect(recPoint);

    }

import com.iver.cit.gvsig.fmap.layers.FBitSet;

	void createSpatialIndex();

	FBitSet queryByRect(Rectangle2D rect) throws DriverException;