335 lines
10 KiB
C++
335 lines
10 KiB
C++
/*
|
|
* Copyright 2011, Blender Foundation.
|
|
*
|
|
* This program is free software; you can redistribute it and/or
|
|
* modify it under the terms of the GNU General Public License
|
|
* as published by the Free Software Foundation; either version 2
|
|
* of the License, or (at your option) any later version.
|
|
*
|
|
* This program is distributed in the hope that it will be useful,
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
* GNU General Public License for more details.
|
|
*
|
|
* You should have received a copy of the GNU General Public License
|
|
* along with this program; if not, write to the Free Software Foundation,
|
|
* Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
|
|
*
|
|
* Contributor:
|
|
* Jeroen Bakker
|
|
* Monique Dewanchand
|
|
*/
|
|
|
|
#include "COM_ScreenLensDistortionOperation.h"
|
|
#include "BLI_math.h"
|
|
#include "BLI_utildefines.h"
|
|
extern "C" {
|
|
#include "BLI_rand.h"
|
|
}
|
|
|
|
ScreenLensDistortionOperation::ScreenLensDistortionOperation() : NodeOperation()
|
|
{
|
|
this->addInputSocket(COM_DT_COLOR);
|
|
this->addInputSocket(COM_DT_VALUE);
|
|
this->addInputSocket(COM_DT_VALUE);
|
|
this->addOutputSocket(COM_DT_COLOR);
|
|
this->setComplex(true);
|
|
this->m_inputProgram = NULL;
|
|
this->m_valuesAvailable = false;
|
|
this->m_dispersion = 0.0f;
|
|
this->m_distortion = 0.0f;
|
|
}
|
|
void ScreenLensDistortionOperation::initExecution()
|
|
{
|
|
this->m_inputProgram = this->getInputSocketReader(0);
|
|
this->initMutex();
|
|
this->m_cx = 0.5f * (float)getWidth();
|
|
this->m_cy = 0.5f * (float)getHeight();
|
|
|
|
}
|
|
|
|
void *ScreenLensDistortionOperation::initializeTileData(rcti *rect)
|
|
{
|
|
void *buffer = this->m_inputProgram->initializeTileData(NULL);
|
|
updateDispersionAndDistortion();
|
|
return buffer;
|
|
}
|
|
|
|
void ScreenLensDistortionOperation::executePixel(float *outputColor, int x, int y, void *data)
|
|
{
|
|
const float height = this->getHeight();
|
|
const float width = this->getWidth();
|
|
MemoryBuffer *buffer = (MemoryBuffer *)data;
|
|
|
|
int dr = 0, dg = 0, db = 0;
|
|
float d, t, ln[6] = {0, 0, 0, 0, 0, 0};
|
|
float tc[4] = {0, 0, 0, 0};
|
|
const float v = this->m_sc * ((y + 0.5f) - this->m_cy) / this->m_cy;
|
|
const float u = this->m_sc * ((x + 0.5f) - this->m_cx) / this->m_cx;
|
|
const float uv_dot = u * u + v * v;
|
|
int sta = 0, mid = 0, end = 0;
|
|
|
|
if ((t = 1.f - this->m_kr4 * uv_dot) >= 0.f) {
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
ln[0] = (u * d + 0.5f) * width - 0.5f, ln[1] = (v * d + 0.5f) * height - 0.5f;
|
|
sta = 1;
|
|
}
|
|
if ((t = 1.f - this->m_kg4 * uv_dot) >= 0.f) {
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
ln[2] = (u * d + 0.5f) * width - 0.5f, ln[3] = (v * d + 0.5f) * height - 0.5f;
|
|
mid = 1;
|
|
}
|
|
if ((t = 1.f - this->m_kb4 * uv_dot) >= 0.f) {
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
ln[4] = (u * d + 0.5f) * width - 0.5f, ln[5] = (v * d + 0.5f) * height - 0.5f;
|
|
end = 1;
|
|
}
|
|
|
|
if (sta && mid && end) {
|
|
float jit = this->m_data->jit;
|
|
float z;
|
|
float color[4];
|
|
{
|
|
// RG
|
|
const int dx = ln[2] - ln[0], dy = ln[3] - ln[1];
|
|
const float dsf = sqrtf((float)dx * dx + dy * dy) + 1.f;
|
|
const int ds = (int)(jit ? ((dsf < 4.f) ? 2.f : sqrtf(dsf)) : dsf);
|
|
const float sd = 1.f / (float)ds;
|
|
|
|
for (z = 0; z < ds; ++z) {
|
|
const float tz = ((float)z + (jit ? BLI_frand() : 0.5f)) * sd;
|
|
t = 1.0f - (this->m_kr4 + tz * this->m_drg) * uv_dot;
|
|
d = 1.0f / (1.f + sqrtf(t));
|
|
const float nx = (u * d + 0.5f) * width - 0.5f;
|
|
const float ny = (v * d + 0.5f) * height - 0.5f;
|
|
buffer->readCubic(color, nx, ny);
|
|
tc[0] += (1.f - tz) * color[0], tc[1] += tz * color[1];
|
|
dr++, dg++;
|
|
}
|
|
}
|
|
{
|
|
// GB
|
|
const int dx = ln[4] - ln[2], dy = ln[5] - ln[3];
|
|
const float dsf = sqrtf((float)dx * dx + dy * dy) + 1.f;
|
|
const int ds = (int)(jit ? ((dsf < 4.f) ? 2.f : sqrtf(dsf)) : dsf);
|
|
const float sd = 1.f / (float)ds;
|
|
|
|
for (z = 0; z < ds; ++z) {
|
|
const float tz = ((float)z + (jit ? BLI_frand() : 0.5f)) * sd;
|
|
t = 1.f - (this->m_kg4 + tz * this->m_dgb) * uv_dot;
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
const float nx = (u * d + 0.5f) * width - 0.5f;
|
|
const float ny = (v * d + 0.5f) * height - 0.5f;
|
|
buffer->readCubic(color, nx, ny);
|
|
tc[1] += (1.f - tz) * color[1], tc[2] += tz * color[2];
|
|
dg++, db++;
|
|
}
|
|
|
|
}
|
|
if (dr) outputColor[0] = 2.0f * tc[0] / (float)dr;
|
|
if (dg) outputColor[1] = 2.0f * tc[1] / (float)dg;
|
|
if (db) outputColor[2] = 2.0f * tc[2] / (float)db;
|
|
|
|
/* set alpha */
|
|
outputColor[3] = 1.0f;
|
|
}
|
|
else {
|
|
outputColor[0] = 0.0f;
|
|
outputColor[1] = 0.0f;
|
|
outputColor[2] = 0.0f;
|
|
outputColor[3] = 0.0f;
|
|
}
|
|
}
|
|
|
|
void ScreenLensDistortionOperation::deinitExecution()
|
|
{
|
|
this->deinitMutex();
|
|
this->m_inputProgram = NULL;
|
|
}
|
|
|
|
void ScreenLensDistortionOperation::determineUV(float result[4], float x, float y, float distortion, float dispersion)
|
|
{
|
|
if (!this->m_valuesAvailable) {
|
|
updateVariables(distortion, dispersion);
|
|
}
|
|
determineUV(result, x, y);
|
|
}
|
|
|
|
void ScreenLensDistortionOperation::determineUV(float result[4], float x, float y) const
|
|
{
|
|
const float height = this->getHeight();
|
|
const float width = this->getWidth();
|
|
|
|
float d, t, ln[6] = {0, 0, 0, 0, 0, 0};
|
|
const float v = this->m_sc * ((y + 0.5f) - this->m_cy) / this->m_cy;
|
|
const float u = this->m_sc * ((x + 0.5f) - this->m_cx) / this->m_cx;
|
|
const float uv_dot = u * u + v * v;
|
|
|
|
if ((t = 1.f - this->m_kr4 * uv_dot) >= 0.f) {
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
ln[0] = (u * d + 0.5f) * width - 0.5f, ln[1] = (v * d + 0.5f) * height - 0.5f;
|
|
}
|
|
if ((t = 1.f - this->m_kg4 * uv_dot) >= 0.f) {
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
ln[2] = (u * d + 0.5f) * width - 0.5f, ln[3] = (v * d + 0.5f) * height - 0.5f;
|
|
}
|
|
if ((t = 1.f - this->m_kb4 * uv_dot) >= 0.f) {
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
ln[4] = (u * d + 0.5f) * width - 0.5f, ln[5] = (v * d + 0.5f) * height - 0.5f;
|
|
}
|
|
|
|
float jit = this->m_data->jit;
|
|
float z;
|
|
{
|
|
// RG
|
|
const int dx = ln[2] - ln[0], dy = ln[3] - ln[1];
|
|
const float dsf = sqrtf((float)dx * dx + dy * dy) + 1.f;
|
|
const int ds = (int)(jit ? ((dsf < 4.f) ? 2.f : sqrtf(dsf)) : dsf);
|
|
const float sd = 1.f / (float)ds;
|
|
|
|
z = ds;
|
|
const float tz = ((float)z + (1.0f)) * sd;
|
|
t = 1.0f - (this->m_kr4 + tz * this->m_drg) * uv_dot;
|
|
d = 1.0f / (1.f + sqrtf(t));
|
|
const float nx = (u * d + 0.5f) * width - 0.5f;
|
|
const float ny = (v * d + 0.5f) * height - 0.5f;
|
|
result[0] = nx;
|
|
result[1] = ny;
|
|
}
|
|
{
|
|
// GB
|
|
const int dx = ln[4] - ln[2], dy = ln[5] - ln[3];
|
|
const float dsf = sqrtf((float)dx * dx + dy * dy) + 1.f;
|
|
const int ds = (int)(jit ? ((dsf < 4.f) ? 2.f : sqrtf(dsf)) : dsf);
|
|
const float sd = 1.f / (float)ds;
|
|
|
|
z = ds;
|
|
const float tz = ((float)z + (1.0f)) * sd;
|
|
t = 1.f - (this->m_kg4 + tz * this->m_dgb) * uv_dot;
|
|
d = 1.f / (1.f + sqrtf(t));
|
|
const float nx = (u * d + 0.5f) * width - 0.5f;
|
|
const float ny = (v * d + 0.5f) * height - 0.5f;
|
|
result[2] = nx;
|
|
result[3] = ny;
|
|
}
|
|
}
|
|
|
|
bool ScreenLensDistortionOperation::determineDependingAreaOfInterest(rcti *input, ReadBufferOperation *readOperation, rcti *output)
|
|
{
|
|
rcti newInputValue;
|
|
newInputValue.xmin = 0;
|
|
newInputValue.ymin = 0;
|
|
newInputValue.xmax = 2;
|
|
newInputValue.ymax = 2;
|
|
|
|
NodeOperation *operation = getInputOperation(1);
|
|
if (operation->determineDependingAreaOfInterest(&newInputValue, readOperation, output) ) {
|
|
return true;
|
|
}
|
|
|
|
operation = getInputOperation(2);
|
|
if (operation->determineDependingAreaOfInterest(&newInputValue, readOperation, output) ) {
|
|
return true;
|
|
}
|
|
|
|
#define MARGIN 96
|
|
|
|
#define UPDATE_INPUT \
|
|
newInput.xmin = MIN3(newInput.xmin, coords[0], coords[2]); \
|
|
newInput.ymin = MIN3(newInput.ymin, coords[1], coords[3]); \
|
|
newInput.xmax = MAX3(newInput.xmax, coords[0], coords[2]); \
|
|
newInput.ymax = MAX3(newInput.ymax, coords[1], coords[3]);
|
|
|
|
rcti newInput;
|
|
float margin;
|
|
float coords[4];
|
|
if (m_valuesAvailable) {
|
|
determineUV(coords, input->xmin, input->ymin);
|
|
newInput.xmin = coords[0];
|
|
newInput.ymin = coords[1];
|
|
newInput.xmax = coords[0];
|
|
newInput.ymax = coords[1];
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmin, input->ymax);
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmax, input->ymax);
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmax, input->ymin);
|
|
UPDATE_INPUT;
|
|
margin = (ABS(this->m_distortion) + this->m_dispersion) * MARGIN + 2.0f;
|
|
}
|
|
else
|
|
{
|
|
determineUV(coords, input->xmin, input->ymin, 1.0f, 1.0f);
|
|
newInput.xmin = coords[0];
|
|
newInput.ymin = coords[1];
|
|
newInput.xmax = coords[0];
|
|
newInput.ymax = coords[1];
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmin, input->ymin, -1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
|
|
determineUV(coords, input->xmin, input->ymax, -1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmin, input->ymax, 1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
|
|
determineUV(coords, input->xmax, input->ymax, -1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmax, input->ymax, 1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
|
|
determineUV(coords, input->xmax, input->ymin, -1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
determineUV(coords, input->xmax, input->ymin, 1.0f, 1.0f);
|
|
UPDATE_INPUT;
|
|
margin = MARGIN;
|
|
}
|
|
|
|
#undef UPDATE_INPUT
|
|
newInput.xmin -= margin;
|
|
newInput.ymin -= margin;
|
|
newInput.xmax += margin;
|
|
newInput.ymax += margin;
|
|
|
|
operation = getInputOperation(0);
|
|
if (operation->determineDependingAreaOfInterest(&newInput, readOperation, output) ) {
|
|
return true;
|
|
}
|
|
return false;
|
|
}
|
|
|
|
void ScreenLensDistortionOperation::updateVariables(float distortion, float dispersion)
|
|
{
|
|
this->m_kg = maxf(minf(distortion, 1.0f), -0.999f);
|
|
// smaller dispersion range for somewhat more control
|
|
const float d = 0.25f * maxf(minf(dispersion, 1.0f), 0.0f);
|
|
this->m_kr = maxf(minf((this->m_kg + d), 1.0f), -0.999f);
|
|
this->m_kb = maxf(minf((this->m_kg - d), 1.0f), -0.999f);
|
|
this->m_maxk = MAX3(this->m_kr, this->m_kg, this->m_kb);
|
|
this->m_sc = (this->m_data->fit && (this->m_maxk > 0.f)) ? (1.f / (1.f + 2.f * this->m_maxk)) : (1.f / (1.f + this->m_maxk));
|
|
this->m_drg = 4.f * (this->m_kg - this->m_kr);
|
|
this->m_dgb = 4.f * (this->m_kb - this->m_kg);
|
|
|
|
this->m_kr4 = this->m_kr * 4.0f;
|
|
this->m_kg4 = this->m_kg * 4.0f;
|
|
this->m_kb4 = this->m_kb * 4.0f;
|
|
}
|
|
|
|
void ScreenLensDistortionOperation::updateDispersionAndDistortion()
|
|
{
|
|
if (this->m_valuesAvailable) return;
|
|
|
|
this->lockMutex();
|
|
if (!this->m_valuesAvailable) {
|
|
float result[4];
|
|
this->getInputSocketReader(1)->read(result, 0, 0, COM_PS_NEAREST);
|
|
this->m_distortion = result[0];
|
|
this->getInputSocketReader(2)->read(result, 0, 0, COM_PS_NEAREST);
|
|
this->m_dispersion = result[0];
|
|
updateVariables(this->m_distortion, this->m_dispersion);
|
|
this->m_valuesAvailable = true;
|
|
}
|
|
this->unlockMutex();
|
|
}
|