|
18 | 18 | from PyQt4.QtCore import *
|
19 | 19 | from qgis.core import *
|
20 | 20 | from qgis.analysis import *
|
21 |
| -from utilities import (unitTestDataPath, |
22 |
| - getQgisTestApp, |
23 |
| - TestCase, |
24 |
| - unittest, |
25 |
| - ) |
| 21 | +from utilities import ( |
| 22 | + unitTestDataPath, |
| 23 | + getQgisTestApp, |
| 24 | + TestCase, |
| 25 | + unittest) |
26 | 26 |
|
27 | 27 | QGISAPP, CANVAS, IFACE, PARENT = getQgisTestApp()
|
28 | 28 |
|
| 29 | + |
29 | 30 | class TestQgsZonalStatistics(TestCase):
|
| 31 | + """Tests for zonal stats class.""" |
30 | 32 |
|
31 |
| - def testStatistics(self): |
32 |
| - sep = os.sep |
33 |
| - TEST_DATA_DIR = unitTestDataPath() + sep + "zonalstatistics" + sep |
34 |
| - myTempPath = QDir.tempPath() + sep |
35 |
| - testDir = QDir(TEST_DATA_DIR) |
36 |
| - for f in testDir.entryList(QDir.Files): |
37 |
| - QFile.remove(myTempPath + f) |
38 |
| - QFile.copy(TEST_DATA_DIR + f, myTempPath + f) |
| 33 | + def testStatistics(self): |
| 34 | + """Test zonal stats""" |
| 35 | + sep = os.sep |
| 36 | + TEST_DATA_DIR = unitTestDataPath() + sep + "zonalstatistics" + sep |
| 37 | + myTempPath = QDir.tempPath() + sep |
| 38 | + testDir = QDir(TEST_DATA_DIR) |
| 39 | + for f in testDir.entryList(QDir.Files): |
| 40 | + QFile.remove(myTempPath + f) |
| 41 | + QFile.copy(TEST_DATA_DIR + f, myTempPath + f) |
39 | 42 |
|
40 |
| - myVector = QgsVectorLayer(myTempPath + "polys.shp", "poly", "ogr") |
41 |
| - myRasterPath = myTempPath + "edge_problem.asc" |
42 |
| - zs = QgsZonalStatistics( myVector, myRasterPath, "", 1 ) |
43 |
| - zs.calculateStatistics(None) |
| 43 | + myVector = QgsVectorLayer(myTempPath + "polys.shp", "poly", "ogr") |
| 44 | + myRasterPath = myTempPath + "edge_problem.asc" |
| 45 | + zs = QgsZonalStatistics(myVector, myRasterPath, "", 1) |
| 46 | + zs.calculateStatistics(None) |
44 | 47 |
|
45 |
| - feat = QgsFeature() |
46 |
| - # validate statistics for each feature |
47 |
| - request = QgsFeatureRequest().setFilterFid(0) |
48 |
| - feat = myVector.getFeatures(request).next() |
49 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (12.0, feat[1])) |
50 |
| - assert feat[1] == 12.0, myMessage |
51 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (8.0, feat[2])) |
52 |
| - assert feat[2] == 8.0, myMessage |
53 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (0.666666666666667, feat[3])) |
54 |
| - assert abs(feat[3] - 0.666666666666667 ) < 0.00001, myMessage |
| 48 | + feat = QgsFeature() |
| 49 | + # validate statistics for each feature |
| 50 | + request = QgsFeatureRequest().setFilterFid(0) |
| 51 | + feat = myVector.getFeatures(request).next() |
| 52 | + myMessage = ('Expected: %f\nGot: %f\n' % (12.0, feat[1])) |
| 53 | + assert feat[1] == 12.0, myMessage |
| 54 | + myMessage = ('Expected: %f\nGot: %f\n' % (8.0, feat[2])) |
| 55 | + assert feat[2] == 8.0, myMessage |
| 56 | + myMessage = ('Expected: %f\nGot: %f\n' % (0.666666666666667, feat[3])) |
| 57 | + assert abs(feat[3] - 0.666666666666667) < 0.00001, myMessage |
55 | 58 |
|
56 |
| - request.setFilterFid(1) |
57 |
| - feat = myVector.getFeatures(request).next() |
58 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (9.0, feat[1])) |
59 |
| - assert feat[1] == 9.0, myMessage |
60 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (5.0, feat[2])) |
61 |
| - assert feat[2] == 5.0, myMessage |
62 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (0.555555555555556, feat[3])) |
63 |
| - assert abs( feat[3] - 0.555555555555556) < 0.00001, myMessage |
| 59 | + request.setFilterFid(1) |
| 60 | + feat = myVector.getFeatures(request).next() |
| 61 | + myMessage = ('Expected: %f\nGot: %f\n' % (9.0, feat[1])) |
| 62 | + assert feat[1] == 9.0, myMessage |
| 63 | + myMessage = ('Expected: %f\nGot: %f\n' % (5.0, feat[2])) |
| 64 | + assert feat[2] == 5.0, myMessage |
| 65 | + myMessage = ('Expected: %f\nGot: %f\n' % (0.555555555555556, feat[3])) |
| 66 | + assert abs(feat[3] - 0.555555555555556) < 0.00001, myMessage |
64 | 67 |
|
65 |
| - request.setFilterFid(2) |
66 |
| - feat = myVector.getFeatures(request).next() |
67 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (6.0, feat[1])) |
68 |
| - assert feat[1] == 6.0, myMessage |
69 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (5.0, feat[2])) |
70 |
| - assert feat[2] == 5.0, myMessage |
71 |
| - myMessage = ('Expected: %f\nGot: %f\n' % (0.833333333333333, feat[3])) |
72 |
| - assert abs( feat[3] - 0.833333333333333 ) < 0.00001, myMessage |
| 68 | + request.setFilterFid(2) |
| 69 | + feat = myVector.getFeatures(request).next() |
| 70 | + myMessage = ('Expected: %f\nGot: %f\n' % (6.0, feat[1])) |
| 71 | + assert feat[1] == 6.0, myMessage |
| 72 | + myMessage = ('Expected: %f\nGot: %f\n' % (5.0, feat[2])) |
| 73 | + assert feat[2] == 5.0, myMessage |
| 74 | + myMessage = ('Expected: %f\nGot: %f\n' % (0.833333333333333, feat[3])) |
| 75 | + assert abs(feat[3] - 0.833333333333333) < 0.00001, myMessage |
73 | 76 |
|
74 | 77 | if __name__ == '__main__':
|
75 | 78 | unittest.main()
|
0 commit comments