Бокплоттинг нескольких замаскированных массивов - PullRequest
0 голосов
/ 26 февраля 2019

Я пытаюсь построить несколько массивов, которые замаскировали все значения -9999 (значения неисправности оборудования) в моем наборе данных.Кажется, что np.boxplot игнорирует маску значений и строит их в любом случае.Как я могу получить недоумение, чтобы исключить маскированные значения в каждом наборе данных?

import numpy as np
import matplotlib.pyplot as plt

a = np.ma.array([229.5,374.0,536.5,-9999,-9999,-9999,-9999,-9999,182.0,42.5,49.0])
b = np.ma.array([363.0,118.5,159.0,-9999,311.0,516.0,380.0,338.5,223.0,211.5,128.5])
c = np.ma.array([205.5,277.5,141.5,278.0,302.0,251.0,299.0,250.0,315.5,92.1,211.9])

a = np.ma.masked_less( a, -9000 )
b = np.ma.masked_less( b, -9000 )
c = np.ma.masked_less( c, -9000 )

data_to_plot = [ a , b, c ]
labels = [ 'a', 'c', 'c' ]

fig = plt.figure( 3, figsize = ( 8, 10 ) )
ax = fig.add_subplot( 111, frameon = False )
bp = ax.boxplot( data_to_plot, patch_artist = True, widths = .85, 
                labels = labels, vert = True)
plt.setp( bp[ 'boxes' ], color = 'black', alpha = .8, linewidth = 3 )
plt.setp( bp[ 'whiskers' ], color = 'black', linewidth = 3 )
plt.setp( bp[ 'fliers' ], markeredgecolor = 'black', alpha = 1, markersize = 30, marker = '.' )
plt.setp( bp[ 'medians' ], color = 'orange', linewidth = 3 )
plt.setp( bp[ 'means' ], color = 'black' )
plt.setp( bp[ 'caps' ], color = 'black', linewidth = 4 )

colors = ['r', 'b', 'g', ]
for b in ( bp ):
    for patch, color in zip( bp[ 'boxes' ], colors ):
        patch.set_facecolor( color )
plt.show()

boxplot of masked -9999 data

Ответы [ 2 ]

0 голосов
/ 28 февраля 2019

https://github.com/matplotlib/matplotlib/issues/13533

Обход решения для matplotlib.boxplots, не распознающего маскированные значения.

0 голосов
/ 26 февраля 2019

Ваш код выдает ошибки, поэтому трудно сказать, используется ли этот код.Я не уверен, как matplotlib обрабатывает замаскированный массив, потому что я обычно создаю новый массив только с теми данными, которые мне нужны.Вы можете заменить:

 a = np.ma.masked_less( a, -9000 )
 b = np.ma.masked_less( b, -9000 )
 c = np.ma.masked_less( c, -9000 )

на

 a = a[a>-9000]
 b = b[b>-9000]
 c = c[c>-9000]
...