def main(*args):
# Test pq con original phrase no funciona
store = gvsig.currentLayer().getFeatureStore()
dm = DALLocator.getDataManager()
filterExpression = ExpressionEvaluatorLocator.getManager().createExpression()
filterExpression.setPhrase("")
spatialGeometry = geom.createGeometryFromWKT("POLYGON ((-74.26662451007567 40.497289998, -74.26662451007567 40.517289997999995, -74.24930400199999 40.527289998, -74.2319834939243 40.517289997999995, -74.2319834939243 40.497289998, -74.24930400199999 40.487289997999994, -74.26662451007567 40.497289998))")
spatialGeometry.setProjection(gvsig.currentLayer().getProjection())
## store, exp, spatialGeometry
print "Spatial"
builder = ExpressionEvaluatorLocator.getManager().createExpressionBuilder()
"""
expr1 = builder.and(
builder.custom(exp.getPhrase()),
builder.ST_Intersects(
builder.geometry(spatialGeometry),
builder.column("GEOMETRY")
)
).toString()
"""
expr1 = builder.ST_Intersects(
builder.geometry(spatialGeometry),
builder.column("GEOMETRY")
).toString()
## Creating new expression
eva = ExpressionEvaluatorLocator.getManager().createExpression()
eva.setPhrase(expr1)
evaluator = DALLocator.getDataManager().createExpresion(eva)
# Query
fq = store.createFeatureQuery()
fq.addFilter(evaluator)
fq.retrievesAllAttributes()
fset = store.getFeatureSet(fq)
print "SIZE: ",fset.getSize()
### Reusing expression
cloneExpression = filterExpression.clone()
print "not cloning"
cloneExpression.setPhrase(expr1)
evaluator = DALLocator.getDataManager().createExpresion(cloneExpression)
## Query
fq = store.createFeatureQuery()
fq.addFilter(evaluator)
fq.retrievesAllAttributes()
fset = store.getFeatureSet(fq)
print "SIZE: ",fset.getSize()
def main(*args):
i18n = ToolsLocator.getI18nManager()
layer = currentLayer()
if layer == None:
msgbox(i18n.getTranslation("_Need_a_vector_layer_selected_in_the_ToC"))
return
panel = AddTimeSupport(currentLayer())
winmgr = ToolsSwingLocator.getWindowManager()
dialog = winmgr.createDialog(
panel.asJComponent(), i18n.getTranslation("_Add_basic_time_support"),
i18n.getTranslation("_Add_basic_time_support_to_vectorial_layer"),
winmgr.BUTTONS_OK_CANCEL)
dialog.addActionListener(panel)
dialog.show(winmgr.MODE.WINDOW)
def main(*args):
###
### INFO
###
# DATA HAS FIELD NAMED CMPLNT_F_1 with HH:mm:ss format
# CONVERT FIELD TO DATE
nameFieldHour = "CMPLNT_F_1"
nameFieldDay = "CMPLNT_FR_"
# Transformar de antemano a fecha
patternHour = 'HH:mm:ss'
patternDay = 'yyyy-MM-dd'
rangeHoursParameter = "0-10,13"
rangeDaysParameter = "0, 3-6,7"
expression = ''
expression = "OFNS_DESC = 'BURGLARY'"
xi = 0
yi = 0
proportion = 1
#store = gvsig.currentView().getLayer("20f").getFeatureStore()
store = gvsig.currentLayer().getFeatureStore()
aoristicClock(store,
nameFieldHour,
nameFieldDay,
patternHour,
patternDay,
rangeHoursParameter,
rangeDaysParameter,
expression,
xi,
yi,
proportion)
def getAllValidDocuments(self, documentException=None):
application = ApplicationLocator.getManager()
project = application.getCurrentProject()
mlayer = gvsig.currentLayer()
docs = project.getDocuments()
all = []
try:
for doc in docs:
#print doc, type(doc)
if isinstance(doc, ViewDocument):
for layer in doc.getMapContext().getLayers():
#print "--", layer==mlayer, layer.getName()
if isinstance(layer, FLyrVect):
if documentException != layer:
all.append(layer)
elif isinstance(doc, TableDocument):
#print "--", doc
if doc.getAssociatedLayer() == None:
if documentException != doc:
all.append(doc)
except:
ex = sys.exc_info()[1]
print "Error", ex.__class__.__name__, str(ex)
return all
def main(*args):
fset = gvsig.currentLayer().getFeatureStore().getFeatureSet()
field = "resblob"
formatName = "PNG" #https://docs.oracle.com/javase/tutorial/2d/images/saveimage.html
ft = gvsig.currentLayer().getFeatureStore().getDefaultFeatureType(
).getAttributeDescriptor("id")
outparams = {
"option": 2,
"params": {
"path": File("/home/osc/Working/exportimages/"),
"field": ft
}
}
exp = ExpressionEvaluatorLocator.getManager().createExpression()
exp.setPhrase("CONCAT('/home/osc/Working/exportimages/', id)")
outparams = {"option": 3, "params": {"expression": exp}}
processExportImage(fset, field, formatName, outparams)
def execute(self,actionCommand, *args):
actionCommand = actionCommand.lower()
if actionCommand == "view-time-add-basic-support-to-datasource":
i18n = ToolsLocator.getI18nManager()
layer = currentLayer()
if layer == None:
msgbox(i18n.getTranslation("_Need_a_vector_layer_selected_in_the_ToC"))
return
panel = AddTimeSupport(currentLayer())
winmgr = ToolsSwingLocator.getWindowManager()
dialog = winmgr.createDialog(
panel.asJComponent(),
i18n.getTranslation("_Add_basic_time_support"),
i18n.getTranslation("_Add_basic_time_support_to_vectorial_layer"),
winmgr.BUTTONS_OK_CANCEL
)
dialog.addActionListener(panel)
dialog.show(winmgr.MODE.WINDOW)
def main(*args):
currentView = gvsig.currentView()
currentLayer = gvsig.currentLayer()
area = calcularArea(currentView, currentLayer)
if area != 0:
message = str(
float(area)) + " Km2" # https://es.wikipedia.org/wiki/Uruguay
title = unichr(
193) + "rea" # https://es.wikipedia.org/wiki/ISO/IEC_8859-1
mensajeCartel(message, title, commonsdialog.IDEA)
def execute(self, actionCommand, *args):
actionCommand = actionCommand.lower()
if actionCommand == "settool-reportbypoint":
#print "### QuickinfoExtension.execute(%s)" % repr(actionCommand)
layer = currentLayer()
if not self.isLayerValid(layer):
return
viewPanel = currentView().getWindowOfView()
mapControl = viewPanel.getMapControl()
reportbypoint = ReportByPoint()
reportbypoint.setTool(mapControl)
def execute(self, actionCommand, *args):
actionCommand = actionCommand.lower()
if actionCommand == "settool-quickmapexport":
#print "### QuickinfoExtension.execute(%s)" % repr(actionCommand)
layer = currentLayer()
if not self.isLayerValid(layer):
return
viewPanel = currentView().getWindowOfView()
mapControl = viewPanel.getMapControl()
quickmapexport = QuickMapExport()
quickmapexport.showTool("_Quick_map_export")
def execute(self, actionCommand, *args):
actionCommand = actionCommand.lower()
if actionCommand == "settool-exportmaptiles":
#print "### QuickinfoExtension.execute(%s)" % repr(actionCommand)
layer = currentLayer()
if not self.isLayerValid(layer):
return
viewPanel = currentView().getWindowOfView()
mapControl = viewPanel.getMapControl()
exportmaptiles = ExportMapTilesPanel()
exportmaptiles.showTool("_Export_map_tiles")
def execute(self, actionCommand, *args):
actionCommand = actionCommand.lower()
if actionCommand == "settool-sentineldownloader":
#print "### QuickinfoExtension.execute(%s)" % repr(actionCommand)
layer = currentLayer()
if not self.isLayerValid(layer):
return
viewPanel = currentView().getWindowOfView()
mapControl = viewPanel.getMapControl()
sentineldownloader = SentinelSearchPanel()
sentineldownloader.setToolMapControl(mapControl)
sentineldownloader.showTool("Sentinel Search Panel")
def main(*args):
layer = gvsig.currentLayer()
ls = layer.getLabelingStrategy()
print ls
print ls.getFont()
print ls.getTextField()
print ls.getRotationField()
ds = LabelingFactory().createDefaultStrategy(layer)
print ds
ds.setTextField("LABEL")
ds.setRotationField("ROTATION")
def execute(self, actionCommand, *args):
actionCommand = actionCommand.lower()
if actionCommand == "settool-cadastralparcelsdownloader":
#print "### QuickinfoExtension.execute(%s)" % repr(actionCommand)
layer = currentLayer()
if not self.isLayerValid(layer):
return
viewPanel = currentView().getWindowOfView()
mapControl = viewPanel.getMapControl()
cadastral = CadastralDownloaderPanel()
#cadastral.setTool(mapControl)
cadastral.showTool("_Cadastral_Parcels_Downloader")
def main(*args):
proportionX = 1
proportionY = 1
xi = 0
yi = 0
baseLines = createBaseLayers(proportionX, proportionY)
# Setting coordinates to aoristic clock
nameFieldHour = "HORA"
nameFieldDay = "DIA"
patternHour = '%H:%M:%S'
patternDay = '%Y-%m-%d'
layer = gvsig.currentLayer()
# New points layer
schema = gvsig.createFeatureType(
layer.getFeatureStore().getDefaultFeatureType()) # DefaultFeatureType
newPoints = gvsig.createShape(schema)
# Transform
set = layer.getFeatureStore().getFeatureSet()
newPoints.edit()
store = newPoints.getFeatureStore()
for f in set:
fieldHour = f.get(nameFieldHour)
d = datetime.datetime.strptime(fieldHour, patternHour).time()
totalSecs = float(d.minute * 60 + d.second) / 3600
x = float(d.hour) + float(totalSecs)
x = x * proportionX
fieldDay = f.get(nameFieldDay)
dday = datetime.datetime.strptime(fieldDay, patternDay)
y = dday.weekday()
y = y * proportionY
nf = store.createNewFeature(f)
newGeom = geom.createPoint(geom.D2, x, y)
nf.setDefaultGeometry(newGeom)
store.insert(nf)
newPoints.commit()
gvsig.currentView().addLayer(newPoints)
mp = MapContextLocator.getMapContextManager()
leg = mp.createLegend("HeatmapLegend")
leg.setROI(baseLines.getFullEnvelope().getGeometry())
leg.setUseFixedViz(False)
leg.setCorrectionFixedViz(100)
leg.setDistance(30)
colorTables = SymbologySwingLocator.getSwingManager().createColorTables(
).get(5)
leg.setColorTable(colorTables.getColors())
newPoints.setLegend(leg)
def main(*args):
layer = gvsig.currentLayer()
#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 )
pass
def main(*args):
menu = JPopupMenu()
ep = ToolsLocator.getExtensionPointManager().get("View_TocActions")
layer = gvsig.currentLayer()
tocItem = TocItemLeaf(None, layer.getName(), layer.getShapeType())
for x in ep.iterator():
action = x.create()
print "action", action
continue
if action.isVisible(tocItem, (layer, )):
print action
menu.add(MenuItem(action, layer, tocItem))
else:
print "*** else:", action
def main2(*args):
store = gvsig.currentLayer().getFeatureStore()
dm = DALLocator.getDataManager()
i="True"
evaluator = dm.createExpresion(i)
fq = store.createFeatureQuery()
fq.addFilter(evaluator)
fq.retrievesAllAttributes()
fset = store.getFeatureSet(fq)
count = 0
for f in fset:
count+=1
print count
###
print "Spatial"
spatialGeometry = geom.createGeometryFromWKT("POLYGON ((-74.26662451007567 40.497289998, -74.26662451007567 40.517289997999995, -74.24930400199999 40.527289998, -74.2319834939243 40.517289997999995, -74.2319834939243 40.497289998, -74.24930400199999 40.487289997999994, -74.26662451007567 40.497289998))")
spatialGeometry.setProjection(gvsig.currentLayer().getProjection())
builder = ExpressionEvaluatorLocator.getManager().createExpressionBuilder()
expr1 = builder.and(
builder.custom(i),
builder.ST_Intersects(
builder.geometry(spatialGeometry),
builder.column("GEOMETRY")
)
).toString()
eva = ExpressionEvaluatorLocator.getManager().createExpression()
eva.setPhrase(expr1)
evaluator = DALLocator.getDataManager().createExpresion(eva)
fq = store.createFeatureQuery()
fq.addFilter(evaluator)
fq.retrievesAllAttributes()
fset = store.getFeatureSet(fq)
count = 0
for f in fset:
count+=1
print count
def main(*args): sextantelayer = gvsig.currentLayer() envelope = sextantelayer.getFullEnvelope() distancegrid = 0.1 gridType = GRID_HEXAGON_HORIZONTAL addEmptyGrids = True dm = DALLocator.getDataManager() i="KY_CD=344" valueExpression = dm.createExpresion(i) store = sextantelayer.getFeatureStore() projection = sextantelayer.getProjection() pointDensityGridCreation(None, store, gridType, distancegrid, addEmptyGrids, projection, envelope, valueExpression)
def main(*args):
layer = gvsig.currentLayer()
store = layer.getFeatureStore()
fset = store.getFeatureSet()
for feature in fset:
field = feature.get("HORA")
pattern = "HH:mm:ss"
#field = "23:10:10"
formatter = SimpleDateFormat(pattern)
newDate = formatter.parse(field)
cal = Calendar.getInstance()
cal.setTime(newDate)
hour = cal.get(Calendar.HOUR_OF_DAY)
print hour
def main1(*args):
layer = gvsig.currentLayer()
store = layer.getFeatureStore()
fset = store.getFeatureSelection() #.getFeatureSet()
cal = Calendar.getInstance()
print "First day: ", cal.getFirstDayOfWeek()
print Calendar.MONDAY, Calendar.SUNDAY, Calendar.SATURDAY
cal.setFirstDayOfWeek(Calendar.MONDAY)
print "First day: ", cal.getFirstDayOfWeek()
for feature in fset:
field = feature.get(0)
print type(field), field
cal.setTime(field)
print type(cal)
value = cal.get(Calendar.DAY_OF_WEEK)
print value
def main(*args):
viewDoc = gvsig.currentView()
layer = viewDoc.getLayer("manzanas_pob")
layer = gvsig.currentLayer()
panel = ReportByPointPanel(layer)
panel.setPreferredSize(400, 300)
winmgr = ToolsSwingLocator.getWindowManager()
dialog = winmgr.createDialog(panel.asJComponent(), "ReportByPoint test",
"ReportByPoint information",
winmgr.BUTTONS_OK_CANCEL)
dialog.show(winmgr.MODE.DIALOG)
if dialog.getAction() == winmgr.BUTTON_OK:
panel.save()
print "Ok"
print "Show field: "
else:
print "Cancel"
def main(*args):
proportionX = 1
proportionY = 1
nameFieldHour = "CDATE"
nameFieldDay = "CMPLNT_FR_"
nameFieldHour = "CMPLNT_FR_"
dm = DALLocator.getDataManager()
evaluator = dm.createExpresion("")
exp = ExpressionEvaluatorLocator.getManager().createExpression()
exp.setPhrase('')
print evaluator
expression = exp
xi = 0
yi = 0
store = gvsig.currentLayer().getFeatureStore()
aoristicClockGrid(store,
proportionX,
proportionY,
nameFieldHour,
nameFieldDay,
expression,
0,
0)
def main(*args):
currentProject = gvsig.currentProject() # Proyecto actual
# Vista actual
currentView = gvsig.currentView()
print "Vista actual:", currentView
# Capas
print "Capas de la vista actual:"
if currentView != None:
currentViewLayers = currentView.getLayers()
for currentViewLayer in currentViewLayers:
print " ", currentViewLayer
# Capa actual
currentLayer = gvsig.currentLayer()
print "Sumatory of all of the areas in the polygon layer: ", currentLayer
area = calcularArea(currentView, currentLayer)
if area != 0:
message = str(
float(area)) + " Km2" # https://es.wikipedia.org/wiki/Uruguay
title = unichr(
193) + "rea" # https://es.wikipedia.org/wiki/ISO/IEC_8859-1
messageType = commonsdialog.IDEA
root = None
commonsdialog.msgbox(message, title, messageType, root)
def main(*args):
#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 isEnabled(self):
layer = currentLayer()
#if not self.isLayerValid(layer):
# return False
return True
def main(*args):
gvsig.currentLayer().getDataStore().getFeatureSelection()
def main(*args):
store1 = gvsig.currentLayer().getFeatureStore()
builder = store1.createExpressionBuilder()
def main(*args): layer = gvsig.currentLayer() store = layer.getFeatureStore() print store.getDefaultFeatureType()
def main(*args):
store = gvsig.currentLayer().getFeatureStore()
l = FieldCalculatorTool(store)
l.showTool("Testing")
pass
def main(*args):
field = "MAPA"
des = gvsig.currentLayer().getFeatureStore().getDefaultFeatureType(
).getAttributeDescriptor(field)
print des.getObjectClass(), des.getSubtype()
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 main(*args):
#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")
print gvsig.currentLayer()
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")
#not working v6 = runalg("gvSIG-xyshift", "randompoints", "false", "-250.0", "-250.0", PATH=["C:/gvsig/ran10.shp","C:/gvsig/ran20.shp","C:/gvsig/ran30.shp"])
algHelp("tablebasicstats")
v7 = runalg("gvSIG-buffer", "randompoints", False, 50.0, 0, False, True, 0, 0, PATH="C:/gvsig/buffer_gvsig010.shp")
print "End"
def main(*args):
layer = gvsig.currentLayer()
p = 500000000
print areaTotal(layer)
gvsig.currentView().addLayer(copyLayerFilterByArea(p, layer))
