old-nall/image/interpolation.hpp

63 lines
2.3 KiB
C++
Raw Permalink Normal View History

2016-05-15 10:42:10 +00:00
#pragma once
namespace nall {
auto image::isplit(uint64_t* c, uint64_t color) -> void {
c[0] = (color & _alpha.mask()) >> _alpha.shift();
c[1] = (color & _red.mask() ) >> _red.shift();
c[2] = (color & _green.mask()) >> _green.shift();
c[3] = (color & _blue.mask() ) >> _blue.shift();
}
auto image::imerge(const uint64_t* c) -> uint64_t {
return c[0] << _alpha.shift() | c[1] << _red.shift() | c[2] << _green.shift() | c[3] << _blue.shift();
}
auto image::interpolate1f(uint64_t a, uint64_t b, double x) -> uint64_t {
return a * (1.0 - x) + b * x;
}
auto image::interpolate1f(uint64_t a, uint64_t b, uint64_t c, uint64_t d, double x, double y) -> uint64_t {
return a * (1.0 - x) * (1.0 - y) + b * x * (1.0 - y) + c * (1.0 - x) * y + d * x * y;
}
auto image::interpolate1i(int64_t a, int64_t b, uint32_t x) -> uint64_t {
return a + (((b - a) * x) >> 32); //a + (b - a) * x
}
auto image::interpolate1i(int64_t a, int64_t b, int64_t c, int64_t d, uint32_t x, uint32_t y) -> uint64_t {
a = a + (((b - a) * x) >> 32); //a + (b - a) * x
c = c + (((d - c) * x) >> 32); //c + (d - c) * x
return a + (((c - a) * y) >> 32); //a + (c - a) * y
}
auto image::interpolate4f(uint64_t a, uint64_t b, double x) -> uint64_t {
uint64_t o[4], pa[4], pb[4];
isplit(pa, a), isplit(pb, b);
for(unsigned n = 0; n < 4; n++) o[n] = interpolate1f(pa[n], pb[n], x);
return imerge(o);
}
auto image::interpolate4f(uint64_t a, uint64_t b, uint64_t c, uint64_t d, double x, double y) -> uint64_t {
uint64_t o[4], pa[4], pb[4], pc[4], pd[4];
isplit(pa, a), isplit(pb, b), isplit(pc, c), isplit(pd, d);
for(unsigned n = 0; n < 4; n++) o[n] = interpolate1f(pa[n], pb[n], pc[n], pd[n], x, y);
return imerge(o);
}
auto image::interpolate4i(uint64_t a, uint64_t b, uint32_t x) -> uint64_t {
uint64_t o[4], pa[4], pb[4];
isplit(pa, a), isplit(pb, b);
for(unsigned n = 0; n < 4; n++) o[n] = interpolate1i(pa[n], pb[n], x);
return imerge(o);
}
auto image::interpolate4i(uint64_t a, uint64_t b, uint64_t c, uint64_t d, uint32_t x, uint32_t y) -> uint64_t {
uint64_t o[4], pa[4], pb[4], pc[4], pd[4];
isplit(pa, a), isplit(pb, b), isplit(pc, c), isplit(pd, d);
for(unsigned n = 0; n < 4; n++) o[n] = interpolate1i(pa[n], pb[n], pc[n], pd[n], x, y);
return imerge(o);
}
}