EWA sampling is designed for downsampling images, i.e. scaling down the size of input image pixels, which happens regularly in compositing. While the standard sampling methods (linear, cubic) work reasonably well for linear transformations, they don't yield good results in non-linear cases like perspective projection or arbitrary displacement. EWA sampling is comparable to mipmapping, but avoids problems with discontinuities. To work correctly the EWA algorithm needs partial derivatives of the mapping functions which convert output pixel coordinates back into the input image space (2x2 Jacobian matrix). With these derivatives the EWA algorithm projects ellipses into the input space and accumulates colors over their area. This calculation was not done correctly in the compositor, only the derivatives du/dx and dv/dy were calculation, basically this means it only worked for non-rotated input images. The patch introduces full derivative calculations du/dx, du/dy, dv/dx, dv/dy for the 3 nodes which use EWA sampling currently: PlaneTrackWarp, MapUV and Displace. In addition the calculation of ellipsis area and axis-aligned bounding boxes has been fixed. For the MapUV and Displace nodes the derivatives have to be estimated by evaluating the UV/displacement inputs with 1-pixel offsets, which can still have problems on discontinuities and sub-pixel variations. These potential problems can only be alleviated by more radical design changes in the compositor functions, which are out of scope for now. Basically the values passed to the UV/Displacement inputs would need to be associated with their 1st order derivatives, which requires a general approach to derivatives in all nodes.
57 lines
2.3 KiB
C++
57 lines
2.3 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
|
|
*/
|
|
|
|
#ifndef _COM_ReadBufferOperation_h
|
|
#define _COM_ReadBufferOperation_h
|
|
|
|
#include "COM_NodeOperation.h"
|
|
#include "COM_MemoryProxy.h"
|
|
|
|
class ReadBufferOperation : public NodeOperation {
|
|
private:
|
|
MemoryProxy *m_memoryProxy;
|
|
bool m_single_value; /* single value stored in buffer, copied from associated write operation */
|
|
unsigned int m_offset;
|
|
MemoryBuffer *m_buffer;
|
|
public:
|
|
ReadBufferOperation();
|
|
int isBufferOperation() { return true; }
|
|
void setMemoryProxy(MemoryProxy *memoryProxy) { this->m_memoryProxy = memoryProxy; }
|
|
MemoryProxy *getMemoryProxy() { return this->m_memoryProxy; }
|
|
void determineResolution(unsigned int resolution[2], unsigned int preferredResolution[2]);
|
|
|
|
void *initializeTileData(rcti *rect);
|
|
void executePixelSampled(float output[4], float x, float y, PixelSampler sampler);
|
|
void executePixelExtend(float output[4], float x, float y, PixelSampler sampler,
|
|
MemoryBufferExtend extend_x, MemoryBufferExtend extend_y);
|
|
void executePixelFiltered(float output[4], float x, float y, float dx[2], float dy[2], PixelSampler sampler);
|
|
const bool isReadBufferOperation() const { return true; }
|
|
void setOffset(unsigned int offset) { this->m_offset = offset; }
|
|
unsigned int getOffset() const { return this->m_offset; }
|
|
bool determineDependingAreaOfInterest(rcti *input, ReadBufferOperation *readOperation, rcti *output);
|
|
MemoryBuffer *getInputMemoryBuffer(MemoryBuffer **memoryBuffers) { return memoryBuffers[this->m_offset]; }
|
|
void readResolutionFromWriteBuffer();
|
|
void updateMemoryBuffer();
|
|
};
|
|
|
|
#endif
|