test2/source/blender/compositor/operations/COM_KeyingScreenOperation.cpp

/*
 * Copyright 2012, 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
 *		Sergey Sharybin
 */

#include "COM_KeyingScreenOperation.h"

#include "MEM_guardedalloc.h"

#include "BLI_listbase.h"
#include "BLI_math.h"
#include "BLI_math_color.h"

#include "DNA_scene_types.h"

extern "C" {
	#include "BKE_movieclip.h"
	#include "BKE_tracking.h"

	#include "IMB_imbuf.h"
	#include "IMB_imbuf_types.h"
}

KeyingScreenOperation::KeyingScreenOperation(): NodeOperation()
{
	this->addOutputSocket(COM_DT_COLOR);
	this->movieClip = NULL;
	this->framenumber = 0;
	this->trackingObject[0] = 0;
	setComplex(true);
}

void KeyingScreenOperation::initExecution()
{
	initMutex();
	this->cachedTriangulation = NULL;
}

void KeyingScreenOperation::deinitExecution()
{
	if (this->cachedTriangulation) {
		TriangulationData *triangulation = cachedTriangulation;

		if (triangulation->triangulated_points)
			MEM_freeN(triangulation->triangulated_points);

		if (triangulation->triangles)
			MEM_freeN(triangulation->triangles);

		MEM_freeN(this->cachedTriangulation);

		this->cachedTriangulation = NULL;
	}
}

KeyingScreenOperation::TriangulationData *KeyingScreenOperation::buildVoronoiTriangulation()
{
	MovieClipUser user = {0};
	TriangulationData *triangulation;
	MovieTracking *tracking = &movieClip->tracking;
	MovieTrackingTrack *track;
	VoronoiSite *sites;
	ImBuf *ibuf;
	ListBase *tracksbase;
	ListBase edges = {NULL, NULL};
	int sites_total;
	int i;
	int width = this->getWidth();
	int height = this->getHeight();

	if (this->trackingObject[0]) {
		MovieTrackingObject *object = BKE_tracking_named_object(tracking, this->trackingObject);

		if (!object)
			return NULL;

		tracksbase = BKE_tracking_object_tracks(tracking, object);
	}
	else
		tracksbase = BKE_tracking_get_tracks(tracking);

	sites_total = BLI_countlist(tracksbase);

	if (!sites_total)
		return NULL;

	BKE_movieclip_user_set_frame(&user, framenumber);
	ibuf = BKE_movieclip_get_ibuf(movieClip, &user);

	if (!ibuf)
		return NULL;

	triangulation = (TriangulationData *) MEM_callocN(sizeof(TriangulationData), "keying screen triangulation data");

	sites = (VoronoiSite *) MEM_callocN(sizeof(VoronoiSite) * sites_total, "keyingscreen voronoi sites");
	track = (MovieTrackingTrack *) tracksbase->first;
	i = 0;
	while (track) {
		VoronoiSite *site = &sites[i];
		MovieTrackingMarker *marker = BKE_tracking_get_marker(track, framenumber);
		ImBuf *pattern_ibuf = BKE_tracking_get_pattern_imbuf(ibuf, track, marker, TRUE, FALSE);
		int j;

		zero_v3(site->color);
		for (j = 0; j < pattern_ibuf->x * pattern_ibuf->y; j++) {
			if (pattern_ibuf->rect_float) {
				add_v3_v3(site->color, &pattern_ibuf->rect_float[4 * j]);
			}
			else {
				unsigned char *rrgb = (unsigned char *)pattern_ibuf->rect;

				site->color[0] += srgb_to_linearrgb((float)rrgb[4 * j + 0] / 255.0f);
				site->color[1] += srgb_to_linearrgb((float)rrgb[4 * j + 1] / 255.0f);
				site->color[2] += srgb_to_linearrgb((float)rrgb[4 * j + 2] / 255.0f);
			}
		}

		mul_v3_fl(site->color, 1.0f / (pattern_ibuf->x * pattern_ibuf->y));
		IMB_freeImBuf(pattern_ibuf);

		site->co[0] = marker->pos[0] * width;
		site->co[1] = marker->pos[1] * height;

		track = track->next;
		i++;
	}

	IMB_freeImBuf(ibuf);

	BLI_voronoi_compute(sites, sites_total, width, height, &edges);

	BLI_voronoi_triangulate(sites, sites_total, &edges, width, height,
	                        &triangulation->triangulated_points, &triangulation->triangulated_points_total,
                            &triangulation->triangles, &triangulation->triangles_total);

	MEM_freeN(sites);
	BLI_freelistN(&edges);

	return triangulation;
}

void *KeyingScreenOperation::initializeTileData(rcti *rect, MemoryBuffer **memoryBuffers)
{
	if (this->movieClip == NULL)
		return NULL;

	if (this->cachedTriangulation)
		return this->cachedTriangulation;

	BLI_mutex_lock(getMutex());
	if (this->cachedTriangulation == NULL) {
		this->cachedTriangulation = buildVoronoiTriangulation();
	}
	BLI_mutex_unlock(getMutex());

	return this->cachedTriangulation;
}

void KeyingScreenOperation::determineResolution(unsigned int resolution[], unsigned int preferredResolution[])
{
	resolution[0] = 0;
	resolution[1] = 0;

	if (this->movieClip) {
		MovieClipUser user = {0};
		int width, height;

		BKE_movieclip_user_set_frame(&user, framenumber);
		BKE_movieclip_get_size(this->movieClip, &user, &width, &height);

		resolution[0] = width;
		resolution[1] = height;
	}
}

void KeyingScreenOperation::executePixel(float *color, int x, int y, MemoryBuffer *inputBuffers[], void *data)
{
	color[0] = 0.0f;
	color[1] = 0.0f;
	color[2] = 0.0f;
	color[3] = 1.0f;

	if (this->movieClip && data) {
		TriangulationData *triangulation = (TriangulationData *) data;
		int i;
		for (i = 0; i < triangulation->triangles_total; i++) {
			int *triangle = triangulation->triangles[i];
			VoronoiTriangulationPoint *a = &triangulation->triangulated_points[triangle[0]],
			                          *b = &triangulation->triangulated_points[triangle[1]],
			                          *c = &triangulation->triangulated_points[triangle[2]];
			float co[2] = {(float) x, (float) y}, w[3];

			if (barycentric_coords_v2(a->co, b->co, c->co, co, w)) {
				if (barycentric_inside_triangle_v2(w)) {
					color[0] += a->color[0] * w[0] + b->color[0] * w[1] + c->color[0] * w[2];
					color[1] += a->color[1] * w[0] + b->color[1] * w[1] + c->color[1] * w[2];
					color[2] += a->color[2] * w[0] + b->color[2] * w[1] + c->color[2] * w[2];
				}
			}
		}
	}
}