# -*- coding: utf-8 -*-
#
# File: gvpy.py
# Version: v0.4
#

__author__ = """Oscar Martinez Olmos <masquesig@gmail.com>"""

import gvsig
import geom
import gvsig_raster
from org.gvsig.app import ApplicationLocator
from org.gvsig.fmap.dal import DALLocator
import os.path
from org.gvsig.andami import PluginsLocator
from org.gvsig.fmap.mapcontext import MapContextLocator
import java.awt
import java.awt.event
import java.awt.geom
from java.io import File

def addDependencyWithPlugin(pluginCode):
    pluginsManager = PluginsLocator.getManager()
    scriptingPlugin = pluginsManager.getPlugin("org.gvsig.scripting.app.extension")
    scriptingPlugin.addDependencyWithPlugin(pluginsManager.getPlugin(pluginCode))

addDependencyWithPlugin("org.gvsig.geoprocess.app.mainplugin")

from es.unex.sextante.core import Sextante, OutputFactory, AnalysisExtent
from es.unex.sextante.outputs import FileOutputChannel
from es.unex.sextante.gui.core import SextanteGUI
from es.unex.sextante.dataObjects import IRasterLayer

from org.gvsig.geoprocess.lib.sextante.dataObjects import FlyrVectIVectorLayer, FLyrRasterIRasterLayer, TableDocumentITable
from org.gvsig.fmap.mapcontext.layers import FLayer
from java.awt.geom import RectangularShape, Rectangle2D
from org.gvsig.fmap.mapcontext.layers.vectorial import FLyrVect
from org.gvsig.raster.fmap.layers import DefaultFLyrRaster

#Constant
TYPE_POLYGON = 0
TYPE_LINE = 1
TYPE_POINT = 2

#With True will appear a lot of extra info about all steps through the process
DEV_INFO = False

