1/2
LFP data from Steinmetz et al. 2019
dataset
posted on 2019-08-25, 19:06 authored by Nick SteinmetzNick Steinmetz, Peter Zatka-Haas, Matteo CarandiniMatteo Carandini, Kenneth HarrisKenneth HarrisRaw LFP data during the recording sessions reported in Steinmetz et al. 2019.
Companion dataset to doi:10.6084/m9.figshare.9598406.
Each LFP file has 385 channels of int16 data at 2500Hz. They can be loaded directly into software like Neuroscope, and can be read in matlab like:
>> fid = fopen('thisfile.bin', 'r');
>> dat = fread(fid, [385 Inf], '*int16');
>> fclose(fid);
Note that the file may well be too large to fit into memory so you may not want to read "Inf" samples as in this example.
Each file has different start and end times, and there are drifts in the clocks between probes - so the timestamp files, which give the aligned timing information, must be used. These files are located in the companion dataset (with spiking and behavioral data). The aligned timestamps are specified in a particular way: an Nx2 matrix where the first column is sample numbers and the second column is the time in seconds at which those samples occurred. Since LFP was sampled regularly, N=2, just the first and last sample. So, to get the time at which every sample occurred, use this line in matlab:
>> tsData = readNPY('...lf.timestamps.npy');
>> allTS = interp1(tsData(:,1), tsData(:,2), tsData(1,1):tsData(2,1));
That's just linearly interpolating between the times given for the first and last samples.
