Hi - yeah, sorry - I originally wrote SUSAN as part of a real-time autonomous-guided-vehicle system, hence why the code is organised to be as fast as possible at the expense of readability.... ;-) Note though that we are planning on recoding SUSAN over the next year into C++, making it consistent with the rest of FSL and librafying it. Cheers, Steve. On Wed, 16 Nov 2005, Christopher Bailey wrote: > Hi Lucas, > > You're right, I can do it using system calls. However, it gets a bit > inefficient. Since I have the data in (ML) memory, it would be much > faster to operate on those instead of dumping to disk, calling susan, > and reading back in. > > For testing I've done as you suggested. I still think I'll try to find > the time to write the mex wrapper for it though. The code is quite an > interesting piece of work, all sorts of freaky pointer shifting taking > place... ;) > > -Chris > > On Fri, 2005-11-11 at 11:06 +0100, Lukas Scheef wrote: > > Hallo, Chris! > > > > > > Do you really need a mex wrapper? You could simply call SUSAN within > > a matlab script. I use this approach for 'bet' and it works fine. > > The code might look a bit ugly usig eval etc but it works. I > > > > > > command_string=['dos(','''YOUR ABS PATH comes in here bet "',input,'" > > "',output,'" -f',blanks(1),num2str(f), ' - > > g',blanks(1),num2str(g),options,'''',',','''-echo''',')']; > > > > [status, pipe]=eval(command_string); > > > > > > Best wishes, > > > > > > Lukas > > > > > > > > > > > > > > yoAm 28 Oct 2005 um 13:57 hat Stephen Smith geschrieben: > > > > > > > Hi - there are nifti/analyze reading routines for matlab in > > > $FSLDIR/etc/matlab > > > > > > Wrt SUSAN smoothing, you can write an image out to disk and then > > call the > > > command line program susan_smooth - note that this version only runs > > on > > > 16SI datatype, so you will need to convert to that type. > > > > > > Cheers, Steve. > > > > > > > > > > > > On Fri, 28 Oct 2005, Christopher Bailey wrote: > > > > > > > Dear FSLers, > > > > > > > > I have been toying with the idea of using SUSAN in some Matlab > > scripts > > > > of mine (intermediate smoothing stage in PET kinetic modelling). > > > > > > > > Has anyone implemented mex-wrappers for the SUSAN implementations? > > I > > > > would be very grateful for such code. > > > > > > > > Best regards, > > > > > > > > Chris > > > > > > > > > > > > > > -- > > > Stephen M. Smith DPhil > > > Associate Director, FMRIB and Analysis Research Coordinator > > > > > > Oxford University Centre for Functional MRI of the Brain > > > John Radcliffe Hospital, Headington, Oxford OX3 9DU, UK > > > +44 (0) 1865 222726 (fax 222717) > > > > > > [log in to unmask] http://www.fmrib.ox.ac.uk/~steve > > > > > > > -- Stephen M. Smith DPhil Associate Director, FMRIB and Analysis Research Coordinator Oxford University Centre for Functional MRI of the Brain John Radcliffe Hospital, Headington, Oxford OX3 9DU, UK +44 (0) 1865 222726 (fax 222717) [log in to unmask] http://www.fmrib.ox.ac.uk/~steve