class Geoprocess:
  def __init__(self):
    Sextante.initialize()
    SextanteGUI.initialize()
    self.__FlyrVectIVectorLayer = None
    self.__outputFactory = SextanteGUI.getOutputFactory()
    self.__algorithms = dict()
    self.__defaultAE = None #default AExtension

    keySetIterator = Sextante.getAlgorithms().keySet().iterator()
    while(keySetIterator.hasNext()):
      key = keySetIterator.next()
      algorithms = Sextante.getAlgorithms().get(str(key))
      for name in algorithms.keySet():
        self.__algorithms[str(name)] = algorithms.get(name)

  def __createSextanteLayer(self, layer):
      """ gvsig layer -> SEXTANTE """
      slayer = FlyrVectIVectorLayer()
      slayer.create(layer)
      return slayer

  def __createSextanteRaster(self, layer):
      """ gvsig raster -> SEXTANTE """
      rlayer = FLyrRasterIRasterLayer()
      rlayer.create(layer)
      return rlayer

  def __createSextanteTable(self, layer):
      """ gvsig table -> SEXTANTE """
      table = TableDocumentITable()
      table.create(layer)
      return table

  def getAlgorithms(self):
      return self.__algorithms

  def __defineParameters_str2value(self, i, param, kwparams):
      if param.getParameterName() in kwparams:
          paramValue = kwparams[param.getParameterName()]
      else:
          paramValue = kwparams[i]

      #Input params: Tranform STRING to NUMERIC
      cond1 = (str(param) == "Numerical Value")
      cond2 = (str(param) == "Selection")
      isstr = isinstance(paramValue, str)
      cond4 = (str(param) == "Boolean")
      cond5 = (str(param) == "Table Field")
      if "Multi" in (str(param)):
          cond6 = True
      else:
          cond6 = False
      cond7 = (str(param) == "Point")

      if isstr:
          if cond1:
              paramValue = float(paramValue)
          elif cond2:
              paramValue = int(paramValue)
          elif cond4:
              paramValue = eval(paramValue.capitalize())
          elif cond5:
              paramValue = int(paramValue)
          elif cond6:
              #  <type 'java.util.ArrayList'>
              paramValue = list(paramValue)
          elif cond7: #Point
              paramValue = paramValue.split(',')
              x = float(paramValue[0])
              y = float(paramValue[1])
              paramValue = java.awt.geom.Point2D.Double(x, y)
          else: #is str
              pass
      else:
          if cond7: #punto de gvsig. agregar condicion punto
              x = paramValue.getX()
              y = paramValue.getY()
              paramValue = java.awt.geom.Point2D.Double(x, y)
      if DEV_INFO:
          print "Param Name: ", param.getParameterName()
          print "Param: ", param
          print "Resulta: ", paramValue, type(paramValue)
      return paramValue

  def __defineParameters_vector2sextante(self, param, paramValue):
              if isinstance(paramValue, str):
                  layer = gvsig.currentView().getLayer(paramValue)
                  paramValue = self.__createSextanteLayer(layer())
              else:
                  paramValue = self.__createSextanteLayer(paramValue())

              if DEV_INFO: print "extent add+"
              if self.__defaultAE == None:
                  self.__defaultAE = AnalysisExtent(paramValue)
              else:
                  self.__defaultAE.addExtent(AnalysisExtent(paramValue))
              if DEV_INFO: print "PARAM VALUE 1: ", paramValue
              return paramValue

  def __defineParameters_multi2sextante(self, param, paramValue):
              if DEV_INFO: print "PARAM VALUE PRE", paramValue
              paramValue2 = []
              for i in paramValue:
                  if isinstance(i, str):
                      layer = gvsig.currentView().getLayer(i)
                      i = self.__createSextanteLayer(layer())
                      ii = self.__createSextanteLayer(layer())

                  else:
                      if DEV_INFO: print type(i)
                      if isinstance(i, tuple): #si dentro de la lista viene una tupla identificamos que es un raster y su banda
                          from es.unex.sextante.parameters import RasterLayerAndBand
                          rlb = RasterLayerAndBand(self.__createSextanteRaster(i[0]),i[1])  #puesta numero de banda

                          ii = self.__createSextanteRaster(i[0])
                          i = rlb
                          if DEV_INFO: print "*************************** RASTERLAYERBAND", rlb

                      else:
                          if str(type(i)) == "<type 'org.gvsig.raster.fmap.layers.DefaultFLyrRaster'>":
                              ii = self.__createSextanteRaster(i)
                              i = self.__createSextanteRaster(i)
                          else:
                              ii = self.__createSextanteLayer(i())
                              i = self.__createSextanteLayer(i())

                      if DEV_INFO: print "createStextenatnte Layer i :", i
                      paramValue2.append(i)

                  if self.__defaultAE == None:
                          try:
                              self.__defaultAE = AnalysisExtent(AnalysisExtent(ii))
                              self.__defaultAE.setCellSize(0.49) ### test
                              self.__defaultAE.setCellSizeZ(0.49)  ### test
                              #self.__defaultAE.addExtent( )
                              if DEV_INFO: print "----------- ini defaultAE"
                          except:
                              try:
                                  self.__defaultAE = AnalysisExtent(ii())
                              except:
                                  self.__defaultAE = AnalysisExtent(ii)
                  elif self.__defaultAE != None:
                          try:
                              ae = self.__defaultAE
                              #ae.addExtent(AnalysisExtent(ii))
                              if DEV_INFO: print "------------ diferente defaultAE"
                          except:

                              self.__defaultAE.addExtent(AnalysisExtent(ii()))

                  """
                  if self.__defaultAE == None:
                      try:
                          self.__defaultAE = AnalysisExtent(i)
                      except:
                          print "excepttttttttttttttt", i
                          ae = i.getRasterLayer()
                          self.__defaultAE = AnalysisExtent(ae)
                  else:
                      try:

                          self.__defaultAE.addExtent(AnalysisExtent(i))
                      except:
                          print 'exceptttttttttttttt 222222222"', i
                          ae = i.getRasterLayer()
                          print 'ae',ae
                          self.__defaultAE.addExtent(AnalysisExtent(ae))
                          #print self.__defaultAE
                  """

              paramValue = paramValue2
              from java.util import ArrayList
              if DEV_INFO:
                  print "PARAM VALUE 2: ", paramValue
                  print "****************** CONDI6: ", paramValue
                  print "****************** CONDI6: ", ArrayList(paramValue)
              newParamValue = ArrayList()
              for i in paramValue:
                  newParamValue.add(i)
              paramValue = newParamValue
              #print "***************** CONDI6: ", paramValue
              #from es.unex.sextante.parameters import ParameterMultipleInput
              #pmi = ParameterMultipleInput()
              #pmi.setParameterValue(paramValue)
              #paramValue = pmi
              if DEV_INFO: print "#####################################################################"
              return paramValue

  def __defineParameters_raster2sextante(self, param, paramValue):
              if isinstance(paramValue, str):
                  layer = gvsig.currentView().getLayer(paramValue)
                  paramValue = self.__createSextanteRaster(layer)
              else:
                  paramValue = self.__createSextanteRaster(paramValue)

              if DEV_INFO: print "** extent add+"
              if self.__defaultAE == None:
                  self.__defaultAE = AnalysisExtent(paramValue)
              else:
                  self.__defaultAE.addExtent(AnalysisExtent(paramValue))
              return paramValue

  def __defineParameters_table2sextante(self, param, paramValue):
              if isinstance(paramValue, str):
                  layer = gvsig.currentProject().getTable(paramValue)
                  paramValue = self.__createSextanteTable(layer())
              else:
                  paramValue = self.__createSextanteTable(paramValue())
              return paramValue

  def __defineParameters(self, algorithm, kwparams):
      """ Define input parameters """
      params = algorithm.getParameters()
      if DEV_INFO: print self.__defaultAE

      for i in xrange(0,params.getNumberOfParameters()):
          param = params.getParameter(i)
          if DEV_INFO: print "****************************************** PARAMETER *******************************************"
          #print "********************* ", dir(param)
          #print "********************* ", type(param)
          #print "********************* ", param.getParameterName()
          #print "********************* ", param.getParameterDescription()
          #print "********************* ", param.getParameterTypeName()
          #print "********************* ", param.getParameterTooltip()
          #print "********************* ", param.getParameterClass()
          #print "********************* ", param.getClass()
          #print "********************* ", param.parameterName
          #print "********************* ", param.parameterTypeName
          paramValue = self.__defineParameters_str2value( i, param, kwparams)
          #print type(paramValue)

          #Vector to SEXTANTE
          if DEV_INFO: print  "PARAMETER TYPE: ",param.getParameterTypeName()
          if param.getParameterTypeName() == "Vector Layer":
              paramValue = self.__defineParameters_vector2sextante(param, paramValue)

          # Multiple input: Vector and Raster
          elif param.getParameterTypeName() == "Multiple Input": #para vectores y raster
              paramValue = self.__defineParameters_multi2sextante(param, paramValue)
          #Raster to SEXTANTE
          elif param.getParameterTypeName() == "Raster Layer":
              paramValue = self.__defineParameters_raster2sextante(param, paramValue)
          #Table to SEXTANTE
          elif param.getParameterTypeName() == "Table":
              paramValue = self.__defineParameters_table2sextante(param, paramValue)
          #Set parameter value
          #print  "@@@@@@@@@@@@@@@@@@@@@@@@ cluster error: ", paramValue, type(paramValue)
          #print "@@@@@@@@@@@@@@@@@@@@@@@@ ", type(param)

          param.setParameterValue(paramValue)
          #print "@@@@@@@@@@@@@@", dir(param)

          #print "@@@@@@@@@@@@@@", param.isParameterValueCorrect()
          #print "@@@@@@@@@@@@@@ param:", param
          #print "@@@@@@@@@@@@@@ ", param.getParameterValueAsArrayList()


  def __defineExtent(self, algorithm, kwparams):
      """ Define Analysis Extent """
      if self.__defaultAE == None:
          if DEV_INFO: print "-------- cambiar extent"
          change=True
      else:
          if DEV_INFO: print "-------- no cambiar extent"
          change = False
          AExtent = self.__defaultAE
      changeCell = False
      if 'EXTENT' in kwparams.keys() and algorithm.getUserCanDefineAnalysisExtent() :
          changeCell = True
          frame = kwparams['EXTENT']
          if isinstance(frame, str): frame = gvsig.currentView().getLayer(frame)
          #print ("|"+str(frame)+"||"+str(type(frame)))
          if isinstance(frame, str) or isinstance(frame, gvsig.View):
              AExtent = AnalysisExtent()
              if DEV_INFO: print "| EXTENT from VIEW"
              if isinstance(frame, gvsig.View): view = frame
              else: view = gvsig.currentProject().getView(frame)
              envelope = view.getMap().getFullEnvelope()
              xlow = envelope.getLowerCorner().getX()
              ylow = envelope.getLowerCorner().getY()
              xup = envelope.getUpperCorner().getX()
              yup = envelope.getUpperCorner().getY()
              AExtent.setXRange(xlow, xup, False)
              AExtent.setYRange(ylow, yup, False)
              AExtent.setZRange(0, 0, False)
          elif isinstance(frame, DefaultFLyrRaster):
              if DEV_INFO: print "| EXTENT from RASTER"
              layer = self.__createSextanteRaster(frame)
              AExtent = AnalysisExtent(layer)
              changeCell = False
              #AExtent.setZRange(0, 1, True)
              if DEV_INFO: print AExtent
          elif isinstance(frame, list):
              if DEV_INFO: print "| EXTENT from LIST"
              AExtent = AnalysisExtent()
              AExtent.setCellSize(0.49)
              AExtent.setCellSizeZ(0.49)
              xlow, ylow, zlow, xup, yup, zup  = frame[0], frame[1], frame[2], frame[3], frame[4], frame[5]
              AExtent.setXRange(xlow, xup, True)
              AExtent.setYRange(ylow, yup, True)
              AExtent.setZRange(zlow, zup, True)
          elif isinstance(frame, gvsig.Layer):
              if DEV_INFO: print "| EXTENT from Layer"
              layer = self.__createSextanteLayer(frame())
              AExtent = AnalysisExtent(layer)
          else:
              raise Exception("Not Extent Define")
          algorithm.setAnalysisExtent(AExtent)

      elif change == True:
          if DEV_INFO: print ("| Not Extent: No input data")
          AExtent = AnalysisExtent()
          changeCell = True
          if self.__defaultAE != None:
              if DEV_INFO: print "| Extent from Algorithm Layers"
              AExtent = self.__defaultAE
          else:
              try:
                  if DEV_INFO: print "| Extent from View"
                  envelope = gvsig.currentView().getMap().getFullEnvelope()
              except:
                  raise Exception("None open View")

              if DEV_INFO: print "| Setting AExtent: ",
              try: #from view
                  xlow = envelope.getLowerCorner().getX()
                  ylow = envelope.getLowerCorner().getY()
                  zlow = 0
                  xup = envelope.getUpperCorner().getX()
                  yup = envelope.getUpperCorner().getY()
                  zup = 0
                  if DEV_INFO: print "| View: ",
                  if DEV_INFO: print xlow, ylow, xup,yup
              except: # default
                  xlow, ylow, zlow, xup, yup, zup = 0,0,0,100,100,0
                  if DEV_INFO: print "| Default:", xlow, ylow, xup, yup

              frame = Rectangle2D.Double(xlow, ylow, xup, yup)
              AExtent.setXRange(xlow, xup, False)
              AExtent.setYRange(ylow, yup, False)
              AExtent.setZRange(zlow, zup, False)
          algorithm.setAnalysisExtent(AExtent)

      #Set: cellsize

      if 'CELLSIZE' in kwparams.keys():
          AExtent.setCellSize(kwparams['CELLSIZE'])
          if DEV_INFO: print "| New Cellsize: ", kwparams['CELLSIZE'], AExtent.getCellSize()
      elif changeCell == True:
          AExtent.setCellSize(AExtent.getCellSize())
          if DEV_INFO: print "| Cellsize: ", AExtent.getCellSize()

      if 'CELLSIZEZ' in kwparams.keys():
          AExtent.setCellSizeZ(kwparams['CELLSIZEZ'])
          if DEV_INFO: print "| New Cellsize Z: ", kwparams['CELLSIZEZ'], AExtent.getCellSizeZ()
      elif changeCell == True:
          AExtent.setCellSize(AExtent.getCellSizeZ())
          if DEV_INFO: print "| CellsizeZ: ", AExtent.getCellSizeZ()

      if DEV_INFO: print "Å‚Å‚ end extent: ", AExtent
      #AExtent.setZRange(0.0, 0.0, True)
      algorithm.setAnalysisExtent(AExtent)
      if DEV_INFO: print ("| Set Extent")

  def __defineOutput(self, algorithm, kwparams):
    if 'PATH' in kwparams.keys():
        path = kwparams['PATH']
        checkFilesExist(path)
        outputSet = algorithm.getOutputObjects()
        if outputSet.getOutputDataObjectsCount() == 1:
            out1 = outputSet.getOutput(0)
            out1.setOutputChannel(FileOutputChannel("New_name"))
            out1channel = out1.getOutputChannel()
            if isinstance(path, str):
                out1channel.setFilename(path)
            elif isinstance(path, list):
                out1channel.setFilename(path[0])
            else:
                raise Exception("No valid path")
            if DEV_INFO: print "| PATH: Good path"
        elif outputSet.getOutputDataObjectsCount() > 1 and isinstance(path, list):
            for n in xrange(0, outputSet.getOutputDataObjectsCount()):
                out1 = outputSet.getOutput(n)
                out1.setOutputChannel(FileOutputChannel("New_name"))
                out1channel = out1.getOutputChannel()
                out1channel.setFilename(path[n])
                if DEV_INFO: print "| PATH: Good path"
        else:
            raise Exception("Bad path")

    elif algorithm.getOutputObjects().getOutput(0).getOutputChannel():
        output0 = algorithm.getOutputObjects().getOutput(0)
        out0 = output0.getOutputChannel()
        out0.setFilename(None)
        if DEV_INFO: print "| PATH: Without path"

  def __executeAlgorithm(self, algorithm):
    """Execute Algorithm"""
    #Check algorithm
    correctValues = algorithm.hasCorrectParameterValues()
    if DEV_INFO: print "| Parameter values:", correctValues
    if not correctValues: raise Exception("Not correct values")
    try:
        if DEV_INFO: print "| Pre-algorithm:", list(algorithm.algorithmAsCommandLineSentences)
    except:
        if DEV_INFO: print "| Not - algorithm"

    algorithm.execute( None, self.__outputFactory)
    print "| Algorithm:", list(algorithm.algorithmAsCommandLineSentences)

  def __getOutputObjects(self, algorithm):
    """Take outputObjets of the algorithm"""
    oos = algorithm.getOutputObjects()
    ret = dict()
    for i in xrange(0,oos.getOutputObjectsCount()):
      oo = oos.getOutput(i)
      value = oo.getOutputObject()
      if isinstance(value, FlyrVectIVectorLayer):
          if DEV_INFO: print "| Vector"
          store = value.getFeatureStore()
          layer = MapContextLocator.getMapContextManager().createLayer(value.getName(),store)
          store.dispose()
          ret[value.getName()] = layer
          layer.dispose()
      elif isinstance(value, IRasterLayer):
          if DEV_INFO: print "| Raster layer"
          dalManager = gvsig.DALLocator.getDataManager()
          mapContextManager = gvsig.MapContextLocator.getMapContextManager()
          params = dalManager.createStoreParameters("Gdal Store")
          params.setFile(File(value.getFilename()))
          dataStore = dalManager.createStore(params)
          layer = mapContextManager.createLayer(value.getName(), dataStore)
          ret[value.getName()] = layer
      else:
          try:
              ret[value.getName()] = value
          except:
              if not value == None:
                  x = 0
                  while True:
                      field = "value_" + str(x)
                      if any(field in s for s in ret.keys()):
                          x += 1
                      ret[str(x)] = value
                      break
    if DEV_INFO: print "@@@@@@@@ ret:", ret
    return ret

  def __returnOutputObjects(self, r, kwparams):
      if r == None: return
      outList = []
      print "| Output layers: "
      for value in r.values():
        if isinstance(value, unicode):
            outList.append(value.encode("UTF-8"))
        elif isinstance(value, FLyrVect):
            print "|\t Value:", value.getName().encode("UTF-8")
            path = value.getDataStore().getFullName()
            print "|\t\tPath: ", path
            if "OUTVIEW" in kwparams: # other view
                viewName = (kwparams["OUTVIEW"]).decode("UTF-8")
                if "TOCNAME" in kwparams:
                    value = loadShapeFileNew(str(path), view=viewName, tocName=kwparams["TOCNAME"])
                else:
                    value = loadShapeFileNew(str(path), view=viewName)
            else: #current view
                if "TOCNAME" in kwparams:
                    value = loadShapeFileNew(str(path), tocName=kwparams["TOCNAME"])
                else:
                    value = loadShapeFileNew(str(path))


            outList.append(value)
        elif isinstance(value,FLayer):
            print "|\t Value:", value.getName()
            print "|\t\t", value.getFileName()[0]
            #Not yet: Waiting for new loadRasterLayer that can set OUTVIEW
            if "TOCNAME" in kwparams:
                value = loadRasterLayerNew(value.getFileName()[0], tocName=kwparams["TOCNAME"])
            else:
                value = loadRasterLayerNew(value.getFileName()[0])
                
            outList.append(value)
        else:
            print "|\t Non-type"
            print "\tValue: ", value
      print "\n"

      #Return object or list
      if len(r.values()) > 1: return outList
      elif len(r.values()) == 1: return value
      else: return None


  def execute(self, algorithmId, kwparams):
    print "| Algoritmo: ", algorithmId
    algorithm = self.getAlgorithms()[algorithmId]
    #print algorithm.getAlgorithmAsCommandLineSentences()

    #Input params
    self.__defineParameters(algorithm, kwparams)

    #Analisys Extension
    self.__defineExtent(algorithm, kwparams)

    #Output files
    self.__defineOutput(algorithm, kwparams)


    #***Cambiar nombre NO el de las capas
    """
    print "** Define name"
    if 'PATH' in kwparams.keys():
        print "| Nombre: Capa_resultados"
        out0 = java.util.HashMap()
        out0["RESULT"] = "Capa_resultados"
    """

    #Exec algorithm
    self.__executeAlgorithm(algorithm)

    #Output objects
    ret = self.__getOutputObjects(algorithm)
    r = self.__returnOutputObjects(ret, kwparams)

    return r

