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; 
     } 
    } 
}