Package nom.tam.fits

Examples of nom.tam.fits.BinaryTableHDU

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
      double minMagErr = Double.MAX_VALUE;
      double maxMagErr = Double.MIN_VALUE;

      for (BasicHDU hdu : hdus) {
        if (hdu instanceof BinaryTableHDU) {
          BinaryTableHDU tableHDU = (BinaryTableHDU) hdu;

          for (int row = 0; row < tableHDU.getNRows()
              && !wasInterrupted(); row++) {
            try {
              double barytime = ((double[]) tableHDU.getElement(
                  row, 0))[0];
              float ap_corr_flux = ((float[]) tableHDU
                  .getElement(row, 7))[0];
              float ap_corr_err = ((float[]) tableHDU.getElement(
                  row, 8))[0];

              // Include only valid magnitude fluxes.
              // Question: why do we see such values in Kepler
              // data sets?
              if (!Float.isInfinite(ap_corr_flux)
                  && !Float.isInfinite(ap_corr_err)
                  && !Float.isNaN(ap_corr_flux)
                  && !Float.isNaN(ap_corr_err)) {

                double hjd = barytime + 2454833.0;
                double mag = 15.0 - 2.5
                    * Math.log(ap_corr_flux)
                    / Math.log(10.0);

                double magErr = 1.086 * ap_corr_err
                    / ap_corr_flux;

                if (magErr < minMagErr) {
                  minMagErr = magErr;
                } else if (magErr > maxMagErr) {
                  maxMagErr = magErr;
                }

                ValidObservation ob = new ValidObservation();
                ob.setDateInfo(new DateInfo(hjd));
                ob.setMagnitude(new Magnitude(mag, magErr));
                ob.setBand(keplerSeries);
                ob.setRecordNumber(row);
                collectObservation(ob);
              }
            } catch (Exception e) {
              String input = tableHDU.getRow(row).toString();
              String error = e.getLocalizedMessage();
              InvalidObservation ob = new InvalidObservation(
                  input, error);
              ob.setRecordNumber(row);
              addInvalidObservation(ob);
View Full Code Here
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
      double minMagErr = Double.MAX_VALUE;
      double maxMagErr = Double.MIN_VALUE;

      for (BasicHDU hdu : hdus) {
        if (hdu instanceof BinaryTableHDU) {
          BinaryTableHDU tableHDU = (BinaryTableHDU) hdu;

          for (int row = 0; row < tableHDU.getNRows()
              && !wasInterrupted(); row++) {
            try {
              int tmid = ((int[]) tableHDU.getElement(row, 0))[0];
              float flux = ((float[]) tableHDU.getElement(row, 1))[0];
              float fluxErr = ((float[]) tableHDU.getElement(row,
                  2))[0];

              float tamFlux = ((float[]) tableHDU.getElement(row,
                  3))[0];
              float tamFluxErr = ((float[]) tableHDU.getElement(
                  row, 4))[0];
              String imageId = (String) tableHDU.getElement(row,
                  5);
              short ccdX = ((short[]) tableHDU.getElement(row, 6))[0];
              short ccdY = ((short[]) tableHDU.getElement(row, 7))[0];
              // short flag = ((short[]) tableHDU.getElement(row,
              // 8))[0];

              if (flux > 1 && flux - fluxErr > 0) {
                double hjd = tmid / 86400.0 + jdRef;
                double mag = 15.0 - 2.5 * Math.log(flux)
                    / Math.log(10.0);
                double magErr = 1.086 * fluxErr / flux;

                if (magErr < minMagErr) {
                  minMagErr = magErr;
                } else if (magErr > maxMagErr) {
                  maxMagErr = magErr;
                }

                ValidObservation ob = new ValidObservation();
                ob.setDateInfo(new DateInfo(hjd));
                ob.setMagnitude(new Magnitude(mag, magErr));
                ob.setBand(superWaspSeries);
                ob.setRecordNumber(row);
                ob.addDetail("IMAGE_ID", imageId, "Image ID");
                ob.addDetail("CCDX", ccdX + "", "CCD X");
                ob.addDetail("CCDY", ccdY + "", "CCD Y");
                ob.addDetail("FLUX", tamFlux + "", "Flux");
                ob.addDetail("FLUXERR", tamFluxErr + "",
                    "Flux Error");
//                ob.addDetail("FLAG", flag + "", "Flag");
                obs.add(ob);
              }
            } catch (Exception e) {
              String input = tableHDU.getRow(row).toString();
              String error = e.getLocalizedMessage();
              InvalidObservation ob = new InvalidObservation(
                  input, error);
              ob.setRecordNumber(row);
              addInvalidObservation(ob);
View Full Code Here
145
146
147
148
149
150
151
152
153
154
155
156
157
158
        BinaryTable table = new BinaryTable();
        FitsFactory.setUseAsciiTables(false);
        table.addColumn(typeList.toArray(new String[typeList.size()]));
        table.addColumn(coordList.toArray(new String[coordList.size()]));
        table.addColumn(configList.toArray(new String[configList.size()]));
        BinaryTableHDU hdu = (BinaryTableHDU) Fits.makeHDU(table);
        hdu.getHeader().addValue("EXTNAME", extName, "Contains saved JSkyCat graphics");
        hdu.setColumnName(0, "type", null);
        hdu.setColumnName(1, "coords", null);
        hdu.setColumnName(2, "config", null);
        deleteHDU(extName);
        fits.addHDU(hdu);
        imageDisplay.checkExtensions(true);
    }
View Full Code Here
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
        if (name == null) {
            return null;
        }

        // find the "CATINFO" table, if it exists
        BinaryTableHDU hdu = findBinaryTableHDU(fits, CATINFO);
        if (hdu == null) {
            return null;
        }

        int ncols = hdu.getNCols();

        // find the entry for this table
        int rowIndex = findConfigEntryRow(hdu, name);
        if (rowIndex == -1) {
            return null;
        }

        // if we get here, we found the entry
        Object[] row = hdu.getRow(rowIndex);
        Properties properties = new Properties();
        properties.setProperty("serv_type", "local");
        properties.setProperty("long_name", name);
        properties.setProperty("url", getFilename());
        for (int colIndex = 0; colIndex < ncols; colIndex++) {
            String value = (String) row[colIndex];
            if (value.length() != 0) {
                properties.setProperty(hdu.getColumnName(colIndex).toLowerCase(), value);
            }
        }
        return new SkycatConfigEntry(properties);
    }
View Full Code Here
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
        //
        // XXX Note: the Fits class doesn't currently handle empty columns correctly,
        //           so we have to add one row with dummy values ("-").
        int ncols = CATINFO_COLUMNS.length;
        Object[][] data;
        BinaryTableHDU hdu = findBinaryTableHDU(fits, CATINFO);
        int newRowIndex = 0;
        if (hdu != null) {
            // merge data from old table to new table
            int nrows = hdu.getNRows();
            int entryIndex = findConfigEntryRow(hdu, name);
            int dummyIndex = findConfigEntryRow(hdu, "-");
            int n = nrows;
            if (entryIndex == -1) {
                n++;
            }
            if (dummyIndex == -1) {
                n++;
            }
            data = new Object[n][];
            for (int rowIndex = 0; rowIndex < nrows; rowIndex++) {
                if (rowIndex == entryIndex) {
                    continue;
                }
                if (rowIndex == dummyIndex) {
                    continue;
                }
                data[newRowIndex++] = _copyCatInfoRow(hdu, hdu.getRow(rowIndex));
            }
            deleteTable(fits, CATINFO);
        } else {
            data = new Object[2][];
        }

        // add a new row for the given table
        Object[] row = new Object[ncols];
        row[0] = name;  // use the given name rather than the original one
        if (table instanceof SkycatTable) {
            SkycatConfigEntry entry = ((SkycatTable) table).getConfigEntry();
            Properties properties = entry.getProperties();
            for (int colIndex = 1; colIndex < ncols; colIndex++) {
                String value = properties.getProperty(CATINFO_COLUMNS[colIndex].toLowerCase());
                if (value != null) {
                    row[colIndex] = value;
                } else {
                    row[colIndex] = "";
                }
            }
        } else {
            for (int colIndex = 1; colIndex < ncols; colIndex++) {
                row[colIndex] = "";
            }
        }
        data[newRowIndex++] = row;

        // add a dummy row as a workaround to a Fits lib bug with empty columns
        row = new Object[ncols];
        for (int colIndex = 0; colIndex < ncols; colIndex++) {
            row[colIndex] = "-";
        }
        data[newRowIndex] = row;

        BinaryTable binTable = new BinaryTable(data);
        FitsFactory.setUseAsciiTables(false);
        hdu = (BinaryTableHDU) Fits.makeHDU(binTable);
        hdu.getHeader().addValue("EXTNAME", CATINFO, "Contains catalog config info");
        for (int i = 0; i < ncols; i++) {
            hdu.setColumnName(i, CATINFO_COLUMNS[i], null);
        }
        fits.addHDU(hdu);
    }
View Full Code Here
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
     * @param name the name to give the FITS table
     */
    private static void _removeCatInfo(Fits fits, String name)
            throws FitsException, IOException {

        BinaryTableHDU hdu = findBinaryTableHDU(fits, CATINFO);
        if (hdu != null) {
            int rowIndex = findConfigEntryRow(hdu, name);
            if (rowIndex != -1) {
                hdu.deleteRows(rowIndex, 1);
            }
            if (hdu.getNRows() <= 1) {
                deleteTable(fits, CATINFO);
            }
        }
    }
View Full Code Here
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
     * @param filename The name of the FITS file
     * @param fits object to use for FITS I/O
     * @param navigator window object managing the table display and plotting
     */
    public static void plotTables(String filename, Fits fits, Navigator navigator) throws FitsException, IOException {
        BinaryTableHDU hdu = findBinaryTableHDU(fits, CATINFO);
        if (hdu == null) {
            return;
        }
        int nrows = hdu.getNRows();
        for (int rowIndex = 0; rowIndex < nrows; rowIndex++) {
            Object[] row = hdu.getRow(rowIndex);
            String name = (String) row[0];
            _plotTable(filename, fits, navigator, name);
        }
    }
View Full Code Here
640
641
642
643
644
645
646
647
648
649
650
     * @param filename The name of the FITS file
     * @param fits object to use for FITS I/O
     * @param name the name of the table to plot (FITS keyword EXTNAME)
     */
    private static TableQueryResult _findTable(String filename, Fits fits, String name) throws FitsException, IOException {
        BinaryTableHDU hdu = findBinaryTableHDU(fits, name);
        if (hdu == null) {
            return null;
        }

        return getTable(filename, fits, hdu);
View Full Code Here
TOP

Related Classes of nom.tam.fits.BinaryTableHDU

  • jsky.image.graphics.gui.FITSGraphics
  • jsky.navigator.NavigatorFITSTable
  • org.aavso.tools.vstar.external.plugin.KeplerFITSObservationSource$KeplerFITSObservationRetriever
  • org.aavso.tools.vstar.external.plugin.SuperWASPFITSObservationSource$FITSObservationRetriever
    • Copyright © 2018 www.massapicom. All rights reserved.
    All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.