Median filter helper function

src/MATLAB/Medianfilter.m

+function imageOut = Medianfilter(image1,Neighborhood,device,useMatlab)
+%	NeighborhoodX, NeighborhoodY, and NeighborhoodZ are the directions and area to look for a given pixel.
+    
+if (~exist('device','var') || isempty(device))
+    device = 1;
+end
+
+if (~exist('useMatlab','var') || isempty(useMatlab))
+    useMatlab = false;
+end
+
+if (~useMatlab)
+    [numCudaDevices, memoryStats] = CudaMex('DeviceCount');
+    if (numCudaDevices < 1)
+        useMatlab = true;
+    elseif (device > numCudaDevices)
+        device = numCudaDevices;
+    end
+    
+    % do something smart here if there is not enough memory...
+end
+
+%Use same neighborhood for all 3 dimensions - Hack 
+NeighborhoodX = Neighborhood;
+NeighborhoodY = Neighborhood;
+NeighborhoodZ = Neighborhood;
+
+if (useMatlab)
+    imageOut = immedian(image1,NeighborhoodX,NeighborhoodY,NeighborhoodZ);
+else
+    imageOut = CudaMex('MedianFilter',image1,[NeighborhoodX,NeighborhoodY,NeighborhoodZ],device);
+end
+end
\ No newline at end of file