|
@@ -51,19 +51,6 @@ lamino_bp_generic (global float *proj,
|
|
|
const float oldx = newx * param->mat_0 + newy * param->mat_1 + newz * param->mat_2 + param->proj_ox;
|
|
|
const float oldy = newx * param->mat_3 + newy * param->mat_4 + newz * param->mat_5 + param->proj_oy;
|
|
|
|
|
|
- /* bilinear interpolation */
|
|
|
- const float yf_1 = oldy - floor(oldy);
|
|
|
- const float yf_0 = 1.0f - yf_1;
|
|
|
- const float xf_1 = oldx - floor(oldx);
|
|
|
- const float xf_0 = 1.0f - xf_1;
|
|
|
- const int base = ((int) floor(oldx)) + ((int) floor(oldy)) * param->proj_sx;
|
|
|
- float result;
|
|
|
-
|
|
|
- /* TODO: check that out of proj plain */
|
|
|
- result = proj[base ] * xf_0 * yf_0;
|
|
|
- result += proj[base + 1] * xf_1 * yf_0;
|
|
|
- result += proj[base + param->proj_sx ] * xf_0 * yf_1;
|
|
|
- result += proj[base + param->proj_sx + 1] * xf_1 * yf_1;
|
|
|
volume[idx] += interpolate (proj, param, oldx, oldy);
|
|
|
}
|
|
|
|