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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <gset.h>
include <mach.h>
include <imhdr.h>
include <mwset.h>
# T_PROW -- Plot an image row.
procedure t_prow ()
pointer image, wcslab, fmt
pointer im, mw, ct, sp, x_vec, y_vec
int row, ncols, nlines
real zmin, zmax
int clgeti()
pointer immap(), mw_openim(), mw_sctran()
begin
call smark (sp)
call salloc (image, SZ_FNAME, TY_CHAR)
call salloc (wcslab, SZ_LINE, TY_CHAR)
call salloc (fmt, SZ_LINE, TY_CHAR)
# Open image
call clgstr ("image", Memc[image], SZ_FNAME)
im = immap (Memc[image], READ_ONLY, 0)
call clgstr ("wcs", Memc[wcslab], SZ_LINE)
mw = mw_openim (im)
call mw_seti (mw, MW_USEAXMAP, NO)
ct = mw_sctran (mw, "logical", Memc[wcslab], 0)
ncols = IM_LEN(im,1)
nlines = IM_LEN(im,2)
row = clgeti ("row")
if (row < 1 || row > nlines) {
call imunmap (im)
call error (2, "line index references outside image")
}
# Get the requested row from the image.
call malloc (x_vec, ncols, TY_REAL)
call malloc (y_vec, ncols, TY_REAL)
call plt_grows (im, mw, ct, row, row, Memr[x_vec], Memr[y_vec],
zmin, zmax, Memc[wcslab], Memc[fmt], SZ_LINE)
# Now draw the vector to the screen.
call pr_draw_vector (Memc[image], Memr[x_vec], Memr[y_vec], ncols,
zmin, zmax, row, row, Memc[wcslab], Memc[fmt], false)
# Free resources.
call mfree (x_vec, TY_REAL)
call mfree (y_vec, TY_REAL)
call imunmap (im)
call sfree (sp)
end
|
