C++，Bilinear Interpolation来修补位图上的孔Image

Question

为了将位图图像缩放3倍，我实现了这个代码并找到了漏洞。我决定使用双线性插值，上面这些代码修补了像这个代码一样计算重量的最近像素。输出[j * hInfo.biWidth + i] =图像[3/2 * j * hInfo.biWidth + 3/1 * i]。
逻辑很简单，所以我认为它会正常工作，我非常确定，但Image的结果看起来像我感觉的西瓜，请看看代码。任何意见都欢迎。谢谢阅读。

void main() 
{ 
BITMAPFILEHEADER hf; 
BITMAPINFOHEADER hInfo; 
RGBQUAD hRGB[256]; 
FILE *fp; 
fp = fopen("input.bmp", "rb"); 
if(fp == NULL) return; 
fread(&hf, sizeof(BITMAPFILEHEADER), 1, fp); 
fread(&hInfo, sizeof(BITMAPINFOHEADER), 1, fp); 
fread(hRGB, sizeof(RGBQUAD), 256, fp); 
hInfo.biWidth *=3; 
hInfo.biHeight *=3; 
int ImgSize = hInfo.biWidth * hInfo.biHeight; 

BYTE *Image = (BYTE *) malloc(ImgSize); 
BYTE *Output = (BYTE *)malloc(ImgSize); 
fread(Image, sizeof(BYTE), ImgSize, fp); 
fclose(fp); 
double Sx, Sy; 
scanf("%lf", & Sx);// scaling factor x : 3 
scanf("%lf", & Sy);// scaling factor y : 3 

for(int i= 0 ; i < hInfo.biHeight; i++){ 
    for(int j=0; j < hInfo.biWidth; j++){ 
     if(j*Sx < hInfo.biWidth && i*Sy<hInfo.biHeight) 

Output[(int)(i*Sy)*hInfo.biWidth+(int)(j*Sx)] = Image[i*hInfo.biWidth + j]; 
    } 
} 


// horizontally executed first. 

for(int i= 0 ; i < hInfo.biHeight; i++){ 
    for(int j=0; j < hInfo.biWidth; j++){ 
     if(j%3 == 1) 
     { 
      Output[i*hInfo.biWidth+j] = (int)Image[((hInfo.biWidth*(1/3) + j+2) + (hInfo.biWidth* (2/3) + (j-1)))]; 
     } 
     else if(j%3 == 2) 
     { 
      Output[i*hInfo.biWidth+j] = (int)Image[((hInfo.biWidth*(2/3) + j+1) + (hInfo.biWidth* (1/3) + (j-2)))]; 
     } 
     else if( j%3 ==1 && j+2 >= hInfo.biWidth) 
     { 
      Output[i*hInfo.biWidth+j] = Image[i*hInfo.biWidth + (j-2)] ; 
     } 
     else if( j%3 ==2 && j+1 >= hInfo.biWidth) 
     { 
      Output[i*hInfo.biWidth+j] = Image[i*hInfo.biWidth + (j-2)] ; 
     } 

    } 
} 


// vertically executed last 
for(int j= 0 ; j < hInfo.biWidth; j++){ 

    for(int i=0; i < hInfo.biHeight; i++){ 

     if(j%3 == 1) 
     { 
      Output[j*hInfo.biWidth+i] = Image[hInfo.biWidth*(1/3) + (i+2) + hInfo.biWidth* (2/3) + (i-1); 
     } 
     else if(j%3 == 2) 
     { 
      Output[j*hInfo.biWidth+i] = Image[hInfo.biWidth*(2/3) + (i+1) + hInfo.biWidth* (1/3) + (i-2)]; 
     } 
     else if( j%3 ==1 && j+2>=hInfo.biWidth) 
     { 
      Output[j*hInfo.biWidth+i] = Image[j*hInfo.biWidth + (i-2)] ; 
     } 
     else if( j%3 ==2 && j+1>=hInfo.biWidth) 
     { 
      Output[j*hInfo.biWidth+i] = Image[j*hInfo.biWidth + (i-2)] ; 
     } 

    } 
} 

fp = fopen("output.bmp", "wb"); 
fwrite(&hf, sizeof(BYTE), sizeof(BITMAPFILEHEADER), fp); 
fwrite(&hInfo, sizeof(BYTE), sizeof(BITMAPINFOHEADER), fp); 
fwrite(hRGB, sizeof(RGBQUAD), 256, fp); 
fwrite(Output, sizeof(BYTE),ImgSize, fp); 
fclose(fp); 
free(Image); 
free(Output); 
} 

Answer 1

感谢你们所有人，我终于完成了代码。是的，如果已经被for循环条件限制的话，我不需要其他的。

for(int i= 0 ; i < hInfo.biHeight; i++){ 
    for(int j=0; j < hInfo.biWidth; j++){ 
     if(j%3 == 1) 
     { 
      Output[i*hInfo.biWidth+j] = (int)Output[i*hInfo.biWidth+j-1]*2/3+(int)Output[i*hInfo.biWidth+j+2]*1/3; 
     } 
     else if(j%3 == 2) 
     { 
      Output[i*hInfo.biWidth+j] = (int)Output[i*hInfo.biWidth+j-1]*1/3 + (int)Output[hInfo.biWidth*i + (j-2)]*2/3; 
     } 

    } 
} 
for(int i= 0 ; i < hInfo.biHeight; i++){  
    for(int j= 0; j < hInfo.biWidth; j++){  
     if(j%3 == 1) 
     { 
      Output[j*hInfo.biWidth+i] = (int)Output[(j-1)*hInfo.biWidth+i]*2/3+(int)Output[(j+2)*hInfo.biWidth+i]*1/3; 
     } 
     else if(j%3 == 2) 
     { 
      Output[j*hInfo.biWidth+i] = (int)Output[(j-1)*hInfo.biWidth+i]*1/3 + (int)Output[hInfo.biWidth*(j-2) + i]*2/3; 
     } 
    } 
}