Hi Paul,

I'm so sorry that there is a bug in the code. It is fixed now.

For question 2, I would look at the images to choose the window size, and may try different window size using the MATLAB stdfilt function to find a good window size. For more accurate prior estimation, I will try SNR estimation methods to estimate the prior. I think you can find a good SNR estimation method that fits your data. For simplest, you may use the default window size of the stdfilt function ([3 3 3]) or your smoothing window size.

Thanks for your interest!

Xue