next up previous
Next: Relative Surface Orientation Up: Description of Three Dimensional Previous: Absolute Surface Area

Surface Elongation

The elongation of a surface region is also a distinguishing characteristic. This has been a traditional pattern recognition measurement applied to two dimensional regions, but a three dimensional version can also be obtained.

For a planar surface, elongation is the ratio of the longest to the shortest dimension of a planar rectangular box circumscribing the surface, which is similar to the two dimensional definition. The definition can be adapted to curved surfaces by using the maximum and minimum arc lengths of the intersection of the surface and a plane normal to the surface. The arc length should be calculated for all intersections of the surface, but this is time-consuming. Instead, only the cross-section widths about a central point are used.

The four factors involved in the estimation of the surface's dimensions are: the image region's dimensions, the surface slant relative to the viewer, the curvature of the surface, and the distance from the viewer. It is possible to approximate the elongation after separating the effects of these factors.

Figure 6.12: Cross-Section Length Distortions
\begin{figure}\epsfysize =3.45in
\epsfbox{FIGURES/Fig6.12.ps}\end{figure}

Figure 6.12 part (a) shows a sketch of the viewing relationship at one cross-section through the surface. By the discussion in the previous section, the cross-section length $S$ is approximately related to the chord length $L$ as:

\begin{displaymath}
S = L * \theta / (2 * sin(\theta / 2))
\end{displaymath}

Then, if the cross-section is slanted away from the viewer by an angle $ \alpha$ , the observed slanted length $L'$ is approximately related to the chord length $L$ (assuming the viewing distance is large) by:

\begin{displaymath}
L' = L * cos( \alpha )
\end{displaymath}

Finally, the observed image length $I$ for the surface at depth $D$ with conversion factor $G$ (as in here) is:

\begin{displaymath}
I = L' * (G / D)
\end{displaymath}

This analysis is then modified for observed slant compression at angles other than the tilt angle. Figure 6.12 part (b) shows the geometry used for the following analysis. This figure shows a unit circle compressed by a slant angle $\sigma $ in the direction $\tau$ and orthographically projected onto the image plane. Elementary trigonometry and algebra show that the observed length $R$ at the angle $\beta$ is:

\begin{displaymath}
R = 1 / \sqrt{1 + (tan( \sigma ) * cos(\beta))^2}
\end{displaymath}

The computation of the elongation value is then:

Let:  
  $ \vec{N}$ be a nominal point in the center of the surface image region
  $w( \alpha )$ be the image cross-section width at image angle $ \alpha$ about $ \vec{N}$
  $ \theta ( \alpha )$ be the change in surface orientation across
  the cross-section at image angle $ \alpha$
  $(P,Q,-T)$ be the unit surface normal at $ \vec{N}$
  $D$ be the distance from the viewer to $ \vec{N}$
  $G$ be the conversion factor for the number of pixels per centimeter of object length when seen at one centimeter distance
  $\vec{v}$ be the unit vector pointing to the viewer from the point $ \vec{N}$
   

Then, the tilt angle is:

\begin{displaymath}
\tau = arctan(Q / P)
\end{displaymath}

the relative slant direction $\beta$ is:

\begin{displaymath}
\beta = \alpha - \tau
\end{displaymath}

the slant angle $\sigma $ is:

\begin{displaymath}
\sigma = arccos( \vec{v} \circ (P,Q,-T))
\end{displaymath}

the slant correction factor is:

\begin{displaymath}
M = \sqrt{1 + (tan( \sigma ) * cos(\beta))^2}
\end{displaymath}

the projected chord length L' is:

\begin{displaymath}
L'( \alpha ) = w( \alpha ) * (D / G)
\end{displaymath}

the unprojected chord length $L$ is:

\begin{displaymath}
L( \alpha ) = L'( \alpha ) * M
\end{displaymath}

and the estimated three dimensional cross-section is:

\begin{displaymath}
cross( \alpha ) = L( \alpha ) * \theta ( \alpha ) / (2 * sin( \theta ( \alpha ) / 2))
\end{displaymath}

Finally, the elongation is:

\begin{displaymath}
E = max_{ \alpha }(cross( \alpha )) / min_{ \alpha }(cross( \alpha ))
\end{displaymath}

The elongations for all unobscured image regions that directly correspond to model SURFACEs are listed in Table 6.9.

Table 6.9: Summary of Estimated Elongations
IMAGE REGION PLANAR OR CURVED ESTIMATED ELONGATION TRUE ELONGATION
8 C 3.3 2.0
9 C 1.8 1.7
16 C 2.9 1.5
26 P 1.4 1.0
29 C 3.6 3.1

These results show that the estimation process gives approximate results when applied to unobscured regions. In part, small regions should be more affected because single pixel errors are significant, but this is not always the case.

A weakness of the process is that it only estimates the dimensions based on reconstructions about a single point, which produces lower bounds for the maximum and upper bounds for the minimum cross-sections. This should result in an estimate that is lower than the true elongation. However, the key source of error is the small size of the regions coupled with image quantization. Other sources of error are the hand measured surface data, surface interpolation and the curvature correction process assuming uniform curvature along the cross-section path. However, the approximations were felt to be justifiable on practical grounds and the above results show that the approach is acceptable.


next up previous
Next: Relative Surface Orientation Up: Description of Three Dimensional Previous: Absolute Surface Area
Bob Fisher 2004-02-26