To find the major and minor axis for disk sources we can use the raw and central image moments within the effective radius of the source.

More on image moments can be found here and a great walkthrough for application can be found here.

The raw image moments are defined as:

Where $I(x,y)$ is the value at pixel $x,y$.

The central image moments we use are:

To get the major and minor axis of the source we need to construct a covariance matrix for the pixels using the second order central moments.

Where :

The eigenvectors of this matrix correspond to the major and minor axis of the source.


The source we’ll use is Deep 2 10064:


With the segmap


Here are the major(red) and minor(blue) axis as measured from the h band plotted over the segmap with the effective radius highlighted in yellow.