| gvSIG desktop 1

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

 *

 * Copyright (C) 2004 IVER T.I. 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

 *

 *   IVER T.I. S.A

 *   Salamanca 50

 *   46005 Valencia

 *   Spain

 *

 *   +34 963163400

 *   dac@iver.es

 */

// 18/09/2007 fjp

// @author: Fco. Jos? Pe?arrubia        fpenarru@gmail.com

package org.gvsig.graph.solvers;

import java.io.File;

import java.sql.Types;

import java.util.ArrayList;

import java.util.HashMap;

import java.util.HashSet;

import java.util.Iterator;

import java.util.Map;

import java.util.Set;

import org.gvsig.exceptions.BaseException;

import org.gvsig.fmap.algorithm.triangulation.visad.DelaunayWatson;

import org.gvsig.graph.core.GvEdge;

import org.gvsig.graph.core.GvNode;

import org.gvsig.graph.core.IGraph;

import org.gvsig.graph.core.Network;

import org.gvsig.graph.core.NetworkUtils;

import com.hardcode.gdbms.driver.exceptions.InitializeDriverException;

import com.hardcode.gdbms.driver.exceptions.InitializeWriterException;

import com.hardcode.gdbms.driver.exceptions.ReadDriverException;

import com.hardcode.gdbms.engine.values.Value;

import com.hardcode.gdbms.engine.values.ValueFactory;

import com.iver.cit.gvsig.fmap.core.DefaultFeature;

import com.iver.cit.gvsig.fmap.core.FShape;

import com.iver.cit.gvsig.fmap.core.IGeometry;

import com.iver.cit.gvsig.fmap.core.ShapeFactory;

import com.iver.cit.gvsig.fmap.core.v02.FConverter;

import com.iver.cit.gvsig.fmap.drivers.FieldDescription;

import com.iver.cit.gvsig.fmap.drivers.SHPLayerDefinition;

import com.iver.cit.gvsig.fmap.edition.DefaultRowEdited;

import com.iver.cit.gvsig.fmap.edition.IRowEdited;

import com.iver.cit.gvsig.fmap.edition.writers.shp.ShpWriter;

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

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

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

import com.vividsolutions.jts.algorithm.CGAlgorithms;

import com.vividsolutions.jts.algorithm.CentroidPoint;

import com.vividsolutions.jts.geom.Coordinate;

import com.vividsolutions.jts.geom.Geometry;

import com.vividsolutions.jts.geom.GeometryFactory;

import com.vividsolutions.jts.geom.LineString;

import com.vividsolutions.jts.geom.LinearRing;

import com.vividsolutions.jts.geom.Polygon;

/**

 * @author fjp

 * FROM BRANCH

 * 

 * This class can label nodes with distances and costs to a flag. You will

 * obtain a temp shp layer with fields IdArc, IdEdge, CostOrig, DistOrig,

 * CostEnd, DistEnd, IdFlag

 * 

 * La diferencia con ServiceAreaExtractor es que esta versi?n escucha al

 * algoritmo Dijkstra, y va montando el shp de l?neas conforme va siendo

 * explorada la red. La gran ventaja de hacerlo as? es que no dependes del

 * tama?o de la red. Solo recorres los tramos y nodos que exploras, de forma que

 * si limitas el ?rea de servicio a una distancia m?xima, la red solo se explora

 * hasta esa distancia / coste.

 * 

 */

public class ServiceAreaExtractor2_From_Branch implements IDijkstraListener {

        private static String tempDirectoryPath = System

                        .getProperty("java.io.tmpdir");

        static FieldDescription[] fields = new FieldDescription[7];

        static {

                FieldDescription fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("IDARC");

                fieldDesc.setFieldType(Types.INTEGER);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(0);

                fields[0] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("IDEDGE");

                fieldDesc.setFieldType(Types.INTEGER);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(0);

                fields[1] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("COSTORIG");

                fieldDesc.setFieldType(Types.DOUBLE);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fields[2] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("DISTORIG");

                fieldDesc.setFieldType(Types.DOUBLE);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fields[3] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("COSTEND");

                fieldDesc.setFieldType(Types.DOUBLE);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fields[4] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("DISTEND");

                fieldDesc.setFieldType(Types.DOUBLE);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fields[5] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("IDFLAG");

                fieldDesc.setFieldType(Types.INTEGER);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fields[6] = fieldDesc;

        }

