anaglyph_algorithm Namespace Reference

Classes
class	Anaglyph
	This class represents an algorithm to compute two stereo images from gimp layers and to combine them to an anaglyph. More...

Functions
def	interp2 (X, Y, Z)
	interpolates Z at locations X,Y output of function is of shape(X)==shape(Y) More...

Function Documentation

def anaglyph_algorithm.interp2	(		X,
			Y,
			Z
	)

interpolates Z at locations X,Y output of function is of shape(X)==shape(Y)

Definition at line 5 of file anaglyph_algorithm.py.

def interp2(X,Y,Z):
    """! interpolates Z at locations X,Y
    output of function is of shape(X)==shape(Y)
    """
    from numpy import *
    X1 = array(floor(X),dtype=int);
    X2 = array(ceil(X),dtype=int);
    Y1 = array(floor(Y),dtype=int);
    Y2 = array(ceil(Y),dtype=int);

    M1 = logical_or(logical_or(Y1<0,Y1>=Z.shape[0]),logical_or(Y2<0,Y2>=Z.shape[0]))
    M2 = logical_or(logical_or(X1<0,X1>=Z.shape[1]),logical_or(X2<0,X2>=Z.shape[1]))
    M = logical_not(logical_or(M1,M2));

    X1[X1<0]=0;
    X1[X1>=Z.shape[1]]=Z.shape[1]-1;
    X2[X2<0]=0;
    X2[X2>=Z.shape[1]]=Z.shape[1]-1;
    XR = X-X1;
    Y1[Y1<0]=0;
    Y1[Y1>=Z.shape[0]]=Z.shape[0]-1;
    Y2[Y2<0]=0;
    Y2[Y2>=Z.shape[0]]=Z.shape[0]-1;
    YR = Y-Y1;
    ERG1 = Z[Y1,X1]+(Z[Y1,X2]-Z[Y1,X1])*XR;
    ERG2 = Z[Y2,X1]+(Z[Y2,X2]-Z[Y2,X1])*XR;
    ERG  = ERG1+(ERG2-ERG1)*YR;
    return ERG,M;