def unionParameters(params, kwparams):
    for i in range(0,len(params)):
      kwparams[i]=params[i]
    return kwparams

def runalg(algorithmId,*params, **kwparams):
    kwparams = unionParameters(params, kwparams)
    geoprocess = Geoprocess()
    r = geoprocess.execute(algorithmId, kwparams)
    del(geoprocess)
    return r


def loadShapeFileNew(shpFile, CRS='CRS:84', tocName=False, active=False, view=gvsig.currentView()):
    try:
        CRS = gvsig.currentView().getProjectionCode()
    except:
        pass
    layer = loadLayerNew('Shape', shpFile=shpFile, CRS=CRS, tocName=tocName, )
    if isinstance(view,str):
        view = gvsig.currentProject().getView(view)
    else:
        view= gvsig.currentView()
    view.addLayer(layer)
    layer.setActive(active)
    return gvsig.Layer(layer)
    
def loadLayerNew(layerType, tocName=False, **parameters): #
        try:
            application = ApplicationLocator.getManager()
            datamanager =  application.getDataManager()
            mapcontextmanager = application.getMapContextManager()
            store_parameters = datamanager.createStoreParameters(layerType)
            gvsig.copyToDynObject(parameters, store_parameters)
            store = datamanager.openStore(layerType, store_parameters)       
            if tocName!=False:
                nameStore = tocName
            else:
                nameStore = store.getName()
            layer = mapcontextmanager.createLayer(nameStore, store)
        except Throwable, ex:
            raise RuntimeException("Can't load layer, "+ str(ex))  

        return layer

