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
|
# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
include <imhdr.h>
include "imtext.h"
# RT_CONVERT_PIXELS -- Called once for each text file to be converted. All
# pixels in the text file are converted to image pixels.
procedure rt_convert_pixels (tf, im, format, pixels)
int tf # File descriptor of input text file
pointer im # Pointer to image header
int format # Format of text pixels (integer/floating point)
int pixels # Get pixels from input text file?
pointer bufptr, sp, word1, pattern
int stat, nlines, npix, i
long v[IM_MAXDIM], start
int impnll(), impnld(), impnlx()
int fscan(), stridxs(), patmatch(), patmake()
long note()
errchk impnll, impnld, impnlx
errchk rt_get_lineptr, rt_output_line, fscan, seek, amovkl
begin
# Determine if text file pixels were written with an integer, complex
# or floating point format. This information may have been already
# determined from the header. If not, the first pixel is read
# from text file. If it contains a decimal point, the character E,
# or a + or - sign not in the first position, it is a floating point
# number. Complex numbers are assumed to be written as "(r,i)".
if (pixels == YES && format == UNSET) {
call smark (sp)
call salloc (word1, SZ_LINE, TY_CHAR)
call salloc (pattern, SZ_LINE, TY_CHAR)
# Note position so we can return to it
start = note (tf)
stat = fscan (tf)
call gargwrd (Memc[word1], SZ_LINE)
if (patmake ("[DdEe]", Memc[pattern], SZ_LINE) == ERR)
call error (7, "Error creating format pattern")
if (stridxs ("(", Memc[word1]) > 0)
format = CPX_FORM
else if (stridxs (".", Memc[word1]) > 0)
format = FP_FORM
else if (patmatch (Memc[word1], Memc[pattern]) > 0)
format = FP_FORM
else if (stridxs ("+", Memc[word1]) > 1)
format = FP_FORM
else if (stridxs ("-", Memc[word1]) > 1)
format = FP_FORM
else
format = INT_FORM
call sfree (sp)
call seek (tf, start)
}
# Pixel type may not have been set by this point either...
if (IM_PIXTYPE(im) == UNSET) {
switch (format) {
case FP_FORM:
IM_PIXTYPE(im) = TY_REAL
case INT_FORM:
IM_PIXTYPE(im) = TY_INT
case CPX_FORM:
IM_PIXTYPE(im) = TY_COMPLEX
default:
call error (0, "Unrecognized format type")
}
}
nlines = 1
do i = 2, IM_NDIM(im)
nlines = nlines * IM_LEN (im, i)
call amovkl (long(1), v, IM_MAXDIM)
npix = IM_LEN (im, 1)
# Initialize text buffer
call rt_rinit
switch (IM_PIXTYPE(im)) {
case TY_SHORT, TY_INT, TY_USHORT, TY_LONG:
do i = 1, nlines {
stat = impnll (im, bufptr, v)
if (pixels == YES)
call rt_output_linel (tf, format, bufptr, npix)
else
call aclrl (Meml[bufptr], npix)
}
case TY_REAL, TY_DOUBLE:
do i = 1, nlines {
stat = impnld (im, bufptr, v)
if (pixels == YES)
call rt_output_lined (tf, format, bufptr, npix)
else
call aclrd (Memd[bufptr], npix)
}
case TY_COMPLEX:
do i = 1, nlines {
stat = impnlx (im, bufptr, v)
if (pixels == YES)
call rt_output_linex (tf, format, bufptr, npix)
else
call aclrx (Memx[bufptr], npix)
}
default:
call error (0, "Image pixel type unset")
}
end
|