        static FieldDescription[] fieldsPol = new FieldDescription[2];

        static {

                FieldDescription fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("COST");

                fieldDesc.setFieldType(Types.DOUBLE);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fieldsPol[0] = fieldDesc;

                fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("IDFLAG");

                fieldDesc.setFieldType(Types.INTEGER);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fieldsPol[1] = fieldDesc;

        }

        static GeometryFactory gf = new GeometryFactory();

        private class VisitedEdge {

                private GvEdge edge;

                private double percentcost;

                public VisitedEdge(GvEdge edge) {

                        this.edge = edge;

                        IGraph g = net.getGraph();

                        GvNode nOrig = g.getNodeByID(edge.getIdNodeOrig());

                        double maxCost = costs[costs .length-1];

                        double costCalculated = nOrig.getBestCost() + edge.getWeight();

                        if (costCalculated < maxCost)

                                percentcost = 1.0;

                        else

                        {

                                double percentCostCalculated = (maxCost - nOrig.getBestCost())/ edge.getWeight();

                                percentcost = percentCostCalculated;

                        }

                }

                public GvEdge getEdge() {

                        return edge;

                }

                public double getPercentcost() {

                        return percentcost;

                }

                public void setPercentCost(double d) {

                        this.percentcost = d;

                }

        }

        private Network net;

        private ShpWriter shpWriter;

        private ShpWriter shpWriterPol;

        private ShpWriter shpWriterTri;

        private File fTempPol;

        private File fTempTri;

        private SHPLayerDefinition layerDefPol;

        private SHPLayerDefinition layerDefTri;

        private HashMap<String, VisitedEdge> visitedEdges = new HashMap();

        private File fTemp;

        private SHPLayerDefinition layerDef;

        private int idFlag;

        private ReadableVectorial adapter;

//        private double maxCost;

        private Geometry serviceArea;

        private ArrayList <Geometry> serviceAreaPolygons;

        private double[] costs = null;        

        private boolean bDoCompactArea = false;

        private ArrayList borderCoords = new ArrayList();

        private HashSet<Coordinate> nodes;

//        DelaunayFast tri2;

        DelaunayWatson tri2;

        /**

         * @param net

         * @throws Exception 

         * @throws InitializeWriterException

         * @throws ReadDriverException

         * @throws InitializeDriverException

         */

        public ServiceAreaExtractor2_From_Branch(Network net) throws Exception {

                this.net = net;

                int aux = (int) (Math.random() * 1000);

                nodes = new HashSet<Coordinate>();

                String nameLine = "tmpServiceAreaLine" + aux + ".shp";

                String namePol = "tmpServiceAreaPol" + aux + ".shp";

                String nameTri = "tmpTri" + aux + ".shp";

                fTemp = new File(tempDirectoryPath + "/" + nameLine );

                fTempPol = new File(tempDirectoryPath + "/" + namePol );

                fTempTri = new File(tempDirectoryPath + "/" + nameTri );

                layerDef = new SHPLayerDefinition();

                layerDef.setFile(fTemp);

                layerDef.setName(nameLine);                

                layerDef.setFieldsDesc(fields);

                layerDef.setShapeType(FShape.LINE);

                layerDefPol = new SHPLayerDefinition();

                layerDefPol.setFile(fTempPol);

                layerDefPol.setName(namePol);                

                layerDefPol.setFieldsDesc(fieldsPol);

                layerDefPol.setShapeType(FShape.POLYGON);

                layerDefTri = new SHPLayerDefinition();

                layerDefTri.setFile(fTempTri);

                layerDefTri.setName(nameTri);                

                layerDefTri.setFieldsDesc(fieldsPol);

                layerDefTri.setShapeType(FShape.POLYGON);

                shpWriter = new ShpWriter();

                shpWriter.setFile(fTemp);

                shpWriter.initialize(layerDef);

                shpWriterPol = new ShpWriter();

                shpWriterPol.setFile(fTempPol);

                shpWriterPol.initialize(layerDefPol);

                shpWriterTri = new ShpWriter();

                shpWriterTri.setFile(fTempTri);

                shpWriterTri.initialize(layerDefTri);

                shpWriter.preProcess();

                shpWriterPol.preProcess();

                shpWriterTri.preProcess();

                FLyrVect lyr = net.getLayer();

                adapter = lyr.getSource();

                adapter.start();

                serviceAreaPolygons = new ArrayList<Geometry>();

        }

