rfits 0.0.1

data/test/rfits_test.rb

@@ -0,0 +1,310 @@
+require 'test/unit'
+require 'fileutils'
+require 'rubygems' rescue LoadError
+
+require 'rfits'
+
+class RFitsTest < Test::Unit::TestCase
+  def setup
+    @fits = RFits::File.new('!test/readwrite.test.fits', 'rw')
+  end
+  
+  def test_all
+    assert_equal 'test/readwrite.test.fits', @fits.path
+    assert_equal RFits::IO::Proxy::READWRITE, @fits.mode
+    assert_equal 'file://', @fits.url_type
+    
+    hdu()
+    header()
+    image_data()
+    table_data()
+  end
+  
+  def hdu()
+    assert_equal 0, @fits.length
+    im1 = RFits::Image.new(@fits, nil, RFits::IO::Proxy::LONG_IMG, [4, 4],
+      :compression_type => :hcompress,
+      :tile_dim => [4, 4],
+      :hcomp_scale => 1,
+      :hcomp_smooth => 1
+    )
+    btbl = RFits::BinaryTable.new(@fits, nil,
+            ['column_1', 'column_2', 'column_3', 'column_4'],
+            ['I', '10A', 'E', 'C'], 'Binary_Table_Test')
+    im2 = RFits::Image.new(@fits, nil, RFits::IO::Proxy::FLOAT_IMG, [3, 3])
+    im3 = RFits::Image.new(@fits, nil, RFits::IO::Proxy::LONG_IMG, [2, 2],
+      :compression_type => :rice,
+      :tile_dim => [2, 2]
+    )
+
+    @fits[0] = im1
+    @fits[2] = btbl
+    @fits << im2
+    @fits << im3
+    
+    img1 = @fits[1]
+    assert_equal :image, img1.hdu_type
+    assert_kind_of RFits::Image, img1
+    assert_equal 1, img1.position
+    assert_equal [4, 4], img1.naxes
+    assert_equal 2, img1.naxis
+    assert_equal :long, img1.image_type
+    assert_equal RFits::IO::Proxy::LONG_IMG, img1.bitpix
+    assert_equal true, img1.compressed?
+    assert_equal :hcompress, img1.compression_type
+    assert_equal [4, 4], img1.tile_dim
+    assert_equal nil, img1.noise_bits
+    assert_equal 1, img1.hcomp_scale
+    assert_equal 1, img1.hcomp_smooth
+    
+    tbl2 = @fits[2]
+    assert_equal :binary_tbl, tbl2.hdu_type
+    assert_kind_of RFits::BinaryTable, tbl2
+    assert_equal 2, tbl2.position
+    
+    img3 = @fits[3]
+    assert_equal :image, img3.hdu_type
+    assert_kind_of RFits::Image, img3
+    assert_equal 3, img3.position
+    assert_equal [3, 3], img3.naxes
+    assert_equal 2, img3.naxis
+    assert_equal :float, img3.image_type
+    assert_equal RFits::IO::Proxy::FLOAT_IMG, img3.bitpix
+    assert_equal false, img3.compressed?
+    
+    img4 = @fits[4]
+    assert_equal :image, img4.hdu_type
+    assert_kind_of RFits::Image, img4
+    assert_equal 4, img4.position
+    assert_equal [2, 2], img4.naxes
+    assert_equal 2, img4.naxis
+    assert_equal :long, img4.image_type
+    assert_equal RFits::IO::Proxy::LONG_IMG, img4.bitpix
+    assert_equal true, img4.compressed?
+    assert_equal :rice, img4.compression_type
+    assert_equal [2, 2], img4.tile_dim
+    assert_equal nil, img4.noise_bits
+  end
+  
+  def header
+    hdr = @fits[1].header
+    assert_equal 22, hdr.length
+    
+    # integer
+    hdr['TEST'] = [123, 'A test']
+    assert_equal [123, 'A test'], hdr['TEST', true]
+    
+    # floating point
+    hdr['TEST'] = 12.3
+    assert_equal 12.3, hdr['TEST']
+    
+    # string
+    hdr['TEST'] = 'Hello, world!'
+    assert_equal 'Hello, world!', hdr['TEST']
+    
+    # boolean
+    hdr['TEST'] = true
+    assert_equal true, hdr['TEST']
+    
+    # complex
+    hdr['TEST'] = Complex.new(12.3, 3.1)
+    assert_equal Complex.new(12.3, 3.1), hdr['TEST']
+    
+    assert_equal 23, hdr.length
+    
+    hdr.delete('TEST')
+    assert_equal 22, hdr.length
+    
+    hdr.comment = "NOAO test comment"
+    assert_equal "NOAO test comment", hdr.comment
+    
+    hdr.history = "NOAO test history"
+    assert_equal "NOAO test history", hdr.history
+  end
+  
+  def image_data
+    pix = @fits[1].data
+    int_array = [
+      1, 2, 3, 4,
+      4, 3, 2, 1,
+      5, 6, 7, 8,
+      9, 8, 7, 6
+    ]
+    pix[0] = int_array
+    assert_equal 16, pix.length
+    assert_equal int_array, pix[0..16]
+    assert_equal [5, 6, 7, 8], pix.row(2)
+    assert_equal [3, 2, 7, 7], pix.column(2)
+    assert_equal 1, pix.first
+    assert_equal 6, pix.last
+    
+    pix = @fits[3].data
+    float_array = [
+      1.1, 2.2, 3.3,
+      4.4, 5.5, 4.4,
+      3.3, 2.2, 1.1
+    ]
+    pix[0] = float_array
+    assert_equal 9, pix.length
+    float_array.each_with_index do |val, i|
+      assert_in_delta val, pix[i], 0.1
+    end
+       
+    pix = @fits[1].data
+    
+    pix[1] = 7
+    assert_equal 7, pix[1]
+    
+    pix[0..4] = [1, 2, 3, 4]
+    assert_equal [1, 2, 3, 4], pix[0..4]
+    
+    pix[[0, 0]] = [2, 4, 6, 8]
+    assert_equal [2, 4, 6, 8], pix[[0, 0], 4]
+    
+    pix[0, 0] = [3, 4, 5, 6]
+    assert_equal [3, 4, 5, 6], pix[[0, 0], 4]
+    
+    pix[[0, 0], 2] = [2, 2]
+    assert_equal [2, 2], pix[[0, 0], 2]
+    
+    pix[[0, 0], [1, 1]] = [-1, 0, 1, 2]
+    assert_equal [-1, 0, 1, 2], pix[[0, 0], [1, 1]]
+    
+    pix.set_row(1, [5, 5, 5, 5])
+    assert_equal pix.row(1), [5, 5, 5, 5]
+    
+    pix.set_column(2, [9, 9, 9, 9])
+    assert_equal pix.column(2), [9, 9, 9, 9]
+     
+    assert_raise(ArgumentError){ pix['blah'] = 2 }
+    assert_raise(ArgumentError){ pix['blah1', 'blah2'] = 3 }
+  end
+  
+  def table_data
+    tbl = @fits[2]
+    
+    assert_equal 4, tbl.column_information.size
+    
+    cinfo = tbl.column_information[0]
+    assert_equal 'column_1', cinfo.name
+    assert_equal :short, cinfo.data_type
+    assert_equal 1, cinfo.repeat
+    assert_equal 2, cinfo.width
+    assert_equal 6, cinfo.display_width
+    assert_equal 1, cinfo.dim
+    assert_equal [1], cinfo.size
+    
+    cinfo = tbl.column_information[1]
+    assert_equal 'column_2', cinfo.name
+    assert_equal :string, cinfo.data_type
+    assert_equal 10, cinfo.repeat
+    assert_equal 10, cinfo.width
+    assert_equal 10, cinfo.display_width
+    assert_equal 1, cinfo.dim
+    assert_equal [10], cinfo.size
+    
+    cinfo = tbl.column_information[2]
+    assert_equal 'column_3', cinfo.name
+    assert_equal :float, cinfo.data_type
+    assert_equal 1, cinfo.repeat
+    assert_equal 4, cinfo.width
+    assert_equal 14, cinfo.display_width
+    assert_equal 1, cinfo.dim
+    assert_equal [1], cinfo.size
+    
+    cinfo = tbl.column_information[3]
+    assert_equal 'column_4', cinfo.name
+    assert_equal :complex, cinfo.data_type
+    assert_equal 1, cinfo.repeat
+    assert_equal 8, cinfo.width
+    assert_equal 31, cinfo.display_width
+    assert_equal 1, cinfo.dim
+    assert_equal [1], cinfo.size
+    
+    tbl.column_information[2] = {:name => 'column_5', :format => 'A15'}
+    tbl.column_information << {:name => 'column_6', :format => 'I'}
+    assert_equal 6, tbl.column_information.size
+    
+    tbl.column_information.delete(5)
+    tbl.column_information.delete(2)
+    
+    assert_equal 4, tbl.column_information.size
+    
+    data = tbl.data
+    
+    int_array = [1, 2, 3, 4, 5]
+    data.set_column(0, int_array)
+    assert_equal int_array, data.column(0)
+    
+    string_array = ['a1', 'b2', 'c3', 'd4', 'e5']
+    data.set_column(1, string_array)
+    assert_equal string_array, data.column(1)
+    
+    float_array = [1.2, 2.3, 3.4, 4.5, 5.5]
+    data.set_column(2, float_array)
+    data.column(2).each_with_index do |val, i|
+      assert_in_delta float_array[i], val, 0.01
+    end
+    
+    complex_array = [Complex.new(1, 2), Complex.new(2, 3), Complex.new(3, 4), Complex.new(4, 5), Complex.new(5, 6)]
+    data.set_column(3, complex_array)
+    assert_equal complex_array, data.column(3)
+    
+    row = data[0]
+    assert_equal 1, row[0]
+    assert_equal 'a1', row[1]
+    assert_in_delta 1.2, row[2], 0.01
+    assert_equal Complex.new(1, 2), row[3]
+  
+    row = data[2]
+    assert_equal 3, row[0]
+    assert_equal 'c3', row[1]
+    assert_in_delta 3.4, row[2], 0.01
+    assert_equal Complex.new(3, 4), row[3]
+    
+    assert_equal 5, data.size
+    
+    data[5] = [9, 'new1', 2.2, Complex.new(1, 1)]
+    row = data[5]
+    assert_equal 9, row[0]
+    assert_equal 'new1', row[1]
+    assert_in_delta 2.2, row[2], 0.01
+    assert_equal Complex.new(1, 1), row[3]
+    
+    data[5, 1] = 'new3'
+    assert_equal 'new3', data[5,1]
+    
+    assert_equal 6, data.size
+    data.delete(5)
+    assert_equal 5, data.size
+    
+    first = data.first
+    assert_equal int_array[0], first[0]
+    assert_equal string_array[0], first[1]
+    assert_in_delta float_array[0], first[2], 0.01
+    assert_equal complex_array[0], first[3]
+    
+    last = data.last
+    assert_equal int_array[-1], last[0]
+    assert_equal string_array[-1], last[1]
+    assert_in_delta float_array[-1], last[2], 0.01
+    assert_equal complex_array[-1], last[3]
+    
+    data << [9, 'new1', 2.2, Complex.new(1, 1)]
+    last = data.last
+    assert_equal 9, last[0]
+    assert_equal 'new1', last[1]
+    assert_in_delta 2.2, last[2], 0.01
+    assert_equal Complex.new(1, 1), last[3]
+    
+    data.delete(5)
+    
+    assert_equal 1, data[0, 0]
+    assert_equal 'c3', data[2, 1]
+  end
+  
+  def teardown
+    @fits.close
+    FileUtils.rm('test/readwrite.test.fits')
+  end
+end

