http://csharpimage.blog60.fc2.com/b

質問

http://csharpimage.blog60.fc2.com/blog-entry-19.html
をみて、直線補間をC言語で記述したのですが、結果が思う様にいきません。(_____)
プログラムは以下です。間違っているのは、おそらく、座標を示すところだとはおもっているのですが、
そこをどうすればいいのか。。。。レベルが低くてごめんなさい（＞＿＜
どなたかご教授願います。
FILE *fpt;
_wfopen_s(&fpt,Common_Data_Raw->filename,L"rb");
FILE *fpt_output;
int width=Common_Data_Raw->width;
int height=Common_Data_Raw->height;
// 拡大縮小後の画像サイズ
int hxSize=Common_Data_Raw->width_rescale;//拡大、縮小後の幅が入っている
int hySize=Common_Data_Raw->height_rescale;//拡大、縮小後の高さが入っている。
// 拡大縮小用
int xSize=width;
int ySize=height;
double xpos, ypos;
double hokanX = (double)(xSize - 1) / hxSize;
double hokanY = (double)(ySize - 1) / hySize;
int i,j;
unsigned char **layer,**bufdata,**rescale;
//メモリの確保→省略します。

if ((xSize < hxSize)&&(ySize<hySize))
{
// X方向の補間
for (i=0;i<ySize;i++)
{
xpos = 0.0;
for (j = 0; j < hxSize; j+=3)
{
bufdata[i][j] = (unsigned char)(((double)layer[i][(int)xpos + 1] - (double)layer[i][(int)xpos]) *(xpos - (int)xpos) +(double)layer[i][(int)xpos]);
bufdata[i][j+1] = (unsigned char)(((double)layer[i][(int)xpos + 1+1] - (double)layer[i][(int)xpos+1]) *(xpos - (int)xpos) +(double)layer[i][(int)xpos+1]);
bufdata[i][j+2] = (unsigned char)(((double)layer[i][(int)xpos + 1+2] - (double)layer[i][(int)xpos+2]) *(xpos - (int)xpos) +(double)layer[i][(int)xpos+2]);
xpos += hokanX;}
}

// Y方向の補間
for (i = 0; i < hxSize; i+=3)
{ypos = 0.0;for (j = 0; j <hySize; j++)
{rescale[j][i] = (unsigned char)(((double)bufdata[(int)ypos + 1][i] - (double)bufdata[(int)ypos][i]) *(ypos - (int)ypos) +
(double)bufdata[(int)ypos][i]);rescale[j][i+1] = (unsigned char)(((double)bufdata[(int)ypos + 1][i+1] - (double)bufdata[(int)ypos][i+1]) *(ypos - (int)ypos) + (double)bufdata[(int)ypos][i+1]);
rescale[j][i+2] = (unsigned char)(((double)bufdata[(int)ypos + 1][i+2] - (double)bufdata[(int)ypos][i+2]) * (ypos - (int)ypos) + (double)bufdata[(int)ypos][i+2]);
ypos += hokanY;}}
}
else
{
//rescaledata = SimpleRescaleImage(data, hxSize, hySize);

回答No.1

salsberry · ベストアンサー

xposの計算がデータ構造に合っていないのだと思われます。
一つの座標に対して、配列layerは3要素 (RGBでしょうか) を使っていますよね。たとえばlayer[i][12], layer[i][13], layer[i][14]で一組です。
xposとhokanXの計算はこのことを考慮していないように見えるので、xposの値が途中14.2になったりするとlayer[i][14], layer[i][15], layer[i][16]を組にしてアクセスします。これでは思うような結果が得られなくて当然です。
さらに、隣の座標をアクセスするのにlayer[i][(int)xpos+1+α]としている部分はlayer[i][(int)xpos+3+α]でなくてはいけないはずです。

データ構造を変えていいのであれば、layer[i][j], layer[i][j+1], layer[i][j+2]のような形ではなく、構造体の二重配列にしてlayer[i][j].r, layer[i][j].g, layer[i][j].bのようにアクセスすれば簡単になります。bufdataとrescaleも同様に。

お礼

ありがとうございました☆