Application: gvSIG desktop » gvSIG scripting

# encoding: utf-8

import gvsig

from gvsig import currentView

from gvsig import createFeatureType

from gvsig.geom import POINT, D2, createPoint

from gvsig import createShape

import random

def createPointLayer():

  #

  # Vamos a crear los puntos sobre el area representada en la vista.

  # usaremos la funcion random.randint de python, que solo usa enteros

  # y devuelbe enteros, asi que usaremos valores decimales con tres

  # decimanes representandolos como de coma fija en un entero.

  # 

  # 1.123 ===> 1123 

  #

  # Y cuando vayamos a crear los puntos con esos valores los volveremos

  # a transformar en decimales dividienlos por 1000.0

  #

  #

  # Para eso, lo primero obtendremos el Envelope de la vista

  envelope = currentView().getMapContext().getViewPort().getEnvelope()

  #

  # Calculamos los valores maximos y minimos de las coordenadas de

  # ese envelope como valores con tres decimales en coma fija representados

  # en un entero.

  #  

  minx = int(envelope.getMinimum(0) * 1000)

  maxx = int(envelope.getMaximum(0) * 1000)

  miny = int(envelope.getMinimum(1) * 1000)

  maxy = int(envelope.getMaximum(1) * 1000)

  print minx, maxx

  print miny, maxy

  #

  # Creamos un shape con solo un campo geometria de tipo

  # punto 2D

  schema = createFeatureType() # DefaultFeatureType

  schema.append("GEOMETRY", "GEOMETRY")

  schema.get("GEOMETRY").setGeometryType(POINT, D2)

  puntos = createShape(schema, prefixname="puntos")

  # 

  # Ponemos en edicion el nuevo shape vacio

  puntos.edit()

  # Vamos a crear 20 puntos

  for x in range(20):

    #

    # Calculamos la x e y y la dividimos por 1000.0

    # para pasarla a double.

    x = random.randint(minx, maxx)/1000.0

    y = random.randint(miny, maxy)/1000.0

    # Creamos el punto

    point = createPoint(D2,x,y)

    # Lo añadimos a la capa de puntos

    puntos.append( { "GEOMETRY": point } )

  #

  # Cuando hemos añadido todos los puntos terminamos la edicion

  puntos.finishEditing()

  # Por ultimo añadiremos la capa a la vista

  currentView().addLayer(puntos)

  # Y mostraremos por consola la ruta en la que se ha creado la capa

  print puntos.getDataStore().getFullName()

  return puntos

def main(*args):

  createPointLayer()

[script]

enable = true

isolation = 0

[Unit]

type = Script

name = 05_crearPuntosAleatorios

description =

createdBy =

version =

[Script]

main = main

Lang = python

PROJCS["ED50 / UTM zone 30N",GEOGCS["ED50",DATUM["European_Datum_1950",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],TOWGS84[-87,-98,-121,0,0,0,0],AUTHORITY["EPSG","6230"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4230"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","23030"]]

PROJCS["ED50 / UTM zone 30N",GEOGCS["ED50",DATUM["European_Datum_1950",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],TOWGS84[-87,-98,-121,0,0,0,0],AUTHORITY["EPSG","6230"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4230"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","23030"]]

PROJCS["ED50 / UTM zone 30N",GEOGCS["ED50",DATUM["European_Datum_1950",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],TOWGS84[-87,-98,-121,0,0,0,0],AUTHORITY["EPSG","6230"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4230"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","23030"]]