|
|
| Line 1: |
Line 1: |
| The nodal functions used in inverse distance weighted interpolation can be higher degree polynomial functions constrained to pass through the scatter point and approximate the nearby points in a least squares manner. Quadratic polynomials have been found to work well in many cases (Franke & Nielson 1980; Franke 1982). The resulting surface reproduces local variations implicit in the data set, is smooth, and approximates the quadratic nodal functions near the scatter points. The equation used for the quadratic nodal function centered at point k is as follows:
| | #REDIRECT [[WMS:Inverse Distance Weighted Interpolation]] |
| | |
| | |
| [[Image:Idw_eq7.gif]]
| |
| | |
| | |
| To define the function, the six coefficients ak1..ak6 must be found. Since the function is centered at the point k and passes through point k, we know beforehand that ak1=fk where fk is the function value at point k. The equation simplifies to:
| |
| | |
| | |
| [[Image:Idw_eq8.gif]]
| |
| | |
| | |
| Now there are only five unknown coefficients. The coefficients are found by fitting the quadratic to the nearest NQ scatter points using a weighted least squares approach. In order for the matrix equation used to solve for the coefficients to be stable, there should be at least five scatter points in the set.
| |
| | |
| ===3D Interpolation===
| |
| For 3D interpolation, the equation for the quadratic nodal function is:
| |
| | |
| | |
| [[Image:WMSiidw_eq3.jpg]]
| |
| | |
| | |
| [[Image:WMSiidw_eq4.jpg]]
| |
| | |
| | |
| [[Image:WMSiidw_eq5.jpg]]
| |
| | |
| | |
| To define the function, the ten coefficients ak1..ak10 must be found. Since the function is centered on point k, we know that ak1=fk where fk is the data value at point k. The equation simplifies to:
| |
| | |
| | |
| [[Image:WMSiidw_eq6.jpg]]
| |
| | |
| | |
| [[Image:WMSiidw_eq7.jpg]]
| |
| | |
| | |
| [[Image:WMSiidw_eq8.jpg]]
| |
| | |
| | |
| Now there are only nine unknown coefficients. The coefficients are found by fitting the quadratic to a subset of the neighboring scatter points in a weighted least squares fashion. In order for the matrix equation used to be solve for the coefficients to be stable, there should be at least ten non-coplanar scatter point in the set.
| |
| | |
| | |
| ==Related Topics==
| |
| * [[WMS:Inverse Distance Weighted Interpolation|Inverse Distance Weighted Interpolation]]
| |
| * [[WMS:Computation of Interpolation Weights|Computation of Interpolation Weights]]
| |
| * [[WMS:Subset Definition |Subset Definition ]]
| |
| | |
| | |
| {{WMSMain}}
| |
