C ++,Bilinear Interpolation在扩大的位图图像上修补孔

时间:2016-06-07 10:58:49

标签: c++ bmp bilinear-interpolation

为了将位图图像缩放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);
}

1 个答案:

答案 0 :(得分:0)

感谢大家,我终于完成了代码。是的,如果已经受到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;
        }
    }
}