Выбор Python Pandas hdfstore (где = '') возвращает неквалифицированные результаты

Когда я запрашиваю большой файл hdfstore (> 10G) следующим образом:

hdf = pd.HDFStore('raw_sample_storage.h5')
nrows = hdf.get_storer('raw_sample_all').nrows
chunksize = 300000

for i in xrange(nrows//chunksize + 1):
    chunk = hdf.select('raw_sample_all', where=[pd.Term('node_id', '==', 1)], start=i*chunksize, stop=(i+1)*chunksize)
    print chunk.head(2)

Я получил результаты, в которых node_id большинства записей равен 1, но у некоторых записей node_id отличается от 1. Значит, это сбой hdfstore или я сделал что-то не так?

Вот часть результатов, которые вы можете увидеть, есть некоторые записи с node_id, отличным от 1.

                 time   GW_time  node_id      X      Y      Z  status  seq  \
2 2013-10-22 17:20:58  39821888        1  16927  21438  22722       0   34   
6 2013-10-22 17:20:58  39822144        1  16927  21438  22722       0   35   

   rssi  lqi  
2   -46   48  
6   -51   48  
                      time   GW_time  node_id      X      Y      Z  status  \
300002 2013-10-22 17:30:50  59223744        3  19915  20840  22003       0   
300006 2013-10-22 17:30:50  59224000        3  19913  20844  22002       0   

        seq  rssi  lqi  
300002   46   -64   50  
300006   47   -64   48  
                      time   GW_time  node_id      X      Y      Z  status  \
600000 2013-10-22 17:40:55  79050561        1  17612  22536  21198       0   
600004 2013-10-22 17:40:55  79050817        1  17613  22535  21201       0   

        seq  rssi  lqi  
600000   55   -67   46  
600004   56   -67   49  
                      time   GW_time  node_id      X      Y      Z  status  \
900003 2013-10-22 17:50:44  98345217        4  18934  20212  19364       0   
900007 2013-10-22 17:50:44  98345473        4  18935  20212  19359       0   

        seq  rssi  lqi  
900003   32   -60   46  
900007   33   -60   48  
                       time    GW_time  node_id      X      Y      Z  status  \
1200003 2013-10-22 18:00:31  117600065        1  17618  22541  21191       0   
1200007 2013-10-22 18:00:31  117600321        1  17620  22538  21187       0   

         seq  rssi  lqi  
1200003  111   -66   47  
1200007  112   -66   48  

Заметив, что строка 300002 является нежелательным результатом, я пытаюсь выбрать узел 1 вокруг этой конкретной области следующим образом:

chunk = hdf.select('raw_sample_all', start=300002-20, stop=300002+20, 
                       where=[pd.Term('node_id', '==', 1)])                           

В результате возвращается только узел 3:

                   time     GW_time  node_id    X      Y       Z status seq rssi lqi
299982  2013-10-22 17:30:50 59222464    3   19912   20838   22003   0   41  -64 48
299986  2013-10-22 17:30:50 59222720    3   19912   20838   22003   0   42  -64 48
299990  2013-10-22 17:30:50 59222976    3   19913   20840   22007   0   43  -64 50
299994  2013-10-22 17:30:50 59223232    3   19913   20840   22007   0   44  -64 50
299998  2013-10-22 17:30:50 59223488    3   19915   20840   22003   0   45  -64 48
300002  2013-10-22 17:30:50 59223744    3   19915   20840   22003   0   46  -64 50
300006  2013-10-22 17:30:50 59224000    3   19913   20844   22002   0   47  -64 48
300010  2013-10-22 17:30:50 59224256    3   19913   20844   22002   0   48  -64 50
300014  2013-10-22 17:30:50 59224512    3   19914   20844   22010   0   49  -64 49
300018  2013-10-22 17:30:50 59224768    3   19914   20844   22010   0   50  -64 50                         

Затем я пытаюсь использовать index вместо start/stop следующим образом:

chunk = hdf.select('raw_sample_all', 
                   where=[pd.Term('index', '>=', 300002-20),
                          pd.Term('index', '<=', 300002+20),
                          pd.Term('node_id', '==', 1)])

И на этот раз он вернул правильные результаты:

                time        GW_time node_id  X       Y      Z   status seq rssi lqi
299984  2013-10-22 17:30:50 59222593    1   17613   22543   21203   0   42  -80 48
299988  2013-10-22 17:30:50 59222849    1   17613   22543   21203   0   43  -81 48
299992  2013-10-22 17:30:50 59223105    1   17610   22547   21194   0   44  -81 48
299996  2013-10-22 17:30:50 59223361    1   17610   22547   21194   0   45  -81 47
300000  2013-10-22 17:30:50 59223617    1   17609   22545   21190   0   46  -81 45
300004  2013-10-22 17:30:50 59223873    1   17609   22545   21190   0   47  -81 49
300008  2013-10-22 17:30:50 59224129    1   17606   22547   21199   0   48  -81 48
300012  2013-10-22 17:30:50 59224385    1   17606   22547   21199   0   49  -81 48
300016  2013-10-22 17:30:50 59224641    1   17607   22548   21191   0   50  -81 49
300020  2013-10-22 17:30:50 59224897    1   17607   22548   21191   0   51  -80 48  

Я думаю, что я мог бы обойти эту проблему с выбором по индексу, но я не совсем уверен, потому что метод с start/stop также дает правильные результаты большую часть времени, так что даже если метод с index правильно понял, где start/stop потерпел неудачу, он может потерпеть неудачу где-то еще.

И мне бы очень хотелось, чтобы работал метод start/stop, потому что он намного быстрее, а у меня большой набор данных, медленный метод действительно трудоемкий.

Кстати, если вам интересно, я не могу использовать 'chunksize' следующим образом:

df = hdf.select('raw_sample_all',chunksize=300000, where="node_id==1")
for chunk in df:
    print chunk.head(2)

Каждый раз, когда я пытаюсь использовать размер фрагмента, я получаю ошибку MemoryError, например -still-get-memory">это. Борясь со многими проблемами, Pandas действительно сложен для новичка, такого как я. Любая помощь приветствуется.


person Ewan    schedule 12.06.2015    source источник


Ответы (1)


Это была недавно исправленная ошибка в PyTables, см. связанную проблему здесь. В некоторых больших хранилищах индексаторы вычисляются неправильно при использовании where и start/stop.

Вам нужно будет обновиться до PyTables 3.2, затем переписать сам магазин. Вы можете воссоздать его так же, как в первый раз, или использовать ptrepack

person Jeff    schedule 12.06.2015
comment
Я использую Anaconda и не могу обновиться с помощью conda update pytables , он говорит, что он уже установлен. Затем я попробовал pip install --upgrade tables, там написано LINK : fatal error LNK1181: cannot open input file 'hdf5dll.lib' ... ERROR:: Could not find a local HDF5 installation. Я искал на своем жестком диске и не смог найти файл hdf5dll.lib. - person Ewan; 15.06.2015