|
| 1 | +from sextante.core.GeoAlgorithm import GeoAlgorithm |
| 2 | +from sextante.outputs.OutputVector import OutputVector |
| 3 | +from sextante.parameters.ParameterVector import ParameterVector |
| 4 | +from qgis.core import * |
| 5 | +from PyQt4.QtCore import * |
| 6 | +from PyQt4.QtGui import * |
| 7 | +from sextante.parameters.ParameterString import ParameterString |
| 8 | +from sextante.core.QGisLayers import QGisLayers |
| 9 | +import os |
| 10 | +from PyQt4 import QtGui |
| 11 | +from sextante.core.GeoAlgorithmExecutionException import GeoAlgorithmExecutionException |
| 12 | +from sextante.parameters.ParameterBoolean import ParameterBoolean |
| 13 | +import sys |
| 14 | + |
| 15 | + |
| 16 | +class FieldsPyculator(GeoAlgorithm): |
| 17 | + |
| 18 | + INPUT_LAYER = "INPUT_LAYER" |
| 19 | + USE_SELECTED = "USE_SELECTED" |
| 20 | + FIELD_NAME = "FIELD_NAME" |
| 21 | + #USE_GLOBAL = "USE_GLOBAL" |
| 22 | + GLOBAL = "GLOBAL" |
| 23 | + FORMULA = "FORMULA" |
| 24 | + OUTPUT_LAYER ="OUTPUT_LAYER" |
| 25 | + RESULT_VAR_NAME = "value" |
| 26 | + |
| 27 | + def getIcon(self): |
| 28 | + return QtGui.QIcon(os.path.dirname(__file__) + "/../images/toolbox.png") |
| 29 | + |
| 30 | + def defineCharacteristics(self): |
| 31 | + self.name = "Field Pyculator" |
| 32 | + self.group = "Algorithms for vector layers" |
| 33 | + self.addParameter(ParameterVector(self.INPUT_LAYER, "Input layer", ParameterVector.VECTOR_TYPE_ANY, False)) |
| 34 | + self.addParameter(ParameterBoolean(self.USE_SELECTED, "Use only selected features", False)) |
| 35 | + self.addParameter(ParameterString(self.FIELD_NAME, "Result field name", "NewField")) |
| 36 | + #self.addParameter(ParameterBoolean(self.USE_GLOBAL, "Use global expression", False)) |
| 37 | + self.addParameter(ParameterString(self.GLOBAL, "Global expression", multiline = True)) |
| 38 | + self.addParameter(ParameterString(self.FORMULA, "Formula", "value = ", multiline = True)) |
| 39 | + self.addOutput(OutputVector(self.OUTPUT_LAYER, "Output layer")) |
| 40 | + |
| 41 | + |
| 42 | + def processAlgorithm(self, progress): |
| 43 | + fieldname = self.getParameterValue(self.FIELD_NAME) |
| 44 | + code = self.getParameterValue(self.FORMULA) |
| 45 | + globalExpression = self.getParameterValue(self.GLOBAL) |
| 46 | + #useGlobal = self.getParameterValue(self.USE_GLOBAL) |
| 47 | + useSelected = self.getParameterValue(self.USE_SELECTED) |
| 48 | + settings = QSettings() |
| 49 | + systemEncoding = settings.value( "/UI/encoding", "System" ).toString() |
| 50 | + output = self.getOutputValue(self.OUTPUT_LAYER) |
| 51 | + layer = QGisLayers.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER)) |
| 52 | + vprovider = layer.dataProvider() |
| 53 | + allAttrs = vprovider.attributeIndexes() |
| 54 | + vprovider.select( allAttrs ) |
| 55 | + fields = vprovider.fields() |
| 56 | + fields[len(fields)] = QgsField(fieldname, QVariant.Double) |
| 57 | + writer = QgsVectorFileWriter( output, systemEncoding,fields, vprovider.geometryType(), vprovider.crs() ) |
| 58 | + outFeat = QgsFeature() |
| 59 | + nFeatures = vprovider.featureCount() |
| 60 | + nElement = 0 |
| 61 | + new_ns = {} |
| 62 | + |
| 63 | + #run global code |
| 64 | + if globalExpression.strip() != "": |
| 65 | + try: |
| 66 | + bytecode = compile(globalExpression, '<string>', 'exec') |
| 67 | + exec bytecode in new_ns |
| 68 | + except: |
| 69 | + raise GeoAlgorithmExecutionException("FieldPyculator code execute error\n" + |
| 70 | + "Global code block can't be executed!%s \n %s" % |
| 71 | + (unicode(sys.exc_info()[0].__name__), unicode(sys.exc_info()[1]))) |
| 72 | + |
| 73 | + #replace all fields tags |
| 74 | + field_map = vprovider.fields() |
| 75 | + for num, field in field_map.iteritems(): |
| 76 | + field_name = unicode(field.name()) |
| 77 | + replval = '__attr[' + str(num) + ']' |
| 78 | + code = code.replace("<"+field_name+">",replval) |
| 79 | + |
| 80 | + #replace all special vars |
| 81 | + code = code.replace('$id','__id') |
| 82 | + code = code.replace('$geom','__geom') |
| 83 | + need_id = code.find("__id") != -1 |
| 84 | + need_geom = code.find("__geom") != -1 |
| 85 | + need_attrs = code.find("__attr") != -1 |
| 86 | + |
| 87 | + |
| 88 | + #compile |
| 89 | + try: |
| 90 | + bytecode = compile(code, '<string>', 'exec') |
| 91 | + except: |
| 92 | + raise GeoAlgorithmExecutionException("FieldPyculator code execute error\n"+ |
| 93 | + "Field code block can't be executed! %s \n %s" |
| 94 | + (unicode(sys.exc_info()[0].__name__), unicode(sys.exc_info()[1]))) |
| 95 | + |
| 96 | + |
| 97 | + #run |
| 98 | + if not useSelected: |
| 99 | + feat = QgsFeature() |
| 100 | + if need_attrs: |
| 101 | + attr_ind = vprovider.attributeIndexes() |
| 102 | + else: |
| 103 | + attr_ind = [] |
| 104 | + vprovider.select(attr_ind, QgsRectangle(), True) |
| 105 | + |
| 106 | + while vprovider.nextFeature(feat): |
| 107 | + progress.setPercentage(int((100 * nElement)/nFeatures)) |
| 108 | + attrMap = feat.attributeMap() |
| 109 | + feat_id = feat.id() |
| 110 | + |
| 111 | + #add needed vars |
| 112 | + if need_id: |
| 113 | + new_ns['__id'] = feat_id |
| 114 | + |
| 115 | + if need_geom: |
| 116 | + geom = feat.geometry() |
| 117 | + new_ns['__geom'] = geom |
| 118 | + |
| 119 | + if need_attrs: |
| 120 | + attr = [] |
| 121 | + for num,a in attrMap.iteritems(): |
| 122 | + attr.append(self.Qvar2py(a)) |
| 123 | + new_ns['__attr'] = attr |
| 124 | + |
| 125 | + #clear old result |
| 126 | + if new_ns.has_key(self.RESULT_VAR_NAME): |
| 127 | + del new_ns[self.RESULT_VAR_NAME] |
| 128 | + |
| 129 | + |
| 130 | + #exec |
| 131 | + #try: |
| 132 | + exec bytecode in new_ns |
| 133 | + #except: |
| 134 | + # raise e |
| 135 | + #=============================================================== |
| 136 | + # GeoAlgorithmExecutionException("FieldPyculator code execute error\n"+ |
| 137 | + # "Field code block can't be executed for feature %s\n%s\n%s" % |
| 138 | + # (unicode(sys.exc_info()[0].__name__), |
| 139 | + # unicode(sys.exc_info()[1]), |
| 140 | + # unicode(feat_id))) |
| 141 | + #=============================================================== |
| 142 | + |
| 143 | + #check result |
| 144 | + if not new_ns.has_key(self.RESULT_VAR_NAME): |
| 145 | + raise GeoAlgorithmExecutionException("FieldPyculator code execute error\n" + |
| 146 | + "Field code block does not return '%s1' variable! Please declare this variable in your code!" % |
| 147 | + self.RESULT_VAR_NAME) |
| 148 | + |
| 149 | + |
| 150 | + #write feature |
| 151 | + nElement += 1 |
| 152 | + outFeat.setGeometry( feat.geometry() ) |
| 153 | + outFeat.setAttributeMap( attrMap ) |
| 154 | + outFeat.addAttribute(len(vprovider.fields()), QVariant(new_ns[self.RESULT_VAR_NAME])) |
| 155 | + writer.addFeature(outFeat) |
| 156 | + |
| 157 | + else: |
| 158 | + features = layer.selectedFeatures() |
| 159 | + nFeatures = len(features) |
| 160 | + for feat in features: |
| 161 | + progress.setPercentage(int((100 * nElement)/nFeatures)) |
| 162 | + attrMap = feat.attributeMap() |
| 163 | + feat_id = feat.id() |
| 164 | + |
| 165 | + #add needed vars |
| 166 | + if need_id: |
| 167 | + new_ns['__id'] = feat_id |
| 168 | + |
| 169 | + if need_geom: |
| 170 | + geom = feat.geometry() |
| 171 | + new_ns['__geom'] = geom |
| 172 | + |
| 173 | + if need_attrs: |
| 174 | + attrMap = feat.attributeMap() |
| 175 | + attr = [] |
| 176 | + for num,a in attrMap.iteritems(): |
| 177 | + attr.append(self.Qvar2py(a)) |
| 178 | + new_ns['__attr'] = attr |
| 179 | + |
| 180 | + #clear old result |
| 181 | + if new_ns.has_key(self.RESULT_VAR_NAME): |
| 182 | + del new_ns[self.RESULT_VAR_NAME] |
| 183 | + |
| 184 | + #exec |
| 185 | + exec bytecode in new_ns |
| 186 | + |
| 187 | + #check result |
| 188 | + if not new_ns.has_key(self.RESULT_VAR_NAME): |
| 189 | + raise GeoAlgorithmExecutionException("FieldPyculator code execute error\n" + |
| 190 | + "Field code block does not return '%s1' variable! Please declare this variable in your code!" % |
| 191 | + self.RESULT_VAR_NAME) |
| 192 | + |
| 193 | + #write feature |
| 194 | + nElement += 1 |
| 195 | + outFeat.setGeometry( feat.geometry() ) |
| 196 | + outFeat.setAttributeMap( attrMap ) |
| 197 | + outFeat.addAttribute(len(vprovider.fields()), QVariant(new_ns[self.RESULT_VAR_NAME])) |
| 198 | + writer.addFeature(outFeat) |
| 199 | + |
| 200 | + del writer |
| 201 | + |
| 202 | + |
| 203 | + def Qvar2py(self,qv): |
| 204 | + if qv.type() == 2: |
| 205 | + return qv.toInt()[0] |
| 206 | + if qv.type() == 10: |
| 207 | + return unicode(qv.toString()) |
| 208 | + if qv.type() == 6: |
| 209 | + return qv.toDouble()[0] |
| 210 | + return None |
| 211 | + |
| 212 | + |
| 213 | + def checkParameterValuesBeforeExecuting(self): |
| 214 | + ##TODO check that formula is correct and fields exist |
| 215 | + pass |
| 216 | + |
| 217 | + |
0 commit comments