        /**

         * Devuelve el ?ndice del intervalo m?s alto que contiene a ese valor.

         * @param bestCost

         * @param costs

         * @return

         */

        private int getCostInterval(double bestCost, double[] costs) {

                int ret = 0;

                if (bestCost > costs[costs.length-1])

                        return -1;

                for (int i=costs.length-1; i>=0; i--) {

                        if (bestCost > costs[i])

                        {

                                ret = i+1;

                                break;

                        }

                }

                return ret;

        }

        /**

         * We process each edge and prepare a list of polygons, classified by

         * cost

         * @param edge

         * @param nodeOrig

         * @param nodeEnd

         * @param geom

         * @param costs

         */

        private void processEdgeForPolygon(GvEdge edge, GvNode nodeOrig, GvNode nodeEnd, IGeometry geom, double[] costs) {

                if (nodeEnd.getBestCost() > nodeOrig.getBestCost())

                {                

                        // miramos en qu? pol?gono cae ese edge POR COMPLETO 

                        // El coste de su punto final es menor que uno de los costes.

                        int indexInterval = getCostInterval(nodeEnd.getBestCost(), costs);

                        // Un pol?gono por cada zona

                        Geometry jtsGeom = geom.toJTSGeometry();

                        if (indexInterval != -1)

                        {

                                for (int i=costs.length-1; i >= indexInterval; i--) {

                                        calculateConvexHull(jtsGeom, i);

                                }

                        }

                        double maxCost = costs[costs.length-1];

                        // Es -1 si caso l?mite externo

                        if (indexInterval < costs.length-1)

                        {

                                // Caso l?mite externo

                                if ((nodeEnd.getBestCost() > maxCost) &&                                                

                                                (nodeOrig.getBestCost() < maxCost))

                                {

                                        double pct = (maxCost - nodeOrig.getBestCost())/ edge.getWeight();

                                        LineString partial = NetworkUtils.getPartialLineString(jtsGeom, pct, edge.getDirec());

                                        calculateConvexHull(partial, costs.length-1);

                                        return;

                                }

                                // Parcial interno

                                maxCost = costs[indexInterval+1];

                                if ((nodeOrig.getBestCost() < maxCost) &&

                                                (nodeEnd.getBestCost() > maxCost)) 

                                {

                                        // A ese tramo hemos llegado parcialmente

                                        double pct = (maxCost - nodeOrig.getBestCost())/ edge.getWeight();

                                        try {

                                                LineString partial = NetworkUtils.getPartialLineString(jtsGeom, pct, edge.getDirec());

                                                calculateConvexHull(partial, indexInterval+1);                                                        

                                        }

                                        catch (Exception e)

                                        {

                                                e.printStackTrace();

                                        }

                                }

                        }

                } 

        }

        /**

         * @param jtsGeom

         * @param i

         */

        private void calculateConvexHull(Geometry jtsGeom, int i) {

                if (serviceAreaPolygons.size() <= i) { // se crea por primera vez

                        Geometry gIni = jtsGeom; 

                        serviceAreaPolygons.add(i, gIni);

                }

                else

                {

                        Geometry antG = serviceAreaPolygons.get(i);

                        if (antG == null)

                                antG = jtsGeom;

                        else

                        {

                                antG = antG.union(jtsGeom);                                

                        }

                        antG = antG.convexHull();

                        serviceAreaPolygons.set(i, antG);

                }

        }

