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
|
include <imhdr.h>
include "ms.h"
# T_NEW_EXTRACTION -- Create a new extraction database.
#
# This is the first step in using the MULTISPEC package. The new database
# may be created from scratch or intialized from an template image.
procedure t_new_extraction ()
# Task parameters:
char image[SZ_FNAME] # Multi-spectra image
char template[SZ_FNAME] # Template image
int samples[3, MAX_RANGES] # Sample line range array
char comment[SZ_LINE]
char database[SZ_FNAME], old_database[SZ_FNAME]
pointer im, ms
bool strne()
int clgranges(), expand_ranges()
pointer immap(), msmap()
begin
# Get database and image name. Map the image and check that
# it is two dimensional.
call clgstr ("image", image, SZ_FNAME)
im = immap (image, READ_ONLY, 0)
if (IM_NDIM(im) != 2)
call error (MS_ERROR, "Image file must be two dimensional.")
# Get the template image name.
call clgstr ("template", template, SZ_FNAME)
if (strne (template, "")) {
# If a template is given then map the template and check
# that the new image dimensions agree with the old dimensions.
ms = msmap (template, READ_ONLY, 0)
if ((MS_LEN(ms, 1) != IM_LEN(im, 1)) ||
(MS_LEN(ms, 2) != IM_LEN(im, 2)))
call error (MS_ERROR,
"New image size does not agree with the old image size.")
call msunmap (ms)
# Copy the old database. Map the new database and clear the
# the old comments before adding a new comment.
call sprintf (database, SZ_FNAME, "%s.db")
call pargstr (image)
call sprintf (old_database, SZ_FNAME, "%s.db")
call pargstr (template)
call fcopy (old_database, database)
ms = msmap (image, READ_WRITE, 0)
COMMENT(ms, 1) = EOS
call sprintf (comment, SZ_LINE,
"Database initialized from the template image %s.")
call pargstr (template)
call history (ms, comment)
} else {
# For a new database initialize the header parameters.
ms = msmap (image, NEW_FILE, MS_DB_ENTRIES)
MS_LEN(ms, 1) = IM_LEN(im, 1)
MS_LEN(ms, 2) = IM_LEN(im, 2)
MS_NSPECTRA(ms) = 0
# Get the sample line ranges and set the number of sample lines
# in the database header.
MS_NSAMPLES(ms) = clgranges ("sample_lines", 1, MS_LEN (ms, 2),
samples, MAX_RANGES)
# Make an entry in the database for the sample lines and then
# access the entry in order to allocate memory for the sample
# line array.
call dbenter (MS_DB(ms), NAME(ms, SAMPLE), MS_NSAMPLES(ms)*SZ_INT,1)
call msgsample (ms)
# Expand the sample line range array into the sample line array.
# Then put the sample line array in the database.
MS_NSAMPLES(ms) = expand_ranges (samples, LINE(ms, 1),
MS_NSAMPLES(ms))
call mspsample (ms)
# Add a history line.
call history (ms, "New MULTISPEC database created.")
}
# Set the image name and image title in the database.
call strcpy (image, MS_IMAGE(ms), SZ_MS_IMAGE)
call strcpy (IM_TITLE(im), MS_TITLE(ms), SZ_MS_TITLE)
# Close image and database. Write the database header record before
# closing the database.
call imunmap (im)
call msphdr (ms)
call msunmap (ms)
end
|
