In this paper, NL-PCA is used as a first denoising stage to estimate spatially variant noise and to guide patch selection on the second stage (PRI-NLM3D). So basically they work in tandem, with the resulting two-stage algorithm called PRI-NL-PCA. 

Thanks Mauro. Do you have code for this paper implementation?

You`re welcome, I found a demo by the authors. It works!

http://personales.upv.es/jmanjon/denoising/nlpca.html

I also have this code. Can you tell how .p file runs?

I am facing some problems in running it.

Which kind of problems are you facing with the code?

This is my file listing:

•  demo.m
•  cPRI_NL_PCA.mexw64
• ssim_index3d.m 
• debias.p 
• t1_icbm_normal_1mm_pn0_rf0.rawb 
• NL_PCA.p demo.m 
• NoiseMapping.p 
• cNL_PCA.mexw64 medfilt3.m

I put all of them in a separate directory. Open MATLAB, change the working directory and run "demo.m". This works fine. 

If you want to run the routines on your own images, you should convert them to double precision first. Just type "my_image = double(my_image)" in the command line. Otherwise the mexw64 compiled code will crash. I found execution times to be rather long on my I5 workstation, so you should first try denoising small volumes. 

Best,

Mauro

Which kind of problems are you facing with the code?

This is my file listing:

•  demo.m
•  cPRI_NL_PCA.mexw64
• ssim_index3d.m 
• debias.p 
• t1_icbm_normal_1mm_pn0_rf0.rawb 
• NL_PCA.p demo.m 
• NoiseMapping.p 
• cNL_PCA.mexw64 medfilt3.m

I put all of them in a separate directory. Open MATLAB, change the working directory and run "demo.m". This works fine. 

If you want to run the routines on your own images, you should convert them to double precision first. Just type "my_image = double(my_image)" in the command line. Otherwise the mexw64 compiled code will crash. I found execution times to be rather long on my I5 workstation, so you should first try denoising small volumes. 

Best,

Mauro

Dear Mauro,

I am getting this error

Dear Laraib,

Could you type the "dir *.mexw64" command and paste the results? Are you using a 64 bit Windows version?

Dear Mauro,

I am using 32- bit operating system. Output is attached below.

Dear Laraib,

That's the problem, you need a 64bit windows to run mex64 compiled files.

Thanks Mauro. The author have protected the source code in NL_PCA.p ,  cPRI_NL_PCA.mexw64 and i want to understand the steps of source code. Is there any way to understand the code?