Glacier surface motion computation from digital image sequences

A technique for computing the field of short-term glacier surface motion on a local scale is presented. Time-lapsed image negatives, digitized to a high resolution, provide the raw data for the three-stage technique, First, cross-correlation is used to establish a number of candidate displacement vectors for a series of regularly spaced templates. A relaxation-labeling routine is then applied to select the most appropriate candidate vectors, according to the local flow. Novel aspects of the relaxation algorithm include a new, efficient form of the support function and the absence of a null-match category. A new development is the application of a past filter to the resultant flow field? providing suitable displacement vectors for templates that were originally unmatched and correcting vectors that are still inconsistent with the local flow. Results from an image sequence from New Zealand's Mount Cook National Park show the superiority of the technique over the maximum cross-correlation method and demonstrate the effectiveness of the post filter in improving correlation-relaxation labeling.