Large Area Distortion Calculation - IDMS V1.0 Section 13.3.4

Compute large area distortion per the International Committee for Display Metrology (ICDM) Information Display Measurements Standard (IDMS)

The function computes large area distortion based on section 13.3.4. in the ICDM IMDS publication. The published section has a few errors and limitations that have been addressed in this MATLAB function:
1) The published equations for pincushion distortion are missing a factor of 2 needed to average rather than sum the corner heights
2) The published equations assume the centerlines (horizontal and vertical) of the input coordinates are effectively the x and y axes of a cartesian coordinate system with the center of the image being at (0,0). The function input will automatically be transformed such that the center image point is at (0,0) for the purposes of the calculations.
3) In MATLAB, images have a y-axis value of zero at the top of the image and positive y in the downward direction, but the equations assume y is positive in the upward direction. Input data is assumed to follow the MATLAB convention and is flipped for the purposes of the calculations.
4) In the MATLAB function, the distances used for the quadratic distortions are calculated from the center vertical and horizontal lines determined from the linear distortion calculations and NOT from the default x and y axes. This removes the effect of any rotation between the test image distortion pattern and the pixel grid from the quadratic distortion calculation.
5) The published example has two issues. 1) It uses different coordinate labels than the published figures and equations. The MATLAB function inputs match the pictures and equations. 2) Some of the results in the published sample are calculated incorrectly. The MATLAB header includes the corrected results.