def loadRasterLayerNew(rasterfile, mode = "r", tocName=False ):
    ## Load a Raster file in a Layer
    gvsig_raster.sourceFileName[0]=rasterfile
    if not isinstance (rasterfile,File):
        rasterfile = File(rasterfile)
    
    name, ext = gvsig_raster.splitext(rasterfile.getName())
        
    if tocName!=False:
        name = tocName
    else:
        name = name


    view = gvsig.currentView()
    
    # Get the manager to use
    dalManager = DALLocator.getDataManager()
    mapContextManager = MapContextLocator.getMapContextManager()

    if ext.lower() == ".ecw" or ext.lower() == ".jp2" :
        # FIXME
        pass
    elif ext.lower() == ".mrsid":
        # FIXME
        pass
    else:
        # Create the parameters to open the raster store based in GDAL
        params = dalManager.createStoreParameters("Gdal Store")
        params.setFile(rasterfile)

    # Create the raster store
    dataStore = dalManager.createStore(params)

    # Create a raster layer based in this raster store
    layer = mapContextManager.createLayer(name, dataStore);

    view.addLayer(layer)
    return layer
    
def algHelp(geoalgorithmId):
    geoprocess = Geoprocess()
    for algorithmId, algorithm in geoprocess.getAlgorithms().items():
        if algorithmId.encode('UTF-8') == geoalgorithmId.encode('UTF-8') or geoalgorithmId == "All": 
            pass
        else: 
            continue
        print "* Algorithm help: ", algorithm.getName().encode('UTF-8')
        print "*", algorithm.commandLineHelp.encode('UTF-8')
    del(geoprocess)

