Loading large binary files in Matlab, quickly

39 次查看(过去 30 天)
I have some pretty massive data files (256 channels, on the order of 75-100 million samples) in int16 format. It is written in flat binary format, so the structure is something like: CH1S1,CH2S1,CH3S1 ... CH256S1,CH1S2,CH2S2,...
I need to read in each channel separately, filter and offset correct it, then save. My current bottleneck is loading each channel, which takes about 7-8 minutes... scale that up 256 times, and I'm looking at nearly 30 hours just to load the data! I am trying to intelligently use fread, to skip bytes as I read each channel; I have the following code in a loop over all 256 channels to do this:
offset = i - 1;
fseek(fid,offset*2,'bof');
dat = fread(fid,[1,nSampsTotal],'*int16',(nChan-1)*2);
Reading around, this is typically the fastest way to load parts of a large binary file, but is the file simply too large to do this any faster? Any suggestions would be much appreciated!
System details: MATLAB 2017a, Windows 7, 64bit
  4 个评论
dpb
dpb 2018-8-20
How much RAM do you actually have? Sounds like the performance hit is probably that you're running into actually being swapped in/out of virtual memory; fread is pretty quick for straight data transfer to/from memory.
Is the processing required dependent upon having the whole timeseries in memory or can you do it piecewise on each channel?
You may just have a system limitation here...
Josh
Josh 2018-8-20
编辑:Josh 2018-8-20
We have 32GB installed, but it is a shared computer, so the availability varies.
For the processing, I only need one channel at a time, but for the filtering and offset correction I'm doing it's necessary to have the entire timeseries per channel, to avoid filtering artifacts that might arise from splitting the timeseries.
I'm a bit confused about the RAM allocation, though. As I'm only trying to load in a subset of the of data (using the "skip" parameter in fread), it should definitely be doable from a RAM standpoint... (for the 37GB file I'm testing now, 1 channel out of the 256 should only be 149MB). Unless the 'skip' function of fread allocates memory in a way that I don't know of?

请先登录,再进行评论。

采纳的回答

OCDER
OCDER 2018-8-20
Seems like you have to use stream processing. Essentially load N frames of data for 256 channels, do the processing, save the frame, and repeat until done. Trying to do channel by channel by skipping 256 channel x 2 bytes seems slow. Here are some example for how to set that up.
The other option is to buy >64 GB RAM.
  7 个评论
OCDER
OCDER 2018-9-12
Hi Livio, to get an answer for your problem, please create a separate Question post instead of responding to this thread that is closed (answer is accepted).
Also, in your new Question post, format your code by selecting your code and pushing the {}Code button.
this is how to format code
for j = 1
end

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Data Import and Analysis 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by