Thanks for sharing this tool! works very well and fast, but i'm experiencing some problems:
I'm trying to filter a 3D MRI volume of size 227x227x147. Depending on the value of 'rho' I use I get a variable number of slices in the third dimension that are all 'NaN'. The bigger rho, the more NaN slices are in the output image volume. For every value of T used the number of NaN-slices accumulates, so that after several timesteps the whole image volume consists of NaNs.
What could be the reason for the NaNs in the output?
My original image is single and normalized to a range of 0 to 1 and does not contain any NaNs.
thx, Elmar
Hi,
there's a little typo in compile_c_files.m.
Last file to be compiled has not the .c extension.
This: "diffusion_scheme_3D_novel_getUpdate.c"
Thanks.
Hi
I am getting the allowing error…. any help is greatly appreciated. I am using matlab 2012b on Mac os. thanks
Sourin
derivatives.c: In function 'gradient2Dx':
derivatives.c:24: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient2Dy':
derivatives.c:57: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient3Dx_double':
derivatives.c:110: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient3Dx_float':
derivatives.c:195: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient3Dy_float':
derivatives.c:280: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient3Dz_float':
derivatives.c:364: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient3Dy_double':
derivatives.c:449: warning: incompatible implicit declaration of built-in function 'malloc'
derivatives.c: In function 'gradient3Dz_double':
derivatives.c:533: warning: incompatible implicit declaration of built-in function 'malloc'
mex: compile of ' "derivatives.c"' failed.
Error using mex (line 206)
Unable to complete successfully.
Error in compile_c_files (line 8)
mex(filename,'-v');
When filtering 3D images I noticed that there is not option to specify the voxel resolution. To me this implies that in creating this function the author assumed that the voxels of the input image have isotropic resolution. Unfortunately this is usually not the case and must be accounted for when discretizing the diffusion equation.
Aside from this, I think this is a great submission. Many thanks!.
The original stencel can be found here:
http://www.scribd.com/doc/52387266/Jahne-B-Handbook-of-Computer-Vision-and-Applications-Vol-2-Signal-Processing-and-Pattern-Recognition
Comment only