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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <syserr.h>
include "mwcs.h"
# MW_SWSAMPD -- Set the sampled WCS curve for an axis.
procedure mw_swsampd (mw, axis, pv, wv, npts)
pointer mw #I pointer to MWCS descriptor
int axis #I axis which gets the wsamp vector
double pv[ARB] #I physical coordinates of points
double wv[ARB] #I world coordinates of points
int npts #I number of data point in curve
pointer wp
int mw_allocd()
errchk syserrs, mw_allocd
begin
# Get the current WCS.
wp = MI_WCS(mw)
if (wp == NULL)
call syserrs (SYS_MWNOWCS, "mw_swsampd")
# Overwrite the current curve, if any, else allocate new storage.
if (WCS_PV(wp,axis) == NULL || WCS_NPTS(wp,axis) < npts)
WCS_PV(wp,axis) = mw_allocd (mw, npts)
call amovd (pv, D(mw,WCS_PV(wp,axis)), npts)
if (WCS_WV(wp,axis) == NULL || WCS_NPTS(wp,axis) < npts)
WCS_WV(wp,axis) = mw_allocd (mw, npts)
call amovd (wv, D(mw,WCS_WV(wp,axis)), npts)
WCS_NPTS(wp,axis) = npts
end
|
