Search Comments and Ratings

go

   
Date File Comment by Comment Rating
27 Aug 2014 Robust smoothing for 1-D to N-D data (easy version of SMOOTHN) EZSMOOTHN allows unsupervised robust smoothing and can deal with missing values Author: Damien Garcia Fei

17 Jul 2014 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia ogadol

Hi Damien

Im applying smoothn in order to smooth bathmetric data sets and it works perfectly. I'm using a smooth factor of 10^9 on a 1601X1601 grid. The spacing of the points in realty is 12.5 m. can you please tell me what is the spatial meaning of the smooth factor in my case ?

Thanks

Ogadol

28 May 2014 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Rogier Westerhoff

Works fine for me!

09 May 2014 Label connected components in 2-D array LABEL is a generalization of BWLABEL Author: Damien Garcia Georg Stillfried

@asan masraf: You need to provide a matrix as input.
For example:
data = [
1 1 1 NaN NaN NaN NaN NaN
1 1 1 NaN 1 1 NaN NaN
1 1 1 NaN 1 1 NaN NaN
1 1 1 NaN NaN NaN 1 NaN
1 1 1 NaN NaN NaN 1 NaN
1 1 1 NaN NaN NaN 1 NaN
1 1 1 NaN NaN 1 1 NaN
1 1 1 NaN NaN NaN NaN NaN];
l4 = label(data,4);
l8 = label(data,8);

30 Apr 2014 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia John Booker

29 Apr 2014 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Matteo Niccoli

29 Apr 2014 Robust smoothing for 1-D to N-D data (easy version of SMOOTHN) EZSMOOTHN allows unsupervised robust smoothing and can deal with missing values Author: Damien Garcia Matteo Niccoli

19 Apr 2014 Label connected components in 2-D array LABEL is a generalization of BWLABEL Author: Damien Garcia asan masraf

I tried to run but I couldn't
it gave me this error:

Error using label (line 80)
Not enough input arguments.

24 Mar 2014 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Jospeh Chang

28 Feb 2014 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia arnold

hey Damien,

is there a way to suppress filling of missing values (nans to be specific)? it'd be really useful if this were an option.

Thanks
Arnold

22 Feb 2014 Robust smoothing for 1-D to N-D data (easy version of SMOOTHN) EZSMOOTHN allows unsupervised robust smoothing and can deal with missing values Author: Damien Garcia Warwick

04 Feb 2014 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Keith

Great work, been running well for me until I ran into an issue today. I received the error code:

??? Subscript indices must either be real positive integers or logicals.

Error in ==> inpaintn>InitialGuess at 179
z(~I) = y(L(~I));

Error in ==> inpaintn at 142
[y,s0] = InitialGuess(x,isfinite(x));

Error in ==> test at 32
betas_paint2(:,k) = inpaintn(betas_paint(:,k));

Was trying to run inpaintn through an external loop feeding it vectors. It made it through the first 8 then gave me an error. Any help is appreciated.

31 Jan 2014 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Ryan Muir

This is a great program! It is easy to use, and is really impressively fast.

I was using it for 3D and 4D interpolation of rather large problems, and was taking hours per image to process. I found that I could offload much of the problem to the GPUs, and get an order of magnitude more speed with a couple of lines of extra code! It's a bit sloppy, but putting y = gpuArray(y) near the beginning of the cdtn() and icdtn() functions bought me this performance enhancement. There are likely more efficient and elegant ways of sending it to the GPU, though this way worked great in a pinch. Small problems (including the example problem) run a bit slower with this modification, so this really is for large problems only.

19 Dec 2013 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Tete Tevi

14 Oct 2013 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia francisco fuentes

Really nice tool, thank you.

12 Sep 2013 Simpson's rule for numerical integration The Simpson's rule uses parabolic arcs instead of the straight lines used in the trapezoidal rule Author: Damien Garcia fernando

its very simple,and now i understand the trapz code but then the simps seems to not work on me ..it says
"Undefined function or method 'simps' for input arguments of
type 'double'." how can i fix this?

