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
|