        private void writePartialEdge(int i, IGeometry geom, GvEdge edge, GvNode nodeOrig, GvNode nodeEnd, int idFlag, double maxCost) throws BaseException {

                Geometry jtsGeom = geom.toJTSGeometry();

                double pct = (maxCost - nodeOrig.getBestCost())/ edge.getWeight();

                if (edge.getDirec() == 0) // Sentido inverso

                        pct = 1-pct;

                LineString partial = NetworkUtils.getPartialLineString(jtsGeom, pct, edge.getDirec());

                if (serviceArea == null)

                        serviceArea = partial;

                else

                {

                        serviceArea = serviceArea.union(partial);

                        serviceArea = serviceArea.convexHull();

                }

                IGeometry newGeom = FConverter.jts_to_igeometry(partial);

                Value[] values = new Value[7];

                values[0] = ValueFactory.createValue(i);

                values[1] = ValueFactory.createValue(edge.getIdEdge());

                values[2] = ValueFactory.createValue(nodeOrig.getBestCost());

                values[3] = ValueFactory.createValue(nodeOrig.getAccumulatedLength());

                values[4] = ValueFactory.createValue(maxCost);

                values[5] = ValueFactory.createValue(nodeOrig.getAccumulatedLength() + edge.getDistance()*pct);

                values[6] = ValueFactory.createValue(idFlag);

                DefaultFeature feat = new DefaultFeature(newGeom, values);

                IRowEdited row = new DefaultRowEdited(feat, DefaultRowEdited.STATUS_ADDED, i);

                shpWriter.process(row);

//                // TODO: TROZO DEL RESTO

//                int direc = 0;

//                pct = pct+0.02;

//                if (edge.getDirec() == 0)

//                {

//                        direc = 1;

//                        pct = pct - 0.04;

//                }

//                LineString partial2 = NetworkUtils.getPartialLineString(jtsGeom, pct, direc);                

//                IGeometry newGeom2 = FConverter.jts_to_igeometry(partial2);

//                values[6] = ValueFactory.createValue(-2);

//                DefaultFeature feat2 = new DefaultFeature(newGeom2, values);

//                IRowEdited row2 = new DefaultRowEdited(feat2, DefaultRowEdited.STATUS_ADDED, i);

//                shpWriter.process(row2);

                borderCoords.add(partial.getCoordinateN(partial.getNumPoints()-1)); 

                if (bDoCompactArea) {

//                        processCompact(partial);

                        Coordinate cLimit = null;

                        if (edge.getDirec() == 0) // Sentido inverso

                                cLimit = partial.getCoordinateN(0);

                        else

                                cLimit = partial.getCoordinateN(partial.getNumPoints()-1);

                        processCompact(cLimit.x, cLimit.y);

//                        processCompact(nodeEnd.getX(), nodeEnd.getY());

                }

        }

        public FLyrVect getBorderPoints() throws Exception {

                Value[] values = new Value[1];

                values[0] = ValueFactory.createValue(costs[costs.length-1]);

                File fTempPoints = new File(tempDirectoryPath + "/borderPoints.shp");

                FieldDescription[] fieldsPoints = new FieldDescription[1];

                FieldDescription fieldDesc = new FieldDescription();

                fieldDesc.setFieldName("COST");

                fieldDesc.setFieldType(Types.DOUBLE);

                fieldDesc.setFieldLength(20);

                fieldDesc.setFieldDecimalCount(5);

                fieldsPoints[0] = fieldDesc;

                SHPLayerDefinition layerDef = new SHPLayerDefinition();

                layerDef.setFile(fTempPoints);

                layerDef.setName("BorderPoints");                

                layerDef.setFieldsDesc(fieldsPoints);

                layerDef.setShapeType(FShape.POINT);

                ShpWriter shpWriter = new ShpWriter();

                shpWriter.setFile(fTempPoints);

                shpWriter.initialize(layerDef);

                int i=0;

                for (Iterator it = borderCoords.iterator(); it.hasNext();) {

                        Coordinate c = (Coordinate) it.next();

                        IGeometry geom = ShapeFactory.createPoint2D(c.x, c.y);

                        DefaultFeature feat = new DefaultFeature(geom, values);

                        IRowEdited row = new DefaultRowEdited(feat,

                                        DefaultRowEdited.STATUS_ADDED, i++);

                        shpWriter.process(row);

                }

                shpWriter.postProcess();

                FLyrVect lyr = (FLyrVect) LayerFactory.createLayer(layerDef.getName(), "gvSIG shp driver", 

                                layerDef.getFile(), null);

                return lyr;

        }