18 Aug 2013 Resize N-D arrays and images Y = RESIZE(X,newsize) resizes input array X using a discrete cosine transform. Author: Damien Garcia William Warriner

Thank you for this. Useful in my work for voxel array sampling for further use with Floyd-Warshall algorithm. Large voxel arrays take obscene quantities of time with Floyd-Warshall, so periodic sampling is useful.

Thank you also for packaging the two required files.

09 Aug 2013 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Jeff

I'm incorrect. Apologies.

Although I had downloaded Damien's new version of inpaintn, I had earlier used the older version (that I had neglected to delete), in a folder that put it ahead of the new version.

08 Aug 2013 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Jeff

The call to test4DCTNANDIDCTN can be eliminated because the functions dctn and idctn are now included in the inpaintn.m file.

Otherwise, works well. Results in my application compare reasonably with John D'Errico's inpaint_nans3, and is faster.

31 Jul 2013 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Marc Crapeau

Very good submission. I was working on my own version of this problem when a google research send me here and save me several hours of work. Thanks!

21 Jul 2013 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Ben

Hi Damien,
What if I have non-uniformly sampled array? Could it be possible to us this function? Currently, I found the smoothed curve did not maintain the start and end points in 3D. The smoothed curve looks being shrunk.

24 May 2013 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Martin Lechner

Super tool for smoothing!

19 May 2013 Simpson's rule for numerical integration The Simpson's rule uses parabolic arcs instead of the straight lines used in the trapezoidal rule Author: Damien Garcia Matthew

18 May 2013 Simpson's rule for numerical integration The Simpson's rule uses parabolic arcs instead of the straight lines used in the trapezoidal rule Author: Damien Garcia Angelos

That is a very good example for the for understand the Simpson's rule!! Excellent work!

10 May 2013 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Shrinivas

