Description:
This routine takes data with a variance array and x y positions and regrids it onto a rectangular grid
using a spline interpolation algorithm (from PDA). Each integration is regridded separately and
added into a mean output image.
Invocation
CALL SCULIB_SPLINE_REGRID( METHOD, SFACTOR,
MAX_MAPS, N_MAPS, N_BOLS, N_INTS, DIAMETER, WAVELENGTH, PXSIZE, NX_OUT, NY_OUT,
I_CENTRE, J_CENTRE, WEIGHT, INT_PTR, DATA_PTR, VAR_PTR, XPOS_PTR,YPOS_PTR,
OUT_DATA, OUT_VARIANCE, OUT_QUALITY, OUT_WEIGHT, STATUS )
Arguments
METHOD =
CHARACTER (Given)
Spline method to use.
SFACTOR = REAL (Given)
Smoothing factor to
be used for PDA_SURFIT
MAX_MAPS = INTEGER (Given)
Max allowed number of
input maps (eg FILES in rebin)
N_MAPS = INTEGER (Given)
Number of input maps
N_BOLS( N_MAPS ) = INTEGER (Given)
Number of bolometers in each map
MAX_INTS =
INTEGER (Given)
Maximum number of separate integrations
N_INTS( N_MAPS )
= INTEGER (Given)
Number of integrations in each input map
DIAMETER = REAL
(Given)
Diameter of dish in metres
WAVELENGTH = REAL (Given)
Wavelength of
observation
PXSIZE = REAL (Given)
Pixel size of output map (radians)
NX_OUT =
INTEGER (Given)
Number of X pixels in output map
NY_OUT = INTEGER (Given)
Number of Y pixels in output map
I_CENTRE = INTEGER (Given)
Reference pixel
in X
J_CENTRE = INTEGER (Given)
Reference pixel in Y
WEIGHT( N_MAPS ) =
REAL (Given)
Relative weight of each input map
INT_LIST( MAX_MAPS, MAX_INTS
1 ) =
INTEGER (Given)
Position of each integration in data (cf DEM_PNTR for integrations)
DATA_PTR(
N_MAPS ) = INTEGER (Given)
Array of pointers to input data
VAR_PTR( N_MAPS ) = INTEGER
(Given)
Array of pointers to input data variance
XPOS_PTR( N_MAPS ) = INTEGER (Given)
Array of pointers to input data X positions
YPOS_PTR( N_MAPS ) = INTEGER (Given)
Array of
pointers to input data Y positions
OUT_DATA( NX_OUT, NY_OUT ) = REAL (Returned)
Regular
gridded data
OUT_VARIANCE( NX_OUT, NY_OUT ) = REAL (Returned)
Variance of
output data
OUT_DATA( NX_OUT, NY_OUT ) = BYTE (Returned)
Quality of output
map
OUT_WEIGHT( NX_OUT, NY_OUT) = REAL (Returned)
Total weight of each
pixel
STATUS = _INTEGER (Given & Returned)
Global status
Copyright
Copyright
©1995,1996,1997,1998,1999 Particle Physics and Astronomy Research Council. All Rights Reserved.