        private void processCompact(LineString partial) {

                Coordinate cIni = partial.getCoordinateN(0);

                Coordinate cEnd = partial.getCoordinateN(partial.getNumPoints()-1);

//                System.out.println("PARTIAL c1=" + cIni + " cEnd=" + cEnd);

                processCompact(cIni.x, cIni.y);

                processCompact(cEnd.x, cEnd.y);

        }

        private void writeTotalEdge(int i, IGeometry geom, GvEdge edge,

                        GvNode nodeOrig, GvNode nodeEnd, int idFlag)

                        throws Exception {

                Value[] values = new Value[7];

                values[0] = ValueFactory.createValue(i);

                values[1] = ValueFactory.createValue(edge.getIdEdge());

                values[2] = ValueFactory.createValue(nodeOrig.getBestCost());

                values[3] = ValueFactory.createValue(nodeOrig.getAccumulatedLength());

                values[4] = ValueFactory.createValue(nodeEnd.getBestCost());

                values[5] = ValueFactory.createValue(nodeEnd.getAccumulatedLength());

                values[6] = ValueFactory.createValue(idFlag);

                if (bDoCompactArea) {

//                        System.out.println(" c1=" + cIni + " cEnd=" + cEnd);

//                        processCompact(nodeOrig.getX(), nodeOrig.getY());

//                        processCompact(nodeEnd.getX(), nodeEnd.getY());

                }

                DefaultFeature feat = new DefaultFeature(geom, values);

                IRowEdited row = new DefaultRowEdited(feat,

                                DefaultRowEdited.STATUS_ADDED, i);

                shpWriter.process(row);

        }

        private void processCompact(double x, double y) {

//                FPoint2D p = new FPoint2D(x,y);

                Coordinate c = new Coordinate(x,y);

//                System.out.println("PARTIAL c1=" + cIni + " cEnd=" + cEnd);

                if (!nodes.contains(c))

                        nodes.add(c);

                else

                {

                        System.out.print("Nodo ya contenido");

                }

        }

        public boolean adjacentEdgeVisited(GvNode fromNode, GvEdge edge) {

                insertVisitedEdge(edge);

                return false;

        }

        /**

         * Si el coste m?nimo del edge > costemax, salimos del m?todo.

         * Miramos si edge est? ya en la lista.

         * Si no est?, lo a?adimos.

         * Si est?, hay que mirar el porcentaje recorrido sobre ese tramo.

         * Casos posibles:

         * EdgeA al 100 % => No se a?ade este.

         * EdgeA.percentCost < 1.0. 

         *                 Comprobamos el percent de el nuevo. Si entre los dos suman > 1.0

         *                 marcamos el antiguo al 1.0 para que se escriba el tramo completo.

         * 

         *                 Si no suman 1.0, hay que a?adir este nuevo Edge, con el porcentaje

         *                 correspondiente.

         * @param edge

         */

        private void insertVisitedEdge(GvEdge edge) {

                IGraph g = net.getGraph();

                GvNode n1 = g.getNodeByID(edge.getIdNodeOrig());

                double maxCost = costs[costs .length-1];

                if (n1.getBestCost() > maxCost)

                        return; // edge outside service area.

                String key = "" + edge.getIdArc();

                if (!visitedEdges.containsKey(key))        

                {

                        visitedEdges.put(key, new VisitedEdge(edge));

//                        System.out.println("idEdge adjacent= " + edge.getIdEdge());

                }

                else

                {

                        VisitedEdge edgeAnt = visitedEdges.get(key);

//                        GvEdge savedEdge = edgeAnt.getEdge();

                        if (edgeAnt.getPercentcost() == 1.0)

                                return; // Ya est? completo, no a?adimos nada.

                        double percentCostCalculated = (maxCost - n1.getBestCost())/ edge.getWeight();

                        if ((percentCostCalculated + edgeAnt.getPercentcost()) >= 1.0)

                                edgeAnt.setPercentCost(1.0);

                        else

                        {

                                visitedEdges.put(key + "_", new VisitedEdge(edge));

                        }

                }

        }

