-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmain_program.py
52 lines (44 loc) · 2.32 KB
/
main_program.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
# -*- coding: UTF-8 -*-
import os
import csv
import imp
# ============ impression des images à partir d'Autocad ================================
#foo=imp.load_source('pilot_acad4','c:/georef/test8/pilot_acad4.py')
#chemin_acad='"C:\\Program Files\\Autodesk\\AutoCAD 2013\\acad.exe"'
#chemin_dvb='C:\\georef\\test8\\DGtoPNG.dvb'
#chemin_dwg='C:\\georef\\test8\\plan_install.dwg'
#nb_tot_process=4
#foo.impression(chemin_acad,chemin_dvb,chemin_dwg,nb_tot_process)
# ============= Géoréférencement des images =============================================
path_prog='/home/cvial/data/tuilage/'
path_drawings='/home/cvial/data/drawings/'
# alternative a l'import : sys.path.append()
foo = imp.load_source('georef', path_prog + 'georef.py')
nom_dwg = 'calepinage.dwg'
dossier_in = path_drawings + nom_dwg[:-4] + '/2_printed_images'
dossier_out = path_drawings + nom_dwg[:-4] + '/3_georef_images'
srs1 = '"+proj=utm +zone=23 +south +ellps=aust_SA +towgs84=-166.65,100.10,52.88,0,0 +no_defs"'
#srs2='"+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs"'
srs2 = '"+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs +over"'
foo.georef(dossier_in, dossier_out, srs1, srs2)
# ============= Création des tuiles =========================================
# ---- on commence par créer le fichier xml qui définit la carte ----
foo=imp.load_source('xml_carte', path_prog + 'xml_carte.py')
dossier_in=dossier_out
mapfile=dossier_in + '/calepinage.xml'
points_in=path_drawings + nom_dwg[:-4] + '/3_georef_images/pts_georef_finaux.csv'
# pour mapnik, il faut que le repertoire courant soit celui au dessus du fichier xml
os.chdir(path_drawings + nom_dwg[:-4] + '/')
foo.xmlage(dossier_in, srs2, mapfile, points_in)
# ---- puis on crée les tuiles ----
foo=imp.load_source('__main__', path_prog + 'generate_tiles_multinodes.py')
tile_dir=dossier_in.replace('3_georef_images', '4_tiles') + '/'
minZoom=19
maxZoom=23
foo.tuilage(mapfile, tile_dir, minZoom, maxZoom)
# ---- les tuiles qui n'ont pas été crées, on les récupère d'OSM ----
# !!!!!!!!!!! à créer !!!!!!!!!!!
# ============= Création d'un fichier avec les temps de création ===============
foo=imp.load_source('stats', path_prog + 'stats.py')
dossier = path_drawings + nom_dwg[:-4] + '/'
foo.stats(dossier)