OtsuThresholdFilter implementation

src/MATLAB/OtsuThresholdFilter.m

+function imageOut = OtsuThresholdFilter(image1,device,useMatlab)
+%	This will set each pixel/voxel to the min value of the neighborhood defined by the given kernel.
+    
+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
+
+if (useMatlab)
+    threshold = graythresh(image1);
+    imageOut = image1 > threshold;
+    imageOut = im2uint8(imageOut);
+else
+    imageOut = CudaMex('OtsuThresholdFilter',image1,device);
+end
+end
\ No newline at end of file