Using sin and cosine of the radians can produce the north-southness
and east-westness but how do you handle flat values which are depicted
as -1's in a GIS raster? One could interpolate the edges of flat
areas to a sloped center bottom cell based on the dem but not sure if
there is a routine available. Perhaps this question motivates the
search for an aml or other rotune that will perform the calculations.
I've treated N, E, S, W, and flat as a categorical variable and
produced meaningful results.
A. Yost