请问高分辨率图像这里的坐标怎么修改？

fanzhenjie

我用POCS实现超分辨率图像重建，但是当低分辨率图像上点映射到高分辨率图像上时，这里的Wr始终为1，这不是想要的结果，我认为是upref(n1,n2)这里有问题，请哪位高手帮忙修改，谢谢！

1.    mcX = X + u;       %运动补偿到高分辨率图像的像素的坐标
   
2.        mcY=Y+ v;
   
3.     BW=edge(upref,'canny'); %边缘检测
   
4.     %%% Loop over entire (low-res) frame
   
5.     for m2 = 1:size(frame,2),
   
6.       for m1 = 1:size(frame,1),
   
7. 
   
8.         %%% Get high-resolution coordinates
   
9.         n1 = 2*m1;   %行坐标
   
10.      n2 = 2*m2;   %列坐标
    
11. 
    
12.     %%% Get coordinates of the motion compensated pixel
    
13.         N2 = mcX(n1,n2);
    
14.         N1 = mcY(n1,n2);
    
15. 
    
16.         %%% If not a border pixel
    
17.         if ( N1>3 & N1<size(upref,1)-2 & N2>3 & N2<size(upref,2)-2 )
    
18.       
    
19.         %%% Find center of the window where the PSF will be applied
    
20.         rN1 = round(N1);  %PSF窗口中心值
    
21.     rN2 = round(N2);
    
22.       
    
23.         %%% Calculate the effective window         5*5
    
24.         windowX = Y(rN1-2:rN1+2,rN2-2:rN2+2);
    
25.         windowY = X(rN1-2:rN1+2,rN2-2:rN2+2);
    
26.         
    
27.         weights = exp(-((N1-windowX).^2+(N2-windowY).^2)./2);   %计算高斯PSF权值，假设高斯方差为1，得到5x5的一个区域
    
28.        if  (BW(rN1,rN2)==1)   %如果PSF窗口中心在边缘上，就对weights进行修改
    
29.        [color=Red] Wr = exp(-(upref(n1,n2)-upref(rN1,rN2)).^2./2*(0.02).^2); [/color]        weights=weights.*Wr;
    
30.      
    
31.         end

复制代码

fanzhenjie

Wr = exp(-(upref(n1,n2)-upref(rN1,rN2)).^2./2*(0.02).^2);
我认为是这里的有问题，upref(rN1,rN2)表示PSF窗口中心对应的高分辨率图像那点的像素值，upref(n1,n2)应该表示为低分辨率图像映射到高分辨率图像对应点的像素值，但是这里坐标如果是(n1,n2)就导致upref(n1,n2)-upref(rN1,rN2)总是0，请问这里怎么修改？