data/test/rfitsio_test.rb

@@ -0,0 +1,315 @@
+require 'test/unit'
+require 'fileutils'
+
+require 'rfits/ext/rfitsio'
+include RFits::IO
+
+class RFitsCompressionTest < Test::Unit::TestCase
+  def setup
+    @fits = Proxy.fits_create_file('!test/test.1.working.fits')
+  end
+  
+  def test_compression
+    Proxy.fits_set_compression_type(@fits, Proxy::RICE_1)
+    Proxy.fits_insert_img(@fits, Proxy::FLOAT_IMG, 2, [5, 5])
+    Proxy.fits_write_subset(@fits, Proxy::TFLOAT, [1, 1], [5, 5], [
+      1.0, 2.1, 3.2, 4.3, 5.4,
+      6.5, 7.6, 8.7, 9.8, 10.9,
+      11.0, 12.1, 13.2, 14.3, 15.4,
+      16.5, 17.6, 18.7, 19.8, 20.9,
+      21.0, 22.1, 23.2, 24.3, 25.4
+    ])
+    assert_equal 1, Proxy.fits_is_compressed_image(@fits)
+    
+    Proxy.fits_set_compression_type(@fits, 0)
+    Proxy.fits_insert_img(@fits, Proxy::LONG_IMG, 2, [4, 4])
+    Proxy.fits_write_subset(@fits, Proxy::TLONG, [1, 1], [4, 4], [
+     1, 2, 3, 4,
+     4, 3, 2, 1,
+     2, 2, 2, 2,
+     3, 3, 3, 3
+    ])
+    assert_equal 0, Proxy.fits_is_compressed_image(@fits)
+    
+    Proxy.fits_set_compression_type(@fits, Proxy::GZIP_1)
+    Proxy.fits_insert_img(@fits, Proxy::LONG_IMG, 2, [3, 3])
+    Proxy.fits_write_subset(@fits, Proxy::TLONG, [1, 1], [3, 3], [
+      1, 5, 9,
+      2, 6, 8,
+      3, 7, 7
+    ])
+    assert_equal 1, Proxy.fits_is_compressed_image(@fits)
+    
+    assert_equal 4, Proxy.fits_get_num_hdus(@fits)
+  end
+  
+  def teardown
+    Proxy.fits_close_file(@fits)
+    FileUtils.rm('test/test.1.working.fits')
+  end
+end
+
+class RFitsIOWriteTest < Test::Unit::TestCase
+  def setup
+    @fits = Proxy.fits_create_file('!test/test.3.working.fits')
+  end
+  
+  def test_rfitsio
+    image_writing()
+    table_writing()
+    keyword_writing()
+  end
+  
+  def image_writing
+    ipixels = [
+      1, 2, 3, 4, 5,
+      2, 3, 4, 5, 1,
+      3, 4, 5, 1, 2,
+      4, 5, 1, 2, 3,
+      5, 1, 2, 3, 4
+    ]
+    
+    fpixels = [
+      1.0, 2.1, 3.2, 4.3, 5.4,
+      6.5, 7.6, 8.7, 9.8, 10.9,
+      11.0, 12.1, 13.2, 14.3, 15.4,
+      16.5, 17.6, 18.7, 19.8, 20.9,
+      21.0, 22.1, 23.2, 24.3, 25.4
+    ]
+    
+    Proxy.fits_create_img(@fits, Proxy::LONG_IMG, 2, [5, 5])
+    Proxy.fits_write_subset(@fits, Proxy::TLONG, [1, 1], [5, 5], ipixels)
+    assert_equal [ipixels, 0], Proxy.fits_read_subset(@fits, Proxy::TLONG, [1, 1], [5, 5], [1, 1], 0)
+    
+    Proxy.fits_create_img(@fits, Proxy::SHORT_IMG, 2, [5, 5])
+    Proxy.fits_write_img(@fits, Proxy::TSHORT, 1, 25, ipixels)
+    assert_equal [ipixels, 0], Proxy.fits_read_img(@fits, Proxy::TLONG, 1, 25, nil)
+    Proxy.fits_delete_hdu(@fits)
+    
+    Proxy.fits_set_compression_type(@fits, Proxy::HCOMPRESS_1)
+    assert_equal Proxy::HCOMPRESS_1, Proxy.fits_get_compression_type(@fits)
+    Proxy.fits_set_tile_dim(@fits, 2, [1, 1])
+    assert_equal [1, 1], Proxy.fits_get_tile_dim(@fits, 2)
+    Proxy.fits_set_noise_bits(@fits, 5)
+    assert_equal 5, Proxy.fits_get_noise_bits(@fits)
+    Proxy.fits_set_hcomp_scale(@fits, 2)
+    assert_equal 2, Proxy.fits_get_hcomp_scale(@fits)
+    Proxy.fits_set_hcomp_smooth(@fits, 2)
+    assert_equal 2, Proxy.fits_get_hcomp_smooth(@fits)
+    
+    Proxy.fits_create_img(@fits, Proxy::FLOAT_IMG, 2, [5, 5])
+    Proxy.fits_write_pixnull(@fits, Proxy::TFLOAT, [1, 1], 25, fpixels, 0.0)
+    Proxy.fits_read_pix(@fits, Proxy::TFLOAT, [1, 1], 25, nil).first.each_with_index do |p, i|
+      assert_in_delta fpixels[i], p, 0.0001
+    end
+    Proxy.fits_delete_hdu(@fits)
+  end
+  
+  def table_writing
+    Proxy.fits_create_tbl(@fits, Proxy::BINARY_TBL, 0, 5,
+      ['column_a', 'column_b', 'column_c', 'column_d', 'column_e'],
+      ['I', '10A', 'E', 'C', 'D'],
+      nil,
+      "Binary table test")
+      
+    assert_equal 5, Proxy.fits_get_num_cols(@fits)
+    assert_equal 3, Proxy.fits_get_colnum(@fits, Proxy::CASEINSEN, 'column_c')
+    assert_equal ['column_d', 4], Proxy.fits_get_colname(@fits, Proxy::CASEINSEN, 'column_d')
+    
+    assert_equal [Proxy::TSHORT, 1, 2], Proxy.fits_get_coltype(@fits, 1)
+    assert_equal [Proxy::TSTRING, 10, 10], Proxy.fits_get_coltype(@fits, 2)
+    assert_equal [Proxy::TFLOAT, 1, 4], Proxy.fits_get_coltype(@fits, 3)
+    assert_equal [Proxy::TCOMPLEX, 1, 8], Proxy.fits_get_coltype(@fits, 4)
+    assert_equal [Proxy::TDOUBLE, 1, 8], Proxy.fits_get_coltype(@fits, 5)
+    
+    assert_equal [Proxy::TSHORT, 1, 2], Proxy.fits_get_eqcoltype(@fits, 1)
+    assert_equal [Proxy::TSTRING, 10, 10], Proxy.fits_get_eqcoltype(@fits, 2)
+    assert_equal [Proxy::TFLOAT, 1, 4], Proxy.fits_get_eqcoltype(@fits, 3)
+    assert_equal [Proxy::TCOMPLEX, 1, 8], Proxy.fits_get_eqcoltype(@fits, 4)
+    assert_equal [Proxy::TDOUBLE, 1, 8], Proxy.fits_get_eqcoltype(@fits, 5)
+    
+    assert_equal 6, Proxy.fits_get_col_display_width(@fits, 1)
+    assert_equal 10, Proxy.fits_get_col_display_width(@fits, 2)
+    assert_equal 14, Proxy.fits_get_col_display_width(@fits, 3)
+    assert_equal 31, Proxy.fits_get_col_display_width(@fits, 4)
+    assert_equal 23, Proxy.fits_get_col_display_width(@fits, 5)
+    
+    assert_equal [1, [10]], Proxy.fits_read_tdim(@fits, 2, 1)
+    
+    Proxy.fits_insert_rows(@fits, 0, 5)
+    assert_equal 5, Proxy.fits_get_num_rows(@fits)
+    Proxy.fits_delete_rows(@fits, 2, 2)
+    assert_equal 3, Proxy.fits_get_num_rows(@fits)
+    Proxy.fits_insert_rows(@fits, 2, 2);
+    Proxy.fits_delete_rowrange(@fits, '2-4')
+    assert_equal 2, Proxy.fits_get_num_rows(@fits)
+    Proxy.fits_insert_rows(@fits, 2, 3);
+    Proxy.fits_delete_rowlist(@fits, [1, 3, 5], 3)
+    assert_equal 2, Proxy.fits_get_num_rows(@fits)
+    Proxy.fits_delete_rows(@fits, 1, 2)
+    
+    Proxy.fits_insert_cols(@fits, 6, 2, ['column_f', 'column_g'], ['I', 'I'])
+    assert_equal 7, Proxy.fits_get_num_cols(@fits)
+    Proxy.fits_delete_col(@fits, 7)
+    Proxy.fits_delete_col(@fits, 6)
+    assert_equal 5, Proxy.fits_get_num_cols(@fits)
+    
+    Proxy.fits_modify_vector_len(@fits, 2, 15)
+    assert_equal [Proxy::TSTRING, 15, 15], Proxy.fits_get_eqcoltype(@fits, 2)
+    
+    short_list = [5, 4, 3, 2, 1]
+    Proxy.fits_write_colnull(@fits, Proxy::TSHORT, 1, 1, 1, 5, short_list, nil)
+    assert_equal [short_list, 0], Proxy.fits_read_col(@fits, Proxy::TSHORT, 1, 1, 1, 5, nil)
+    
+    string_list = ['a1', 'b2', 'c3', 'd4', 'e5']
+    Proxy.fits_write_colnull(@fits, Proxy::TSTRING, 2, 1, 1, 5, string_list, nil)
+    assert_equal [string_list, 0], Proxy.fits_read_col(@fits, Proxy::TSTRING, 2, 1, 1, 5, nil)
+    
+    float_list = [1.1, 2.2, 3.3, 4.4, 5.5]
+    Proxy.fits_write_colnull(@fits, Proxy::TFLOAT, 3, 1, 1, 5, float_list, nil)
+    Proxy.fits_read_col(@fits, Proxy::TFLOAT, 3, 1, 1, 5, nil).first.each_with_index do |f, i|
+      assert_in_delta float_list[i], f, 0.001
+    end
+    
+    complex_list = [Complex.new(1.1, 2.2), Complex.new(2.2, 3.3), Complex.new(3.3, 4.4), Complex.new(4.4, 5.5), Complex.new(5.5, 6.5)]
+    Proxy.fits_write_colnull(@fits, Proxy::TCOMPLEX, 4, 1, 1, 5, complex_list, nil)
+    Proxy.fits_read_col(@fits, Proxy::TCOMPLEX, 4, 1, 1, 5, nil).first.each_with_index do |rc, i|
+      assert_in_delta complex_list[i].real, rc.real, 0.001
+      assert_in_delta complex_list[i].image, rc.image, 0.001
+    end
+    
+    double_list = [10.5, 9.4, 8.3, 7.2, 6.1]
+    Proxy.fits_write_colnull(@fits, Proxy::TDOUBLE, 5, 1, 1, 5, double_list, nil)
+    Proxy.fits_read_col(@fits, Proxy::TDOUBLE, 5, 1, 1, 5, nil).first.each_with_index do |d, i|
+      assert_in_delta double_list[i], d, 0.001
+    end
+  end
+
+  def keyword_writing
+    Proxy.fits_movabs_hdu(@fits, 1)
+    
+    Proxy.fits_update_key(@fits, Proxy::TSTRING, 'TEST1', 'Hello!', 'A testing string')
+    assert_equal ["'Hello!  '", 'A testing string'], Proxy.fits_read_keyword(@fits, 'TEST1')
+
+    Proxy.fits_update_key(@fits, Proxy::TLOGICAL, 'TEST2', true, 'A testing logical');
+    assert_equal ['T', 'A testing logical'], Proxy.fits_read_keyword(@fits, 'TEST2')
+    
+    Proxy.fits_update_key(@fits, Proxy::TLOGICAL, 'TEST3', false, 'A testing logical');
+    assert_equal ['F', 'A testing logical'], Proxy.fits_read_keyword(@fits, 'TEST3')
+
+    Proxy.fits_update_key(@fits, Proxy::TINT, 'TEST4', 999999, 'A testing integer')
+    assert_equal ['999999', 'A testing integer'], Proxy.fits_read_keyword(@fits, 'TEST4')
+
+    Proxy.fits_update_key(@fits, Proxy::TDOUBLE, 'TEST5', 345.78, 'A testing double')
+    assert_equal ['345.78', 'A testing double'], Proxy.fits_read_keyword(@fits, 'TEST5')
+
+    Proxy.fits_update_key(@fits, Proxy::TCOMPLEX, 'TEST6', Complex.new(4.1, 5.2), 'A testing complex')
+    assert_equal ['(4.1, 5.2)', 'A testing complex'], Proxy.fits_read_keyword(@fits, 'TEST6')
+
+    Proxy.fits_update_key_null(@fits, 'TEST7', 'A testing null')
+    assert_equal ['', 'A testing null'], Proxy.fits_read_keyword(@fits, 'TEST7')
+  end
+  
+  def teardown
+    Proxy.fits_close_file(@fits)
+    FileUtils.rm('test/test.3.working.fits')
+  end
+end
+
+class RFitsIOReadTest < Test::Unit::TestCase
+  def setup
+    FileUtils.cp('test/test.2.fits', 'test/test.2.working.fits')
+    @fits = Proxy.fits_open_file('test/test.2.working.fits', Proxy::READWRITE)
+  end
+  
+  def test_rfitsio
+    file_access()
+    hdu_access()
+    keyword_reading()
+    image_routines()
+    compression()
+  end
+  
+  def file_access
+    assert_equal 'test/test.2.working.fits', Proxy.fits_file_name(@fits)
+    assert_equal Proxy::READWRITE, Proxy.fits_file_mode(@fits)
+    assert_equal 'file://', Proxy.fits_url_type(@fits)
+  end
+  
+  def hdu_access
+    assert_equal 9, Proxy.fits_get_num_hdus(@fits)
+    assert_equal Proxy::IMAGE_HDU, Proxy.fits_movabs_hdu(@fits, 3)
+    assert_equal 3, Proxy.fits_get_hdu_num(@fits)
+    assert_equal Proxy::IMAGE_HDU, Proxy.fits_get_hdu_type(@fits)
+  end
+  
+  def keyword_reading
+    Proxy.fits_movabs_hdu(@fits, 1)
+    
+    assert_equal 62, Proxy.fits_get_hdrspace(@fits).first
+    
+    ra_obj = Proxy.fits_read_keyword(@fits, 'RA_OBJ')
+    assert_equal ['182.635454000001', 'R.A. of the object (degrees)'], ra_obj
+    assert_equal [182.635454000001, 'R.A. of the object (degrees)'], Proxy.fits_read_key(@fits, Proxy::TDOUBLE, 'RA_OBJ')
+    assert_equal 'F', Proxy.fits_get_keytype(ra_obj.first)
+    assert_equal ['RA_OBJ', '182.635454000001', 'R.A. of the object (degrees)'], Proxy.fits_read_keyn(@fits, 11)
+    
+    object = Proxy.fits_read_keyword(@fits, 'OBJECT')
+    assert_equal ["'NGC 4151'", 'Name of observed object'], object
+    assert_equal ["NGC 4151", 'Name of observed object'], Proxy.fits_read_key(@fits, Proxy::TSTRING, 'OBJECT')
+    assert_equal 'C', Proxy.fits_get_keytype(object.first)
+    assert_equal ['OBJECT', "'NGC 4151'", 'Name of observed object'], Proxy.fits_read_keyn(@fits, 10)
+    
+    off_axis = Proxy.fits_read_keyword(@fits, 'OFF-AXIS')
+    assert_equal ['T', 'Was this pointing done off-axis'], off_axis
+    assert_equal [true, 'Was this pointing done off-axis'], Proxy.fits_read_key(@fits, Proxy::TLOGICAL, 'OFF-AXIS')
+    assert_equal 'L', Proxy.fits_get_keytype(off_axis.first)
+    assert_equal ['OFF-AXIS', 'T', 'Was this pointing done off-axis'], Proxy.fits_read_keyn(@fits, 25)
+    
+    bitpix = Proxy.fits_read_keyword(@fits, 'BITPIX')
+    assert_equal ['8', 'Character information'], bitpix
+    assert_equal [8, 'Character information'], Proxy.fits_read_key(@fits, Proxy::TINT, 'BITPIX')
+    assert_equal 'I', Proxy.fits_get_keytype(bitpix.first)
+    assert_equal ['BITPIX', '8', 'Character information'], Proxy.fits_read_keyn(@fits, 2)
+  end
+  
+  def image_routines
+    Proxy.fits_movabs_hdu(@fits, 3)
+    
+    assert_equal Proxy::SHORT_IMG, Proxy.fits_get_img_type(@fits)
+    assert_equal Proxy::SHORT_IMG, Proxy.fits_get_img_equivtype(@fits)
+    assert_equal 2, Proxy.fits_get_img_dim(@fits)
+    assert_equal [2048, 300], Proxy.fits_get_img_size(@fits, 2)
+    assert_equal [Proxy::SHORT_IMG, 2, [2048, 300]], Proxy.fits_get_img_param(@fits, 2)
+    
+    assert_equal [[17.0, 15.0, 7.0, 15.0, 5.0, 4.0, 5.0, 5.0, 11.0, 4.0], 0], Proxy.fits_read_img(@fits, Proxy::TDOUBLE, 112, 10, 0)
+    assert_equal [[17, 15, 7, 15, 5, 4, 5, 5, 11, 4], 0], Proxy.fits_read_img(@fits, Proxy::TLONG, 112, 10, 0)
+    
+    assert_equal [[1.0, 1.0, 3.0, 0.0, 0.0, 2.0, 0.0, 1.0, 0.0], 0],
+      Proxy.fits_read_subset(@fits, Proxy::TDOUBLE, [902, 122], [910, 122], [1, 1], 0)
+    assert_equal [[1, 1, 3, 0, 0, 2, 0, 1, 0], 0],
+      Proxy.fits_read_subset(@fits, Proxy::TLONG, [902, 122], [910, 122], [1, 1], 0)
+      
+    assert_equal [[1.0, 1.0, 3.0, 0.0, 0.0, 2.0, 0.0, 1.0, 0.0], 0],
+      Proxy.fits_read_pix(@fits, Proxy::TDOUBLE, [902, 122], 9, 0)
+    assert_equal [[1, 1, 3, 0, 0, 2, 0, 1, 0], 0],
+      Proxy.fits_read_pix(@fits, Proxy::TLONG, [902, 122], 9, 0)
+  end
+  
+  def compression
+    Proxy.fits_movabs_hdu(@fits, 3)
+    
+    # Image is not compressed, so this is sort of boring.
+    assert_equal 0, Proxy.fits_is_compressed_image(@fits)
+    assert_equal 0, Proxy.fits_get_compression_type(@fits)
+    assert_equal [0, 0], Proxy.fits_get_tile_dim(@fits, 2)
+    assert_equal 0, Proxy.fits_get_noise_bits(@fits)
+    assert_equal 0, Proxy.fits_get_hcomp_scale(@fits)
+    assert_equal 0, Proxy.fits_get_hcomp_smooth(@fits)
+  end
+  
+  def teardown
+    Proxy.fits_close_file(@fits)
+    FileUtils.rm('test/test.2.working.fits')
+  end
+end