1 For every two consecutive video frames I1,I2 do
2 D(x,y,t)=I1(x,y,t)-I2(x,y,t-1)
3 A(x,y,t)=D(x,y,t) \\ with enhanced Local Motion as mentioned in Algorithm 1.
4 F(v,w,t)=FT(A(x,y,t)) \\ FT denotes 2D fourier transform.
5 Apply Radon transform, R(d,θ)=Radon(f(v,w,t))
6 Normalize, P(θ)=max(R(d,θ))
7 K=Kurtosis(R(d,θ))
8 M=Mean(R(d,θ))
9 EndFor
10 Return Histogram(P, nbins), Histogram(K,nbins), Histogram(M,nbins).