自适应滑动窗口过滤器

Question

我想创建一个在稀疏矩阵（Image）上运行的自适应滑动窗口滤镜。我想计算每个具有非零值的物镜像素在9x9窗口内具有非零值的所有像素的平均值。我尝试使用以下代码

fun = @(x) mean(x(x~=0));

B = nlfilter(A,[9 9],fun);

A是要过滤的稀疏图像。代码运行但对未经过滤的稀疏图像中具有零像素值的元素返回非零值。我还想确保计算目标像素的平均值，窗口内至少存在5个非零样本。如果这是不可能的，我想继续扩展滑动窗口，直到满足为止。

非常感谢任何建议。