Oncologic therapy of laryngeal cancer may necessitate a total excision of the larynx which results in loss of voice. Voice rehabilitation can be achieved using mucosal tissue vibrations at the upper part of the esophagus which serves as substitute voice generating element (PE segment). The quality of the substitute voice is closely related to vibratory characteristics of the PE segment. By means of a high-speed camera the dynamics of the PE segment can be recorded in real-time. Using image processing the deformations of the PE segment are extracted from the image series as deforming contours. Commonly, the characterization of PE dynamics bases on the spectral analysis of the time varying contour area. However, this constitutes an integral approach which masks most of the specific dynamics of PE deformations.
We present an algorithm that automatically registers one segmented contour in a frame of the video sequence to the contour in the next frame to derive discrete 2-D trajectories of PE vibrations. By concatenation of the obtained transformations this approach provides a total registration of PE segment contours. We suggest a mixed-integer programming formulation for the problem that combines an advanced outlier and deformation handling with the introduction of dummy points in regions that newly open up, and that includes normal information in the objective function to avoid unwanted deformations. Numerical experiments show that the implemented alternate convex search algorithm produces robust results which is demonstrated at the example of five high-speed recordings of laryngectomee subjects.