Application: gvSIG desktop » gvSIG scripting

# -*- coding: utf-8 -*-

#

# File: gvsig.py

#

# Copyright (c) 2011 by Model Driven Development sl and Antonio Carrasco Valero

#

# GNU General Public License (GPL)

#

# 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., 51 Franklin Street, Fifth Floor, Boston, MA

# 02110-1301, USA.

#

#

__author__ = """Antonio Carrasco Valero

Model Driven Development sl and Antonio Carrasco Valero

<carrasco@modeldd.org>

Victor Acevedo Royer <vacevedor@gmail.com>

"""

__docformat__ = 'plaintext'

from org.gvsig.fmap.crs import CRSFactory

from org.gvsig.fmap.geom import Geometry, GeometryLocator

from org.gvsig.fmap.geom.primitive import Point

from org.gvsig.tools import ToolsLocator

#GeometryTypes

AGGREGATE = Geometry.TYPES.AGGREGATE

ARC = Geometry.TYPES.ARC

CIRCLE = Geometry.TYPES.CIRCLE

CURVE = Geometry.TYPES.CURVE

ELLIPSE = Geometry.TYPES.ELLIPSE

ELLIPTICARC = Geometry.TYPES.ELLIPTICARC

GEOMETRY = Geometry.TYPES.GEOMETRY

MULTICURVE = Geometry.TYPES.MULTICURVE

MULTIPOINT = Geometry.TYPES.MULTIPOINT

MULTISOLID = Geometry.TYPES.MULTISOLID

MULTISURFACE = Geometry.TYPES.MULTISURFACE

NULL = Geometry.TYPES.NULL

POINT = Geometry.TYPES.POINT

SOLID =  Geometry.TYPES.SOLID

SPLINE = Geometry.TYPES.SPLINE

SURFACE = Geometry.TYPES.SURFACE

# Common named geometry types

POLIGON = Geometry.TYPES.SURFACE

LINE = Geometry.TYPES.CURVE

MULTILINE = Geometry.TYPES.MULTICURVE

MULTIPOLIGON = Geometry.TYPES.MULTISURFACE

# geometrySubTypes 

D2= Geometry.SUBTYPES.GEOM2D

DM2= Geometry.SUBTYPES.GEOM2DM

D3= Geometry.SUBTYPES.GEOM3D

DM3= Geometry.SUBTYPES.GEOM3DM

UNKNOWN= Geometry.SUBTYPES.UNKNOWN

def createGeometry(type, subtype=D2):

  """

  Create a new geometry with a concrete type and subtype 

    :param type: geometry type

    :type type: string

    :param subtype: geometry subtype

    :type type: string

    :return: geometry

    :rtype: geometry

  """

  #type = getGeometryType(type)

  #subtype = getGeometrySubType(subtype)

  geometryManager = GeometryLocator.getGeometryManager()

  geometry = geometryManager.create(type, subtype)

  return geometry

def createPoint(x=0, y=0, subtype=D2):

  """

  Create a new point with a subtype and sets the value for the X and the Y

    :param x: X coordinate value 

    :param y: Y coordinate value

    :type x: double

    :type y: double

    :return: Point

    :rtype: Point  

  """

  geometryManager = GeometryLocator.getGeometryManager()

  point = geometryManager.createPoint(x, y, subtype)

  return point

def createMultiPoint(points=None, subtype=D2):

  """

  Create a new multipoint with a subtype from a list of Points. Also values of 

  X and Y tuples like ([x1, y1], [x2, y2], ...., [xn, yn]) are allowed

    :param points: list of points

    :type points: list

    :return: multipoint

    :rtype: multipoint

  """

  multipoint = createGeometry(MULTIPOINT, subtype)

  if points == None:

    return multipoint

  for point in points: 

    if not isinstance(point, Point):

      point = createPoint(point[0], point[1], subtype)

    multipoint.add(point)

  return multipoint 

def createEnvelope(pointMax=None, pointMin=None, dimension=2):

  """

  Create envelope as a minimum bounding box or rectangle. This envelope is 

  equivalent to the GM_Envelope specified in ISO 19107. 

    :param pointMax: 

    :type pointMax: geometry POINT

    :param pointMin: 

    :type pointMin: geometry POINT

  """

  geometryManager = GeometryLocator.getGeometryManager()

  if None in (pointMax, pointMin):

     envelope = geometryManager.createEnvelope()    

  else:

     envelope = geometryManager.createEnvelope(pointMax, pointMin)

  return envelope

def getCRS(crs):

  return CRSFactory.getCRS(crs)

def main():

  point = createPoint(1,1)