1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include "im2interpdef.h"
include <math/iminterp.h>
# MSIGRID -- Procedure to evaluate the interpolant on a rectangular
# grid. The procedure assumes that 1 <= x <= nxpix and 1 <= y <= nypix.
# The x and y vectors must be ordered such that x[i] < x[i+1] and
# y[i] < y[i+1].
procedure msigrid (msi, x, y, zfit, nx, ny, len_zfit)
pointer msi # pointer to interpolant descriptor structure
real x[ARB] # array of x values
real y[ARB] # array of y values
real zfit[len_zfit,ARB] # array of fitted values
int nx # number of x points
int ny # number of y points
int len_zfit # row length of zfit
errchk ii_grnearest, ii_grlinear, ii_grpoly3, ii_grpoly5, ii_grspline3
errchk ii_grsinc, ii_grlsinc, ii_grdirz
begin
switch (MSI_TYPE(msi)) {
case II_BINEAREST:
call ii_grnearest (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), x, y, zfit, nx, ny, len_zfit)
case II_BILINEAR:
call ii_grlinear (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), x, y, zfit, nx, ny, len_zfit)
case II_BIPOLY3:
call ii_grpoly3 (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), x, y, zfit, nx, ny, len_zfit)
case II_BIPOLY5:
call ii_grpoly5 (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), x, y, zfit, nx, ny, len_zfit)
case II_BISPLINE3:
call ii_grspline3 (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), x, y, zfit, nx, ny, len_zfit)
case II_BISINC:
call ii_grsinc (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), MSI_NYCOEFF(msi), x, y, zfit, nx, ny, len_zfit,
MSI_NSINC(msi), DX, DY)
case II_BILSINC:
call ii_grlsinc (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), MSI_NYCOEFF(msi), x, y, zfit, nx, ny, len_zfit,
LTABLE(MSI_LTABLE(msi)), 2 * MSI_NSINC(msi) + 1, MSI_NXINCR(msi),
MSI_NYINCR(msi), DX, DY)
case II_BIDRIZZLE:
call ii_grdriz (COEFF(MSI_COEFF(msi)), MSI_FSTPNT(msi),
MSI_NXCOEFF(msi), MSI_NYCOEFF(msi), x, y, zfit, nx, ny,
len_zfit, MSI_XPIXFRAC(msi), MSI_YPIXFRAC(msi),
MSI_BADVAL(msi))
}
end
|
