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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <mach.h>
include <mwset.h>
# T_RIMCURSOR -- Read the image cursor list until EOF is seen on the list,
# transforming the coordinates to the desired system and printing the
# transformed cursor reads on the standard output.
procedure t_rimcursor()
double px, py, wx, wy
pointer gp, ct, mw, sp, wcs, rest, format[2], fmt[2]
int axis, frame, newframe, ntokens, wcscode
pointer gopen(), rim_getctran()
int clscan(), nscan()
errchk gopen, rim_getctran
begin
# Allocate working space.
call smark (sp)
call salloc (wcs, SZ_FNAME, TY_CHAR)
call salloc (rest, SZ_LINE, TY_CHAR)
do axis = 1, 2 {
call salloc (format[axis], SZ_FNAME, TY_CHAR)
call salloc (fmt[axis], SZ_FNAME, TY_CHAR)
}
# Open graphics context (doesn't work currently for stdimage).
iferr (gp = gopen ("stdimage", APPEND, STDIMAGE))
gp = NULL
# Initialize.
frame = 1
mw = NULL
ct = NULL
call clgstr ("wxformat", Memc[format[1]], SZ_FNAME)
call clgstr ("wyformat", Memc[format[2]], SZ_FNAME)
# Read the cursor repeatedly until EOF is seen.
while (clscan ("cursor") != EOF) {
# Get cursor value.
call gargd (px)
call gargd (py)
call gargi (wcscode)
newframe = wcscode / 100
call gargstr (Memc[rest], SZ_LINE)
# Get coordinate transformation.
ntokens = nscan()
if (ntokens < 2)
call error (1, "bad cursor read")
else if (mw == NULL || (ntokens >= 3 && newframe != frame)) {
if (mw != NULL)
call mw_close (mw)
ct = rim_getctran (newframe, Memc[wcs], mw)
if (mw != NULL)
call mw_ssystem (mw, Memc[wcs])
frame = newframe
do axis = 1, 2 {
if (Memc[format[axis]] != EOS)
call strcpy (Memc[format[axis]], Memc[fmt[axis]],
SZ_FNAME)
else if (mw != NULL) {
iferr (call mw_gwattrs (mw, axis, "format",
Memc[fmt[axis]], SZ_FNAME))
call strcpy ("%0.15g", Memc[fmt[axis]], SZ_FNAME)
} else
call strcpy ("%0.15g", Memc[fmt[axis]], SZ_FNAME)
}
}
# Transform coordinates.
if (ct != NULL)
call mw_c2trand (ct, px,py, wx,wy)
else {
wx = px
wy = py
}
# Always output transformed coordinates.
do axis = 1, 2 {
# Output the transformed value.
if (axis == 1) {
call fprintf (STDOUT, Memc[fmt[axis]])
call pargd (wx)
} else {
call fprintf (STDOUT, Memc[fmt[axis]])
call pargd (wy)
}
call putci (STDOUT, ' ')
}
# Output WCS field if present in input.
if (ntokens > 2) {
call fprintf (STDOUT, "%d")
call pargi (wcscode)
}
# Output rest of cursor value if present in input.
if (ntokens > 3)
call putline (STDOUT, Memc[rest])
call putci (STDOUT, '\n')
}
# Shutdown.
if (mw != NULL)
call mw_close (mw)
if (gp != NULL)
call gclose (gp)
call sfree (sp)
end
# RIM_GETCTRAN -- Get the WCS context to be used for coordinate transforms.
# Determine the reference image, load the WCS, and compile a transform from
# logical image coordinates to the desired coordinate system.
pointer procedure rim_getctran (frame, wcs, mw)
int frame #I image frame, if display is used
char wcs[ARB] #O name of requested wcs
pointer mw #O MWCS descriptor
int status
bool use_display
pointer sp, imname, ds, iw, im, ct
pointer imd_mapframe(), iw_open(), immap(), mw_sctran(), mw_openim()
errchk imd_mapframe, iw_open, immap, mw_sctran, mw_openim
int envfind(), clgeti()
bool streq()
begin
call smark (sp)
call salloc (imname, SZ_LINE, TY_CHAR)
# Access image display to get name of reference image?
if (clgeti ("$nargs") > 0)
use_display = false
else if (envfind ("stdimcur", Memc[imname], SZ_LINE) > 0)
use_display = streq (Memc[imname], "stdimage")
else
use_display = false
# Get the name of the reference image.
if (use_display) {
ds = imd_mapframe (frame, READ_ONLY, NO)
iw = iw_open (ds, frame, Memc[imname], SZ_LINE, status)
call iw_close (iw)
call imunmap (ds)
} else
call clgstr ("image", Memc[imname], SZ_LINE)
# Map the image if one was specified.
if (Memc[imname] == EOS)
im = NULL
else iferr (im = immap (Memc[imname], READ_ONLY, 0))
im = NULL
# Get WCS if image was given.
if (im == NULL) {
mw = NULL
ct = NULL
wcs[1] = EOS
} else {
ifnoerr (mw = mw_openim (im)) {
call clgstr ("wcs", wcs, SZ_FNAME)
if (wcs[1] == EOS)
call strcpy ("logical", wcs, SZ_FNAME)
ct = mw_sctran (mw, "logical", wcs, 03B)
} else {
mw = NULL
ct = NULL
}
# We don't need the image once the WCS is loaded.
call imunmap (im)
if (mw != NULL)
call mw_seti (mw, MW_REFIM, NULL)
}
call sfree (sp)
return (ct)
end
|
