Back to the main page.
Bug 1215 - problems plotting continuous resting data at sensor level using ft_miltiplotER
| Status | CLOSED FIXED |
| Reported | 2011-12-05 17:35:00 +0100 |
| Modified | 2012-03-21 15:43:50 +0100 |
| Product: | FieldTrip |
| Component: | plotting |
| Version: | unspecified |
| Hardware: | All |
| Operating System: | All |
| Importance: | P3 normal |
| Assigned to: | Jan-Mathijs Schoffelen |
| URL: | |
| Tags: | |
| Depends on: | |
| Blocks: | |
| See also: |
Francesco Di Pompeo - 2011-12-05 17:35:06 +0100
I used in the past ft_multiplotER in order to plotting raw data to have a quickly idea of broken channels. I'm analyzing 4d resting state data. With recente versions of fieldtrip it doesn't work and I visualize an empty head scheme and no errors. I have the same problem using both Windows and Mac. In the attached file I resampled the data only to minimize the size (generally I use no resampled raw data) This is the code I used to create the file cfg = []; cfg.dataset = '0'; cfg.channel = {'MEG'}; cfg.bpfilter = 'yes'; cfg.bpfreq = [1 150]; data_org = ft_preprocessing(cfg); cfg=[]; cfg.resamplefs = 200; cfg.detrend = 'no' ; cfg.demean = 'no' ; cfg.feedback = 'text'; cfg.trials = 'all' ; [data_org] = ft_resampledata(cfg, data_org); save preproc_data_dec data_org This piece of code can be used to reproduce the bug load preproc_data_dec cfg = []; sens=ft_read_sens('0'); cfg.grad=sens; cfg.channel = 'MEG'; clf; ft_multiplotER(cfg,data_org);
Boris Reuderink - 2011-12-07 12:02:32 +0100
Dear Francesco, Could you provide a dataset for your example? When I run the script I get the complaint that the file '0' is not found.
Francesco Di Pompeo - 2011-12-13 11:36:42 +0100
(In reply to comment #1) Hi Boris, sorry for the late reply...
Francesco Di Pompeo - 2011-12-16 16:19:49 +0100
you can download the raw data here.. http://dl.dropbox.com/u/45423174/0 and use the following script cfg = []; cfg.dataset = '0'; cfg.channel = {'MEG'}; cfg.bpfilter = 'yes'; %'no' or 'yes' bandpass filter (default = 'no') cfg.bpfreq = [1 150]; % bandpass frequency range, specified as [low high] in Hz data_org = ft_preprocessing(cfg); cfg = []; cfg.xlim = [-1. 360]; cfg.ylim = [-5e-12 5e-12]; cfg.channel = 'MEG'; cfg.showlabels = 'yes'; clf; cfg.xparam = 'time'; ft_multiplotER(cfg,data_org);
Boris Reuderink - 2011-12-22 15:33:43 +0100
Would you also have a config file for me?
Francesco Di Pompeo - 2011-12-22 15:40:43 +0100
Ok, it is there now...
Boris Reuderink - 2012-01-03 11:53:34 +0100
Sorry, my mistake. I get the following error with your script: ??? Error using ==> read_4d_hdr at 236 Cannot open config file Error in ==> ft_read_header at 171 orig = read_4d_hdr(datafile); Error in ==> ft_preprocessing at 338 hdr = ft_read_header(cfg.headerfile, 'headerformat', cfg.headerformat); Error in ==> bug1215 at 8 data_org = ft_preprocessing(cfg); Could you perhaps provide a complete example in a zip file?
Francesco Di Pompeo - 2012-01-04 10:51:22 +0100
Sorry, the correct link for the config file is http://dl.dropbox.com/u/45423174/config But note in previous version the visualization was working..
Boris Reuderink - 2012-01-10 11:23:49 +0100
Thank you. The script is running now.
Boris Reuderink - 2012-01-10 11:56:21 +0100
Created attachment 210 Pathological plot.
Boris Reuderink - 2012-01-10 11:57:58 +0100
Confirmed, an empty plot is displayed indeed. Verified with SVN revision 5104. Note: I have stored the files to reproduce this problem.
Jan-Mathijs Schoffelen - 2012-01-27 08:43:58 +0100
assigned bug to get all bugs lined up before the grand bug binge
Jan-Mathijs Schoffelen - 2012-02-01 15:09:05 +0100
Hi Francesco, Please specify cfg.layout = '4D248.lay'; Then it works