def algSearch(strSearch):
    print "\nInicio de busqueda.."
    geoprocess = Geoprocess()
    search = strSearch.lower().encode('ASCII','ignore')
    for algorithmId, algorithm in geoprocess.getAlgorithms().items():
        name = (algorithm.getName()).lower().encode('ASCII','ignore')
        group = (algorithm.getGroup()).lower().encode('ASCII','ignore')
        con1 = str(name).find(search) >= 0
        con2 = str(group).find(search) >= 0
        con3 = algorithmId.encode('ASCII').find(search) >= 0
        if con1 or con2 or con3:
            if con1 or con2:
                print "ID: ", algorithmId, " || GROUP: ", algorithm.getGroup().encode('UTF-8'), " || NAME: ", algorithm.getName().encode('UTF-8')
            else:
                print "*", algorithm.commandLineHelp.encode('UTF-8')
    print "..Busqueda finalizada\n"
    del(geoprocess)

def currentRaster():
    """Return first raster active layer on the View"""
    layers = gvsig.currentView().getLayers()
    lyrlist = [ layers[i] for i in range(0,layers.__len__())]
    for i in lyrlist:
      if i.isActive() and isinstance(i, DefaultFLyrRaster):  return i
    return None

def firstRaster():
    """Return first raster active layer on the View"""
    layers = gvsig.currentView().getLayers()
    lyrlist = [ layers[i] for i in range(0,layers.__len__())]
    for i in lyrlist:
      if isinstance(i, DefaultFLyrRaster):  return i
    return None

def sRaster(n):
    """Return first raster active layer on the View"""
    layers = gvsig.currentView().getLayers()
    lyrlist = [ layers[i] for i in range(0,layers.__len__())]
    count = 0
    for i in lyrlist:
      if isinstance(i, DefaultFLyrRaster):
          if count == n:
              return i
          else:
             count+=1
    return None
def currentActive():
    """Return first active layer on the View"""
    layers = gvsig.currentView().getLayers()
    lyrlist = [ layers[i] for i in range(0, layers.__len__())]
    for i in lyrlist:
        if i.isActive(): return i
    return None

def getProjectLayer(view,layer):
    """Get vector layer or raster"""

    try:
        if isinstance(view, str): view = gvsig.currentProject().getView(view)
        if isinstance(layer, str):
            return view.getLayer(layer)
        else:
            return layer
    except: #Bug: Raster problem with getLayer
        for i in gvsig.currentProject().getView(view).getLayers():
            if i.name == layer: return i

def getProjectTable(table):
    """Get table"""
    return gvsig.currentProject().getTable(table)

def checkFilesExist(files):
    """Path or Paths of files"""
    #raise a exception

    if isinstance(files, str):
        if os.path.isfile(files): raise Exception("File already exist" + files)
    elif isinstance(files, list):
        for fname in files:
            if os.path.isfile(fname): raise Exception("File already exist" + fname)

