MEASUREMENT PRECISION AND ACCURACY OF THE CENTRE LOCATION OF AN ELLIPSE BY WEIGHTED CENTROID METHOD

• Main Author: R. Matsuoka
• Format: Article
• Language: English
• Published: Copernicus Publications 2015-03-01
• Series: ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences
• Online Access: http://www.isprs-ann-photogramm-remote-sens-spatial-inf-sci.net/II-3-W4/111/2015/isprsannals-II-3-W4-111-2015.pdf
• ISSN: 2194-9042; 2194-9050

Circular targets are often utilized in photogrammetry, and a circle on a plane is projected as an ellipse onto an oblique image. This paper reports an analysis conducted in order to investigate the measurement precision and accuracy of the centre location of an ellipse on a digital image by an intensity-weighted centroid method. An ellipse with a semi-major axis <i>a</i>, a semi-minor axis <i>b</i>, and a rotation angle &theta; of the major axis is investigated. In the study an equivalent radius <i>r</i> = (<i>a</i>²cos²&theta; + <i>b</i>²sin²&theta;)^1/2 is adopted as a measure of the dimension of an ellipse. First an analytical expression representing a measurement error (&varepsilon;_<i>x</i>, &varepsilon;_<i>y</i>,) is obtained. Then variances <i>V_x</i> of &varepsilon;_<i>x</i> are obtained at 1/256 pixel intervals from 0.5 to 100 pixels in <i>r</i> by numerical integration, because a formula representing <i>V_x</i> is unable to be obtained analytically when <i>r</i> > 0.5. The results of the numerical integration indicate that <i>V_x</i>would oscillate in a 0.5 pixel cycle in <i>r</i> and <i>V_x</i> excluding the oscillation component would be inversely proportional to the cube of <i>r</i>. Finally an effective approximate formula of <i>V_x</i> from 0.5 to 100 pixels in <i>r</i> is obtained by least squares adjustment. The obtained formula is a fractional expression of which numerator is a fifth-degree polynomial of {<i>r</i>−0.5×int(2<i>r</i>)} expressing the oscillation component and denominator is the cube of <i>r</i>. Here int(<i>x</i>) is the function to return the integer part of the value <i>x</i>. Coefficients of the fifth-degree polynomial of the numerator can be expressed by a quadratic polynomial of {0.5×int(2<i>r</i>)+0.25}.