spatial smoothing functions added

both Gaussian and Median filtering options added to eliminate background noise from erratically behaving voxels
I think this should replace each plane with the result of the filtering, and then return the volume, right? As in

```
for z in arange(0, dims[2]):
    im[:,:,z] = gaussian_filter(im[:,:,z], sigma)
return im
```