directionsPerPixel = 4;
dirM = [ cos((2 * pi) / directionsPerPixel) sin((2 * pi) / directionsPerPixel); - sin((2 * pi) / directionsPerPixel) cos((2 * pi) / directionsPerPixel) ]

lPos = [8 5 3]
pPos = [10 3 6]

