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
|
#!perl
# Testing CFITSIO read/write of fits headers
use strict;
use Test::More;
BEGIN {
eval "use Astro::FITS::CFITSIO qw / :longnames /; use Astro::FITS::CFITSIO qw/ :constants /;";
if ($@) {
plan skip_all => "Astro::FITS::CFITSIO module not available";
exit;
} else {
plan tests => 41;
}
}
use File::Spec;
use File::Copy;
# ----------------------------------------------------------------------------
# Copy the test data to a temporary location
BEGIN{ copy( File::Spec->catfile("t", "cfitsio.fit"),
File::Spec->catfile("t", "test.fit" ) )
or croak("Unable to copy test data file to temporary location"); };
# Delete the temporary copy of the test data
END{ unlink File::Spec->catfile("t", "test.fit" )
or croak("Unable to delete test data"); };
# ----------------------------------------------------------------------------
require_ok( "Astro::FITS::Header::CFITSIO" );
# Read from the __DATA__ block
my @cards = <DATA>;
chomp(@cards);
# Version test, need cfitsio library version > 2.1
ok( fits_get_version(my $version) > 2.1, "Check Astro::FITS::CFITSIO version number");
# get the path to the test data file
my $file = File::Spec->catfile("t", "test.fit");
# create a header
my $header = new Astro::FITS::Header::CFITSIO( File => $file );
# overwrite the header in memory with cards stored in __DATA__
# its a cheap and cheerful kludge and not the OO way to do it
$header->configure( Cards => \@cards );
# write our modified header out again
$header->writehdr( File => $file );
# open the header back up again reading directly from the file
my $compare = new Astro::FITS::Header::CFITSIO( File => $file );
# test the header against the raw data
for my $i (0 .. $#cards) {
my @items = $compare->item($i);
is( "$items[0]", $cards[$i], "Compare item $i");
}
exit;
__DATA__
SIMPLE = T / file does conform to FITS standard
BITPIX = -32 / number of bits per data pixel
NAXIS = 2 / number of data axes
NAXIS1 = 59 / length of data axis 1
NAXIS2 = 110 / length of data axis 2
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.
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
CRVAL1 = -0.03125 / Axis 1 reference value
CRPIX1 = 29.5 / Axis 1 pixel value
CTYPE1 = 'LINEAR ' / Quantity represented by axis 1
CRVAL2 = -0.03125 / Axis 2 reference value
CRPIX2 = 55.0 / Axis 2 pixel value
CTYPE2 = 'LINEAR ' / Quantity represented by axis 2
CD1_1 = 0.0625 / Axis rotation and scaling matrix
CD1_2 = 0.0 / Axis rotation and scaling matrix
CD2_1 = 0.0 / Axis rotation and scaling matrix
CD2_2 = 0.0625 / Axis rotation and scaling matrix
OBJECT = ' ' / Title of the dataset
DATE = '2000-12-12T20:28:35'/ file creation date (YYYY-MM-DDThh:mm:ss UTC)
ORIGIN = 'Starlink Project, U.K.'/ Origin of this FITS file
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
LINEAR = F / Linearisation disabled
FILTER = 'B1 ' / Combined filter name
FILTERS = 'B1 ' / Combined filter name
LAMP = 'off ' / Name of calibration lamp
END
|
