2013-02-16 00:23:56 +01:00
|
|
|
# -*- coding: utf-8 -*-
|
|
|
|
|
|
|
|
|
|
"""
|
|
|
|
|
***************************************************************************
|
|
|
|
|
MeanAndStdDevPlot.py
|
|
|
|
|
---------------------
|
|
|
|
|
Date : January 2013
|
|
|
|
|
Copyright : (C) 2013 by Victor Olaya
|
|
|
|
|
Email : volayaf 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__ = 'Victor Olaya'
|
|
|
|
|
__date__ = 'January 2013'
|
|
|
|
|
__copyright__ = '(C) 2013, Victor Olaya'
|
2013-10-01 20:52:22 +03:00
|
|
|
|
2013-02-16 00:23:56 +01:00
|
|
|
# This will get replaced with a git SHA1 when you do a git archive
|
2013-10-01 20:52:22 +03:00
|
|
|
|
2013-02-16 00:23:56 +01:00
|
|
|
__revision__ = '$Format:%H$'
|
|
|
|
|
|
|
|
|
|
import matplotlib.pyplot as plt
|
|
|
|
|
import matplotlib.pylab as lab
|
|
|
|
|
import numpy as np
|
2015-01-17 14:57:21 +02:00
|
|
|
|
2013-10-01 20:52:22 +03:00
|
|
|
from processing.core.GeoAlgorithm import GeoAlgorithm
|
2014-07-14 14:19:09 +02:00
|
|
|
from processing.core.parameters import ParameterTable
|
|
|
|
|
from processing.core.parameters import ParameterTableField
|
|
|
|
|
from processing.core.outputs import OutputHTML
|
2015-01-17 14:57:21 +02:00
|
|
|
|
|
|
|
|
from processing.tools import vector
|
|
|
|
|
from processing.tools import dataobjects
|
2013-02-16 00:23:56 +01:00
|
|
|
|
|
|
|
|
|
2013-10-01 20:52:22 +03:00
|
|
|
class MeanAndStdDevPlot(GeoAlgorithm):
|
2013-02-16 00:23:56 +01:00
|
|
|
|
2013-10-01 20:52:22 +03:00
|
|
|
INPUT = 'INPUT'
|
|
|
|
|
OUTPUT = 'OUTPUT'
|
|
|
|
|
NAME_FIELD = 'NAME_FIELD'
|
|
|
|
|
MEAN_FIELD = 'MEAN_FIELD'
|
|
|
|
|
STDDEV_FIELD = 'STDDEV_FIELD'
|
2013-02-16 00:23:56 +01:00
|
|
|
|
2015-01-17 14:57:21 +02:00
|
|
|
def defineCharacteristics(self):
|
2015-07-26 03:48:27 +02:00
|
|
|
self.name, self.i18n_name = self.trAlgorithm('Mean and standard deviation plot')
|
|
|
|
|
self.group, self.i18n_group = self.trAlgorithm('Graphics')
|
2015-01-17 14:57:21 +02:00
|
|
|
|
|
|
|
|
self.addParameter(ParameterTable(self.INPUT,
|
2015-08-22 14:29:41 +02:00
|
|
|
self.tr('Input table')))
|
2015-01-17 14:57:21 +02:00
|
|
|
self.addParameter(ParameterTableField(self.NAME_FIELD,
|
2015-08-22 14:29:41 +02:00
|
|
|
self.tr('Category name field'), self.INPUT,
|
|
|
|
|
ParameterTableField.DATA_TYPE_ANY))
|
2015-01-17 14:57:21 +02:00
|
|
|
self.addParameter(ParameterTableField(self.MEAN_FIELD,
|
2015-08-22 14:29:41 +02:00
|
|
|
self.tr('Mean field'), self.INPUT))
|
2015-01-17 14:57:21 +02:00
|
|
|
self.addParameter(ParameterTableField(self.STDDEV_FIELD,
|
2015-08-22 14:29:41 +02:00
|
|
|
self.tr('StdDev field'), self.INPUT))
|
2015-01-17 14:57:21 +02:00
|
|
|
|
2015-05-22 15:56:45 +02:00
|
|
|
self.addOutput(OutputHTML(self.OUTPUT, self.tr('Plot')))
|
2015-01-17 14:57:21 +02:00
|
|
|
|
2013-02-16 00:23:56 +01:00
|
|
|
def processAlgorithm(self, progress):
|
2015-01-17 14:57:21 +02:00
|
|
|
layer = dataobjects.getObjectFromUri(
|
|
|
|
|
self.getParameterValue(self.INPUT))
|
2013-02-16 00:23:56 +01:00
|
|
|
namefieldname = self.getParameterValue(self.NAME_FIELD)
|
|
|
|
|
meanfieldname = self.getParameterValue(self.MEAN_FIELD)
|
|
|
|
|
stddevfieldname = self.getParameterValue(self.STDDEV_FIELD)
|
2015-01-17 14:57:21 +02:00
|
|
|
|
2013-02-16 00:23:56 +01:00
|
|
|
output = self.getOutputValue(self.OUTPUT)
|
2013-02-28 22:08:32 +01:00
|
|
|
|
2015-01-17 14:57:21 +02:00
|
|
|
values = vector.values(layer, namefieldname, meanfieldname, stddevfieldname)
|
|
|
|
|
plt.close()
|
2013-02-16 00:23:56 +01:00
|
|
|
ind = np.arange(len(values[namefieldname]))
|
|
|
|
|
width = 0.8
|
2013-10-01 20:52:22 +03:00
|
|
|
plt.bar(ind, values[meanfieldname], width, color='r',
|
|
|
|
|
yerr=values[stddevfieldname],
|
|
|
|
|
error_kw=dict(ecolor='yellow'),
|
fix python pep8 warnings and fix some revealed errors
pep8 --ignore=E111,E128,E201,E202,E203,E211,E221,E222,E225,E226,E227,E231,E241,E261,E265,E272,E302,E303,E501,E701 \
--exclude="ui_*.py,debian/*,python/ext-libs/*" \
.
2015-02-01 14:15:42 +01:00
|
|
|
)
|
2013-02-16 00:23:56 +01:00
|
|
|
|
2013-10-01 20:52:22 +03:00
|
|
|
plt.xticks(ind, values[namefieldname], rotation=45)
|
|
|
|
|
plotFilename = output + '.png'
|
2013-02-28 22:08:32 +01:00
|
|
|
lab.savefig(plotFilename)
|
2013-10-01 20:52:22 +03:00
|
|
|
f = open(output, 'w')
|
2015-04-21 10:31:14 +02:00
|
|
|
f.write('<html><img src="' + plotFilename + '"/></html>')
|
2013-02-16 00:23:56 +01:00
|
|
|
f.close()
|
