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
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
|
# Astro::FITS::Header test harness -*-perl-*-
# strict
use strict;
#load test
use Test::More tests => 164;
# load modules
require_ok("Astro::FITS::Header");
require_ok("Astro::FITS::Header::Item");
# T E S T H A R N E S S --------------------------------------------------
# read header from DATA block
my @raw = <DATA>;
chomp(@raw);
# build header array
my $header = new Astro::FITS::Header( Cards => \@raw );
# test the header
for my $i (0 .. $#raw) {
my $card = $header->item($i);
$card->card( undef ); # clear cache
is( "$card", $raw[$i], "Compare card for keyword ". $card->keyword);
}
# See how many items we have of INT type
my @integers = $header->itembytype( "INT" );
is( scalar(@integers), 46, "Count number of INT keywords");
# build a test card
my $int_card = new Astro::FITS::Header::Item(
Keyword => 'LIFE',
Value => 42,
Comment => 'Life the Universe and everything',
Type => 'INT' );
# build another
my $string_card = new Astro::FITS::Header::Item(
Keyword => 'STUFF',
Value => 'Blah Blah Blah',
Comment => 'So long and thanks for all the fish',
Type => 'STRING' );
# and another
my $another_card = new Astro::FITS::Header::Item(
Keyword => 'VALUE',
Value => 34.5678,
Comment => 'A floating point number',
Type => 'FLOAT' );
# and one that contains embedded quotes
my $quote_card = new Astro::FITS::Header::Item(
Keyword=>'STRSTR',
Value => "She said 'Foo!' (\"really?\")",
Comment=> "It was 'foobar'.",
Type => 'STRING');
# Check quoting
my $qcstr = $quote_card->card;
my $qtstr1 = "STRSTR = 'She said ''Foo!'' (\"really?\")'";
my $qtstr2 = "/ It was 'foobar'.";
is(substr($qcstr,0,length($qtstr1)), $qtstr1,"Quote check 1");
is(substr($qcstr,index($qcstr,'/',length($qtstr1)),length($qtstr2)), $qtstr2,
"Quote check 2");
# insert
$header->insert(1, $int_card);
# value
my @test_value = $header->value('LIFE');
is($test_value[0], 42, "Value of LIFE");
# itembyname
my @itembyname = $header->itembyname('LIFE');
is("$int_card","$itembyname[0]","Check LIFE card");
# item
my @item = $header->item(1);
is("$int_card","$itembyname[0]", "Check item 1 is int_card");
# splice
my @cards = $header->splice( 0, 6, $string_card);
my @comp = ( $raw[0], $int_card, $raw[1], $raw[2], $raw[3], $raw[4] );
for my $i (0 .. $#cards) {
is( "$cards[$i]", "$comp[$i]","Splice removal");
}
my $first = $header->item(0);
is( "$first", "$string_card", "Check item 0");
$first = $header->splice(0,1);
is( "$first", "$string_card", "Check removed item 0");
# item
my $test_item = $header->item(1);
is( "$test_item", $raw[6], "Check item 1" );
# itembyname
my @comments = $header->itembyname('COMMENT');
is( scalar(@comments), 4, "Count number of comments");
for my $j (0 .. $#comments) {
is( "$comments[$j]", "$raw[$j+7]", "Compare comment $j");
}
# index
my @index = $header->index('COMMENT');
my @actual = (2,3,4,5);
for my $k (0 .. $#index ) {
is( $index[$k], $actual[$k], "Compare comment position $k" );
}
# insert
$header->insert(5, $string_card);
# comment
my @comment = $header->comment('STUFF');
is( "$comment[0]", "So long and thanks for all the fish", "Check comment");
# replacebyname
my $replacebyname = $header->replacebyname('STUFF', $int_card);
is("$string_card","$replacebyname","Replaced STUFF by name");
# replace
my $replace = $header->replace(5, $another_card);
is("$int_card","$replace", "Replacement by index");
# value
my @floating = $header->value('VALUE');
is($floating[0],34.5678, "Got first VALUE");
# remove
my $remove = $header->remove(5);
is("$another_card","$remove","Removed by position");
@floating = $header->value('VALUE');
is($floating[0],undef,"Got no VALUE");
# insert
$header->insert(5, $string_card);
# removebyname
my $removebyname = $header->removebyname('STUFF');
is("$string_card","$removebyname","Was STUFF removed");
# check regular expressions
@index = $header->index( qr/CLOCK\d/ );
@actual = (53..59);
is( scalar @index, scalar @actual, "Count expected number of CLOCK matches" );
while( @index )
{
is( shift @index, shift @actual, "Compare CLOCK keyword location" );
}
exit;
__DATA__
SIMPLE = T / file does conform to FITS standard
BITPIX = -32 / number of bits per data pixel
NAXIS = 3 / number of data axes
NAXIS1 = 25 / length of data axis 1
NAXIS2 = 36 / length of data axis 2
NAXIS3 = 252 / length of data axis 3
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format defined in Astronomy and
COMMENT Astrophysics Supplement Series v44/p363, v44/p371, v73/p359, v73/p365.
COMMENT Contact the NASA Science Office of Standards and Technology for the
COMMENT FITS Definition document #100 and other FITS information.
CRVAL1 = -0.07249999791383749 / Axis 1 reference value
CRPIX1 = 12.5 / Axis 1 pixel value
CTYPE1 = 'a1 ' / LINEAR
CRVAL2 = -0.07249999791383743 / Axis 2 reference value
CRPIX2 = 18.0 / Axis 2 pixel value
CTYPE2 = 'a2 ' / LINEAR
CRVAL3 = 1.27557086671004E-6 / Axis 3 reference value
CRPIX3 = 126.0 / Axis 3 pixel value
CTYPE3 = 'a3 ' / LAMBDA
OBJECT = 'galaxy ' / Title of the dataset
DATE = '2000-12-13T22:44:53' / file creation date (YYYY-MM-DDThh:mm:ss UTC)
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 1999' / FITS file originator
BSCALE = 1.0 / True_value = BSCALE * FITS_value + BZERO
BZERO = 0.0 / True_value = BSCALE * FITS_value + BZERO
HDUCLAS1= 'NDF ' / Starlink NDF (hierarchical n-dim format)
HDUCLAS2= 'DATA ' / Array component subclass
IRAF-TLM= '23:07:26 (27/02/2000)' / Time of last modification
TELESCOP= 'UKIRT, Mauna Kea, HI' / Telescope name
INSTRUME= 'CGS4 ' / Instrument
OBSERVER= 'SMIRF ' / Observer name(s)
OBSREF = '? ' / Observer reference
DETECTOR= 'fpa046 ' / Detector array used
OBSTYPE = 'OBJECT ' / Type of observation
INTTYPE = 'STARE+NDR' / Type of integration
MODE = 'ND_STARE' / Observing mode
GRPNUM = 0 / Number of observation group
RUN = 54 / Number of run
EXPOSED = 180 / Total exposure time for integration
OBJCLASS= 0 / Class of observed object
CD1_1 = 0.144999980926513672 / Axis rotation and scaling matrix
CD1_2 = 0.0 / Axis rotation and scaling matrix
CD1_3 = 0.0 / Axis rotation and scaling matrix
CD2_1 = 0.0 / Axis rotation and scaling matrix
CD2_2 = 0.144999980926513672 / Axis rotation and scaling matrix
CD2_3 = 0.0 / Axis rotation and scaling matrix
CD3_1 = 0.0 / Axis rotation and scaling matrix
CD3_2 = 0.0 / Axis rotation and scaling matrix
CD3_3 = 2.07933226192836E-10 / Axis rotation and scaling matrix
MEANRA = 10.34629999999999939 / Object RA at equinox (hrs)
MEANDEC = 20.1186000000000007 / Object Dec at equinox (deg)
RABASE = 10.34629999999999939 / Offset zero-point RA at equinox (hrs)
DECBASE = 20.1186000000000007 / Offset zero-point Dec at equinox (deg)
RAOFF = 0 / Offset RA at equinox (arcsec)
DECOFF = 0 / Offset Dec at equinox (arcsec)
DROWS = 178 / No of det. in readout row
DCOLUMNS= 256 / No of det. in readout column
DEPERDN = 6 / Electrons per data number
CLOCK0 = -6.20000000000000018 / ALICE CLOCK0 voltage
CLOCK1 = -3 / ALICE CLOCK1 voltage
CLOCK2 = -7.5 / ALICE CLOCK2 voltage
CLOCK3 = -2.79999999999999982 / ALICE CLOCK3 voltage
CLOCK4 = -6 / ALICE CLOCK4 voltage
CLOCK5 = -2 / ALICE CLOCK5 voltage
CLOCK6 = -7.5 / ALICE CLOCK6 voltage
VSLEW = 4 / ALICE VSLEW voltage
VDET = -3.02000000000000002 / ALICE VDET voltage
DET_BIAS= 0.57999999999999996 / ALICE DET_BIAS voltage
VDDUC = -3.60000000000000009 / ALICE VDDUC voltage
VDETGATE= -4.5 / ALICE VDETGATE voltage
VGG_A = -1.60000000000000009 / ALICE VGG_ACTIVE voltage
VGG_INA = -1.30000000000000004 / ALICE VGG_INACTIVE voltage
VDDOUT = -1 / ALICE VDDOUT voltage
V3 = -2.79999999999999982 / ALICE V3 voltage
VLCLR = -3 / ALICE VLCLR voltage
VLD_A = 4 / ALICE VLOAD_ACTIVE voltage
VLD_INA = 4 / ALICE VLOAD_INACTIVE voltage
WFREQ = 1 / ALICE waveform state freq. (MHz)
RESET_DL= 0.200000000000000011 / NDR reset delay (seconds)
CHOP_DEL= 0.029999998999999999 / Chop delay (seconds)
READ_INT= 5 / NDR read interval (seconds)
NEXP_PH = 0 / Exposures in each chop phase
DEXPTIME= 180 / Exposure time (seconds)
RDOUT_X1= 1 / Start column of array readout
RDOUT_X2= 256 / End column of array readout
RDOUT_Y1= 45 / Start row of array readout
RDOUT_Y2= 222 / End row of array readout
CHOPDIFF= T / Main-offset beam value stored
IF_SHARP= F / Shift & add disabled
LINEAR = F / Linearisation disabled
FILTER = 'B1 ' / Combined filter name
FILTERS = 'B1 ' / Combined filter name
DETINCR = 1 / Increment (pixels) betw scan positions
DETNINCR= 2 / Number of scan positions in scan
WPLANGLE= 0 / IRPOL waveplate angle
SANGLE = -2.19303900000000018 / Angle of slit
SLIT = '0ew ' / Name of slit
SLENGTH = 18 / Length of slit
SWIDTH = 4 / Width of slit
DENCBASE= 800 / Zeropoint (steps) of detector translation
DFOCUS = 1.819309999999999983 / Detector focus position
GRATING = '150_lpmm' / Name of grating
GLAMBDA = 1.274947000000000052 / Grating wavelength
GANGLE = 17.09262000000000015 / Grating wavelength
GORDER = 3 / Grating order
GDISP = 0.00020796522 / Grating dispersion
CNFINDEX= 75488 / Index increments when h/w config changes
CVF = 'open ' / Name of CVF
CLAMBDA = 0 / CVF wavelength
IRTANGLE= 6.396519999999999762 / Image rotator angle
LAMP = 'off ' / Name of calibration lamp
BBTEMP = 0 / Black body temperature
CALAPER = 0 / Aperture of tungsten-halogen lamp (%)
THLEVEL = 0 / Level of tungsten-halogen lamp
IDATE = 19980217 / Date as integer
OBSNUM = 54 / Number of observation
NEXP = 1 / Exposures in integration
AMSTART = 1.334643999999999942 / Airmass at start of obs
AMEND = 1.320149999999999935 / Airmass at end of obs
RUTSTART= 8.000171999999999173 / Start time of obs (hrs)
RUTEND = 8.101883000000000834 / End time of obs (hrs)
NBADPIX = 32
END
|
