-
Notifications
You must be signed in to change notification settings - Fork 2
/
cob_ter_sobrp.py
124 lines (114 loc) · 4.5 KB
/
cob_ter_sobrp.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
"""
/***************************************************************************
LEOXINGU
-------------------
begin : 2017-05-24
copyright : (C) 2017 by Leandro Franca - Cartographic Engineer
email : geoleandro.franca@gmail.com
***************************************************************************/
/***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation. *
* *
***************************************************************************/
"""
# Verificacao de cobertura terrestre (SOBREPOSICAO)
##10. Cob Ter Sobreposicao=name
##LF02) Revisao=group
##Camada_de_Sobreposicao=output vector
saida = Camada_de_Sobreposicao
tol = 0.1
from PyQt4.QtCore import *
from qgis.gui import QgsMessageBar
from qgis.utils import iface
from qgis.core import *
import processing
import time
CT = [ 'hid_massa_dagua_a',
'hid_trecho_massa_dagua_a',
'hid_barragem_a',
'loc_area_edificada_a',
'rel_terreno_exposto_a',
'rel_alter_fisiog_antropica_a',
'rel_rocha_a',
'tra_pista_ponto_pouso_a',
'veg_campo_a',
'veg_floresta_a',
'veg_estepe_a',
'veg_brejo_pantano_a',
'veg_caatinga_a',
'veg_campinarana_a',
'veg_cerrado_cerradao_a',
'veg_macega_chavascal_a',
'veg_mangue_a',
'veg_veg_area_contato_a',
'veg_veg_restinga_a',
'veg_vegetacao_a',
'veg_veg_cultivada_a']
# Pegando o SRC da camada trecho de drenagem
layerList = QgsMapLayerRegistry.instance().mapLayersByName('hid_trecho_drenagem_l')
if layerList:
layer = layerList[0]
SRC = layer.crs()
# Varrer camadas e pegar as geometrias das classes que compoe a cobertura terrestre
# Carregar cada feicao na camada temporaria
lista = []
cont = 1
for layer in QgsMapLayerRegistry.instance().mapLayers().values():
if layer.name() in CT:
for feature in layer.getFeatures():
geom = feature.geometry()
if geom:
coord = geom.asPolygon()
if coord == []:
COORD = geom.asMultiPolygon()
for coord in COORD:
if coord != []:
lista += [(coord, layer.name())]
cont +=1
else:
if coord != []:
lista += [(coord, layer.name())]
cont +=1
# Cirando camada de sobreposicao
fields = QgsFields()
fields.append(QgsField('id', QVariant.Int))
fields.append(QgsField('camada1', QVariant.String))
fields.append(QgsField('camada2', QVariant.String))
encoding = 'utf-8'
formato = 'ESRI Shapefile'
writer = QgsVectorFileWriter(saida, encoding, fields, QGis.WKBPolygon, SRC, formato)
# Verificar Sobreposicao
fet =QgsFeature()
tam = len(lista)
cont = 0
for i in range(tam-1):
for j in range(i+1, tam):
A = QgsGeometry.fromPolygon(lista[i][0])
B = QgsGeometry.fromPolygon(lista[j][0])
if A.intersects(B) or A.equals(B):
C = A.intersection(B)
if C.asPolygon() != []:
att = [cont, lista[i][1], lista[j][1]]
if C.area() > tol:
cont +=1
fet.setGeometry(C)
fet.setAttributes(att)
writer.addFeature(fet)
elif C.asMultiPolygon() != []:
for pol in C.asMultiPolygon():
att = [cont, lista[i][1], lista[j][1]]
geom = QgsGeometry.fromPolygon(pol)
if geom.area()>tol:
cont +=1
fet.setGeometry(geom)
fet.setAttributes(att)
writer.addFeature(fet)
progress.setPercentage(int(((i+1)/float(tam))*100))
del writer
progress.setInfo('<b>Operacao concluida!</b><br/><br/>')
progress.setInfo('<b>Leandro França - Eng Cart</b><br/>')
time.sleep(3)
iface.messageBar().pushMessage(u'Situacao', "Operacao Concluida com Sucesso!", level=QgsMessageBar.INFO, duration=5)