-
Notifications
You must be signed in to change notification settings - Fork 5
/
spectral_angle.py
195 lines (157 loc) · 7.78 KB
/
spectral_angle.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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
# -*- coding: utf-8 -*-
"""
/***************************************************************************
QGIS_Edu
A QGIS plugin
Segmentation using OTB application
-------------------
begin : 2014-05-06
copyright : (C) 2014 by CNES
email : [email protected]
***************************************************************************/
/***************************************************************************
* *
* 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. *
* *
***************************************************************************/
"""
import os
from qgis.gui import QgsMapToolPan, QgsMapTool, QgsMessageBar
from qgis.core import QgsPoint, QgsMapLayerRegistry
from PyQt4 import QtCore, QtGui
# import project lib
# from ptmaptool import ProfiletoolMapTool
import terre_image_processing
# imports for display
import manage_QGIS
# import loggin for debug messages
import logging
logging.basicConfig()
# create logger
logger = logging.getLogger('TerreImage_SpectralAngle')
logger.setLevel(logging.DEBUG)
class SpectralAngleMapTool(QgsMapTool):
__pyqtSignals__ = ("canvas_clicked(PyQt_PyObject)")
def __init__(self, canvas):
# print "SpectralAngleMapTool"
QgsMapTool.__init__(self, canvas)
self.canvas = canvas
self.cursor = QtGui.QCursor(QtCore.Qt.CrossCursor)
def canvasReleaseEvent(self, event):
# print "canvasReleaseEvent before"
# print "event", event.pos()
self.emit(QtCore.SIGNAL("canvas_clicked"), {'x': event.pos().x(), 'y': event.pos().y()})
# print "canvasReleaseEvent after"
def activate(self):
QgsMapTool.activate(self)
self.canvas.setCursor(self.cursor)
def deactivate(self):
# self.emit( SIGNAL("deactivate") )
# self.canvas.setCursor( QCursor(Qt.ArrowCursor))
QgsMapTool.deactivate(self)
def setCursor(self, cursor):
self.cursor = QtGui.QCursor(cursor)
class SpectralAngle(QtCore.QObject):
__pyqtSignals__ = ("anglesComputed(PyQt_PyObject)")
def __init__(self, iface, working_dir, layer, mirror_map, rubberband):
QtCore.QObject.__init__(self)
self.iface = iface
self.canvas = self.iface.mapCanvas()
self.mirrormap_tool = mirror_map
self.working_directory = working_dir
if layer:
self.layer = layer
# self.tool = ProfiletoolMapTool(self.iface.mapCanvas()) #the mouselistener
# logger.debug( "self.tool" + str(self.tool))
self.pointstoDraw = None # Polyline in mapcanvas CRS analysed
self.maptool = self.canvas.mapTool()
# self.get_point_for_angles()
self.rubberband = rubberband
# self.polygon = False
# self.rubberband = QgsRubberBand(self.canvas, self.polygon)
# self.rubberband.setWidth(10)
# self.rubberband.setColor(QtGui.QColor(QtCore.Qt.yellow))
def angles(self, x, y):
widget = self.iface.messageBar().createMessage("Terre Image", "Travail en cours...")
self.iface.messageBar().pushWidget(widget, QgsMessageBar.INFO)
self.iface.mainWindow().statusBar().showMessage("Terre Image : Travail en cours...")
if self.layer:
image_output = terre_image_processing.angles(self.layer, self.working_directory, self.iface, x, y)
# self.tool.deactivate()
# self.deactivate()
self.canvas.setMapTool(self.maptool)
if image_output:
# self.rubberband.addPoint(QgsPoint(mapPos.x(),mapPos.y()))
self.emit(QtCore.SIGNAL("anglesComputed(PyQt_PyObject)"), image_output)
# self.display(image_output)
def freezeCanvas(self, setFreeze):
if setFreeze:
if not self.iface.mapCanvas().isFrozen():
self.iface.mapCanvas().freeze(True)
else:
if self.iface.mapCanvas().isFrozen():
self.iface.mapCanvas().freeze(False)
self.iface.mapCanvas().refresh()
def display(self, output_filename):
self.freezeCanvas(True)
# result_layer = manage_QGIS.get_raster_layer( output_filename, os.path.basename(os.path.splitext(self.layer.source_file)[0]) + "_" + self.processing_name )
result_layer = manage_QGIS.addRasterLayerToQGIS(output_filename, os.path.basename(os.path.splitext(output_filename)[0]) , self.iface)
manage_QGIS.histogram_stretching(result_layer, self.iface.mapCanvas())
self.output_layer = result_layer
# 2 ouvrir une nouvelle vue
self.mirror = self.mirrormap_tool.runDockableMirror("Angle Spectral")
logger.debug(self.mirror)
self.mirror.mainWidget.addLayer(result_layer.id())
# 1 mettre image en queue
ifaceLegend = self.iface.legendInterface()
ifaceLayers = QgsMapLayerRegistry.instance().mapLayers()
logger.debug("ifacelayers" + str(ifaceLayers))
id_layer = result_layer.id()
logger.debug("id_layer" + str(id_layer))
logger.debug("result layer" + str(result_layer))
# QgsMapLayerRegistry.instance().mapLayers()
# {u'QB_1_ortho20140521141641682': <qgis.core.QgsRasterLayer object at 0x6592b00>, u'QB_1_ortho_bande_bleue20140521141927295': <qgis.core.QgsRasterLayer object at 0x6592950>}
ifaceLegend.setLayerVisible(result_layer, False)
self.iface.mapCanvas().refresh()
logger.debug(ifaceLegend.isLayerVisible(result_layer))
# thaw the canvas
self.freezeCanvas(False)
def get_point_for_angles(self, layer):
# print "get point for angles"
self.tool = SpectralAngleMapTool(self.iface.mapCanvas()) # the mouselistener
# print "self.tool", self.tool
logger.debug("self.tool" + str(self.tool))
self.layer = layer
logger.debug("get point for angles")
# init the mouse listener comportement and save the classic to restore it on quit
self.canvas.setMapTool(self.tool)
QtCore.QObject.connect(self.tool, QtCore.SIGNAL("canvas_clicked"), self.onRightClick)
logger.debug("listener set")
# init the temp layer where the polyline is draw
# self.polygon = False
# self.rubberband = QgsRubberBand(self.canvas, self.polygon)
# self.rubberband.setWidth(10)
# self.rubberband.setColor(QtGui.QColor(QtCore.Qt.yellow))
# init the table where is saved the poyline
self.pointstoDraw = []
self.pointstoCal = []
def deactivate(self): # enable clean exit of the plugin
QtCore.QObject.disconnect(self.tool, QtCore.SIGNAL("canvas_clicked"), self.onRightClick)
self.rubberband.reset()
def onRightClick(self, position): # used to quit the current action
# print "right click", position
mapPos = self.canvas.getCoordinateTransform().toMapCoordinates(position["x"], position["y"])
newPoints = [[mapPos.x(), mapPos.y()]]
# if newPoints == self.lastClicked: return # sometimes a strange "double click" is given
# self.rubberband.reset(self.polygon)
# self.rubberband.reset(QGis.Point)
# print "setting rubberband", mapPos.x(),mapPos.y()
self.rubberband.addPoint(QgsPoint(mapPos.x(), mapPos.y()))
self.toolPan = QgsMapToolPan(self.canvas)
self.canvas.setMapTool(self.toolPan)
# create new vlayer ???
self.angles(mapPos.x(), mapPos.y())
self.tool = None