Application: gvSIG desktop

Este algoritmo tiene varios problemas sobre los que deberíamos tomar algunas decisiones:

1. El algoritmo no hace lo que dice la descripción (al menos en español, he mirado en otros idiomas como inglés, italiano,... y allí la descripción en más vaga):

En la descripción dice que va a generar el Bounding Box de un grupo de geometrías pero el algoritmo lo que crea es el bounding box de cada una de las geometrías (sin unirlas en una única feature) copiando todos los atributos de la capa original.
Deberíamos decidir si:

1. Arreglamos el algoritmo para que haga lo que dice la descripción
2. Arreglamos la descripción para que describa lo que hace el algoritmo. Aunque debería parecer que tiene más sentido arreglar el algoritmo, me decanto por arreglar la descripción. Si optamos por 1, perderíamos la posibilidad de crear una capa de envelopes de geometrías, mientras que si optamos por 2 siempre podemos generar una unión de esos envelopes para generar un envelope global.
3. O, incluso, creamos otro algoritmo que haga lo que decidimos no hacer aquí, pero eso ya sería una Feature Request.

1. El algoritmo siempre crea como capa de salida una capa de polígonos, pero el envelope de un punto es un punto y al convertirlo a las geometrías de gvSIG se transforman también en un punto. Después, cuando se intenta guardar ese punto en la capa de polígonos se produce un error y no se añade la feature ==> capa vacía.

Deberíamos decidir, cuando la capa de origen es de puntos, si:

1. en lugar de crear una capa de polígonos, creamos una capa de puntos, con lo que prácticamente en la capa de salida tendríamos lo mismo que en la capa de entrada
2. o si restringimos este algoritmo para que no esté activo con una capa de puntos.

Hi Francisco,
I'm sorry but I didn't check the algorithm behaviour for lines and polygons, when I've opened this ticket. I have read its help and thought that it was not working fine only for point layers.
I think that a good solution could be to simply add a checkbox in the geoalgorithm gui, asking the user if he/she wants the global or local bbox. Maybe this is a feature request.

This how it could be implemented in a Jython script in gvSIG 2.3:

# encoding: utf-8

from gvsig.libs.toolbox import *

from es.unex.sextante.additionalInfo import AdditionalInfoVectorLayer
from es.unex.sextante.dataObjects.vectorFilters import BoundingBoxFilter
from com.vividsolutions.jts.operation.union import CascadedPolygonUnion

from java.lang import Integer

class BoundingBox(ToolboxProcess):

    def defineCharacteristics(self):

        self.setName("Bounding Box")
        self.setGroup("My algorithms")

        params = self.getParameters()

        params.addInputVectorLayer("INPUT",           # name
            Sextante.getText("Vector_layer"),         # description
            AdditionalInfoVectorLayer.SHAPE_TYPE_ANY, # shape type
            1)                                        # mandatory (0=false, 1=true)

        params.addBoolean("DISSOLVE", Sextante.getText("Dissolve_entities"), 0) # 0 = false

        self.addOutputVectorLayer("RESULT",           # name
            Sextante.getText("Bounding_Box")          # description
            )                       

    def processAlgorithm(self):

        params = self.getParameters()
        layer = params.getParameterValueAsVectorLayer("INPUT")
        if not self.isAutoExtent():
            layer.addFilter(BoundingBoxFilter(self.getAnalysisExtent()))
        dissolve = params.getParameterValueAsBoolean("DISSOLVE")

        if dissolve: 
            outputFieldTypes = []
            outputFieldTypes.append(Integer)
            outputFieldNames = []
            outputFieldNames.append("ID")
            output = self.getNewVectorLayer("RESULT", 
                Sextante.getText("Bounding Box"),
                SHAPE_TYPE_POLYGON,
                outputFieldTypes,
                outputFieldNames)
            geometries = []
            outputValues = []
        else:
            output = self.getNewVectorLayer("RESULT", 
                Sextante.getText("Bounding Box"),
                SHAPE_TYPE_POLYGON,
                layer.getFieldTypes(),
                layer.getFieldNames())

        iShapeCount = layer.getShapesCount()
        iterator = layer.iterator()

        self.setRangeOfValues(0, iShapeCount)
        while iterator.hasNext():
            if self.isCanceled():
                break
            feature = iterator.next()
            geom = feature.getGeometry()

            if dissolve:
                geometries.append(geom)
            else:
                outputValues = feature.getRecord().getValues()
                output.addFeature(geom.getEnvelope(), outputValues)

            self.next()

        if dissolve:
            cpu = CascadedPolygonUnion(geometries)
            geom = cpu.union()
            outputValues.append(Integer(0))
            output.addFeature(geom.getEnvelope(), outputValues)

        return not(self.isCanceled())

def main(*args):

    # create the geoprocess
    process = BoundingBox()
    # register the geoprocess between those available in the "Scripting" group
    process.selfregister("Scripting")
    # update the Geoprocess toolbox
    process.updateToolbox()

    # inform the user that the geoprocess is added to the toolbox
    from org.gvsig.app import ApplicationLocator
    from javax.swing import JOptionPane
    application = ApplicationLocator.getManager()
    application.message("The '%s/%s/%s' geoprocess has been added to the Toolbox." % (
        "Scripting",
        process.getGroup(),
        process.getName()
        ),    
        JOptionPane.INFORMATION_MESSAGE)