Compute the dominant direction in the patch.

Hello,i am trying to register CT and MR images of spine. I have applied harris corner points algorithm to extract corner points. now i want to find dominant direction in each image patch centered around the corner point. how can i do that?