I am getting error
??? Error: File: inpaintn.m Line: 176 Column: 7
Expression or statement is incorrect--possibly unbalanced (, {, or [.

09 May 2013 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Scott Richardson

Hi, thank you for sharing this.

Do you happen to have onion.png, and wind.mat to try the other examples?

04 May 2013 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Ben

Thanks for sharing!! This is really a nice pice of code.

25 Apr 2013 Image segmentation using Otsu thresholding OTSU(I,N) segments the image I into N classes by means of Otsu's N-thresholding method. Author: Damien Garcia Venkat

08 Nov 2012 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Michael

What is with the "DCTN and iDCTN are required"? I was able to access them outside the MatLab browser, but not from within it.

30 Oct 2012 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia yo

hey everyone, the looks great, but unfotunately i've been getting the following error:

Error in ==> inpaintn at
140
y =
RF*idctn(Gamma.*dctn(W.*(x-y)+y))
+ (1-RF)*y;

since i don't feel capable enpugh changing the code on my own, can someone tell me what went wrong?

Thanks alot,

yoav

10 Oct 2012 Hybrid median filtering HMF performs hybrid median filtering of a 2-D array or an RGB image Author: Damien Garcia NCC A.

09 Jul 2012 Image segmentation using Otsu thresholding OTSU(I,N) segments the image I into N classes by means of Otsu's N-thresholding method. Author: Damien Garcia Jing

Thanks a lot for the code!
I just have a question: in the n=2 case, what's the meaning for the following code? I tried to derive the equation, but could not make it. Would you please explain a bit? Thanks!!

sigma2B =...
(mu(end)*w(2:end-1)-mu(2:end-1)).^2./w(2:end-1)./(1-w(2:end-1));

07 Jun 2012 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia none

09 May 2012 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Jakub Wosyka

Hello Damien,
very useful and nice function.
I tried to add the non-uniform grid (x values) according to your article but then it only interpolates the data. When I raise the value of smooting parameter, new the function will fall to zero (zero mean). (...probably my fault...). Do you plan to extend this nice function for non-uniform grids? Thank you

01 Mar 2012 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Kyle

Very useful script. I will never use matlab extrapolation routines ever again.

01 Mar 2012 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Kyle

23 Feb 2012 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Erdal Bizkevelci

21 Jan 2012 Image segmentation using Otsu thresholding OTSU(I,N) segments the image I into N classes by means of Otsu's N-thresholding method. Author: Damien Garcia Jamuna Jayashri Mahalingam

I am geting this error...How do I correct it?

Undefined function or method 'otsu' for input arguments of type 'uint8'.

09 Dec 2011 DCT and DST (+ inverse) in arbitrary dimension Discrete cosine/sine transforms and their inverses Author: Damien Garcia G. Sken

these function when used in >resize3d< allow faster volume resizing of large volumes with less memory usage

08 Nov 2011 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia tudor dima

Great work, extremely useful, it works out of the box.
Thanks a lot for posting it!

Would there be a simple way to restrict calculations inside a "gating mask"?

06 Oct 2011 Resize N-D arrays and images Y = RESIZE(X,newsize) resizes input array X using a discrete cosine transform. Author: Damien Garcia Lili

Nice job! Thanks!!

20 Sep 2011 Label connected components in 2-D array LABEL is a generalization of BWLABEL Author: Damien Garcia James Mure-Dubois

Nice work ! I used this function as an alternative to bwlabel. The only required modification in my code was to take into account the syntax difference : label puts 0 binary values into a connected object. As per the author's 'Important remark', zeros should be replaced by NaNs in this use case.

30 Aug 2011 Simpson's rule for numerical integration The Simpson's rule uses parabolic arcs instead of the straight lines used in the trapezoidal rule Author: Damien Garcia Danny Sale

This is a great extension of Simpson's rule. I find it most valuable that this file works correctly for both even and odd length vectors, and that it can correctly handle arbitrary spacing. Also, the method is fully vectorized so it is very fast. Thank you and excellent work!

23 May 2011 Label connected components in 2-D array LABEL is a generalization of BWLABEL Author: Damien Garcia tacoyu

18 May 2011 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Erik

Really nice method!

I does seem to me though that there is a bug when using this in combination with your pppiv-function.

In the function RobustWeights, there is a line:
MAD = median(abs(r(I)-median(r(I))));

However, the values of r are complex when using it together with pppiv. That causes the median function to sort the values by the real component and then take the middle value. I think you probably want to add
r = abs(r);
before taking the median to sort the residuals by magnitude instead. Is that correct?

11 Feb 2011 Image segmentation using Otsu thresholding OTSU(I,N) segments the image I into N classes by means of Otsu's N-thresholding method. Author: Damien Garcia Samuele

works smoothly...thx ;)

08 Dec 2010 Fast 'n easy smoothing SMOOTHN allows automatized and robust smoothing in arbitrary dimension w/wo missing values Author: Damien Garcia Richard

Well done! well documented. very easy to use.
Richard

28 Oct 2010 Label connected components in 2-D array LABEL is a generalization of BWLABEL Author: Damien Garcia Damien Garcia

to Irwin Donis: Because LABEL works with arrays of any class, the 0s are NOT considered as the background. Only NaN values are considered as background. Transform your background into NaNs before using LABEL.

13 Oct 2010 Image segmentation using Otsu thresholding OTSU(I,N) segments the image I into N classes by means of Otsu's N-thresholding method. Author: Damien Garcia Matteo Niccoli

It has worked for me very nicely. I like to run it in a while loop that breaks when sep starts decreasing.

12 Oct 2010 Inpaint over missing data in 1-D, 2-D, 3-D,... N-D arrays Y = INPAINTN(X) computes the missing data in the N-D array X. Author: Damien Garcia Bohuslav Ruzek

I had to interpolate missing 3D geophysical data in order to enable standard filtration process. Inpaint code worked very well for this, the results look pretty well. Using inpaint is easy.
I recommend using this script also for other interested users.

Contact us