Re: [Freesurfer] MEG time series to surface overlay

I see. We don't really support w files anymore. Can you write it out as a .mgz/.mgh instead? Have sourcedata be nvertices x 1 x 1 x ntimepoints and use MRIwrite or save_mgh

On Wed, 15 Mar 2017, Dorothy Sincasto wrote:

Hi Bruce, first I do:

[vertex_coords, faces, magic] = read_surf('rh.13475.midthickness.decimated.05.surf')

I have a time series variable sourcedata =  6149X 149646. Is a double and the numbers are in the range of -1.181711152525857e-14

then i try to write it:

write_wfile(/'path/rh.13475.midthickness.decimated.05.MEG.mgz',sourcedata, vertex_coords),

and I get this error:

Error using bitshift

Double inputs must have integer values in the range of ASSUMEDTYPE.

Error in fwrite3 (line 26)

b1 = bitand(bitshift(val, -16), 255) ;

Error in write_wfile (line 58)

fwrite3(fid, v(i)) ;          % vertex number (0-based)

Thanks

Dorothy

On Tue, Mar 14, 2017 at 7:28 PM, Bruce Fischl fischl@nmr.mgh.harvard.edu wrote: Can you send us the full command you ran and the full screen output?

  Thanks
  Bruce

  > On Mar 14, 2017, at 9:33 PM, Dorothy Sincasto
  <dsincasto@gmail.com> wrote:
  >
  > Hi,
  >
  > I have MEG time series for every vertex in a decimated surface
  (~6000 vertices). I want to create a file (i think is called
  overlay) in which the timeseries are saved, like an fMRI mgz
  overlay. How can I do that with matlab?
  >
  > I have tried read_surf for reading the midthickness surface
  file, and then using the vertex coordinates and the MEG time
  series variable as input for write_wfile. however i get the
  error:
  >
  > Error using bitshift
  > Double inputs must have integer values in the range of
  ASSUMEDTYPE.
  >
  > Error in fwrite3 (line 26)
  > b1 = bitand(bitshift(val, -16), 255) ;
  >
  > Error in write_wfile (line 58)
  >   fwrite3(fid, v(i)) ;          % vertex number (0-based)
  >
  > It seems that's not the function I need to use.
  >
  > Thanks
  > Dorothy

...

---

Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

---

Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.