1 |
# Copyright (c) 2003 by Intevation GmbH |
# Copyright (c) 2003-2004 by Intevation GmbH |
2 |
# Authors: |
# Authors: |
3 |
# Jonathan Coles <[email protected]> |
# Jan-Oliver Wagner <[email protected]> (2004) |
4 |
|
# Bernhard Herzog <[email protected]> (2003) |
5 |
|
# Thomas K�ster <[email protected]> (2003) |
6 |
|
# Jonathan Coles <[email protected]> (2003) |
7 |
# |
# |
8 |
# This program is free software under the GPL (>=v2) |
# This program is free software under the GPL (>=v2) |
9 |
# Read the file COPYING coming with Thuban for details. |
# Read the file COPYING coming with Thuban for details. |
29 |
The resulting classification will consist of one group for each |
The resulting classification will consist of one group for each |
30 |
item in _list whose properties ramp between 'prop1' and 'prop2'. |
item in _list whose properties ramp between 'prop1' and 'prop2'. |
31 |
|
|
32 |
_list -- any object that implements the iterator interface |
_list -- a list of values for each singleton |
33 |
|
|
34 |
ramp -- an object which implements the CustomRamp interface |
ramp -- an object which implements the CustomRamp interface |
35 |
""" |
""" |
38 |
|
|
39 |
i = 0 |
i = 0 |
40 |
maxValue = float(len(_list) - 1) |
maxValue = float(len(_list) - 1) |
41 |
|
if maxValue < 1: maxValue = 1 |
42 |
|
|
43 |
for value in _list: |
for value in _list: |
44 |
prop = ramp.GetProperties(i / maxValue) |
prop = ramp.GetProperties(i / maxValue) |
45 |
clazz.AppendGroup(ClassGroupSingleton(value, prop)) |
clazz.AppendGroup(ClassGroupSingleton(value, prop)) |
63 |
|
|
64 |
end = "[" |
end = "[" |
65 |
maxValue = float(numGroups - 1) |
maxValue = float(numGroups - 1) |
66 |
|
if maxValue < 1: maxValue = 1 |
67 |
|
|
68 |
for i in range(1, numGroups + 1): |
for i in range(1, numGroups + 1): |
69 |
|
|
70 |
prop = ramp.GetProperties(float(i-1) / maxValue) |
prop = ramp.GetProperties(float(i-1) / maxValue) |
71 |
|
|
72 |
if intStep: |
if intStep: |
73 |
cur_max = min + int(round((i * (max - min + 1)) / maxValue)) |
cur_max = min + int(round((i * (max - min + 1)) / float(numGroups))) |
74 |
else: |
else: |
75 |
cur_max = min + (i * (max - min)) / maxValue |
cur_max = min + (i * (max - min)) / float(numGroups) |
76 |
|
|
77 |
if i == numGroups: |
if i == numGroups: |
78 |
cur_max = max |
cur_max = max |
132 |
end = "]" |
end = "]" |
133 |
|
|
134 |
maxValue = float(numGroups - 1) |
maxValue = float(numGroups - 1) |
135 |
|
if maxValue < 1: maxValue = 1 |
136 |
for (q, p) in quantiles[3]: |
for (q, p) in quantiles[3]: |
137 |
|
|
138 |
prop = ramp.GetProperties(float(i-1) / maxValue) |
prop = ramp.GetProperties(float(i-1) / maxValue) |
144 |
max = _list[q] |
max = _list[q] |
145 |
|
|
146 |
group = ClassGroupRange(Range((start, min, max, end)), prop) |
group = ClassGroupRange(Range((start, min, max, end)), prop) |
|
|
|
|
group.SetLabel("%s%% - %s%%" % (round(oldp*100, 2), |
|
|
round(p*100, 2))) |
|
|
oldp = p |
|
|
start = "]" |
|
|
min = max |
|
|
clazz.AppendGroup(group) |
|
|
i += 1 |
|
|
|
|
|
return (adjusted, clazz) |
|
147 |
|
|
|
def GenQuantiles0(_list, percents, ramp, _range): |
|
|
"""Same as GenQuantiles, but the first class won't be added to |
|
|
the classification. |
|
|
|
|
|
Returns a tuple (adjusted, Classification, upper_class0) where |
|
|
upper_class0 is the highest value inside the first class. |
|
|
|
|
|
_list -- a sort list of values |
|
|
|
|
|
percents -- a sorted list of floats in the range 0.0-1.0 which |
|
|
represent the upper bound of each quantile. the |
|
|
union of all percentiles should be the entire |
|
|
range from 0.0-1.0 |
|
|
|
|
|
ramp -- an object which implements the CustomRamp interface |
|
|
|
|
|
_range -- a Range object |
|
|
|
|
|
Raises a Value Error if 'percents' has fewer than two items, or |
|
|
does not cover the entire range. |
|
|
""" |
|
|
|
|
|
clazz = Classification() |
|
|
quantiles = calculate_quantiles(_list, percents, _range) |
|
|
adjusted = True |
|
|
|
|
|
if quantiles is not None: |
|
|
|
|
|
numGroups = len(quantiles[3]) - 1 |
|
|
|
|
|
if numGroups > 0: |
|
|
adjusted = quantiles[0] |
|
|
|
|
|
start, min, endMax, right = _range.GetRange() |
|
|
|
|
|
class0 = quantiles[3][0] |
|
|
min = _list[class0[0]] |
|
|
oldp = class0[1] |
|
|
i = 1 |
|
|
end = "]" |
|
|
|
|
|
maxValue = float(numGroups - 1) |
|
|
for (q, p) in quantiles[3][1:]: |
|
|
prop = ramp.GetProperties(float(i-1) / maxValue) |
|
|
|
|
|
if i == numGroups: |
|
|
max = endMax |
|
|
end = right |
|
|
else: |
|
|
max = _list[q] |
|
|
|
|
|
group = ClassGroupRange(Range((start, min, max, end)), prop) |
|
|
|
|
148 |
group.SetLabel("%s%% - %s%%" % (round(oldp*100, 2), |
group.SetLabel("%s%% - %s%%" % (round(oldp*100, 2), |
149 |
round(p*100, 2))) |
round(p*100, 2))) |
150 |
oldp = p |
oldp = p |
153 |
clazz.AppendGroup(group) |
clazz.AppendGroup(group) |
154 |
i += 1 |
i += 1 |
155 |
|
|
156 |
return (adjusted, clazz, _list[class0[0]]) |
return (adjusted, clazz) |
157 |
|
|
158 |
|
|
159 |
def calculate_quantiles(_list, percents, _range): |
def calculate_quantiles(_list, percents, _range): |
160 |
"""Calculate quantiles for the given _list of percents from the |
"""Calculate quantiles for the given _list of percents from the |
161 |
sorted list of values that are in range. |
sorted list of values that are in range. |
162 |
|
|
163 |
This may not actually generate len(percents) quantiles if |
This may not actually generate len(percents) quantiles if |
164 |
many of the values that fall on quantile borders are the same. |
many of the values that fall on quantile borders are the same. |
165 |
|
|
316 |
class CustomRamp: |
class CustomRamp: |
317 |
|
|
318 |
def __init__(self, prop1, prop2): |
def __init__(self, prop1, prop2): |
319 |
|
"""Create a ramp between prop1 and prop2.""" |
320 |
self.prop1 = prop1 |
self.prop1 = prop1 |
321 |
self.prop2 = prop2 |
self.prop2 = prop2 |
322 |
|
|
323 |
def GetRamp(self): |
def GetRamp(self): |
324 |
|
"""Return this ramp.""" |
325 |
return self |
return self |
326 |
|
|
327 |
def GetProperties(self, index): |
def GetProperties(self, index): |
337 |
|
|
338 |
newProps = ClassGroupProperties() |
newProps = ClassGroupProperties() |
339 |
|
|
340 |
color1 = self.prop1.GetLineColor() |
self.__SetProperty(self.prop1.GetLineColor(), |
341 |
color2 = self.prop2.GetLineColor() |
self.prop2.GetLineColor(), |
342 |
|
index, newProps.SetLineColor) |
343 |
self.__SetProperty(color1, color2, index, newProps.SetLineColor) |
self.__SetProperty(self.prop1.GetFill(), self.prop2.GetFill(), |
344 |
self.__SetProperty(color1, color2, index, newProps.SetFill) |
index, newProps.SetFill) |
345 |
|
|
346 |
w = (self.prop2.GetLineWidth() - self.prop1.GetLineWidth()) \ |
w = (self.prop2.GetLineWidth() - self.prop1.GetLineWidth()) \ |
347 |
* index \ |
* index \ |
348 |
+ self.prop1.GetLineWidth() |
+ self.prop1.GetLineWidth() |
|
|
|
349 |
newProps.SetLineWidth(int(round(w))) |
newProps.SetLineWidth(int(round(w))) |
350 |
|
|
351 |
|
s = (self.prop2.GetSize() - self.prop1.GetSize()) \ |
352 |
|
* index \ |
353 |
|
+ self.prop1.GetSize() |
354 |
|
newProps.SetSize(int(round(s))) |
355 |
|
|
356 |
return newProps |
return newProps |
357 |
|
|
358 |
def __SetProperty(self, color1, color2, index, setf): |
def __SetProperty(self, color1, color2, index, setf): |
359 |
|
"""Use setf to set the appropriate property for the point |
360 |
|
index percent between color1 and color2. setf is a function |
361 |
|
to call that accepts a Color object or Transparent. |
362 |
|
""" |
363 |
|
|
364 |
if color1 is Transparent and color2 is Transparent: |
if color1 is Transparent and color2 is Transparent: |
365 |
setf(Transparent) |
setf(Transparent) |
380 |
(color2.blue - color1.blue) * index + color1.blue)) |
(color2.blue - color1.blue) * index + color1.blue)) |
381 |
|
|
382 |
class MonochromaticRamp(CustomRamp): |
class MonochromaticRamp(CustomRamp): |
383 |
|
"""Helper class to make ramps between two colors.""" |
384 |
|
|
385 |
def __init__(self, start, end): |
def __init__(self, start, end): |
386 |
|
"""Create a Monochromatic Ramp. |
387 |
|
|
388 |
|
start -- starting Color |
389 |
|
|
390 |
|
end -- ending Color |
391 |
|
""" |
392 |
sp = ClassGroupProperties() |
sp = ClassGroupProperties() |
393 |
sp.SetLineColor(start) |
sp.SetLineColor(start) |
394 |
sp.SetFill(start) |
sp.SetFill(start) |
399 |
|
|
400 |
CustomRamp.__init__(self, sp, ep) |
CustomRamp.__init__(self, sp, ep) |
401 |
|
|
402 |
GreyRamp = MonochromaticRamp(Color(1, 1, 1), Color(0, 0, 0)) |
grey_ramp = MonochromaticRamp(Color(1, 1, 1), Color(0, 0, 0)) |
403 |
RedRamp = MonochromaticRamp(Color(1, 1, 1), Color(.8, 0, 0)) |
red_ramp = MonochromaticRamp(Color(1, 1, 1), Color(.8, 0, 0)) |
404 |
GreenRamp = MonochromaticRamp(Color(1, 1, 1), Color(0, .8, 0)) |
green_ramp = MonochromaticRamp(Color(1, 1, 1), Color(0, .8, 0)) |
405 |
BlueRamp = MonochromaticRamp(Color(1, 1, 1), Color(0, 0, .8)) |
blue_ramp = MonochromaticRamp(Color(1, 1, 1), Color(0, 0, .8)) |
406 |
GreenToRedRamp = MonochromaticRamp(Color(1, .8, 1), Color(1, 0, 0)) |
green_to_red_ramp = MonochromaticRamp(Color(0, .8, 0), Color(1, 0, 0)) |
407 |
|
|
408 |
class HotToColdRamp: |
class HotToColdRamp: |
409 |
|
"""A ramp that generates properties with colors ranging from |
410 |
|
'hot' colors (e.g. red, orange) to 'cold' colors (e.g. green, blue) |
411 |
|
""" |
412 |
|
|
413 |
def GetRamp(self): |
def GetRamp(self): |
414 |
|
"""Return this ramp.""" |
415 |
return self |
return self |
416 |
|
|
417 |
def GetProperties(self, index): |
def GetProperties(self, index): |
442 |
|
|
443 |
return prop |
return prop |
444 |
|
|
445 |
|
class FixedRamp: |
446 |
|
"""FixedRamp allows particular properties of a ramp to be |
447 |
|
held constant over the ramp. |
448 |
|
""" |
449 |
|
|
450 |
|
def __init__(self, ramp, fixes): |
451 |
|
""" |
452 |
|
ramp -- a source ramp to get the default properties |
453 |
|
|
454 |
|
fixes -- a tuple (lineColor, lineWidth, fillColor) such that |
455 |
|
if any item is not None, the appropriate property will |
456 |
|
be fixed to that item value. |
457 |
|
""" |
458 |
|
|
459 |
|
self.fixes = fixes |
460 |
|
self.ramp = ramp |
461 |
|
|
462 |
|
def GetRamp(self): |
463 |
|
"""Return this ramp.""" |
464 |
|
return self |
465 |
|
|
466 |
|
def GetProperties(self, index): |
467 |
|
"""Return a ClassGroupProperties object whose properties |
468 |
|
represent a point at 'index' between the properties in |
469 |
|
the ramp that initialized this FixedRamp. |
470 |
|
|
471 |
|
index -- a value such that 0 <= index <= 1 |
472 |
|
""" |
473 |
|
|
474 |
|
props = self.ramp.GetProperties(index) |
475 |
|
if self.fixes[0] is not None: props.SetLineColor(self.fixes[0]) |
476 |
|
if self.fixes[1] is not None: props.SetLineWidth(self.fixes[1]) |
477 |
|
if self.fixes[2] is not None: props.SetFill(self.fixes[2]) |
478 |
|
|
479 |
|
return props |