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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <imhdr.h>
include <imset.h>
include <plset.h>
include <plio.h>
# PL_LOADIM -- Load a mask stored as a conventional image, i.e., convert an
# image to a mask.
procedure pl_loadim (pl, imname, title, maxch)
pointer pl #I mask descriptor
char imname[ARB] #I image name or section
char title[ARB] #O mask title
int maxch #I max chars out
bool sampling
pointer im, px, im_pl, bp
long vs_l[PL_MAXDIM], vs_p[PL_MAXDIM]
long ve_l[PL_MAXDIM], ve_p[PL_MAXDIM]
int npix, naxes, maxdim, maxval, depth, sz_buf, i
long v_in[PL_MAXDIM], v_out[PL_MAXDIM], vn[PL_MAXDIM]
pointer immap()
int imgnli(), imstati(), im_pmsvhdr()
errchk immap, imgnli
begin
# Open the input image.
im = immap (imname, READ_ONLY, 0)
# Encode and output the image header.
bp = NULL; i = im_pmsvhdr (im, bp, sz_buf)
call strcpy (Memc[bp], title, maxch)
call mfree (bp, TY_CHAR)
# Determine the mask depth in bits.
maxval = IM_MAX(im)
if (maxval <= 0)
maxval = I_PVMAX
do i = 0, ARB
if (2**i > min (I_PVMAX, maxval)) {
depth = i
break
}
# Initialize the mask to the size of the image section.
npix = IM_LEN(im,1)
naxes = IM_NDIM(im)
maxdim = min (IM_MAXDIM, PL_MAXDIM)
call amovl (IM_LEN(im,1), vn, maxdim)
call pl_ssize (pl, naxes, vn, depth)
# If the image is already a mask internally, check whether any
# subsampling, axis flipping, or axis mapping is in effect.
# If so we can't use PLIO to copy the mask section.
im_pl = imstati (im, IM_PLDES)
sampling = false
if (im_pl != NULL) {
call amovkl (long(1), vs_l, maxdim)
call amovl (IM_LEN(im,1), ve_l, maxdim)
call imaplv (im, vs_l, vs_p, maxdim)
call imaplv (im, ve_l, ve_p, maxdim)
do i = 1, maxdim {
vn[i] = ve_l[i] - vs_l[i] + 1
if (vn[i] != ve_p[i] - vs_p[i] + 1) {
sampling = true
break
}
}
}
# If the source image is already a mask internally and no axis
# geometry is in effect in the image section (if any), then we can
# use a PLIO rasterop to efficiently copy the mask subsection.
if (im_pl != NULL && !sampling) {
# Copy a mask subsection (or entire image if no section).
call pl_rop (im_pl, vs_p, pl, vs_l, vn, PIX_SRC)
} else {
# Copy image pixels. Initialize the vector loop indices.
call amovkl (long(1), v_in, maxdim)
call amovkl (long(1), v_out, maxdim)
# Copy the image.
while (imgnli (im, px, v_in) != EOF) {
call pl_plpi (pl, v_out, Memi[px], 0, npix, PIX_SRC)
call amovl (v_in, v_out, maxdim)
}
}
call pl_compress (pl)
call imunmap (im)
end
|