def main(*args):
    #checkFilesExist(["C:/gvsig/ran10.shp"])

    #geoprocessSearch(" ")
    #geoprocessHelp("closegapsnn")
    #geoprocessHelp("perturbatepointslayer")
    #r = runalg("perturbatepointslayer", LAYER = gvsig.currentLayer(),MEAN = 10, STDDEV = 10 )
    #r = runalg("perturbatepointslayer", EXTENT = "VIEW", LAYER = currentLayer(),MEAN = 10, STDDEV = 10 )
    #r = runalg("perturbatepointslayer", EXTENT = [0,0,500,500], LAYER = currentLayer(), MEAN = 10, STDDEV = 10 )
    #r = runalg("perturbatepointslayer", PATH = "C://gvsig//perturbatepoints028.shp", LAYER = gvsig.currentLayer(),MEAN = 5, STDDEV = 5 )
    #layer = gvsig.currentView().getLayer("data_test_lines.shp")
    #r = runalg("linestoequispacedpoints", LINES=layer,DISTANCE=2)

    #for i in range(10):
    #r = runalg("perturbatepointslayer", LAYER = r[0],PATH = "C://gvsig//perturbatepoints028_" + str(i) + ".shp",MEAN =0.5, STDDEV = 0.5 )
    #r = runalg("fixeddistancebuffer", LAYER = r[0], DISTANCE=1, TYPES="", RINGS=3, NOTROUNDED=False)
    #r = runalg("randomvector", COUNT=20, TYPE=1, EXTENT=gvsig.currentView())
    #RASTER
    #r1 = runalg("generaterandomnormal", EXTENT = [0,0,500,500], PATH = "C://gvsig//perturbatepoints030.tif", MEAN =0.5, STDDEV = 0.5)
    #layer = gvsig.currentView().getLayer("perturbatepoints030")

    #raster = currentRaster()
    #print raster

    #layer = gvsig_raster.loadRasterLayer('c:/gvsig/test_low.tif')
    #r1 = runalg("gradientlines",INPUT = layer, MIN=1, MAX=10, SKIP=1)
    #r = runalg("gridorientation",INPUT=layer,METHOD=0)
    #r = runalg("gridorientation",INPUT=layer,METHOD=0, PATH = "C://gvsig//perturbatepoints010.tif")
    #r = runalg("gridorientation", INPUT=raster, METHOD=0, PATH = "C://gvsig//perturbatepoints011.tif")
    #r = runalg("gradientlines", INPUT = layer, MIN=1, MAX=10, SKIP=1, PATH = "C://gvsig//perturbatepoints012.tif")
    #r = runalg("generaterandomnormal", EXTENT = [0,0,500,500], PATH = "C://gvsig//perturbatepoints013.tif", MEAN =0.5, STDDEV = 0.5)
    #geoprocessHelp("randomvector")
    #r = runalg("randomvector", COUNT=20, TYPE=2, EXTENT=gvsig.currentLayer())
    #r = runalg("randomvector", 200, TYPE_POINT, EXTENT=gvsig.currentLayer(), PATH="C://gvsig//test_puntos_sm01.shp")
    #r = runalg("randomvector", COUNT=20, TYPE=1, EXTENT=gvsig.currentView())
    #r = runalg("randomvector", COUNT=20, TYPE=1, EXTENT="VIEW")
    #r = runalg("randomvector", COUNT=20, TYPE=1, EXTENT=currentRaster())
    #r = runalg("gvSIG-convexhull", LAYER="Puntos_de_interes_01.shp", CHECK=True, PATH = "C://gvsig//gvsigconvexhull_001.shp")
    #r = runalg("generaterandomnormal", PATH = "C://gvsig//per.tif", EXTENT=gvsig.currentLayer(), CELLSIZE = 100, PATH = "C://gvsig//perturbatepoints014.tif", MEAN =5, STDDEV = 5)
    #geoprocessHelp("tablebasicstats")
    #r =runalg("tablebasicstats",TABLE=gvsig.currentTable(), FIELD=0)

    #Without parameters label

    #layers = gvsig.currentView().getLayers()
    ##r = runalg("gridorientation", layer, 0, PATH = "C://gvsig//perturbatepoints012.tif")

    #r = runalg("gridorientation", layer, 0)
    #r = runalg("gradientlines", layers[0], 1, 10, 1, PATH = "C://gvsig//perturbatepoints012.shp")
    #r = runalg("cva", "test_low", "test_low", "test_low", "test_low")
    #r = runalg("cva", currentRaster(), currentRaster(), currentRaster(), currentRaster(),PATH=["C:/gvsig/1.tif","C:/gvsig/2.tif"])

    #layerRaster = gvsig_raster.loadRasterLayer('c:/gvsig/test_low.tif')
    #r = runalg("gridorientation",layerRaster,0, PATH = "C://gvsig//Grid_orientation.tif",EXTENT=layerRaster, CELLSIZE=1, CELLSIZEZ=10)
    #r2 = runalg("cva", r, r, r, r, PATH=["C:/gvsig/1.tif","C:/gvsig/2.tif"])
    #print r2[0], r2[1]
    #print layerRaster
    #print getProjectLayer("Vista1", "test_low")


    #layer = getProjectLayer("Vista1", "as.shp")
    #extent = getProjectLayer("Vista1", "analisis_extent")
    #runalg("difference", "as.shp", "vista2_testeo.shp", PATH="C:/gvsig/recorte_extent.shp", EXTENT=[100, 100, 0, 540, 500, 0])
    #runalg("difference", "vista2_testeo.shp", layer, PATH="C:/gvsig/recorte_extent_2.shp", EXTENT=layer, OUTVIEW="Nueva")
    #r = runalg("tablebasicstats", "species", 0)
    #print r.encode("UTF-8")

    #algHelp("generaterandomnormal")
    """
    r = runalg("generaterandomnormal", 100,100, CELLSIZE=100, EXTENT=[250,250,0,500,500,0])
    r = runalg("generaterandomnormal", 10, 10, CELLSIZE=50, EXTENT=[500,500,0, 1000,1000,0])
    r = runalg("generaterandombernoulli", 50.0, CELLSIZE=25, EXTENT=[1000,1000,0, 1250,1250,0])
    r = runalg("gradientlines", r, 1, 100, 1)

    v1 = runalg("randomvector",10, TYPE_POLYGON, EXTENT=[0,0,0,500,500,0])
    v2 = runalg("randomvector", 5, TYPE_POLYGON, EXTENT=v1)
    v3 = runalg("difference", v1, v2, PATH="C:/gvsig/Diferencia.shp")
    v4 = runalg("randomvector", 5, 0, PATH="C:/gvsig/randomvector.shp", EXTENT=v3)
    v5 = runalg("randomvector", 100, 2, PATH="C:/gvsig/randompoints.shp", EXTENT="randomvector", OUTVIEW="Nueva")
    v6 = runalg("gvSIG-xyshift", "randompoints", "false", "-250.0", "-250.0", PATH=["C:/gvsig/ran1010.shp","C:/gvsig/ran2020.shp","C:/gvsig/ran3030.shp"])
    algHelp("tablebasicstats")
    v7 = runalg("gvSIG-buffer", "randompoints", False, 50.0, 0, False, True, 0, 0, PATH="C:/gvsig/buffer_gvsig0138.shp")
    v5 = runalg("randomvector", 100, 2, EXTENT=[0,0,0,500,500,0])

    #Test for model sextante
    v5 = runalg("randomvector", 10, 2, EXTENT=[0,0,0,500,500,0], OUTVIEW="Nueva")
    v5 = runalg("randomvector", "5", 2, EXTENT=[0,0,0,500,500,0], OUTVIEW="Nueva")

    r = runalg("generaterandomnormal", 100,100, CELLSIZE=1)

    vista= (gvsig.currentView().name).encode('UTF-8')
    """


    """
    vista= "Sin título"
    #vista = vista.encode("UTF-8")

    print gvsig.currentProject().getView(vista)


    endid = "_104.shp"
    v = runalg("randomvector", 20, 0)
    v2 = runalg("randomvector", 20, 0,PATH=["C://gvsig//i"+endid])
    #fitpoints = runalg("fitnpointsinpolygon",v, "2", "74", 1, 0)
    intersection = runalg("gvSIG-intersection", v, getProjectLayer(gvsig.currentView(),"i_102"), True, True,PATH=["C://gvsig//inter_pol"+endid,"C://gvsig//inter_line"+endid,"C://gvsig//inter_point"+endid])
    """
    #algSearch("nodelines")
    #ID:  extractendpointsoflines  || GROUP:  Topolog????a  || NAME:  Obtener extremos de l??neas

    #v7 = runalg("gvSIG-buffer", v, "False", 50.0, 0, "false", "true", 0, 0)
    #print "**Addfield: ", gvpy.addField(v, "ID5")
    #import gvpy
    #gvpy.addField(v, "ID89")
    #runalg("vectoraddfield", v, "Campo1", 0, 1, 0, 1, PATH="C:\\Users\\Oscar\\Desktop\\testi.shp")
    #layer = runalg("randomvector", 20, TYPE_LINE, EXTENT=[0,0,0,500,500,0], OUTVIEW="Po1")

    #2015-1-2


    #runalg("graticulebuilder", "10.0", "10.0", 1, PATH="/home/oscm/temp/reticula10.shp")
    #layer1 = gvsig.currentView().getLayer("1")
    #layer2 = gvsig.currentView().getLayer("2")
    #print layer1, layer2
    #
    #algHelp("cluster")
    #runalg("merge", layer1, [layer2])

    #layer1 = sRaster(0)
    #layer2 = sRaster(1)
    #print "Layer1: ", layer1.name, layer1, layer1.getCellSize()
    #print "Layer2: ", layer2.name, layer2, layer2.getCellSize()


    #runalg("cluster" , INPUT=[(layer1, 0),(layer2,0)] , NUMCLASS=3)#, CELLSIZE=layer1.getCellSize(), CELLSIZEZ=1)#,EXTENT=layer1)
    #runalg("mergegrids", [layer1, layer2], "0")
    v1 = runalg("randomvector",10, TYPE_POLYGON, EXTENT=[0,0,0,500,500,0], TOCNAME="Random vector")
    r = runalg("generaterandomnormal", EXTENT = [0,0,0,500,500,0], MEAN =0.5, STDDEV = 0.5, TOCNAME="Raster Generate Random normal")
    print "end"