        public boolean minimumCostNodeSelected(GvNode node) {

//                IGraph g = net.getGraph();

//                int idEdge = node.getFromLink();

//                if (idEdge == -1) 

//                        return false;

//                GvEdge edge = g.getEdgeByID(idEdge);

//                insertVisitedEdge(edge);

                return false; // true if we want to stop Dijkstra

        }

        public void setIdFlag(int idFlag) {

                this.idFlag = idFlag;

        }

        /**

         * Write edges and polygons associated with active flag and costs

         * @param costs

         * @throws BaseException

         */

        public void writeServiceArea() throws Exception {

                Set<Map.Entry<String, VisitedEdge>> keySet = visitedEdges.entrySet();

                GvEdge edge;

                IGraph g = net.getGraph();

//                Integer idEdge;

                double maxCost = costs[costs .length-1];

                serviceAreaPolygons = new ArrayList<Geometry>(costs.length);

                for (int i=0; i < costs.length-1; i++)

                        serviceAreaPolygons.add(null);

                for (Map.Entry<String, VisitedEdge> entry : keySet) {

//                        idEdge = entry.getKey();

                        VisitedEdge visitedEdge = entry.getValue(); 

                        edge = visitedEdge.getEdge();

                        GvNode nodeEnd = g.getNodeByID(edge.getIdNodeEnd());

                        GvNode nodeOrig = g.getNodeByID(edge.getIdNodeOrig());

                        IGeometry geom;

                        try {

                                geom = adapter.getShape(edge.getIdArc());

                                processEdgeForPolygon(edge, nodeOrig, nodeEnd, geom, costs);

                                double costAux = nodeOrig.getBestCost() + edge.getWeight();

//                                if (nodeEnd.getBestCost() > nodeOrig.getBestCost())

                                {

                                        // TODO: REVISAR EL CASO EN EL QUE HAS LLEGADO POR LOS

                                        // 2 LADOS PERO HAY UN TRAMO INALCANZABLE ENMEDIO

                                        if (visitedEdge.getPercentcost() == 1.0) {

                                                // A ese tramo hemos llegado por completo

                                                // Recuperamos su distancia y etiquetamos.

                                                writeTotalEdge(edge.getIdArc(), geom, edge, nodeOrig, nodeEnd, idFlag);        

                                        }

                                        else

                                        {

                                                if (nodeOrig.getBestCost() < maxCost) {

                                                        // A ese tramo hemos llegado parcialmente

                                                        // Recuperamos su distancia y etiquetamos.

                                                        writePartialEdge(edge.getIdArc(), geom, edge, nodeOrig, nodeEnd, idFlag, maxCost);        

                                                }

                                        } // else

                                } // if nodeEnd > nodeOrig

                        } catch (Exception e) {

                                e.printStackTrace();

                                throw new RuntimeException(e);

                        }

                } // for

                for (int j=serviceAreaPolygons.size()-1; j>=0; j--) {

                        Geometry jtsGeom = serviceAreaPolygons.get(j);

                        writePolygon(idFlag, costs[j], jtsGeom);

                }

                if (bDoCompactArea) {

                        calculateTriangulation();

                }

        }

        private void calculateTriangulation() {

                // Centroide

                CentroidPoint centroid = new CentroidPoint();

                for (Iterator it = nodes.iterator(); it.hasNext();)

                {

                        Coordinate node = (Coordinate) it.next();

                        centroid.add(node);

                }

                nodes.add(centroid.getCentroid());

                int numPoints = nodes.size();

            double[][] samples = new double[2][numPoints];

            double[] samp0 = samples[0];

            double[] samp1 = samples[1];

            Iterator it = nodes.iterator();

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

                    Coordinate node = (Coordinate) it.next();

                    samp0[i] = node.x;

                    samp1[i] = node.y;

            }

