In [2]:
import numpy as np
import os
rec_volume = np.memmap(os.path.join('.','test.tmp'),
                                dtype='float32', mode='w+', 
                                shape=(2029,2096, 4000)
                           )
rec_volume[:] = 0
del rec_volume

rec_volume = np.memmap(os.path.join('.','test.tmp'),
                                dtype='float32', mode='r+', 
                                shape=(2029,2096, 4000)
                           )
rec_volume[:] = 1
del rec_volume

In [5]:
!rm test.tmp

In [ ]: