Newer
Older
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
import sys
from shapely.ops import unary_union
from utils import *
def computeData(df):
"""
Specific treatments
"""
# Init
ENV_targetProj = os.getenv('TARGET_PROJ')
# Reproj
df = checkAndReproj(df, ENV_targetProj)
# Init timer
bufferTimer = startTimerLog('Buffer Timer')
# Buffer 2m
dfBuffered = makeBufferFromGDF(df, 2)
# End timer
endTimerLog(bufferTimer)
# Union timer
unionTimer = startTimerLog('Union Timer')
# Regroup
unionGeom = unary_union(dfBuffered.geometry)
dataUnion = {'geometry': unionGeom}
unionFacadeGDF = gp.GeoDataFrame(dataUnion, crs=ENV_targetProj)
#TODO: + Clean and repair geom ??
# End union timer
endTimerLog(unionTimer)
# Init timer
simplifyTimer = startTimerLog('Simplify Timer')
# Simplify
dfSimplify = unionFacadeGDF.simplify(1)
# Reproj
unionBatiDF = checkAndReproj(dfSimplify, ENV_targetProj)
# End timer
endTimerLog(simplifyTimer)
return unionBatiDF
def wrongArguments():
print(style.RED + '/!\ Subscript error : Please specify a source and result filename as arguments to launch this script \n', style.RESET)
sys.exit(1)
if __name__ == "__main__":
# Init timer
subTimer = startTimerLog('Proximité façade subscript process')
# Get data with temp filename in argv
argv = sys.argv[1:]
firstArgv = None
secArgv = None
# Argv exist ?
if argv:
if len(sys.argv[1:]) > 0:
firstArgv = sys.argv[1:][0]
else:
wrongArguments()
if len(sys.argv[1:]) > 1:
secArgv = sys.argv[1:][1]
else:
wrongArguments()
# Load file Data (geoJSON)
currentGDF = createGDFfromGeoJSON(firstArgv)
# Log & Launch treatment
currentGDF = computeData(currentGDF)
# Write Result in temp file
currentGDF.to_file(secArgv)
# End timer
endTimerLog(subTimer)
else:
wrongArguments()