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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <syserr.h>
include <plset.h>
include <imhdr.h>
include <imio.h>
# IMGGSC -- Get a general section, any datatype (called by one of the typed
# procedures, which subsequently convert the datatype of the pixels returned
# by this routine). The mapping of the subraster in the input buffer to the
# imagefile is described by the section descriptor vectors VS and VE. Images
# of up to IM_MAXDIM dimensions are permitted, and each dimension may be
# accessed in either the forward or reverse direction.
pointer procedure imggsc (im, vs, ve, ndim, dtype, totpix)
pointer im # image descriptor
long vs[ARB], ve[ARB] # logical coords of corners of section
int ndim # dimensionality of section
int dtype # datatype of pixels desired
long totpix # total pixels in section (output)
bool rlio
pointer sp, px, bp, line, rl_high
long v[IM_MAXDIM], vinc[IM_MAXDIM]
long pvs[IM_MAXDIM], pve[IM_MAXDIM]
int sz_pixel, inbounds, npix, xstep, n
pointer imgibf()
int imsinb(), imloop(), pl_p2ri(), sizeof()
errchk imgibf, imrdpx, imrbpx
include <szpixtype.inc>
begin
#sz_pixel = sizeof(IM_PIXTYPE(im))
#sz_pixel = max ( sizeof(dtype), sizeof(IM_PIXTYPE(im)) )
#sz_pixel = pix_size[IM_PIXTYPE(im)]
sz_pixel = sizeof(IM_PIXTYPE(im))
rlio = (and (IM_PLFLAGS(im), PL_RLIO+PL_FAST) == PL_RLIO)
# Check that the section does not extend out of bounds.
inbounds = imsinb (im, vs, ve, ndim)
if (inbounds == ERR)
call imerr (IM_NAME(im), SYS_IMREFOOB)
# Get an (input) buffer to put the pixels into. Map the logical
# section into a physical section. Prepare the section descriptor
# do-loop index and increment vectors V, VINC.
bp = imgibf (im, vs, ve, ndim, dtype)
call imaplv (im, vs, pvs, ndim)
call imaplv (im, ve, pve, ndim)
call imsslv (im, pvs, pve, v, vinc, npix)
# A temporary pixel buffer is required for RLIO conversions.
if (rlio) {
call smark (sp)
call salloc (px, npix, TY_INT)
}
line = bp
totpix = 0
rl_high = bp - 1
# Read the section into the input buffer, line segment by line segment,
# advancing through the dimensions in storage order (leftmost subscript
# varies fastest).
repeat {
xstep = vinc[1]
# Convert the pixel array to a range list? (image masks). This is
# done more efficiently at a lower level if no complex geometric
# transformations are required (due to sections or OOB references).
if (rlio) {
if (inbounds == YES)
call imrdpx (im, Memi[px], npix, v, xstep)
else
call imrbpx (im, Memi[px], npix, v, xstep)
if (rl_high >= line)
call imerr (IM_NAME(im), SYS_IMRLOVFL)
else {
n = pl_p2ri (Memi[px], 1, Memc[line], npix)
rl_high = line + (n * RL_LENELEM * sz_pixel) - 1
}
} else {
if (inbounds == YES)
call imrdpx (im, Memc[line], npix, v, xstep)
else
call imrbpx (im, Memc[line], npix, v, xstep)
}
line = line + (npix * sz_pixel)
totpix = totpix + npix
} until (imloop (v, pvs, pve, vinc, IM_NPHYSDIM(im)) == LOOP_DONE)
if (rlio)
call sfree (sp)
return ((bp - 1) / sizeof(dtype) + 1)
end
|
