Matlab 如何使用反褶积VWNR或其他反褶积方法对反褶积的输出信号施加约束？

我想用点扩散函数去卷积一个信号，以获取地面真值信号。我知道基本事实大致上是高斯的，应该是7x7像素

当我对原始信号进行反褶积时，这就起作用了。然而，当我用插值改变信号，然后尝试去卷积时，我得到了这个奇怪的周期性半高斯信号，它的长度与点扩散函数47像素相同

是否可以将输出的大小限制为7x7映像？或者约束它，使输出具有高斯分布？ 这是代码。变量位于dropbox链接中：

    deconvwnr(temp,sptint,0)

编辑：如果维纳滤波器反褶积不能做到这一点，是否有任何反褶积可以做到？或者，是否有一些数学方法来恢复7x7高斯信号使用我所拥有的

编辑2：我附上了一张逗号分隔的表格：

要生成高斯分布，我需要

w=116.4736
px=106.6667
truth  = fspecial('gaussian',7,w/px);

那么要生成sptint，我需要

load 'Alexa Fluor 647.csv';
spt = Alexa_Fluor_647(:,3);
sptint=interp1(300:900,spt,600:5:800);

然后生成temp，我会这样做

sptimg = conv2(truth,spt');
sptimg(:,end-(3-1):end)=[];
sptimg(:,1:3)=[];
sptimg2 = interp1(300:900,sptimg',600:5:800)';
temp=zeros(7,47)
temp(:,4:44)=sptimg2(:,:,1);

---

可以考虑加权平均数。 使用滑动窗口，将一个7乘7的窗口从temp的左侧移动到右侧，通过该窗口的权重sptint将数字叠加，并迭代添加结果

startAvg = 1:47-6;
avg = zeros(7);
for ii = 1:length(startAvg)
    avg = avg + temp(:,startAvg(ii):startAvg(ii)+6)*sptint(ii);
end
avg = avg/ii

---

过滤图像的输出不能小于原始图像。这就像是试图拍摄一辆汽车的图像，由于长时间曝光而沿路被弄脏，然后尝试获取原始的汽车尺寸。去毛刺是不行的，我明白了。这是真的。但知道原始图像是一辆汽车，知道汽车大小我知道高斯函数的sigma，如果你使用迭代反褶积方法，你不能约束它，使已知汽车大小之外的像素强度值为零吗？你能添加用于生成两个输出的代码吗？@bla我添加了代码和csv谢谢你的回答。我认为这是一个好办法。当我做conv2ptint时，平均答案与temp不太匹配。我想这很接近，我会玩一玩，接受你的回答，一旦我弄明白了，我就给你赏金。我搞不明白。模式接近正确，但值偏离了几个数量级。我想需要一些标准化因子。我也不太明白这其中的逻辑。这个方法有名字吗？