def mainLibrary():
    #showFields(gvsig.currentLayer())
    #removeField(layer, "campo3")
    #removeField(layer, "campo5")
    #addField(layer,"campo1")
    #addField(layer,"campo2")
    #renameField(layer, "campo5", "campo50")

    #print list2geompoly([[1,2],[3,10],[5,30]])
    #addFeature(layer, "test01", "test2", [15,200])
    #addFeature(layer, "test11", "test3", geom.createPoint(300, 301))
    #addFeature(layer, campo1 = "kwparam", campo2 = "kwparam2", GEOMETRY = geom.createPoint(300,300))
    #addFeature(layer, "linea", "01", [[1,2],[3,10],[5,30]])
    #addFeature(layer, "pol", "02", [[50,80],[150,50],[100,10],[0,10],[50,80]])

    """
    layer = newLayer()
    addFeature(layer, "pol", "01", [[50,80],[150,50],[100,10],[0,10],[50,80]])
    addFeature(layer, "pol", "02", [[0,0],[10,5],[10,10],[0,10],[5,5]])
    addFeature(layer, "pol", "03", [[-50, -34],[0,0], [-14,30]])
    addField(layer,"campo3")
    modifyFeatures(layer, "campo3", "nuevo poligono")
    """
    """
    #Create shapes
    layer1 = newLayer(layer,"C:/gvsig/point_shape.shp", 1)
    layer2 = newLayer(layer,"C:/gvsig/line_shape.shp", 2)
    layer3 = newLayer(layer,"C:/gvsig/polygon_shape", 3)

    #Add features
    addFeature(layer1, "point", "01", [50,80])
    addFeature(layer1, "point", "02",[150,50])
    addFeature(layer1, "point", "03",[100,10])
    addFeature(layer1, "point", "04",[0,10])

    addFeature(layer3, "polygon", "01", [[50,80],[150,50],[100,10],[0,10],[50,80]])

    addFeature(layer2, "line", "01", [[50,80],[150,50],[100,10],[0,10],[50,80]])

    #Modify all values in one column
    modifyFeatures(layer1, "campo1", "Points_gsoc")
    modifyFeatures(layer2, "campo1", "Lines_gsoc")

    #Modify schema
    addField(layer1,"Name")
    removeField(layer1,"Surname")

    ##Execute SEXTANTE
    #r = geoprocess("perturbatepointslayer", LAYER = currentLayer(),MEAN = 5, STDDEV = 5 )
    """
    #layer = gvsig.currentView().getLayer("line_04.shp")
    #newLayer(layer, "C:/gvsig/gvpy_test006.shp")
    #layer2 = copyLayer(layer, "C:/gvsig/gvpy_copylayer_012.shp")
    #v = copyLayer(layer, "C:/gvpy_copylayer_new_06.shp")
    #layer = gvsig.currentView().getLayer("gvpy_copylayer_new_06")
    #addFeature(v, "Camino", [[50,00],[50,50],[10,10],[0,1],[50,18]])

    #Basics field
    """
    addField(v, "Direccion")
    modifyFeatures(v, "Direccion", "Av")
    removeField(v, "Direccion")
    """

    """
    #EJEMPLO 1
    #New shapes
    layer = gvpy.runalg("randomvector", 20, gvpy.TYPE_LINE, EXTENT=[0,0,0,500,500,0], OUTVIEW="Po1")
    #layer = gvsig.currentLayer()
    #Advanced field
    removeField(layer, "ID")
    removeField(layer, "Distance")
    addField(layer, "ID") #Add fields
    addField(layer, "Distance", "STRING")
    addField(layer, "Long", "LONG")
    removeField(layer, "Long") #Remove field
    modifyFeatures(layer, "ID", "90") #Modify all features
    modifyField(layer, "ID", "LONG") #Modify type of field
    modifyField(layer, "Distance", "FLOAT")
    addFeature(layer, 1, 0, [[50,0],[1000,0]]) #Add new feature with geometry line
    for feature in layer.features():
        perimeter = feature.geometry().perimeter()
        modifyFeature(layer, feature, "Distance", perimeter) #Modify each feature

    pass
    #FIN EJEMPLO 1
    """

    #modifyFeatures(gvsig.currentLayer(), "ID", 100, "Distance == 90")
    layer = gvsig.currentLayer()
    for feature in gvsig.currentLayer().features():
        value = feature.Distance
        #modifyFeature(layer, feature, "ID", distance)
        feature.edit()
        feature.set("ID", value)
        layer.update(feature)
    #model2script("C://gsoc//test02.model", "C://gsoc//nuevoScript.py")
    print "END TEST"

def copyLayerFeatures2Layer(layer1, layer2):
    for i in layer1.features():
        layer2.append(i.getValues())
    layer2.commit()

def copyLayer(layer, path):
    output = newLayer(layer, path)
    copyLayerFeatures2Layer(layer, output)
    #addLayerView(output)
    return output

def newLayer(layer, path, geometryType=None):
    CRS = layer.getProjectionCode()
    schema = gvsig.createSchema(layer.getSchema())
    if geometryType==None: geometryType = layer.getTypeVectorLayer().getType()
    output = gvsig.createShape( schema, path, CRS=CRS, geometryType=geometryType )
    gvsig.currentView().addLayer(output)
    return output

def addLayerView(layer):
    gvsig.currentView().addLayer(layer)

def addFeature(layer, *params, **kwparams):
    #IN: layer, feature params + geometry
    typeLayer = layer.getTypeVectorLayer().name
    #if kwparams != {}:
    #    layer.append(kwparams)
    #    layer.commit()
    #    return

    if "COMMIT" in kwparams:
        pass
    else:
        COMMIT=1
    if params != ():
        schValues = layer.getSchema().getAttrNames()
        values = {}
        itera = iter((list(params)))
        value = itera.next()
        for sch in schValues:
            #Si el campo a modificar es una geometria
            #print "Comprobación:", sch, isinstance(value, list)
            #re comprobacion si es campo geometry
            #bug: Comprobar si es lista o objeto geom  en primer if
            #... sch == "Geometry" and ES UNA LISTA
            #... sino copia el valor directamente: caso de pasar geometrias
            if sch == "GEOMETRY":
                print typeLayer
                if typeLayer == "Point2D":
                    if isinstance(value, list):
                        values[sch] = geom.createPoint(value[0],value[1])
                elif typeLayer == "MultiCurve2D":
                    if isinstance(value, list):
                        values[sch] = list2geomcurve(value)
                elif typeLayer == "MultiSurface2D":
                    if isinstance(value, list):
                        values[sch] = list2geompoly(value)
                else: #Si son geometrias
                    values[sch] = value
            else:
                values[sch] = value
            try:
                value = itera.next()
            except:
                break
        layer.append(values)
        if COMMIT==1: layer.commit()
    print "Add feature ", params, " to ", layer

