【文章內容簡介】 反向變換得到的浮點坐標為(i+u,j+v)，其中i、j均為非負整數(shù)，u、v為[0,1)區(qū)間的浮點數(shù)，則這個像素得值 f(i+u,j+v) 可由原圖像中坐標為 (i,j)、(i+1,j)、(i,j+1)、(i+1,j+1)所對應的周圍四個像素的值決定，即：f(i+u,j+v) = (1u)(1v)f(i,j) + (1u)vf(i,j+1) + u(1v)f(i+1,j) + uvf(i+1,j+1)其中f(i,j)表示源圖像(i,j)處的的像素值，以此類推這就是雙線性內插值法。雙線性內插值法計算量大，但縮放后圖像質量高，不會出現(xiàn)像素值不連續(xù)的的情況。由于雙線性插值具有低通濾波器的性質，使高頻分量受損，所以可能會使圖像輪廓在一定程度上變得模糊。三次卷積法能夠克服以上兩種算法的不足，計算精度高，但計算量大，他考慮一個浮點坐標(i+u,j+v)周圍的16個鄰點，目的像素值f(i+u,j+v)可由如下插值公式得到：f(i+u,j+v) = [A] * [B] * [C][A]=[ S(u + 1) S(u + 0) S(u 1) S(u 2) ]┏ f(i1, j1) f(i1, j+0) f(i1, j+1) f(i1, j+2) ┓[B]=┃ f(i+0, j1) f(i+0, j+0) f(i+0, j+1) f(i+0, j+2) ┃┃ f(i+1, j1) f(i+1, j+0) f(i+1, j+1) f(i+1, j+2) ┃┗ f(i+2, j1) f(i+2, j+0) f(i+2, j+1) f(i+2, j+2) ┛┏ S(v + 1) ┓[C]=┃ S(v + 0) ┃┃ S(v 1) ┃┗ S(v 2) ┛┏ 12*Abs(x)^2+Abs(x)^3 , 0=Abs(x)1S(x)=｛ 48*Abs(x)+5*Abs(x)^2Abs(x)^3 , 1=Abs(x)2 ┗ 0 , Abs(x)=2S(x)是對 Sin(x*Pi)/x 的逼近（Pi是圓周率——π）最鄰近插值（近鄰取樣法）、雙線性內插值、三次卷積法 等插值算法對于旋轉變換、錯切變換、一般線性變換 和 非線性變換 都適用。5 程序主要源代碼分析 圖像數(shù)據(jù)存儲的主程序如下:FILE*file=fopen(,wb)。writebmpfileheader(file)。/*寫文件頭函數(shù)*/writebmpinfoheader(file)。/* 寫信息頭函數(shù)*/writeimagedata(file)。/* 寫圖像數(shù)據(jù)函數(shù)*/fclose(fi le)。寫文件頭和寫信息頭函數(shù)比較簡單，這里bfSize為230454。b iWidth,biHeight為320,240。bi Sizelmage為230400。biXPelsPerMeter,bi YPelsPerMeter，biClrUsed,biCIrImportant均為Oa以下給出 writeimagedata(file)II 數(shù)的原型:voidwriteimagedata(FILE* file){Int datapointer。Int mydata 。for (int i= 0。 i = 239。 i++)/*共240行*/datapointer data+(239i )*3 20 。/*data為采集映射時圖像數(shù)據(jù)首地址，計算每一要寫入的行在內存中的首地址*/for(int j= 0。j = 31 9。j + + )/* 共 320列*/{mydata=*( datapointer+j)。fwrite(amp。mydata,1, 1, fil e)。 /*寫入一行數(shù)據(jù)*/ } }} 圖像平滑模塊Typedef struct{unsigned short bfType。Unsigned bfSize。Unsigned short bfReserved 1。Unsigned short bfReserved 2。Unsigned bfOffBits。}BMPHEADER。BMPHEADER header。Typedef struct{unsigned biSize。unsigned biWidth。unsigned biHeight。unsigned short biplanes。unsigned short biBitCount。unsigned biCompression。unsigned biSizelmage。unsigned biXpelsPerMeter。unsigned biYpelsPerMeter。unsigned biClrUsed。unsigned biClrUsed。unsigned biClrImportant。 }BMPINFOHEADER。BMPINFOHEADER infobeader。在進行處理之前，先讀取文件信息，；FILE*fp。Fp=fopen(“”,“rb”)。fread(amp。(),2,1,fp)。fread(amp。(),4,1,fp)。fread(amp。(),2,1,fp)。fread(amp。(),4,1,fp)。fread(amp。((BMPINFOHEADER),1,fp)。接下來，將圖像數(shù)據(jù)存入我們的內存中，并將圖像數(shù)據(jù)灰度化：Unsigned char*prodata。Unsigned char*middledata。Int size=。Prodata=malloc(size)。Fread(prodata,1,size,fp)。Fclose(fp)。For(z=0,z。z++){middledata[z]=(unsigned char)(*prodata+z+2)*3 +(unsigned char )(*(prodata+z+1)* + (unsigned char)(*prodata +z)*)。Z++。Middledata[z]=middledata[z1]。Z++。Middledata[z]=middledata[z1]。} Memcpy(prodata,middledata,)。以下是高斯平滑的主程序：Double Template[9]。/*運算模版*/Double Coef=(double)()。/*模版均值*//*模版?zhèn)€稀疏聲明*/Template[0]=。Template[1]=。Template[2]=。 Template[3]=。Template[4]=。 Template[5]=。Template[6]=。 Template[7]=。Template[8]=。For(x=1。xXSaveDim1。x++){for(y=1。yYSaveDim1。y++){/*XsaveDim,YsaveDim為圖像存儲尺寸，x，y從1開始并且以XsaveDim1，YsaveDim1結束，是為了去除邊界*/for（j=1。j3。j++）{for(k=0。k3。k++){/*j，k為模版內部循環(huán)變量*/Result+=(*(prodata+GetPixelOffset(x=j1,y+k1)))*Template[3*j+k]。}}/*GetPixelOffset為我們自己定義的函數(shù)，返回對應x,y坐標的像素在數(shù)據(jù)空間中的存儲位置，返回值為(XSaveDim*y+x)*()*/Result*=Coed。If(Result255){Result=255。}Z=GetPixelOffset(x,y)。Middledata[z]=(unsigned char)(Result+)。Z++。Middledata[z]=middledata[z1]。Z++。Middledata[z]=middledata[z2]。}} 圖像平移主要代碼int Main(int argc, char **argv){ ARMTargetInit()。 OSInit()。