-
Notifications
You must be signed in to change notification settings - Fork 5
/
Datum73MilToETR89PTTM06_Raster.py
100 lines (82 loc) · 4.12 KB
/
Datum73MilToETR89PTTM06_Raster.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
# -*- coding: latin-1 -*-
"""
***************************************************************************
Datum73MilToETR89PTTM06_Raster.py
---------------------
Date : July 2014
Copyright : (C) 2014 by Pedro Venâncio, Giovanni Manghi
Email : pedrongvenancio at gmail dot com
giovanni dot manghi 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__ = 'Pedro Venâncio, Giovanni Manghi'
__date__ = 'July 2014'
__copyright__ = '(C) 2014, Pedro Venâncio, Giovanni Manghi'
# 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.ParameterRaster import ParameterRaster
from processing.parameters.ParameterSelection import ParameterSelection
from processing.outputs.OutputRaster import OutputRaster
except:
from processing.core.parameters import ParameterRaster
from processing.core.parameters import ParameterSelection
from processing.core.outputs import OutputRaster
from processing.core.GeoAlgorithm import GeoAlgorithm
from processing.algs.gdal.GdalUtils import GdalUtils
class Datum73MilToETR89PTTM06_Raster(GeoAlgorithm):
INPUT = 'INPUT'
OUTPUT = 'OUTPUT'
GRID = 'GRELHAS'
GRID_OPTIONS = ['José Alberto Gonçalves',
'Direção-Geral do Territorio']
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 Datum 73 Militar para PT-TM06/ETRS89 Raster'
self.group = 'Transformação de Datum em Rasters'
self.addParameter(ParameterRaster(self.INPUT, 'Ficheiro de Entrada', False))
self.addParameter(ParameterSelection(self.GRID, 'Grelha NTv2 a usar (Fonte)',
self.GRID_OPTIONS))
self.addOutput(OutputRaster(self.OUTPUT, 'Ficheiro de Saída'))
def processAlgorithm(self, progress):
arguments = ['-s_srs']
if self.getParameterValue(self.GRID) == 0:
# Jose Alberto Goncalves
arguments.append('+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=200180.598 +y_0=299913.01 +ellps=intl +nadgrids=' + os.path.dirname(__file__) + '/grids/pt73_e89.gsb +wktext +units=m +no_defs')
else:
# Direccao Geral do Territorio
arguments.append('+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=200180.598 +y_0=299913.01 +ellps=intl +nadgrids=' + os.path.dirname(__file__) + '/grids/D73_ETRS89_geo.gsb +wktext +units=m +no_defs')
arguments.append('-t_srs')
arguments.append('EPSG:3763')
arguments.append('-r')
arguments.append('bilinear')
arguments.append('-dstnodata')
arguments.append('nan')
arguments.append('-of')
out = self.getOutputValue(self.OUTPUT)
arguments.append(GdalUtils.getFormatShortNameFromFilename(out))
arguments.append(self.getParameterValue(self.INPUT))
arguments.append(out)
GdalUtils.runGdal(['gdalwarp', GdalUtils.escapeAndJoin(arguments)],
progress)