def list2geompoly(listPoints):
    #IN: list[[x,y],...]
    #OUT: geometry polygon
    geometry = geom.createGeometry(3)
    for point in listPoints:
        geometry.addVertex(geom.createPoint(point[0],point[1]))
    if listPoints[0] != listPoints[len(listPoints)-1]:
        geometry.addVertex(geom.createPoint(listPoints[0][0], listPoints[0][1]))
    return geometry

def list2geomcurve(listPoints):
    #IN: list [[x,y],...]
    #OUT: geometry line
    geometry = geom.createGeometry(2)
    for point in listPoints:
        geometry.addVertex(geom.createPoint(point[0],point[1]))
    return geometry

def modifyFeatures(layer, field, value, COMMIT=1,FILTER=None):
    #IN: layer, field, new value
    if FILTER == None:
        features = layer.features()
    else:
        features = layer.features(FILTER)

    for feature in features:
        feature.edit()
        feature.set(field, value)
        layer.update(feature)
    if COMMIT==1: layer.commit()
    print "Modify feature ", layer.name, field, value

def showFields(layer):
    #IN: layer
    #OUT: str(atributos)
    print layer.getSchema().getAttrNames()

def addField(layer,field, sType = "STRING",iSize=20):
    #IN: layer, field, *sType, *iSize)
    #OUT: layer
    #addField(layer, "nombre")
    schema = gvsig.createSchema(layer.getSchema())
    schema.modify()
    if isinstance(field,str): schema.append(field,sType,iSize)
    layer.edit()
    layer.updateSchema(schema)
    layer.commit()
    print "Add field ", field, " to ", layer.name
    return layer

def modifyField(layer, field, iType="STRING", iSize=20):
    temp = []
    for i in layer.features():
        temp.append(i.get(field))
    removeField(layer, field)
    addField(layer, field, iType, iSize)
    n = 0
    for i in layer.features():
        modifyFeature(layer, i, field, temp[n],COMMIT=0)
        n += 1
    layer.commit()
    print "Modify field type to: ", field, " in ", layer.name

def modifyFeature(layer, feature, field, value, COMMIT=1):
    feature.edit()
    feature.set(field, value)
    layer.update(feature)
    print "Modify Feature field: ", field , " to ", value
    if COMMIT==1: layer.commit()

def removeField(layer, field):
    #IN: layer, field
    #OUT: layer
    #removeField(layer, "apellido")
    print "Campo %s eliminado" % (str(field))
    schema = layer.getSchema()
    schema.modify()
    if isinstance(field,str): schema.remove(field)
    else: return
    layer.edit()
    layer.updateSchema(schema)
    layer.commit()
    return layer


def model2script(pathXML, pathFile):
    #eliminar la ultima linea
    #pathXML = commonsdialog.openFileDialog("Selecciona archivo", 'C:/gsoc/')
    #pathXML = str(pathXML[0])
    #pathXML = 'C:/gsoc/test02.model'
    #pathFILE = 'C:/gsoc/script0002.py'
    import os.path
    if os.path.isfile(pathXML)==True:
        pass
    else:
        print "No valid model file path"
        return
    fileFile = open(pathXML, 'r')
    document = fileFile.read()
    import xml.dom.minidom
    root = xml.dom.minidom.parseString(document)
    #root
    inputObject = {}
    inputObjectParams = {}
    dataObject = {}
    algorithms = {}
    tab = "    "
    gvpyFile = open(pathFile, "w")
    #Cargamos los parametros
    print "\nData object"
    for child in root.getElementsByTagName("data_object"):
        #data_object - Parametros
        if "INNER" in child.getAttribute("key"):
            inputObjectParams[child.getAttribute("key")] = child.getAttribute("value")
        #data_object - result of algorithms
        else:
            inputObject[child.getAttribute("key")]=[child.getAttribute("value"), child.getAttribute("description").encode("UTF-8")]

    print "\n Attribute"
    for child in root.getElementsByTagName("input"):
        for i in child.getElementsByTagName("attribute"):
                if i.getAttribute("name")=="default": dataObject[child.getAttribute("name")] = i.getAttribute("value")

    print "\n Algorithm"
    order = 1
    for child in reversed(root.getElementsByTagName("algorithm")):
        print "Algoritmo: ", child
        keyAlgorithm = child.getAttribute("key")
        algorithmParams = {}
        algorithmParams["alg_cmd_line_name"]=child.getAttribute("alg_cmd_line_name")
        for i in child.getElementsByTagName("assignment"):
            algorithmParams[i.getAttribute("key")] = i.getAttribute("assigned_to")
        algorithmParams["result_of_algorithm"] = keyAlgorithm
        algorithms[order] = algorithmParams
        order +=1

    print "\n\n******* RESULTADO *******"
    print "inputObject: ", inputObject
    print "inputObjectParams: ", inputObjectParams
    print "algorithms: ", algorithms
    print "data object: ", dataObject

    #Writing script .py
    print "\nTransform to gvpy"
    for i in root.getElementsByTagName("model"):
        modelName = i.getAttribute("name")
    gvpyFile.write("# Modelo de SEXTANTE: " + modelName)
    gvpyFile.write(
"""
import gvpy
import gvsig
import geom
import gvsig_raster

def main(*args):
""")

    print "gvpy - data_object"
    listInputObject = []
    for n in reversed(inputObject.keys()):
        listInputObject.append([n,inputObject[n][1]])

    print "gvpy - inputObjectParams"
    for n in (inputObjectParams.keys()):
        gvpyFile.write( tab + n + ' = "' + inputObjectParams[n] + '"\n' )

    print "gvpy - vars"
    for n in (dataObject.keys()):
        gvpyFile.write( tab + n +' = "' + dataObject[n] + '"\n\n' )


    print "gvpy - algorithms"
    #inputObject list of result algorithms names
    for n in reversed(sorted(algorithms.keys())): #reversed(algorithms.keys()):
        gvpy= ""
        alg = algorithms[n]
        #prefijo: buscar en los data_object el nombre que debe de llevar la capa resultado
        for i in listInputObject:
            if alg["result_of_algorithm"] in i[0]:
                prefix = i[0]
                description = i[1]
        #Escribimos el codigo del algoritmo
        gvpyFile.write( tab + '# '+ description + '\n')
        gvpy += prefix + '= gvpy.runalg("'+alg["alg_cmd_line_name"]+'"'
        for i in alg:
            if i == "alg_cmd_line_name" or i == "result_of_algorithm": continue
            gvpy += ', '+i+'='+ alg[i] + ''
        gvpy += ')'
        gvpyFile.write( tab + gvpy + "\n\n" )

    gvpyFile.close()