forked from qgispt/qgis_pt_processing_pttransform
-
Notifications
You must be signed in to change notification settings - Fork 0
/
UTM29NED50ToETR89PTTM06.py
93 lines (72 loc) · 3.24 KB
/
UTM29NED50ToETR89PTTM06.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
# -*- coding: latin-1 -*-
"""
***************************************************************************
UTM29NED50ToETR89PTTM06.py
---------------------
Date : July 2014
Copyright : (C) 2014 by Alexander Bruy
Email : alexander dot bruy at gmail dot 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; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************
"""
__author__ = 'Alexander Bruy'
__date__ = 'July 2014'
__copyright__ = '(C) 2014, Alexander Bruy'
# This will get replaced with a git SHA1 when you do a git archive
__revision__ = '$Format:%H$'
import inspect
import os
from PyQt4.QtGui import *
from qgis.core import *
from processing.gui.Help2Html import getHtmlFromRstFile
try:
from processing.parameters.ParameterVector import ParameterVector
from processing.outputs.OutputVector import OutputVector
except:
from processing.core.parameters import ParameterVector
from processing.core.outputs import OutputVector
from processing.core.GeoAlgorithm import GeoAlgorithm
from processing.algs.gdal.GdalUtils import GdalUtils
from processing.tools.vector import ogrConnectionString
from processing.tools.system import *
class UTM29NED50ToETR89PTTM06(GeoAlgorithm):
INPUT = 'INPUT'
OUTPUT = 'OUTPUT'
def getIcon(self):
return QIcon(os.path.dirname(__file__) + '/icons/pttransform.svg')
def help(self):
name = self.commandLineName().split(':')[1].lower()
filename = os.path.join(os.path.dirname(inspect.getfile(self.__class__)), 'help', name + '.rst')
try:
html = getHtmlFromRstFile(filename)
return True, html
except:
return False, None
def defineCharacteristics(self):
self.name = 'De UTM 29N ED50 para PT-TM06/ETRS89'
self.group = 'Transformação de Datum em Vectores'
self.addParameter(ParameterVector(self.INPUT, 'Ficheiro de Entrada',
[ParameterVector.VECTOR_TYPE_ANY]))
self.addOutput(OutputVector(self.OUTPUT, 'Ficheiro de Saída'))
def processAlgorithm(self, progress):
inLayer = self.getParameterValue(self.INPUT)
conn = ogrConnectionString(inLayer)
output = self.getOutputFromName(self.OUTPUT)
outFile = output.value
arguments = ['-f',
'ESRI Shapefile',
'-s_srs',
'+proj=utm +zone=29 +ellps=intl +nadgrids=' + os.path.dirname(__file__) + '/grids/ptED_e89.gsb +wktext +units=m +no_defs',
'-t_srs',
'EPSG:3763',
outFile,
conn
]
commands = ['ogr2ogr', GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)