Matlab code I wrote for an orientational sensitive texture filter textureFilter.m computes the local motion orientation based on the Jacobian matrix and displays it with streaks, and a directional color-coded map.