peakfind

Finds peaks in 3d slices of 4d images. Adapted from peak_4dfp by Avi
Snyder and peak_nii by Tony Wilson.

Copyright (c) 2012 Washington University
Author: Kevin A. Archie <karchie@wustl.edu>

All Fortran subroutines in the package were written by Avi Snyder.

This package produces two MATLAB functions (as mex files):

------------------------------------------------------------

[peaks, roi] = peakfind(image, mmppixr, centerr, orad, mask, ...)

Produces a list of peaks and an ROI mask image.

image : single-precision 3D image data in which peaks will be identified
mmppixr : single-precision 3-vector of voxel dimensions in mm
centerr : single-precision 3-vector atlas space offset of image
orad : radius of peak ROIs, in mm
mask : mask applied before producing ROI

optional arguments can be specified like this:

[p, r] = peakfind(i, m, c, o, mask, 'vtpos', 0.5, 'vtneg', 0.1);

valid optional arguments are:

vtpos : positive threshold peak value (default 0)
vtneg : negative threshold peak value (default 0)
ntop : maximum number of output peaks (default 1000)
minvox : minimum output ROI voxel count (default 1)
dthresh : maximum distance for merging peaks (default 0)
ctpos : peak curvature threshold for positive peaks (default 0)
ctneg : peak curvature threshold for negative peaks (default 0)

return values are:

peaks: array of peak structures; each peak structure contains fields:
  i : floating point index (i,j,k) of peak center (rederived from x)
  x : floating point location (x,y,z) of peak center in atlas space
  v : value at peak center
  del2v : Laplacian at peak center
  weight : number of peaks consolidated into this peak
  voxels : size of this peak's ROI sphere, in voxels

------------------------------------------------------------

blurimg = sphereblur(image, mmppixr, radius)

Blurs the provided image by convolving with a spherical kernel of
specified radius.

image : single-precision 3D image to be blurred
mmppixr : single-precision 3-vector of voxel dimensions in mm
radius : sphere radius, in mm