Nearest Neighbour Orientation Distribution

Hello, I want to find the distribution of the orientation of the undirected line connecting the typical point with its nearest neighbour. I have found something similar in R: https://rdrr.io/spatstat.core/man/orient.html and I am wondering how this could be done in Matlab. I have the pattern as an image, and I can extract the points and find their nearest neighbours using knnsearch .

[KSSV]

REad about knnsearch. Get the nearest points and do what you want.

[Image Analyst]

You forgot to attach the image. When you have an image, each pixel has 8 neighbors, so what do you consider as the "undirected line connecting the typical point with its nearest neighbour"??? There are 8 nearest neighbors, or at least 4. You imply there is only one.