def processAlgorithm(self, progress):
GdalUtils.runGdal(self.getConsoleCommands(), progress)
output = self.getOutputValue(self.OUTPUT)
with open(output, 'w') as f:
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(s)
f.write('</pre>')
def processAlgorithm(self, parameters, context, feedback):
GdalUtils.runGdal(self.getConsoleCommands(parameters), feedback)
output = self.getOutputValue(self.OUTPUT)
with open(output, 'w') as f:
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(s)
f.write('</pre>')
def processAlgorithm(self, feedback):
GdalUtils.runGdal(self.getConsoleCommands(), feedback)
output = self.getOutputValue(information.OUTPUT)
with open(output, 'w') as f:
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(str(s))
f.write('</pre>')
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = ogrConnectionString(inLayerA)[1:-1]
layernameA = ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = ogrConnectionString(inLayerB)[1:-1]
layernameB = ogrLayerName(inLayerB)
fieldA = unicode(self.getParameterValue(self.FIELD_A))
fieldB = unicode(self.getParameterValue(self.FIELD_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerA = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER_A))
geomTypeA = layerA.geometryType()
wkbTypeA = layerA.wkbType()
sridA = layerA.crs().postgisSrid()
layerB = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER_B))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
multi = self.getParameterValue(self.MULTI)
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
if multi:
sqlstring = "-sql \"WITH temp_table AS (SELECT ST_Union(" + geomColumnB + ") AS geom FROM " + layernameB + ") SELECT ST_ShortestLine(g1." + geomColumnA + ",g2.geom) AS geom, ST_Distance(g1." + geomColumnA + ",g2.geom) AS distance, g1." + fieldA + " AS id_from FROM " + layernameA + " AS g1, temp_table AS g2\" -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom -nlt LINESTRING --config PG_USE_COPY YES -a_srs EPSG:" + str(sridA) + ""
else:
sqlstring = "-sql \"SELECT ST_ShortestLine(g1." + geomColumnA + ",g2." + geomColumnB + ") AS geom, ST_Distance(g1." + geomColumnA + ",g2." + geomColumnB + ") AS distance, g1." + fieldA + " AS id_from, g2." + fieldB + " AS id_to FROM " + layernameA + " AS g1, " + layernameB + " AS g2\" -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom -nlt LINESTRING --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
#print table
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
GdalUtils.runGdal(self.getConsoleCommands(), progress)
output = self.getOutputValue(information.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, parameters, context, feedback):
GdalUtils.runGdal(self.getConsoleCommands(parameters, context, feedback), feedback)
output = self.parameterAsFileOutput(parameters, self.OUTPUT, context)
with open(output, 'w') as f:
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(str(s))
f.write('</pre>')
return {self.OUTPUT: output}
def processAlgorithm(self, parameters, context, feedback):
GdalUtils.runGdal(self.getConsoleCommands(parameters, context, feedback), feedback)
output = self.parameterAsFileOutput(parameters, self.OUTPUT, context)
with open(output, 'w') as f:
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(str(s))
f.write('</pre>')
return {self.OUTPUT: output}
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = ogrConnectionString(inLayerA)[1:-1]
layernameA = ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = ogrConnectionString(inLayerB)[1:-1]
layernameB = ogrLayerName(inLayerB)
fieldsA = unicode(self.getParameterValue(self.FIELDS_A))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
single = self.getParameterValue(self.SINGLE)
if len(fieldsA) > 0:
fieldstring = fieldsA.replace(",", ", g1.")
fieldstring = ", g1." + fieldstring
else:
fieldstring = ""
if single:
sqlstring = "-sql \"SELECT (ST_Dump(ST_Difference(g1." + geomColumnA + ",ST_Union(g2." + geomColumnB + ")))).geom::geometry(Polygon) AS geom" + fieldstring + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 GROUP BY g1." + geomColumnA + fieldstring + "\""" -nln " + schema + "." + table + " -lco FID=gid -nlt POLYGON -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (ST_Multi(ST_CollectionExtract(ST_Difference(g1." + geomColumnA + ",ST_Union(g2." + geomColumnB + ")),3)))::geometry(MultiPolygon) AS geom" + fieldstring + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 GROUP BY g1." + geomColumnA + fieldstring + "\""" -nln " + schema + "." + table + " -lco FID=gid -nlt MULTIPOLYGON -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = ogrConnectionString(inLayer)[1:-1]
layername = ogrLayerName(inLayer)
fields = unicode(self.getParameterValue(self.FIELDS))
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
layer = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER))
geomType = layer.geometryType()
wkbType = layer.wkbType()
srid = layer.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
if len(fields) > 0:
fieldstring = "," + fields
else:
fieldstring = ""
if wkbType == 3:
layertype = "POLYGON"
sqlstring = "-sql \"SELECT (ST_Dump(ST_MakeValid(g1." + geomColumn + "))).geom::geometry(" + layertype + "," + str(srid) + ") AS geom" + fieldstring + " FROM " + layername + " AS g1\" -nlt " + layertype + " -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
layertype = "MULTIPOLYGON"
sqlstring = "-sql \"SELECT (ST_MakeValid(g1." + geomColumn + "))::geometry(" + layertype + "," + str(srid) + ") AS geom" + fieldstring + " FROM " + layername + " AS g1\" -nlt " + layertype + " -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = ogrConnectionString(inLayer)[1:-1]
layername = ogrLayerName(inLayer)
fields = unicode(self.getParameterValue(self.FIELDS))
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
layer = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER))
geomType = layer.geometryType()
wkbType = layer.wkbType()
srid = layer.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
if len(fields) > 0:
fieldstring = "," + fields
else:
fieldstring = ""
if wkbType == 3:
layertype = "POLYGON"
sqlstring = "-sql \"SELECT (ST_Dump(g1." + geomColumn + ")).geom::geometry(" + layertype + "," + str(srid) + ") AS geom, ST_IsValidReason(" + geomColumn + ") AS invalid_reason" + fieldstring + " FROM " + layername + " AS g1 WHERE NOT ST_IsValid(" + geomColumn + ")\" -nlt " + layertype + " -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
layertype = "MULTIPOLYGON"
sqlstring = "-sql \"SELECT (g1." + geomColumn + ")::geometry(" + layertype + "," + str(srid) + ") AS geom, ST_IsValidReason(" + geomColumn + ") AS invalid_reason" + fieldstring + " FROM " + layername + " AS g1 WHERE NOT ST_IsValid(" + geomColumn + ")\" -nlt " + layertype + " -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
arguments = []
if self.getParameterValue(information.NOGCP):
arguments.append('-nogcp')
if self.getParameterValue(information.NOMETADATA):
arguments.append('-nomd')
arguments.append(self.getParameterValue(information.INPUT))
GdalUtils.runGdal(['gdalinfo', GdalUtils.escapeAndJoin(arguments)],
progress)
output = self.getOutputValue(information.OUTPUT)
f = open(output, 'w')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write('<p>' + str(s) + '</p>')
f.close()
def processAlgorithm(self, progress):
arguments = []
arguments.append("-al")
arguments.append("-so")
layer = self.getParameterValue(self.INPUT)
conn = self.ogrConnectionString(layer)
arguments.append(conn)
GdalUtils.runGdal(["ogrinfo", GdalUtils.escapeAndJoin(arguments)], progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, "w")
f.write("<pre>")
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write("</pre>")
f.close()
def processAlgorithm(self, progress):
arguments = []
if self.getParameterValue(information.NOGCP):
arguments.append("-nogcp")
if self.getParameterValue(information.NOMETADATA):
arguments.append("-nomd")
arguments.append(self.getParameterValue(information.INPUT))
GdalUtils.runGdal(["gdalinfo", GdalUtils.escapeAndJoin(arguments)], progress)
output = self.getOutputValue(information.OUTPUT)
f = open(output, "w")
f.write("<pre>")
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write("</pre>")
f.close()
def processAlgorithm(self, progress):
arguments = []
arguments.append('-al')
arguments.append('-so')
layer = self.getParameterValue(self.INPUT)
conn = self.ogrConnectionString(layer)
arguments.append(conn)
GdalUtils.runGdal(
['ogrinfo', GdalUtils.escapeAndJoin(arguments)], progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = ogrConnectionString(inLayerA)[1:-1]
layernameA = ogrLayerName(inLayerA)
layerA = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER_A))
fieldA = unicode(self.getParameterValue(self.FIELD_A))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
geomTypeA = layerA.geometryType()
wkbTypeA = layerA.wkbType()
sridA = layerA.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
sqlstring = "-sql \"SELECT ST_Makeline(g1." + geomColumnA + ",g2." + geomColumnA + ") AS geom, ST_Distance(g1." + geomColumnA + ",g2." + geomColumnA + ") AS distance, g1." + fieldA + " AS id_from, g2." + fieldA + " AS id_to FROM " + layernameA + " AS g1, " + layernameA + " AS g2 WHERE g1." + fieldA + " > " + "g2." + fieldA + "\" -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom -nlt LINESTRING --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
#print table
commands = []
if isWindows():
commands = [
'cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)
]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = self.ogrConnectionString(inLayer)
layername = self.ogrLayerName(inLayer)
field = unicode(self.getParameterValue(self.FIELD))
fields = unicode(self.getParameterValue(self.FIELDS))
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
if len(fields) > 0:
fieldstring = "," + fields
else:
fieldstring = ""
sqlstring = "-sql \"SELECT ST_MakeValid(g1." + geomColumn + ") AS geom, g1." + field + " AS id_valid" + fieldstring + " FROM " + layername + " AS g1\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = ogrConnectionString(inLayerA)[1:-1]
layernameA = ogrLayerName(inLayerA)
fieldsA = unicode(self.getParameterValue(self.FIELDS_A))
fieldA = unicode(self.getParameterValue(self.FIELD_A))
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = ogrConnectionString(inLayerB)[1:-1]
layernameB = ogrLayerName(inLayerB)
fieldB = unicode(self.getParameterValue(self.FIELD_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerB = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER_B))
layerA = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER_A))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
sridA = layerA.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
single = self.getParameterValue(self.SINGLE)
if len(fieldsA) > 0:
fieldstringA = fieldsA.replace(",", ", g1.")
fieldstringA = ", g1." + fieldstringA
else:
fieldstringA = ""
type = "POINT"
if single:
multistring = "-nlt " + type
caststring = type
st_function = "ST_Dump"
castgeom = ".geom"
else:
multistring = "-nlt MULTI" + type
caststring = "MULTI" + type
st_function = "ST_Multi"
castgeom = ""
sqlstring = "-sql \"SELECT (" + st_function + "(g1." + geomColumnA + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridA
) + ") AS geom, g2." + fieldB + " AS " + fieldA + fieldstringA + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g2." + geomColumnB + ",g1." + geomColumnA + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
print geomTypeB
commands = []
if isWindows():
commands = [
'cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)
]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = ogrConnectionString(inLayer)[1:-1]
layername = ogrLayerName(inLayer)
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
layer = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER))
geomType = layer.geometryType()
wkbType = layer.wkbType()
srid = layer.crs().postgisSrid()
field = unicode(self.getParameterValue(self.FIELD))
statsatt = unicode(self.getParameterValue(self.STATSATT))
stats = self.getParameterValue(self.STATS)
area = self.getParameterValue(self.AREA)
single = self.getParameterValue(self.SINGLE)
count = self.getParameterValue(self.COUNT)
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
dissolveall = self.getParameterValue(self.DISSOLVEALL)
options = unicode(self.getParameterValue(self.OPTIONS))
if single:
layertype = "POLYGON"
else:
layertype = "MULTIPOLYGON"
if dissolveall:
fieldstring = ""
else:
fieldstring = "," + field
if single:
querystart = '-sql "SELECT (ST_Dump(ST_Union(' + geomColumn + '))).geom::geometry(POLYGON,' + str(
srid) + ')' + fieldstring
else:
querystart = '-sql "SELECT (ST_Multi(ST_Union(' + geomColumn + ')))::geometry(MULTIPOLYGON,' + str(
srid) + ')' + fieldstring
if dissolveall:
queryend = ' FROM ' + layername + '"' + " -nln " + schema + "." + table + " -nlt " + layertype + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
queryend = ' FROM ' + layername + ' GROUP BY ' + field + '"' + " -nln " + schema + "." + table + " -nlt " + layertype + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
#if fields:
# queryfields = ",*"
#else:
# queryfields = "," + field
if count:
querycount = ", COUNT(" + geomColumn + ") AS count"
else:
querycount = ""
if stats:
querystats = ", SUM(" + statsatt + ") AS sum_dissolved, MIN(" + statsatt + ") AS min_dissolved, MAX(" + statsatt + ") AS max_dissolved, AVG(" + statsatt + ") AS avg_dissolved"
else:
querystats = ""
if area:
queryarea = ", SUM(ST_area(" + geomColumn + ")) AS area_dissolved, ST_perimeter(ST_union(" + geomColumn + ")) AS perimeter_dissolved"
else:
queryarea = ""
query = querystart + querystats + queryarea + querycount + queryend
#query = querystart + queryfields + querycount + querystats + queryarea + queryend
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(query)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = [
'cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)
]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = ogrConnectionString(inLayerA)[1:-1]
layernameA = ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = ogrConnectionString(inLayerB)[1:-1]
layernameB = ogrLayerName(inLayerB)
fieldsB = unicode(self.getParameterValue(self.FIELDS_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerB = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER_B))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
single = self.getParameterValue(self.SINGLE)
keep = self.getParameterValue(self.KEEP)
keepc = self.getParameterValue(self.KEEPC)
if len(fieldsB) > 0:
fieldstringB = fieldsB.replace(",", ", g2.")
fieldstringB = ", g2." + fieldstringB
else:
fieldstringB = ""
if geomTypeB == 0:
type = "POINT"
elif geomTypeB == 1:
type = "LINESTRING"
else:
type = "POLYGON"
if single:
multistring = "-nlt " + type
caststring = type
st_function = "ST_Dump"
castgeom = ".geom"
else:
multistring = "-nlt MULTI" + type
caststring = "MULTI" + type
st_function = "ST_Multi"
castgeom = ""
if geomTypeB == 0:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"WITH temp_table AS (SELECT ST_Union(" + geomColumnA + ") AS geom FROM " + layernameA + ") SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM temp_table AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES -a_srs EPSG:" + str(
sridB) + ""
elif geomTypeB == 1:
if keepc:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_ContainsProperly(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true AND ST_Touches(g1." + geomColumnA + ",g2." + geomColumnB + ") is false\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
if keepc:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_ContainsProperly(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true OR ST_Contains(g2." + geomColumnB + ",g1." + geomColumnA + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(
sridB
) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true OR ST_Overlaps(g1." + geomColumnA + ",g2." + geomColumnB + ") is true OR ST_Contains(g2." + geomColumnB + ",g1." + geomColumnA + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
print geomTypeB
commands = []
if isWindows():
commands = [
'cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)
]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = ogrConnectionString(inLayer)[1:-1]
layername = ogrLayerName(inLayer)
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
layer = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER))
geomType = layer.geometryType()
wkbType = layer.wkbType()
srid = layer.crs().postgisSrid()
fields = unicode(self.getParameterValue(self.FIELDS))
dissfield = unicode(self.getParameterValue(self.DISSFIELD))
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
options = unicode(self.getParameterValue(self.OPTIONS))
distance = unicode(self.getParameterValue(self.DISTANCE))
dissolveall = self.getParameterValue(self.DISSOLVEALL)
single = self.getParameterValue(self.SINGLE)
if len(fields) > 0:
fieldstring = fields + ","
else:
fieldstring = ""
if dissolveall:
if single:
query = '-sql "SELECT (ST_Dump(ST_Multi(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + '))))).geom::geometry(POLYGON,' + str(
srid
) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt POLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
query = '-sql "SELECT (ST_Multi(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + '))))::geometry(MULTIPOLYGON,' + str(
srid
) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt MULTIPOLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
if dissfield != 'None':
if single:
query = '-sql "SELECT ' + dissfield + ',(ST_Dump(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + ')))).geom::geometry(POLYGON,' + str(
srid
) + ') AS geom FROM ' + layername + ' GROUP BY ' + dissfield + '" -nln ' + schema + '.' + table + ' -nlt POLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
query = '-sql "SELECT ' + dissfield + ',(ST_Multi(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + '))))::geometry(MULTIPOLYGON,' + str(
srid
) + ') AS geom FROM ' + layername + ' GROUP BY ' + dissfield + '" -nln ' + schema + '.' + table + ' -nlt MULTIPOLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
if single:
query = '-sql "SELECT ' + fieldstring + '(ST_Dump(ST_Multi(ST_Buffer(' + geomColumn + ',' + distance + ')))).geom::geometry(POLYGON,' + str(
srid
) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt POLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
query = '-sql "SELECT ' + fieldstring + '(ST_Multi(ST_Buffer(' + geomColumn + ',' + distance + ')))::geometry(MULTIPOLYGON,' + str(
srid
) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt MULTIPOLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
print dissfield
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(query)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = [
'cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)
]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = self.ogrConnectionString(inLayerA)[1:-1]
layernameA = self.ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = self.ogrConnectionString(inLayerB)[1:-1]
layernameB = self.ogrLayerName(inLayerB)
fieldsB = unicode(self.getParameterValue(self.FIELDS_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerB = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER_B))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
single = self.getParameterValue(self.SINGLE)
keep = self.getParameterValue(self.KEEP)
bufferdist = self.getParameterValue(self.BUFFER)
if len(fieldsB) > 0:
fieldstringB = "," + fieldsB
else:
fieldstringB = ""
if geomTypeB == 0:
type = "POINT"
elif geomTypeB == 1:
type = "LINESTRING"
else:
type = "POLYGON"
if single:
multistring = "-nlt " + type
caststring = type
st_function = "ST_Dump"
castgeom = ".geom"
else:
multistring = "-nlt MULTI" + type
caststring = "MULTI" + type
st_function = "ST_Multi"
castgeom = ""
if geomTypeB == 0:
if bufferdist == "0" or bufferdist == "":
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE ST_Intersects(g2."
+ geomColumnB
+ ",g1."
+ geomColumnA
+ ') is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
if keep:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE (ST_Intersects(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ")) is true OR (ST_Touches(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ')) is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE (ST_Intersects(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ")) is true AND (ST_Touches(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ')) is false" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
elif geomTypeB == 1:
if bufferdist == "0" or bufferdist == "":
if keep:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE ST_Intersects(g2."
+ geomColumnB
+ ",g1."
+ geomColumnA
+ ') is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE ST_Intersects(g2."
+ geomColumnB
+ ",g1."
+ geomColumnA
+ ") is true AND ST_Touches(g1."
+ geomColumnA
+ ",g2."
+ geomColumnB
+ ') is false" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
if keep:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE (ST_Intersects(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ")) is true OR (ST_Touches(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ')) is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE (ST_Intersects(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ")) is true AND (ST_Touches(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ')) is false" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
if bufferdist == "0" or bufferdist == "":
if keep:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE ST_Intersects(g2."
+ geomColumnB
+ ",g1."
+ geomColumnA
+ ') is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE ST_Intersects(g2."
+ geomColumnB
+ ",g1."
+ geomColumnA
+ ") is true AND ST_Touches(g2."
+ geomColumnB
+ ",g1."
+ geomColumnA
+ ') is false" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
if keep:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE (ST_Intersects(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ")) is true OR ST_Contains(g2."
+ geomColumnB
+ ",ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ ')) is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
else:
sqlstring = (
'-sql "SELECT ('
+ st_function
+ "(g2."
+ geomColumnB
+ "))"
+ castgeom
+ "::geometry("
+ caststring
+ ","
+ str(sridB)
+ ") AS geom"
+ fieldstringB
+ " FROM "
+ layernameA
+ " AS g1, "
+ layernameB
+ " AS g2 WHERE (ST_Contains(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ")) is true OR ST_Overlaps(ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ "),g2."
+ geomColumnB
+ ") is true OR ST_Contains(g2."
+ geomColumnB
+ ",ST_Buffer(g1."
+ geomColumnA
+ ","
+ bufferdist
+ ')) is true" -nln '
+ schema
+ "."
+ table
+ " -lco FID=gid "
+ multistring
+ " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
)
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append("-f")
arguments.append("PostgreSQL")
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append("-overwrite")
if len(options) > 0:
arguments.append(options)
print geomTypeB
commands = []
if isWindows():
commands = ["cmd.exe", "/C ", "ogr2ogr.exe", GdalUtils.escapeAndJoin(arguments)]
else:
commands = ["ogr2ogr", GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, "w")
f.write("<pre>")
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write("</pre>")
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = ogrConnectionString(inLayer)[1:-1]
layername = ogrLayerName(inLayer)
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
layer = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER))
geomType = layer.geometryType()
wkbType = layer.wkbType()
srid = layer.crs().postgisSrid()
field = unicode(self.getParameterValue(self.FIELD))
statsatt = unicode(self.getParameterValue(self.STATSATT))
stats = self.getParameterValue(self.STATS)
area = self.getParameterValue(self.AREA)
single = self.getParameterValue(self.SINGLE)
count = self.getParameterValue(self.COUNT)
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
dissolveall = self.getParameterValue(self.DISSOLVEALL)
options = unicode(self.getParameterValue(self.OPTIONS))
if single:
layertype = "POLYGON"
else:
layertype = "MULTIPOLYGON"
if dissolveall:
fieldstring = ""
else:
fieldstring = "," + field
if single:
querystart = '-sql "SELECT (ST_Dump(ST_Union(' + geomColumn + '))).geom::geometry(POLYGON,' + str(srid) + ')' + fieldstring
else:
querystart = '-sql "SELECT (ST_Multi(ST_Union(' + geomColumn + ')))::geometry(MULTIPOLYGON,' + str(srid) + ')' + fieldstring
if dissolveall:
queryend = ' FROM ' + layername + '"' + " -nln " + schema + "." + table + " -nlt " + layertype + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
queryend = ' FROM ' + layername + ' GROUP BY ' + field + '"' + " -nln " + schema + "." + table + " -nlt " + layertype + " -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
#if fields:
# queryfields = ",*"
#else:
# queryfields = "," + field
if count:
querycount = ", COUNT(" + geomColumn + ") AS count"
else:
querycount = ""
if stats:
querystats = ", SUM(" + statsatt + ") AS sum_dissolved, MIN(" + statsatt + ") AS min_dissolved, MAX(" + statsatt + ") AS max_dissolved, AVG(" + statsatt + ") AS avg_dissolved"
else:
querystats = ""
if area:
queryarea = ", SUM(ST_area(" + geomColumn + ")) AS area_dissolved, ST_perimeter(ST_union(" + geomColumn + ")) AS perimeter_dissolved"
else:
queryarea = ""
query = querystart + querystats + queryarea + querycount + queryend
#query = querystart + queryfields + querycount + querystats + queryarea + queryend
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(query)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
connection = self.DB_CONNECTIONS[self.getParameterValue(self.DATABASE)]
settings = QSettings()
mySettings = '/PostgreSQL/connections/' + connection
dbname = settings.value(mySettings + '/database')
user = settings.value(mySettings + '/username')
host = settings.value(mySettings + '/host')
port = settings.value(mySettings + '/port')
password = settings.value(mySettings + '/password')
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = self.ogrConnectionString(inLayer)
ssrs = unicode(self.getParameterValue(self.S_SRS))
tsrs = unicode(self.getParameterValue(self.T_SRS))
schema = unicode(self.getParameterValue(self.SCHEMA))
schemastring = "-lco SCHEMA="+schema
table = unicode(self.getParameterValue(self.TABLE))
pk = unicode(self.getParameterValue(self.PK))
pkstring = "-lco FID="+pk
geocolumn = unicode(self.getParameterValue(self.GEOCOLUMN))
geocolumnstring = "-lco GEOMETRY_NAME="+geocolumn
dim = self.DIMLIST[self.getParameterValue(self.DIM)]
dimstring = "-lco DIM="+dim
simplify = unicode(self.getParameterValue(self.SIMPLIFY))
segmentize = unicode(self.getParameterValue(self.SEGMENTIZE))
spat = self.getParameterValue(self.SPAT)
ogrspat = self.ogrConnectionString(spat)
clip = self.getParameterValue(self.CLIP)
where = unicode(self.getParameterValue(self.WHERE))
wherestring = "-where '"+where+"'"
gt = unicode(self.getParameterValue(self.GT))
overwrite = self.getParameterValue(self.OVERWRITE)
append = self.getParameterValue(self.APPEND)
addfields = self.getParameterValue(self.ADDFIELDS)
launder = self.getParameterValue(self.LAUNDER)
launderstring = "-lco LAUNDER=NO"
index = self.getParameterValue(self.INDEX)
indexstring = "-lco SPATIAL_INDEX=OFF"
skipfailures = self.getParameterValue(self.SKIPFAILURES)
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-progress')
arguments.append('--config PG_USE_COPY YES')
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append('PG:"host=')
arguments.append(host)
arguments.append('port=')
arguments.append(port)
arguments.append('user='******'dbname=')
arguments.append(dbname)
arguments.append('password='******'"')
arguments.append(dimstring)
arguments.append(ogrLayer)
if index:
arguments.append(indexstring)
if launder:
arguments.append(launderstring)
if append:
arguments.append('-append')
if addfields:
arguments.append('-addfields')
if overwrite:
arguments.append('-overwrite')
if len(self.GEOMTYPE[self.getParameterValue(self.GTYPE)]) > 0:
arguments.append('-nlt')
arguments.append(self.GEOMTYPE[self.getParameterValue(self.GTYPE)])
if len(schema) > 0:
arguments.append(schemastring)
if len(geocolumn) > 0:
arguments.append(geocolumnstring)
if len(pk) > 0:
arguments.append(pkstring)
if len(table) > 0:
arguments.append('-nln')
arguments.append(table)
if len(ssrs) > 0:
arguments.append('-s_srs')
arguments.append(ssrs)
if len(tsrs) > 0:
arguments.append('-t_srs')
arguments.append(tsrs)
if len(spat) > 0:
regionCoords = ogrspat.split(',')
arguments.append('-spat')
arguments.append(regionCoords[0])
arguments.append(regionCoords[2])
arguments.append(regionCoords[1])
arguments.append(regionCoords[3])
if clip:
arguments.append('-clipsrc spat_extent')
if skipfailures:
arguments.append('-skipfailures')
if where:
arguments.append(wherestring)
if len(simplify) > 0:
arguments.append('-simplify')
arguments.append(simplify)
if len(segmentize) > 0:
arguments.append('-segmentize')
arguments.append(segmentize)
if len(gt) > 0:
arguments.append('-gt')
arguments.append(gt)
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = self.ogrConnectionString(inLayerA)[1:-1]
layernameA = self.ogrLayerName(inLayerA)
fieldsA = unicode(self.getParameterValue(self.FIELDS_A))
fieldA = unicode(self.getParameterValue(self.FIELD_A))
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = self.ogrConnectionString(inLayerB)[1:-1]
layernameB = self.ogrLayerName(inLayerB)
fieldB = unicode(self.getParameterValue(self.FIELD_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerB = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER_B))
layerA = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER_A))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
sridA = layerA.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
single = self.getParameterValue(self.SINGLE)
if len(fieldsA) > 0:
fieldstringA = "," + fieldsA
else:
fieldstringA = ""
type = "POINT"
if single:
multistring = "-nlt " + type
caststring = type
st_function = "ST_Dump"
castgeom = ".geom"
else:
multistring = "-nlt MULTI" + type
caststring = "MULTI" + type
st_function = "ST_Multi"
castgeom = ""
sqlstring = "-sql \"SELECT (" + st_function + "(g1." + geomColumnA + "))" + castgeom + "::geometry(" + caststring + "," + str(sridA) + ") AS geom, g2." + fieldB + " AS " + fieldA + fieldstringA + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g2." + geomColumnB + ",g1." + geomColumnA + ") is true\" -nln " + schema + "." + table + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
print geomTypeB
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = ogrConnectionString(inLayerA)[1:-1]
layernameA = ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = ogrConnectionString(inLayerB)[1:-1]
layernameB = ogrLayerName(inLayerB)
fieldA = unicode(self.getParameterValue(self.FIELD_A))
fieldB = unicode(self.getParameterValue(self.FIELD_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerA = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER_A))
geomTypeA = layerA.geometryType()
wkbTypeA = layerA.wkbType()
sridA = layerA.crs().postgisSrid()
layerB = dataobjects.getObjectFromUri(
self.getParameterValue(self.INPUT_LAYER_B))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
fields = unicode(self.getParameterValue(self.FIELDS))
multi = self.getParameterValue(self.MULTI)
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
if len(fields) > 0:
fieldstring = fields.replace(",", ", g1.")
fieldstring = ", g1." + fieldstring
else:
fieldstring = ""
if multi:
sqlstring = "-sql \"WITH temp_table AS (SELECT ST_Union(" + geomColumnB + ") AS geom FROM " + layernameB + ") SELECT (ST_ClosestPoint(g2.geom, g1." + geomColumnA + "))::geometry(POINT," + str(
sridA
) + ") AS geom, ST_Distance(g1." + geomColumnA + ",g2.geom) AS distance, g1." + fieldA + " AS id_from" + fieldstring + " FROM " + layernameA + " AS g1, temp_table AS g2\" -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom -nlt POINT --config PG_USE_COPY YES -a_srs EPSG:" + str(
sridA) + ""
else:
sqlstring = "-sql \"SELECT (ST_ClosestPoint(g2." + geomColumnB + ",g1." + geomColumnA + "))::geometry(Point," + str(
sridA
) + ") AS geom, ST_Distance(g1." + geomColumnA + ",g2." + geomColumnB + ") AS distance, g1." + fieldA + " AS id_from" + fieldstring + ", g2." + fieldB + " AS id_to FROM " + layernameA + " AS g1, " + layernameB + " AS g2\" -nln " + schema + "." + table + " -lco FID=gid -lco GEOMETRY_NAME=geom -nlt POINT --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
#print table
commands = []
if isWindows():
commands = [
'cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)
]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = self.ogrConnectionString(inLayerA)[1:-1]
layernameA = self.ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = self.ogrConnectionString(inLayerB)[1:-1]
layernameB = self.ogrLayerName(inLayerB)
fieldsB = unicode(self.getParameterValue(self.FIELDS_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
layerB = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER_B))
geomTypeB = layerB.geometryType()
wkbTypeB = layerB.wkbType()
sridB = layerB.crs().postgisSrid()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
single = self.getParameterValue(self.SINGLE)
keep = self.getParameterValue(self.KEEP)
keepc = self.getParameterValue(self.KEEPC)
if len(fieldsB) > 0:
fieldstringB = "," + fieldsB
else:
fieldstringB = ""
if geomTypeB == 0:
type = "POINT"
elif geomTypeB == 1:
type = "LINESTRING"
else:
type = "POLYGON"
if single:
multistring = "-nlt " + type
caststring = type
st_function = "ST_Dump"
castgeom = ".geom"
else:
multistring = "-nlt MULTI" + type
caststring = "MULTI" + type
st_function = "ST_Multi"
castgeom = ""
if geomTypeB == 0:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"WITH temp_table AS (SELECT ST_Union(" + geomColumnA + ") AS geom FROM " + layernameA + ") SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM temp_table AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES -a_srs EPSG:" + str(sridB) + ""
elif geomTypeB == 1:
if keepc:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_ContainsProperly(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true AND ST_Touches(g1." + geomColumnA + ",g2." + geomColumnB + ") is false\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
if keepc:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_ContainsProperly(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
if keep:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true OR ST_Contains(g2." + geomColumnB + ",g1." + geomColumnA + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
else:
sqlstring = "-sql \"SELECT (" + st_function + "(g2." + geomColumnB + "))" + castgeom + "::geometry(" + caststring + "," + str(sridB) + ") AS geom" + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Contains(g1." + geomColumnA + ",g2." + geomColumnB + ") is true OR ST_Overlaps(g1." + geomColumnA + ",g2." + geomColumnB + ") is true OR ST_Contains(g2." + geomColumnB+ ",g1." + geomColumnA + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
print geomTypeB
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayerA = self.getParameterValue(self.INPUT_LAYER_A)
ogrLayerA = self.ogrConnectionString(inLayerA)
layernameA = self.ogrLayerName(inLayerA)
inLayerB = self.getParameterValue(self.INPUT_LAYER_B)
ogrLayerB = self.ogrConnectionString(inLayerB)
layernameB = self.ogrLayerName(inLayerB)
fieldA = unicode(self.getParameterValue(self.FIELD_A))
fieldB = unicode(self.getParameterValue(self.FIELD_B))
fieldsA = unicode(self.getParameterValue(self.FIELDS_A))
fieldsB = unicode(self.getParameterValue(self.FIELDS_B))
dsUriA = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_A))
geomColumnA = dsUriA.geometryColumn()
dsUriB = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER_B))
geomColumnB = dsUriB.geometryColumn()
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
multi = self.getParameterValue(self.MULTI)
if len(fieldsB) > 0:
fieldstringB = "," + fieldsB
else:
fieldstringB = ""
if len(fieldsA) > 0:
fieldstringA = "," + fieldsA
else:
fieldstringA = ""
if multi:
multistring = "-nlt MULTIPOLYGON"
else:
multistring = "-nlt POLYGON"
sqlstring = "-sql \"SELECT ST_Intersection(g1." + geomColumnA + ",g2." + geomColumnB + ") AS geom, g1. " + fieldA + " AS id_clipper, g2. " + fieldB + " AS id_clipped" + fieldstringA + fieldstringB + " FROM " + layernameA + " AS g1, " + layernameB + " AS g2 WHERE ST_Intersects(g1." + geomColumnA + ",g2." + geomColumnB + ") is true\" -nln " + table + " -lco SCHEMA=" + schema + " -lco FID=gid " + multistring + " -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES"
options = unicode(self.getParameterValue(self.OPTIONS))
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayerA)
arguments.append(ogrLayerA)
arguments.append(sqlstring)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT_LAYER)
ogrLayer = ogrConnectionString(inLayer)[1:-1]
layername = ogrLayerName(inLayer)
dsUri = QgsDataSourceURI(self.getParameterValue(self.INPUT_LAYER))
geomColumn = dsUri.geometryColumn()
layer = dataobjects.getObjectFromUri(self.getParameterValue(self.INPUT_LAYER))
geomType = layer.geometryType()
wkbType = layer.wkbType()
srid = layer.crs().postgisSrid()
fields = unicode(self.getParameterValue(self.FIELDS))
dissfield = unicode(self.getParameterValue(self.DISSFIELD))
schema = unicode(self.getParameterValue(self.SCHEMA))
table = unicode(self.getParameterValue(self.TABLE))
options = unicode(self.getParameterValue(self.OPTIONS))
distance = unicode(self.getParameterValue(self.DISTANCE))
dissolveall = self.getParameterValue(self.DISSOLVEALL)
single = self.getParameterValue(self.SINGLE)
if len(fields) > 0:
fieldstring = fields + ","
else:
fieldstring = ""
if dissolveall:
if single:
query = '-sql "SELECT (ST_Dump(ST_Multi(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + '))))).geom::geometry(POLYGON,' + str(srid) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt POLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
query = '-sql "SELECT (ST_Multi(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + '))))::geometry(MULTIPOLYGON,' + str(srid) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt MULTIPOLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
if dissfield != 'None':
if single:
query = '-sql "SELECT ' + dissfield + ',(ST_Dump(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + ')))).geom::geometry(POLYGON,' + str(srid) + ') AS geom FROM ' + layername + ' GROUP BY ' + dissfield + '" -nln ' + schema + '.' + table + ' -nlt POLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
query = '-sql "SELECT ' + dissfield + ',(ST_Multi(ST_Union(ST_Buffer(' + geomColumn + ',' + distance + '))))::geometry(MULTIPOLYGON,' + str(srid) + ') AS geom FROM ' + layername + ' GROUP BY ' + dissfield + '" -nln ' + schema + '.' + table + ' -nlt MULTIPOLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
if single:
query = '-sql "SELECT ' + fieldstring + '(ST_Dump(ST_Multi(ST_Buffer(' + geomColumn + ',' + distance + ')))).geom::geometry(POLYGON,' + str(srid) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt POLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
else:
query = '-sql "SELECT ' + fieldstring + '(ST_Multi(ST_Buffer(' + geomColumn + ',' + distance + ')))::geometry(MULTIPOLYGON,' + str(srid) + ') AS geom FROM ' + layername + '" -nln ' + schema + '.' + table + ' -nlt MULTIPOLYGON -lco FID=gid -lco GEOMETRY_NAME=geom --config PG_USE_COPY YES'
print dissfield
arguments = []
arguments.append('-f')
arguments.append('PostgreSQL')
arguments.append(ogrLayer)
arguments.append(ogrLayer)
arguments.append(query)
arguments.append('-overwrite')
if len(options) > 0:
arguments.append(options)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'ogr2ogr.exe',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
output = self.getOutputValue(self.OUTPUT)
f = open(output, 'w')
f.write('<pre>')
for s in GdalUtils.getConsoleOutput()[1:]:
f.write(unicode(s))
f.write('</pre>')
f.close()
