MouseLand / Kilosort

Fast spike sorting with drift correction for up to a thousand channels
https://kilosort.readthedocs.io/en/latest/
GNU General Public License v3.0
449 stars 235 forks source link

Channel Map File for Cambridge Neurotech H6 Molex ASSY 236 #539

Closed ryanlash2 closed 5 months ago

ryanlash2 commented 1 year ago

Hello All!

I have been having trouble getting my chanMap config file to translate correctly to kilosort. I'm using the bad channels to emphasize the issue that the gui isn't properly indexing the channels as seen in the image. I'm using a cambridge neurotech h6 molex probe 64channel 2shank(32 each). It correctly recognizes the coordinates just not the channel IDs. Any advice? `Nchannels = 64; connected = true(Nchannels, 1);

%channel map is based on adding one to the intan value for channel IDs since matlab doesn't index from zero chanMap = [12, 13, 8, 30, 32, 22, 15, 6, 27, 20, 5, 26, 4, 21, 9, 1, ... 29, 19, 14, 28, 24, 31, 22, 7, 18, 10, 2, 23, 3, 17, 25, 11, ... 63, 62, 39, 47, 36, 61, 33, 57, 49, 41, 51, 34, 37, 59, 46, 54, ... 40, 52, 56, 60, 43, 38, 55, 42, 45, 64, 50, 35, 53, 44, 58, 48]; chanMap0ind = chanMap - 1; bad_channels = [12,48]; connected(bad_channels) = false;

xcoords = [0ones(16,1); 22.5ones(16,1); 250ones(16,1); 272.5ones(16,1)]; xcoords = xcoords(:); ycoords =[(cumsum(25ones(16,1))-25); (cumsum(25ones(16,1))-12.5); (cumsum(25ones(16,1))-25); (cumsum(25ones(16,1))-12.5);] ycoords = ycoords(:); kcoords = ones(Nchannels,1); % grouping of channels (i.e. tetrode groups)

fs = 20000; % sampling frequency ` probeView

ryanlash2 commented 1 year ago

`Nchannels = 64; connected = true(Nchannels, 1);

%channel map is based on adding one to the intan value for channel IDs since matlab doesn't index from zero chanMap = [12, 13, 8, 30, 32, 16, 15, 6, 27, 20, 5, 26, 4, 21, 9, 1,... 29, 19, 14, 28, 24, 31, 22, 7, 18, 10, 2, 23, 3, 17, 25, 11,... 63, 62, 39, 47, 36, 61, 33, 57, 49, 41, 51, 34, 37, 59, 46, 54,... 40, 52, 56, 60, 43, 38, 55, 42, 45, 64, 50, 35, 53, 44, 58, 48] chanMap0ind = chanMap - 1; bad_channels = [12,48]; connected(bad_channels) = false;

xcoords =[0ones(16,1); 22.5ones(16,1); 250ones(16,1); 272.5ones(16,1)]; xcoords = zeros(Nchannels, 1); xcoords(chanMap) =xcoords; % xcoords =xcoords_d(chanMap) ycoords =[(cumsum(25ones(16,1))-25); (cumsum(25ones(16,1))-12.5); (cumsum(25ones(16,1))-25); (cumsum(25ones(16,1))-12.5)]; ycoords = zeros(Nchannels, 1); ycoords(chanMap) = ycoords; % ycoords =ycoords_d(chanMap) kcoords = ones(Nchannels,1);`