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
195
196
197
198
|
#!/usr/bin/env python
# -*- coding: utf-8 -*-
# Copyright (c) 2009, 2013.
# SMHI,
# Folkborgsvägen 1,
# Norrköping,
# Sweden
# Author(s):
# Martin Raspaud <martin.raspaud@smhi.se>
# Adam Dybbroe <adam.dybbroe@smhi.se>
# This file is part of the mpop.
# mpop 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 3 of the License, or
# (at your option) any later version.
# mpop is distributed in the hope that it will be useful, but
# WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
# General Public License for more details.
# You should have received a copy of the GNU General Public License
# along with mpop. If not, see <http://www.gnu.org/licenses/>.
"""Palette holder module.
"""
def tv_legend():
"""Palette for TV.
"""
legend = []
legend.append(( 0, 0, 0)) # Unprocessed: Black
legend.append(( 0, 120, 0)) # Land
legend.append(( 0, 0, 215)) # Sea: Blue
legend.append(( 0, 120, 0)) # Land (Snow on land)
legend.append(( 0, 0, 215)) # Sea: Blue (Snow/Ice on sea)
for i in range(5, 256):
# All other pixel values are grey according to IR temp.
legend.append((i, i, i))
return convert_palette(legend)
def vv_legend():
"""Palette for Swedish road authorities (Vägverket).
"""
legend = []
legend.append(( 0, 0, 0)) # Unprocessed: Black
legend.append(( 0, 120, 0)) # Land
legend.append(( 0, 0, 215)) # Sea: Blue
# Cloud type values 5 to 8:
legend.append((255, 150, 0)) # Very low cumuliform
legend.append((255, 100, 0)) # Very low
legend.append((255, 220, 0)) # Low cumuliform
legend.append((255, 180, 0)) # Low
for i in range(7, 256):
# All other pixel values are grey according to IR temp.
legend.append((i, i, i))
return convert_palette(legend)
def cms_modified():
"""Palette for regular cloud classification.
"""
return nwcsaf_cloudtype()
def nwcsaf_cloudtype():
"""Palette for regular cloud classification.
"""
legend = []
legend.append((100, 100, 100)) # Unprocessed: Grey
legend.append(( 0, 120, 0))
legend.append(( 0, 0, 0)) # Sea: Black
legend.append((250, 190, 250)) # Snow
legend.append((220, 160, 220)) # Sea-ice
legend.append((255, 150, 0)) # Very low cumuliform
legend.append((255, 100, 0)) # Very low
legend.append((255, 220, 0)) # Low cumuliform
legend.append((255, 180, 0)) # Low
legend.append((255, 255, 140)) # Medium cumuliform
legend.append((240, 240, 0)) # Medium
legend.append((250, 240, 200)) # High cumiliform
legend.append((215, 215, 150)) # High
legend.append((255, 255, 255)) # Very high cumuliform
legend.append((230, 230, 230)) # Very high
legend.append(( 0, 80, 215)) # Semi-transparent thin
legend.append(( 0, 180, 230)) # Semi-transparent medium
legend.append(( 0, 240, 240)) # Semi-transparent thick
legend.append(( 90, 200, 160)) # Semi-transparent above
legend.append((200, 0, 200)) # Broken
legend.append(( 95, 60, 30)) # Undefined: Brown
return convert_palette(legend)
def ctth_height():
"""CTTH height palette.
"""
legend = []
legend.append((0, 0, 0))
legend.append((255, 0, 216)) # 0 meters
legend.append((126, 0, 43))
legend.append((153, 20, 47))
legend.append((178, 51, 0))
legend.append((255, 76, 0))
legend.append((255, 102, 0))
legend.append((255, 164, 0))
legend.append((255, 216, 0))
legend.append((216, 255, 0))
legend.append((178, 255, 0))
legend.append((153, 255, 0))
legend.append((0, 255, 0))
legend.append((0, 140, 48))
legend.append((0, 178, 255))
legend.append((0, 216, 255))
legend.append((0, 255, 255))
legend.append((238, 214, 210))
legend.append((239, 239, 223))
legend.append((255, 255, 255)) # 10,000 meters
for i in range(79):
legend.append((255, 255, 255))
legend.append((224,224,224))
return convert_palette(legend)
def ctth_height_pps():
"""CTTH height palette for NWCSAF/PPS.
Identical to the one found in the hdf5 files.
"""
legend = []
legend.append((255, 0, 216)) # 0 meters
legend.append((255, 0, 216)) # 0 meters
legend.append((255, 0, 216)) # 0 meters
legend.append((126, 0, 43))
legend.append((126, 0, 43))
legend.append((153, 20, 47))
legend.append((153, 20, 47))
legend.append((153, 20, 47))
legend.append((178, 51, 0))
legend.append((178, 51, 0))
legend.append((255, 76, 0))
legend.append((255, 76, 0))
legend.append((255, 76, 0))
legend.append((255, 102, 0))
legend.append((255, 102, 0))
legend.append((255, 164, 0))
legend.append((255, 164, 0))
legend.append((255, 164, 0))
legend.append((255, 216, 0))
legend.append((255, 216, 0))
legend.append((216, 255, 0))
legend.append((216, 255, 0))
legend.append((178, 255, 0))
legend.append((178, 255, 0))
legend.append((178, 255, 0))
legend.append((153, 255, 0))
legend.append((153, 255, 0))
legend.append((0, 255, 0))
legend.append((0, 255, 0))
legend.append((0, 255, 0))
legend.append((0, 140, 48))
legend.append((0, 140, 48))
legend.append((0, 178, 255))
legend.append((0, 178, 255))
legend.append((0, 178, 255))
legend.append((0, 216, 255))
legend.append((0, 216, 255))
legend.append((0, 255, 255))
legend.append((0, 255, 255))
legend.append((0, 255, 255))
legend.append((238, 214, 210))
legend.append((238, 214, 210))
legend.append((239, 239, 223))
legend.append((239, 239, 223))
for idx in range(47, 150):
legend.append((255, 255, 255)) # 10,000 meters
for idx in range(150, 256):
legend.append((0, 0, 0))
return convert_palette(legend)
def convert_palette(palette):
"""Convert palette from [0,255] range to [0,1].
"""
new_palette = []
for i in palette:
new_palette.append((i[0] / 255.0,
i[1] / 255.0,
i[2] / 255.0))
return new_palette
|