            try {

                        tri2 = new DelaunayWatson(samples);

                        tri2.improve(samples, 10);

                        System.out.println("Fin de trayecto. Num. tri?ngulos=" + tri2.Tri.length);

                        for (int i=0; i< tri2.Tri.length; i++) {

                              writeTri(tri2.Tri[i], samples);

                    }

                } catch (Exception e) {

                        // TODO Auto-generated catch block

                        e.printStackTrace();

                }

        }

        private void writeTri(int[] vertex, double[][] points) throws BaseException {

                Value[] values = new Value[2];

                values[0] = ValueFactory.createValue(2.0);

                values[1] = ValueFactory.createValue(1);

                Coordinate c1 = new Coordinate(points[0][vertex[0]], points[1][vertex[0]]);

                Coordinate c2 = new Coordinate(points[0][vertex[1]], points[1][vertex[1]]);

                Coordinate c3 = new Coordinate(points[0][vertex[2]], points[1][vertex[2]]);

                Coordinate[] c = new Coordinate[4];

                c[0] = c1;

                c[1] = c3;

                c[2] = c2;

                c[3] = c1;

                LinearRing linRing = null;

                if (CGAlgorithms.isCCW(c))

                {

                        Coordinate[] ccw = new Coordinate[4];

                        ccw[0] = c1;

                        ccw[1] = c2;

                        ccw[2] = c3;

                        ccw[3] = c1;

                        linRing = gf.createLinearRing(ccw);

                }

                else

                {

                        linRing = gf.createLinearRing(c);

//                        return;

                }

                Polygon pol = gf.createPolygon(linRing, null);

                IGeometry geom = FConverter.jts_to_igeometry(pol);

                DefaultFeature feat = new DefaultFeature(geom, values);

                IRowEdited row = new DefaultRowEdited(feat, DefaultRowEdited.STATUS_ADDED, idFlag);

//                shpWriterPol.process(row);

                shpWriterTri.process(row);

        }

        private void writePolygon(int idFlag, double maxCost, Geometry jtsGeom) throws BaseException {

                Value[] values = new Value[2];

                values[0] = ValueFactory.createValue(maxCost);

                values[1] = ValueFactory.createValue(idFlag);

                IGeometry geom = FConverter.jts_to_igeometry(jtsGeom);

                DefaultFeature feat = new DefaultFeature(geom, values);

                IRowEdited row = new DefaultRowEdited(feat, DefaultRowEdited.STATUS_ADDED, idFlag);

                shpWriterPol.process(row);

        }

        /**

         * Close writers.

         * @throws Exception

         */

        public void closeFiles() throws Exception {

//                for (int j=serviceAreaPolygons.size()-1; j>=0; j--) {

//                        Geometry jtsGeom = serviceAreaPolygons.get(j);

//                        writePolygon(idFlag, costs[j], jtsGeom);

//                }

                shpWriter.postProcess();

                shpWriterPol.postProcess();

                shpWriterTri.postProcess();

                adapter.stop();

        }

        public double[] getCosts() {

                return costs;

        }

        public void setCosts(double[] costs) {

                this.costs = costs;

        }

        public FLyrVect getPolygonLayer() throws Exception {

                FLyrVect lyr = (FLyrVect) LayerFactory.createLayer(layerDefPol.getName(), "gvSIG shp driver", 

                                layerDefPol.getFile(), null);

                return lyr;

        }

        public FLyrVect getLineLayer() throws Exception {

                FLyrVect lyr = (FLyrVect) LayerFactory.createLayer(layerDef.getName(), "gvSIG shp driver", 

                                layerDef.getFile(), null);

                return lyr;

        }

        public boolean isDoCompactArea() {

                return bDoCompactArea;

        }

        public void setDoCompactArea(boolean doCompactArea) {

                bDoCompactArea = doCompactArea;

        }

}