public void TestReadBuffered() { string filename = TestFileSetup.GetTargetFilename("test.fits"); BufferedFile bf = new BufferedFile( filename, FileAccess.Read, FileShare.None); Header h = Header.ReadHeader(bf); long n = h.DataSize; int naxes = h.GetIntValue("NAXIS"); int lastAxis = h.GetIntValue("NAXIS" + naxes); HeaderCard hnew = new HeaderCard("NAXIS", naxes - 1, "this is header card with naxes"); h.AddCard(hnew); float[] line = new float[h.DataSize]; for (int i = 0; i < lastAxis; i += 1) { Console.Out.WriteLine("read"); bf.Read(line); } bf.Close(); bf.Dispose(); }
// Make FITS file public void MakeFile() { Initialize(); Fits f = new Fits(); // Make HDUs of various types. f.AddHDU(Fits.MakeHDU(bimg)); f.AddHDU(Fits.MakeHDU(simg)); f.AddHDU(Fits.MakeHDU(iimg)); f.AddHDU(Fits.MakeHDU(limg)); f.AddHDU(Fits.MakeHDU(fimg)); f.AddHDU(Fits.MakeHDU(dimg)); f.AddHDU(Fits.MakeHDU(img3)); f.AddHDU(Fits.MakeHDU(img1)); Assert.AreEqual(f.NumberOfHDUs, 8); // Write a FITS file. BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("image1.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Flush(); bf.Close(); bf.Dispose(); f.Close(); }
void WriteFile(Fits f, String name) { BufferedFile bf = new BufferedFile(name, FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Flush(); bf.Close(); bf.Dispose(); }
public void TestRowDelete() { Fits f = null; try { f = new Fits( TestFileSetup.GetTargetFilename( TestFileSetup.GetTargetFilename("bt1.fits"))); f.Read(); BinaryTableHDU thdu = (BinaryTableHDU)f.GetHDU(1); Assert.AreEqual(50, thdu.NRows); thdu.DeleteRows(10, 20); Assert.AreEqual(30, thdu.NRows); double[] dbl = (double[])thdu.GetColumn(6); Assert.AreEqual(dbl[9], doubles[9]); Assert.AreEqual(dbl[10], doubles[30]); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("bt1x.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt1x.fits")); f.Read(); thdu = (BinaryTableHDU)f.GetHDU(1); dbl = (double[])thdu.GetColumn(6); Assert.AreEqual(30, thdu.NRows); Assert.AreEqual(9, thdu.NCols); Assert.AreEqual(dbl[9], doubles[9]); Assert.AreEqual(dbl[10], doubles[30]); thdu.DeleteRows(20); Assert.AreEqual(20, thdu.NRows); dbl = (double[])thdu.GetColumn(6); Assert.AreEqual(20, dbl.Length); Assert.AreEqual(dbl[0], doubles[0]); Assert.AreEqual(dbl[19], doubles[39]); } finally { if (f != null) { f.Close(); } } }
public void BuildByColumn() { Fits f = null; try { BinaryTable btab = new BinaryTable(); btab.AddColumn(floats); btab.AddColumn(vf); btab.AddColumn(strings); btab.AddColumn(vbool); btab.AddColumn(ints); f = new Fits(); f.AddHDU(Fits.MakeHDU(btab)); BufferedFile bdos = new BufferedFile( TestFileSetup.GetTargetFilename("bt3.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bdos); bdos.Close(); bdos.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt3.fits")); BinaryTableHDU bhdu = (BinaryTableHDU)f.GetHDU(1); btab = (BinaryTable)bhdu.Data; Assert.AreEqual(true, ArrayFuncs.ArrayEquals(floats, bhdu.GetColumn(0))); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(vf, bhdu.GetColumn(1))); // problem is here only String[] col = (String[])bhdu.GetColumn(2); for (int i = 0; i < col.Length; i += 1) { col[i] = col[i].Trim(); } Assert.AreEqual(true, ArrayFuncs.ArrayEquals(strings, col)); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(vbool, bhdu.GetColumn(3))); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(ints, bhdu.GetColumn(4))); } finally { if (f != null) { f.Close(); } } }
public void Delete() { Fits f = null; try { f = new Fits(TestFileSetup.GetTargetFilename("at1.fits"), FileAccess.ReadWrite); TableHDU th = (TableHDU)f.GetHDU(1); Assert.AreEqual(50, th.NRows); th.DeleteRows(2, 2); Assert.AreEqual(48, th.NRows); BufferedFile bf = new BufferedFile(TestFileSetup.GetTargetFilename("at1y.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("at1y.fits")); th = (TableHDU)f.GetHDU(1); Assert.AreEqual(48, th.NRows); Assert.AreEqual(5, th.NCols); th.DeleteColumnsIndexZero(3, 2); Assert.AreEqual(3, th.NCols); th.DeleteColumnsIndexZero(0, 2); Assert.AreEqual(1, th.NCols); bf = new BufferedFile(TestFileSetup.GetTargetFilename("at1z.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("at1z.fits")); th = (TableHDU)f.GetHDU(1); Assert.AreEqual(1, th.NCols); } finally { if (f != null) { f.Close(); } } }
public void TestDegenerate() { Fits f = null; try { String[] sa = new String[10]; int[,] ia = new int[10, 0]; f = new Fits(); for (int i = 0; i < sa.Length; i += 1) { sa[i] = ""; } Object[] data = new Object[] { sa, ia }; BinaryTableHDU bhdu = (BinaryTableHDU)Fits.MakeHDU(data); Header hdr = bhdu.Header; f.AddHDU(bhdu); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("bt7.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); Assert.AreEqual(2, hdr.GetIntValue("TFIELDS")); Assert.AreEqual(10, hdr.GetIntValue("NAXIS2")); Assert.AreEqual(0, hdr.GetIntValue("NAXIS1")); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt7.fits")); bhdu = (BinaryTableHDU)f.GetHDU(1); hdr = bhdu.Header; Assert.AreEqual(2, hdr.GetIntValue("TFIELDS")); Assert.AreEqual(10, hdr.GetIntValue("NAXIS2")); Assert.AreEqual(0, hdr.GetIntValue("NAXIS1")); } finally { if (f != null) { f.Close(); } } }
public void test() { float[][] data = new float[300][]; for (int i = 0; i < 300; i++) { data[i] = new float[300]; } for (int i = 0; i < 300; i += 1) { for (int j = 0; j < 300; j += 1) { data[i][j] = 1000 * i + j; } } Fits f = null; try { f = new Fits(); BufferedFile bf = new BufferedFile(TestFileSetup.GetTargetFilename("tiler1.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.AddHDU(Fits.MakeHDU(data)); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("tiler1.fits")); ImageHDU h = (ImageHDU)f.ReadHDU(); ImageTiler t = h.Tiler; doTile("t1", data, t, 200, 200, 50, 50); doTile("t2", data, t, 133, 133, 72, 26); Object o = h.Data.Kernel; doTile("t3", data, t, 200, 200, 50, 50); doTile("t4", data, t, 133, 133, 72, 26); } finally { if (f != null) { f.Close(); } } }
public static ThreeAxisFITSCubeFrameStream OpenFile(string fileName, VideoController videoController) { var fitsFile = new Fits(); var bufferedFile = new BufferedFile(fileName, FileAccess.Read, FileShare.ReadWrite); fitsFile.Read(bufferedFile); var imageHDU = fitsFile.GetHDU(0); ThreeAxisFITSCubeFrameStream fitsStream = null; videoController.SetCursor(Cursors.WaitCursor); try { var hasher = new SHA1CryptoServiceProvider(); hasher.Initialize(); byte[] combinedFileNamesBytes = Encoding.UTF8.GetBytes(fileName); var hash = hasher.ComputeHash(combinedFileNamesBytes, 0, combinedFileNamesBytes.Length); var fitsFileHash = Convert.ToBase64String(hash); var frm = new frmDefineFitsCube3D(imageHDU, fitsFileHash, videoController); if (videoController.ShowDialog(frm) == DialogResult.OK) { videoController.SetFlipSettings(frm.FlipVertically, frm.FlipHorizontally); fitsStream = new ThreeAxisFITSCubeFrameStream( fileName, fitsFile, bufferedFile, imageHDU, frm.TimeStampReader, frm.WidthIndex, frm.HeightIndex, frm.FrameIndex, frm.MinPixelValue, frm.MaxPixelValue, frm.BitPix, frm.NegPixCorrection); } } finally { if (fitsStream == null) { bufferedFile.Dispose(); } } return(fitsStream); }
public void TestVar() { Fits f = null; try { Object[] data = new Object[] { floats, vf, vs, vd, shorts, vbool }; f = new Fits(); f.AddHDU(Fits.MakeHDU(data)); BufferedFile bdos = new BufferedFile( TestFileSetup.GetTargetFilename("bt2.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bdos); bdos.Close(); bdos.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt2.fits"), FileAccess.Read); f.Read(); BinaryTableHDU bhdu = (BinaryTableHDU)f.GetHDU(1); Header hdr = bhdu.Header; Assert.AreEqual(true, hdr.GetIntValue("PCOUNT") > 0); Assert.AreEqual(6, hdr.GetIntValue("TFIELDS")); for (int i = 0; i < data.Length; i += 1) { Assert.AreEqual(true, ArrayFuncs.ArrayEquals(data[i], bhdu.GetColumn(i))); } } finally { f.Close(); } }
public void TestFitsCopy() { String file = TestFileSetup.GetTargetFilename("test_dup.fits"); Fits f = new Fits(file); int i = 0; BasicHDU h; do { h = f.ReadHDU(); if (h != null) { if (i == 0) { Console.Out.WriteLine("\n\nPrimary header:\n"); } else { Console.Out.WriteLine("\n\nExtension " + i + ":\n"); } i += 1; h.Info(); } } while (h != null); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("gbfits3.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); }
public void TestSimpleImages() { float[][] img = new float[300][]; for (int i = 0; i < 300; i++) { img[i] = new float[300]; } Fits f = null; try { f = new Fits(); ImageHDU hdu = (ImageHDU)Fits.MakeHDU(img); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("ht1.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.AddHDU(hdu); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("ht1.fits")); hdu = (ImageHDU)f.GetHDU(0); Header hdr = hdu.Header; Assert.AreEqual(2, hdr.GetIntValue("NAXIS")); Assert.AreEqual(300, hdr.GetIntValue("NAXIS1")); Assert.AreEqual(300, hdr.GetIntValue("NAXIS2")); Assert.AreEqual(300, hdr.GetIntValue("NAXIS2", -1)); Assert.AreEqual(-1, hdr.GetIntValue("NAXIS3", -1)); Assert.AreEqual(-32, hdr.GetIntValue("BITPIX")); Cursor c = hdr.GetCursor(); c.MoveNext(); HeaderCard hc = (HeaderCard)((DictionaryEntry)c.Current).Value; Assert.AreEqual("SIMPLE", hc.Key); c.MoveNext(); hc = (HeaderCard)((DictionaryEntry)c.Current).Value; Assert.AreEqual("BITPIX", hc.Key); c.MoveNext(); hc = (HeaderCard)((DictionaryEntry)c.Current).Value; Assert.AreEqual("NAXIS", hc.Key); c.MoveNext(); hc = (HeaderCard)((DictionaryEntry)c.Current).Value; Assert.AreEqual("NAXIS1", hc.Key); c.MoveNext(); hc = (HeaderCard)((DictionaryEntry)c.Current).Value; Assert.AreEqual("NAXIS2", hc.Key); } finally { if (f != null) { f.Close(); } } }
public void TestSimpleIO() { Fits f = null; try { FitsFactory.UseAsciiTables = false; f = new Fits(); Object[] data = new Object[] { bytes, bits, bools, shorts, ints, floats, doubles, longs, strings }; f.AddHDU(Fits.MakeHDU(data)); BinaryTableHDU bhdu = (BinaryTableHDU)f.GetHDU(1); bhdu.SetColumnName(0, "bytes", null); bhdu.SetColumnName(1, "bits", "bits later on"); bhdu.SetColumnName(6, "doubles", null); bhdu.SetColumnName(5, "floats", "4 x 4 array"); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("bt1.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Flush(); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt1.fits")); f.Read(); Assert.AreEqual(2, f.NumberOfHDUs); BinaryTableHDU thdu = (BinaryTableHDU)f.GetHDU(1); Header hdr = thdu.Header; Assert.AreEqual(9, hdr.GetIntValue("TFIELDS")); Assert.AreEqual(2, hdr.GetIntValue("NAXIS")); Assert.AreEqual(8, hdr.GetIntValue("BITPIX")); Assert.AreEqual("BINTABLE", hdr.GetStringValue("XTENSION")); Assert.AreEqual("bytes", hdr.GetStringValue("TTYPE1")); Assert.AreEqual("doubles", hdr.GetStringValue("TTYPE7")); for (int i = 0; i < data.Length; i += 1) { Object col = thdu.GetColumn(i); if (i == 8) { String[] st = (String[])col; for (int j = 0; j < st.Length; j += 1) { st[j] = st[j].Trim(); } } Assert.AreEqual(true, ArrayFuncs.ArrayEquals(data[i], col)); } } finally { if (f != null) { f.Close(); } } }
public void TestObj() { FitsFactory.UseAsciiTables = false; /*** Create a binary table from an Object[][] array */ Object[][] x = new Object[5][]; for (int i = 0; i < 5; i += 1) { x[i] = new Object[3]; x[i][0] = new float[] { i }; string temp = string.Concat("AString", i); x[i][1] = new string[] { temp }; int[][] t = new int[2][]; for (int j = 0; j < 2; j++) { t[j] = new int[2]; t[j][0] = j * i; t[j][1] = (j + 2) * i; } x[i][2] = t; } Fits f = null; try { f = new Fits(); BasicHDU hdu = Fits.MakeHDU(x); f.AddHDU(hdu); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("bt5.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); /* Now get rid of some columns */ BinaryTableHDU xhdu = (BinaryTableHDU)hdu; // First column Assert.AreEqual(3, xhdu.NCols); xhdu.DeleteColumnsIndexOne(1, 1); Assert.AreEqual(2, xhdu.NCols); xhdu.DeleteColumnsIndexZero(1, 1); Assert.AreEqual(1, xhdu.NCols); bf = new BufferedFile( TestFileSetup.GetTargetFilename("bt6.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt6.fits")); xhdu = (BinaryTableHDU)f.GetHDU(1); Assert.AreEqual(1, xhdu.NCols); } finally { if (f != null) { f.Close(); } } }
public void BuildByRow() { Fits f = null; try { f = new Fits( TestFileSetup.GetTargetFilename("bt2.fits"), FileAccess.Read); f.Read(); BinaryTableHDU bhdu = (BinaryTableHDU)f.GetHDU(1); Header hdr = bhdu.Header; BinaryTable btab = (BinaryTable)bhdu.Data; for (int i = 0; i < 50; i += 1) { Object[] row = (Object[])btab.GetRow(i); float[] qx = (float[])row[1]; Array[] p = (Array[])row[0]; float[] pt = (float[])p.GetValue(0); pt[0] = (float)(i * Math.Sin(i)); btab.AddRow(row); } f.Close(); f = new Fits(); f.AddHDU(Fits.MakeHDU(btab)); BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("bt4.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Flush(); bf.Close(); bf.Dispose(); f.Close(); f = new Fits( TestFileSetup.GetTargetFilename("bt4.fits"), FileAccess.Read); btab = (BinaryTable)f.GetHDU(1).Data; Assert.AreEqual(100, btab.NRows); // Try getting data before we Read in the table. Array[] xf = (Array[])btab.GetColumn(0); Array[] xft = (Array[])xf.GetValue(50); float[] xftt = (float[])xft.GetValue(0); Assert.AreEqual((float)0, (float)xftt[0]); xft = (Array[])xf.GetValue(99); xftt = (float[])xft.GetValue(0); Assert.AreEqual((float)(49 * Math.Sin(49)), (float)xftt[0]); for (int i = 0; i < xf.Length; i += 3) { bool[] ba = (bool[])btab.GetElement(i, 5); float[] fx = (float[])btab.GetElement(i, 1); int trow = i % 50; Assert.AreEqual(true, ArrayFuncs.ArrayEquals(ba, vbool[trow])); // prob 1 Assert.AreEqual(true, ArrayFuncs.ArrayEquals(fx, vf[trow])); } // Fill the table. Data data = f.GetHDU(1).Data; xf = (Array[])btab.GetColumn(0); xft = (Array[])xf.GetValue(50); xftt = (float[])xft.GetValue(0); Assert.AreEqual(0F, (float)xftt[0]); xft = (Array[])xf.GetValue(99); xftt = (float[])xft.GetValue(0); Assert.AreEqual((float)(49 * Math.Sin(49)), (float)xftt[0]); for (int i = 0; i < xf.Length; i += 3) { bool[] ba = (bool[])btab.GetElement(i, 5); float[] fx = (float[])btab.GetElement(i, 1); int trow = i % 50; Assert.AreEqual(true, ArrayFuncs.ArrayEquals(ba, vbool[trow])); // prob 2 Assert.AreEqual(true, ArrayFuncs.ArrayEquals(fx, vf[trow])); } } finally { if (f != null) { f.Close(); } } }
public void ModifyTable() { Fits f = null; try { f = new Fits(TestFileSetup.GetTargetFilename("at1.fits"), FileAccess.ReadWrite); Object[] samp = GetSampleCols(); AsciiTableHDU hdu = (AsciiTableHDU)f.GetHDU(1); AsciiTable data = (AsciiTable)hdu.GetData(); float[] f1 = (float[])data.GetColumn(0); float[] f2 = (float[])f1.Clone(); for (int i = 0; i < f2.Length; i += 1) { f2[i] = 2 * f2[i]; } data.SetColumn(0, f2); f1 = new float[] { 3.14159f }; data.SetElement(3, 0, f1); hdu.SetNullString(0, "**INVALID**"); data.SetNull(5, 0, true); data.SetNull(6, 0, true); Object[] row = new Object[5]; row[0] = new float[] { 6.28f }; row[1] = new int[] { 22 }; row[2] = new long[] { 0 }; row[3] = new double[] { -3 }; row[4] = new String[] { "A string" }; data.SetRow(5, row); data.SetElement(4, 2, new long[] { 54321 }); BufferedFile bf = new BufferedFile(TestFileSetup.GetTargetFilename("at1x.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bf); bf.Flush(); bf.Close(); bf.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("at1x.fits"), FileAccess.Read); AsciiTable tab = (AsciiTable)f.GetHDU(1).Data; Object[] kern = (Object[])tab.Kernel; float[] fx = (float[])kern[0]; int[] ix = (int[])kern[1]; long[] lx = (long[])kern[2]; double[] dx = (double[])kern[3]; String[] sx = (String[])kern[4]; float[] fy = (float[])samp[0]; int[] iy = (int[])samp[1]; long[] ly = (long[])samp[2]; double[] dy = (double[])samp[3]; String[] sy = (String[])samp[4]; Assert.AreEqual(true, tab.IsNull(6, 0)); Assert.AreEqual(false, tab.IsNull(5, 0)); for (int i = 0; i < data.NRows; i += 1) { if (i != 5) { if (i != 6) { // Null Assert.AreEqual(1f, f2[i] / fx[i], Math.Pow(10, -6)); } Assert.AreEqual(iy[i], ix[i]); if (i == 4) { Assert.AreEqual(54321L, lx[i]); } else { Assert.AreEqual(ly[i], lx[i]); } Assert.AreEqual(1f, dy[i] / dx[i], Math.Pow(10, -14)); Assert.AreEqual(sy[i], sx[i].Trim()); } } Object[] r5 = (Object[])data.GetRow(5); String[] st = (String[])r5[4]; st[0] = st[0].Trim(); for (int i = 0; i < r5.Length; i++) { Assert.AreEqual(true, ArrayFuncs.ArrayEquals(row[i], r5[i], Math.Pow(10, -6), Math.Pow(10, -14))); } } finally { if (f != null) { f.Close(); } } }
public void TestSet() { Fits f = null; try { f = new Fits(TestFileSetup.GetTargetFilename("bt2.fits"), FileAccess.Read); f.Read(); BinaryTableHDU bhdu = (BinaryTableHDU)f.GetHDU(1); Header hdr = bhdu.Header; // Check the various set methods on variable length data. float[] dta = (float[])bhdu.GetElement(4, 1); dta = new float[] { 22, 21, 20 }; bhdu.SetElement(4, 1, dta); BufferedFile bdos = new BufferedFile( TestFileSetup.GetTargetFilename("bt2a.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bdos); bdos.Close(); bdos.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt2a.fits")); bhdu = (BinaryTableHDU)f.GetHDU(1); float[] xdta = (float[])bhdu.GetElement(4, 1); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(dta, xdta)); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(3, 1), vf[3])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(5, 1), vf[5])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), dta)); float[] tvf = new float[] { 101, 102, 103, 104 }; vf[4] = tvf; bhdu.SetColumn(1, vf); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(3, 1), vf[3])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), vf[4])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(5, 1), vf[5])); bdos = new BufferedFile( TestFileSetup.GetTargetFilename("bt2b.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bdos); bdos.Close(); bdos.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt2b.fits")); bhdu = (BinaryTableHDU)f.GetHDU(1); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(3, 1), vf[3])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), vf[4])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(5, 1), vf[5])); Object[] rw = (Object[])bhdu.GetRow(4); float[] trw = new float[] { -1, -2, -3, -4, -5, -6 }; rw[1] = trw; bhdu.SetRow(4, rw); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(3, 1), vf[3])); Assert.AreEqual(false, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), vf[4])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), trw)); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(5, 1), vf[5])); // bdos = new BufferedDataStream(new FileStream("bt2c.fits",FileMode.Open)); bdos = new BufferedFile( TestFileSetup.GetTargetFilename("bt2c.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); f.Write(bdos); bdos.Close(); bdos.Dispose(); f.Close(); f = new Fits(TestFileSetup.GetTargetFilename("bt2c.fits")); bhdu = (BinaryTableHDU)f.GetHDU(1); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(3, 1), vf[3])); Assert.AreEqual(false, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), vf[4])); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(4, 1), trw)); Assert.AreEqual(true, ArrayFuncs.ArrayEquals(bhdu.GetElement(5, 1), vf[5])); } finally { if (f != null) { f.Close(); } } }
public void TestRandomGroup() { //float[,] fa = new float[20,20]; float[][] fa = new float[20][]; for (int i = 0; i < fa.Length; i++) { fa[i] = new float[20]; } float[] pa = new float[3]; BufferedFile bf = new BufferedFile( TestFileSetup.GetTargetFilename("rg1.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); Object[][] data = new Object[1][]; data[0] = new Object[2]; data[0][0] = pa; data[0][1] = fa; Console.Out.WriteLine("***** Write header ******"); BasicHDU hdu = Fits.MakeHDU(data); Header hdr = hdu.Header; // Change the number of groups hdr.AddValue("GCOUNT", 20, "Number of groups"); hdr.Write(bf); Console.Out.WriteLine("***** Write data group by group ******"); for (int i = 0; i < 20; i += 1) { for (int j = 0; j < pa.Length; j += 1) { pa[j] = i + j; } for (int j = 0; j < fa.GetLength(0); j += 1) { try { // fa[j, j] = i * j; fa[j][j] = i * j; } catch (Exception e) { Console.WriteLine(e.Message); Console.WriteLine("i ,j value:" + i + " " + j); } } // Write a group bf.WriteArray(data); } byte[] padding = new byte[FitsUtil.Padding(20 * ArrayFuncs.ComputeSize(data))]; Console.Out.WriteLine("****** Write padding ******"); bf.Write(padding); bf.Flush(); bf.Close(); bf.Dispose(); Fits f = null; try { Console.Out.WriteLine("****** Read data back in ******"); f = new Fits(TestFileSetup.GetTargetFilename("rg1.fits")); BasicHDU[] hdus = f.Read(); data = (Object[][])hdus[0].Kernel; // data = hdus[0].Kernel; Console.Out.WriteLine("**** Check parameter and data info *****"); for (int i = 0; i < data.Length; i += 1) { pa = (float[])data[i][0]; // fa = (float[,]) data[i,1]; Array[] tfa = (Array[])data[i][1]; for (int j = 0; j < pa.Length; j += 1) { Assert.AreEqual((float)(i + j), pa[j]); } for (int j = 0; j < fa.Length; j += 1) { // Assert.AreEqual("dataTest:" + i + " " + j, (float)(i * j), fa[j,j]); Assert.AreEqual((float)(i * j), ((Array)tfa.GetValue(j)).GetValue(j)); } } f.Close(); Console.Out.WriteLine("**** Create HDU from kernel *****"); f = new Fits(); // Generate a FITS HDU from the kernel. f.AddHDU(Fits.MakeHDU(data)); bf = new BufferedFile( TestFileSetup.GetTargetFilename("rg2.fits"), FileAccess.ReadWrite, FileShare.ReadWrite); Console.Out.WriteLine("**** Write new file *****"); f.Write(bf); bf.Flush(); bf.Close(); bf.Dispose(); f.Close(); Console.Out.WriteLine("**** Read and check *****"); f = new Fits(TestFileSetup.GetTargetFilename("rg2.fits")); data = (Object[][])f.Read()[0].Kernel; for (int i = 0; i < data.Length; i += 1) { pa = (float[])data[i][0]; // fa = (float[,]) data[i,1]; Array[] tfa = (Array[])data[i][1]; for (int j = 0; j < pa.Length; j += 1) { Assert.AreEqual((float)(i + j), pa[j]); } for (int j = 0; j < fa.Length; j += 1) { //Assert.AreEqual("dataTest:" + i + " " + j, (float)(i * j), fa[j,j]); Assert.AreEqual((float)(i * j), ((Array)tfa.GetValue(j)).GetValue(j)); } } } finally { if (f != null) { f.Close(); } } }
public void Dispose() { m_BufferedFile.Dispose(); }