source: branches/alma/external/atnf/PKSIO/NROFITSDataset.cc@ 1774

Last change on this file since 1774 was 1757, checked in by Kana Sugimoto, 14 years ago

New Development: Yes

JIRA Issue: Yes (CAS-2211)

Ready for Test: Yes

Interface Changes: Yes

What Interface Changed: ASAP 3.0.0 interface changes

Test Programs:

Put in Release Notes: Yes

Module(s): all the CASA sd tools and tasks are affected.

Description: Merged ATNF-ASAP 3.0.0 developments to CASA (alma) branch.

Note you also need to update casa/code/atnf.
File size: 74.9 KB
Line 
1//#---------------------------------------------------------------------------
2//# NROFITSDataset.cc: Class for NRO 45m FITS dataset.
3//#---------------------------------------------------------------------------
4//# Copyright (C) 2000-2006
5//# Associated Universities, Inc. Washington DC, USA.
6//#
7//# This library is free software; you can redistribute it and/or modify it
8//# under the terms of the GNU Library General Public License as published by
9//# the Free Software Foundation; either version 2 of the License, or (at your
10//# option) any later version.
11//#
12//# This library is distributed in the hope that it will be useful, but WITHOUT
13//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
14//# FITNESS FOR A PARTICULAR PURPOSE. See the GNU Library General Public
15//# License for more details.
16//#
17//# You should have received a copy of the GNU Library General Public License
18//# along with this library; if not, write to the Free Software Foundation,
19//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
20//#
21//# Correspondence concerning AIPS++ should be addressed as follows:
22//# Internet email: aips2-request@nrao.edu.
23//# Postal address: AIPS++ Project Office
24//# National Radio Astronomy Observatory
25//# 520 Edgemont Road
26//# Charlottesville, VA 22903-2475 USA
27//#
28//# $Id$
29//#---------------------------------------------------------------------------
30//# Original: 2009/02/27, Takeshi Nakazato, NAOJ
31//#---------------------------------------------------------------------------
32
33#include <atnf/PKSIO/NROFITSDataset.h>
34#include <scimath/Mathematics/InterpolateArray1D.h>
35
36#include <iostream>
37#include <fstream>
38#include <casa/math.h>
39#include <casa/iomanip.h>
40
41using namespace std ;
42
43//header size (8*2880bytes)
44#define FITS_HEADER_SIZE 23040
45
46// constructor
47NROFITSDataset::NROFITSDataset( string name )
48 : NRODataset( name )
49{
50 LogIO os( LogOrigin( "NROFITSDataset", "NROFITSDataset()", WHERE ) ) ;
51
52 fp_ = NULL ;
53 dataid_ = -1 ;
54 record_ = new NRODataRecord() ;
55 record_->LDATA = NULL ; // never use LDATA
56
57 // open file
58 if ( open() )
59 os << LogIO::SEVERE << "error while opening file " << filename_ << LogIO::EXCEPTION ;
60
61 // data initialization
62 readHeader( numField_, "TFIELDS", same_ ) ;
63 forms_.resize( numField_ ) ;
64 names_.resize( numField_ ) ;
65 units_.resize( numField_ ) ;
66 getField() ;
67
68 // check endian
69 // FITS file is always BIG_ENDIAN, but it is not true for NRO data
70 vector<int> itmp( 6 ) ;
71 if ( readTable( itmp, "LAVST", true, 0 ) != 0 ) {
72 os << LogIO::WARN << "Error while checking endian." << LogIO::POST ;
73 return ;
74 }
75// os << "itmp = " << itmp[0] << " "
76// << itmp[1] << " " << itmp[2] << " "
77// << itmp[3] << " " << itmp[4] << " "
78// << itmp[5] << LogIO::POST ;
79 // check endian by looking month value in LAVST (start time)
80 if ( itmp[1] > 0 && itmp[1] < 13 ) {
81 same_ = 1 ;
82 os << LogIO::NORMAL << "same endian " << LogIO::POST ;
83 }
84 else {
85 same_ = 0 ;
86 os << LogIO::NORMAL << "different endian " << LogIO::POST ;
87 }
88
89 // memory allocation
90 initialize() ;
91}
92
93// destructor
94NROFITSDataset::~NROFITSDataset()
95{
96 // close file
97 close() ;
98}
99
100// data initialization
101void NROFITSDataset::initialize()
102{
103 LogIO os( LogOrigin( "NROFITSDataset", "initialize()", WHERE ) ) ;
104
105 int status = 0 ;
106 status = readHeader( ARYNM, "ARYNM", same_ ) ;
107 if ( status != 0 )
108 ARYNM = 1 ;
109 readHeader( rowNum_, "NAXIS2", same_ ) ;
110 readHeader( scanLen_, "NAXIS1", same_ ) ;
111 status = 0 ;
112 scanNum_ = rowNum_ / ARYNM ;
113 int nchan = 0 ;
114 if ( readTable( nchan, "NCH", same_ ) != 0 ) {
115 os << LogIO::WARN << "Error while checking maximum channel number." << LogIO::POST ;
116 return ;
117 }
118 chmax_ = nchan ;
119 datasize_ = sizeof( int ) * chmax_ ;
120 //record_->JDATA.resize( chmax_ ) ;
121 JDATA.resize( chmax_ ) ;
122 // zero clear
123 for ( uInt i = 0 ; i < JDATA.size() ; i++ )
124 JDATA[i] = 0 ;
125
126 RX.resize( ARYNM ) ;
127 HPBW.resize( ARYNM ) ;
128 EFFA.resize( ARYNM ) ;
129 EFFB.resize( ARYNM ) ;
130 EFFL.resize( ARYNM ) ;
131 EFSS.resize( ARYNM ) ;
132 GAIN.resize( ARYNM ) ;
133 HORN.resize( ARYNM ) ;
134 POLTP.resize( ARYNM ) ;
135 POLDR.resize( ARYNM ) ;
136 POLAN.resize( ARYNM ) ;
137 DFRQ.resize( ARYNM ) ;
138 SIDBD.resize( ARYNM ) ;
139 REFN.resize( ARYNM ) ;
140 IPINT.resize( ARYNM ) ;
141 MULTN.resize( ARYNM ) ;
142 MLTSCF.resize( ARYNM ) ;
143 LAGWIND.resize( ARYNM ) ;
144 BEBW.resize( ARYNM ) ;
145 BERES.resize( ARYNM ) ;
146 CHWID.resize( ARYNM ) ;
147 ARRY.resize( NRO_FITS_ARYMAX ) ;
148 NFCAL.resize( ARYNM ) ;
149 F0CAL.resize( ARYNM ) ;
150 FQCAL.resize( ARYNM ) ;
151 CHCAL.resize( ARYNM ) ;
152 CWCAL.resize( ARYNM ) ;
153 DSBFC.resize( ARYNM ) ;
154
155 ARYTP.resize( ARYNM ) ;
156 arrayid_.resize( ARYNM ) ;
157 for ( int i = 0 ; i < ARYNM ; i++ )
158 arrayid_[i] = -1 ;
159
160 for ( int i = 0 ; i < ARYNM ; i++ ) {
161 FQCAL[i].resize( 10 ) ;
162 CHCAL[i].resize( 10 ) ;
163 CWCAL[i].resize( 10 ) ;
164 }
165
166 datasize_ += sizeof( char ) * ARYNM * 16 // RX
167 + sizeof( double ) * ARYNM * 6 // HPBW, EFFA, EFFB, EFFL, EFSS GAIN
168 + sizeof( char ) * ARYNM * 4 // HORN
169 + sizeof( char ) * ARYNM * 4 // POLTP
170 + sizeof( double ) * ARYNM * 3 // POLDR, POLAN, DFRQ
171 + sizeof( char ) * ARYNM * 4 // SIDBID
172 + sizeof( int ) * ARYNM * 3 // REFN, IPINT, MULTN
173 + sizeof( double ) * ARYNM // MLTSCF
174 + sizeof( char ) * ARYNM * 8 // LAGWIND
175 + sizeof( double ) * ARYNM * 3 // BEBW, BERES, CHWID
176 + sizeof( int ) * NRO_FITS_ARYMAX // ARRY
177 + sizeof( int ) * ARYNM // NFCAL
178 + sizeof( double ) * ARYNM // F0CAL
179 + sizeof( double ) * ARYNM * 10 * 3 // FQCAL, CHCAL, CWCAL
180 + sizeof( char ) * 180 ; // CDMY1
181}
182
183// fill data header
184int NROFITSDataset::fillHeader()
185{
186 LogIO os( LogOrigin( "NROFITSDataset", "fillHeader()", WHERE ) ) ;
187
188 // open file
189 if ( open() ) {
190 os << LogIO::SEVERE << "Error opening file " << filename_ << "." << LogIO::EXCEPTION ;
191 return -1 ;
192 }
193
194 // fill
195 int status = fillHeader( same_ ) ;
196
197 return status ;
198}
199
200int NROFITSDataset::fillHeader( int sameEndian )
201{
202 LogIO os( LogOrigin( "NROFITSDataset", "fillHeader()", WHERE ) ) ;
203
204 // fill array type
205 fillARYTP() ;
206
207 // read data header
208 float ftmp = 0.0 ;
209 if ( readHeader( LOFIL, "LOFIL" ) != 0 ) {
210 os << LogIO::NORMAL << "LOFIL set to 'FITS'." << LogIO::POST ;
211 LOFIL = "FITS" ;
212 }
213 // DEBUG
214 //cout << "LOFIL = \'" << LOFIL << "\'" << endl ;
215 //
216 if ( readHeader( VER, "VER" ) != 0 ) {
217 if ( readHeader( VER, "HISTORY NEWSTAR VER" ) != 0 ) {
218 os << LogIO::NORMAL << "VER set to 'V000'." << LogIO::POST ;
219 VER = "V000" ;
220 }
221 }
222 // DEBUG
223 //cout << "VER = \'" << VER << "\'" << endl ;
224 //
225 if ( readHeader( GROUP, "GROUP" ) != 0 ) {
226 if ( readHeader( GROUP, "HISTORY NEWSTAR GROUP" ) != 0 ) {
227 os << LogIO::NORMAL << "GROUP set to 'GRP0'." << LogIO::POST ;
228 GROUP = "GROUP0" ;
229 }
230 }
231 // DEBUG
232 //cout << "GROUP = \'" << GROUP << "\'" << endl ;
233 //
234 if ( readHeader( PROJ, "PROJECT" ) != 0 ) {
235 if ( readHeader( PROJ, "HISTORY NEWSTAR PROJECT" ) != 0 ) {
236 os << LogIO::NORMAL << "PROJ set to 'PROJ0'." << LogIO::POST ;
237 PROJ = "PROJECT0" ;
238 }
239 }
240 // DEBUG
241 //cout << "PROJ = \'" << PROJ << "\'" << endl ;
242 //
243 if ( readHeader( SCHED, "SCHED" ) != 0 ) {
244 if ( readHeader( SCHED, "HISTORY NEWSTAR SCHED" ) != 0 ) {
245 os << LogIO::NORMAL << "SCHED set to 'SCHED0'." << LogIO::POST ;
246 SCHED = "SCHED0" ;
247 }
248 }
249 // DEBUG
250 //cout << "SCHED = \'" << SCHED << "\'" << endl ;
251 //
252 if ( readHeader( OBSVR, "OBSERVER" ) != 0 ) {
253 os << LogIO::NORMAL << "OBSVR set to 'SOMEONE'" << LogIO::POST ;
254 OBSVR = "SOMEONE" ;
255 }
256 // DEBUG
257 //cout << "OBSVR = \'" << OBSVR << "\'" << endl ;
258 //
259 if ( readHeader( LOSTM, "STRSC" ) != 0 ) {
260 os << LogIO::WARN << "Error while reading data LOSTM." << LogIO::POST ;
261 return -1 ;
262 }
263 if ( LOSTM[0] == '9' ) {
264 LOSTM = "19" + LOSTM ;
265 }
266 else if ( LOSTM[0] == '0') {
267 LOSTM = "20" + LOSTM ;
268 }
269 // DEBUG
270 //cout << "LOSTM = \'" << LOSTM << "\'" << endl ;
271 //
272 if ( readHeader( LOETM, "STPSC" ) != 0 ) {
273 os << LogIO::WARN << "Error while reading data LOETM." << LogIO::POST ;
274 return -1 ;
275 }
276 if ( LOETM[0] == '9' ) {
277 LOETM = "19" + LOETM ;
278 }
279 else if ( LOETM[0] == '0') {
280 LOETM = "20" + LOETM ;
281 }
282 // DEBUG
283 //cout << "LOETM = \'" << LOETM << "\'" << endl ;
284 //
285 if ( readHeader( NSCAN, "NAXIS2", sameEndian ) != 0 ) {
286 os << LogIO::WARN << "Error while reading data NSCAN." << LogIO::POST ;
287 return -1 ;
288 }
289 // DEBUG
290 //cout << "NSCAN = " << NSCAN << endl ;
291 //
292 string subt ;
293 if ( readHeader( TITLE, "TITLE" ) != 0 ) {
294 int stat1 = readHeader( TITLE, "HISTORY NEWSTAR TITLE1" ) ;
295 int stat2 = readHeader( subt, "HISTORY NEWSTAR TITLE2" ) ;
296 if ( stat1 != 0 && stat2 != 0 ) {
297 os << LogIO::NORMAL << "TITLE set to 'NOTITLE'." << LogIO::POST ;
298 TITLE = "NOTITLE" ;
299 }
300 else {
301 //cout << "TITLE = \'" << TITLE << "\'" << endl ;
302 //cout << "subt = \'" << subt << "\'" << endl ;
303 TITLE = TITLE + subt ;
304 }
305 }
306 // DEBUG
307 //cout << "TITLE = \'" << TITLE << "\'" << endl ;
308 //
309 if ( readHeader( OBJ, "OBJECT" ) != 0 ) {
310 os << LogIO::NORMAL << "OBJ set to 'SOMEWHERE'." << LogIO::POST ;
311 OBJ = "NOOBJ" ;
312 }
313 // DEBUG
314 //cout << "OBJ = \'" << OBJ << "\'" << endl ;
315 //
316 if ( readHeader( ftmp, "EPOCH", sameEndian ) != 0 ) {
317 os << LogIO::WARN << "Error while reading data EPOCH." << LogIO::POST ;
318 return -1 ;
319 }
320 if ( ftmp == 1950.0 )
321 EPOCH = "B1950" ;
322 else if ( ftmp == 2000.0 )
323 EPOCH = "J2000" ;
324 else {
325 os << LogIO::WARN << "Unidentified epoch. set to 'XXXXX'" << LogIO::POST ;
326 EPOCH = "XXXXX" ;
327 }
328 // DEBUG
329 //cout << "EPOCH = \'" << EPOCH << "\'" << endl ;
330 //
331 string stmp ;
332 if ( readHeader( stmp, "RA" ) != 0 ) {
333 os << LogIO::WARN << "Error while reading data RA0." << LogIO::POST ;
334 return -1 ;
335 }
336 RA0 = radRA( stmp.c_str() ) ;
337 // DEBUG
338 //cout << "RA0 = " << RA0 << endl ;
339 //
340 if ( readHeader( stmp, "DEC" ) != 0 ) {
341 os << LogIO::WARN << "Error while reading data DEC0." << LogIO::POST ;
342 return -1 ;
343 }
344 DEC0 = radDEC( stmp.c_str() ) ;
345 // DEBUG
346 //cout << "DEC0 = " << DEC0 << endl ;
347 //
348 if ( readHeader( GLNG0, "GL0", sameEndian ) != 0 ) {
349 os << LogIO::WARN << "Error while reading data GLNG0." << LogIO::POST ;
350 return -1 ;
351 }
352 // DEBUG
353 //cout << "GLNG0 = " << GLNG0 << endl ;
354 //
355 if ( readHeader( GLAT0, "GB0", sameEndian ) != 0 ) {
356 os << LogIO::WARN << "Error while reading data GLAT0." << LogIO::POST ;
357 return -1 ;
358 }
359 // DEBUG
360 //cout << "GLAT0 = " << GLAT0 << endl ;
361 //
362 if ( readHeader( NCALB, "NCALB", sameEndian ) != 0 ) {
363 os << LogIO::WARN << "Error while reading data NCALB." << LogIO::POST ;
364 return -1 ;
365 }
366 // DEBUG
367 //cout << "NCALB = " << NCALB << endl ;
368 //
369 if ( readHeader( SCNCD, "SCNCD", sameEndian ) != 0 ) {
370 os << LogIO::NORMAL << "SCNCD set to 0 (RADEC)." << LogIO::POST ;
371 SCNCD = 0 ;
372 }
373 // DEBUG
374 //cout << "SCNCD = " << SCNCD << endl ;
375 //
376 if ( readHeader( SCMOD, "SCMOD1" ) != 0 ) {
377 os << LogIO::WARN << "Error while reading data SCMOD." << LogIO::POST ;
378 return -1 ;
379 }
380 string::size_type pos = SCMOD.find( ' ' ) ;
381 if ( pos != string::npos ) {
382 SCMOD = SCMOD.substr( 0, pos ) ;
383 SCMOD = SCMOD + " " ;
384 }
385 //cout << "SDMOD1 = \'" << SCMOD << "\'" << endl ;
386 if ( readHeader( stmp, "SCMOD2" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
387 if ( ( pos = stmp.find( ' ' ) ) != string::npos )
388 stmp = stmp.substr( 0, pos ) ;
389 SCMOD = SCMOD + stmp + " " ;
390 //cout << "SCMOD2 = \'" << SCMOD << "\'" << endl ;
391 }
392 if ( readHeader( stmp, "SCMOD3" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
393 if ( ( pos = stmp.find( ' ' ) ) != string::npos )
394 stmp = stmp.substr( 0, pos ) ;
395 SCMOD = SCMOD + stmp + " " ;
396 //cout << "SCMOD3 = \'" << SCMOD << "\'" << endl ;
397 }
398 if ( readHeader( stmp, "SCMOD4" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
399 if ( ( pos = stmp.find( ' ' ) ) != string::npos )
400 stmp = stmp.substr( 0, pos ) ;
401 SCMOD = SCMOD + stmp + " " ;
402 //cout << "SCMOD4 = \'" << SCMOD << "\'" << endl ;
403 }
404 if ( readHeader( stmp, "SCMOD5" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
405 if ( ( pos = stmp.find( ' ' ) ) != string::npos )
406 stmp = stmp.substr( 0, pos ) ;
407 SCMOD = SCMOD + stmp + " " ;
408 //cout << "SCMOD5 = \'" << SCMOD << "\'" << endl ;
409 }
410 if ( readHeader( stmp, "SCMOD6" ) == 0 && stmp.compare( 0, 1, " " ) != 0 ) {
411 if ( ( pos = stmp.find( ' ' ) ) != string::npos )
412 stmp = stmp.substr( 0, pos ) ;
413 SCMOD = SCMOD + stmp + " " ;
414 //cout << "SCMOD6 = \'" << SCMOD << "\'" << endl ;
415 }
416 // DEBUG
417 //cout << "SCMOD = \'" << SCMOD << "\'" << endl ;
418 //
419 if ( readHeader( URVEL, "VEL", sameEndian ) != 0 ) {
420 os << LogIO::WARN << "Error while reading data URVEL." << LogIO::POST ;
421 return -1 ;
422 }
423 // DEBUG
424 //cout << "URVEL = " << URVEL << endl ;
425 //
426 if ( readHeader( VREF, "VREF" ) != 0 ) {
427 os << LogIO::WARN << "Error while reading data VREF." << LogIO::POST ;
428 return -1 ;
429 }
430 // DEBUG
431 //cout << "VREF = \'" << VREF << "\'" << endl ;
432 //
433 if ( readHeader( VDEF, "VDEF" ) != 0 ) {
434 os << LogIO::WARN << "Error while reading data VDEF." << LogIO::POST ;
435 return -1 ;
436 }
437 // DEBUG
438 //cout << "VDEF = \'" << VDEF << "\'" << endl ;
439 //
440 if ( readHeader( SWMOD, "SWMOD" ) != 0 ) {
441 os << LogIO::WARN << "Error while reading data SWMOD." << LogIO::POST ;
442 return -1 ;
443 }
444 // DEBUG
445 //cout << "SWMOD = \'" << SWMOD << "\'" << endl ;
446 //
447 if ( readHeader( FRQSW, "FRQSW", sameEndian ) != 0 ) {
448 os << LogIO::NORMAL << "FRQSW set to 0." << LogIO::POST ;
449 FRQSW = 0.0 ;
450 }
451 // DEBUG
452 //cout << "FRQSW = " << FRQSW << endl ;
453 //
454 if ( readHeader( DBEAM, "DBEAM", sameEndian ) != 0 ) {
455 os << LogIO::WARN << "Error while reading data DBEAM." << LogIO::POST ;
456 return -1 ;
457 }
458 // DEBUG
459 //cout << "DBEAM = " << DBEAM << endl ;
460 //
461 if ( readHeader( MLTOF, "MLTOF", sameEndian ) != 0 ) {
462 os << LogIO::NORMAL << "MLTOF set to 0." << LogIO::POST ;
463 MLTOF = 0.0 ;
464 }
465 // DEBUG
466 //cout << "MLTOF = " << MLTOF << endl ;
467 //
468 if ( readHeader( CMTQ, "CMTQ", sameEndian ) != 0 ) {
469 os << LogIO::WARN << "Error while reading data CMTQ." << LogIO::POST ;
470 return -1 ;
471 }
472 // DEBUG
473 //cout << "CMTQ = " << CMTQ << endl ;
474 //
475 if ( readHeader( CMTE, "CMTE", sameEndian ) != 0 ) {
476 os << LogIO::WARN << "Error while reading data CMTE." << LogIO::POST ;
477 return -1 ;
478 }
479 // DEBUG
480 //cout << "CMTE = " << CMTE << endl ;
481 //
482 if ( readHeader( CMTSOM, "CMTSOM", sameEndian ) != 0 ) {
483 os << LogIO::WARN << "Error while reading data CMTSOM." << LogIO::POST ;
484 return -1 ;
485 }
486 // DEBUG
487 //cout << "CMTSOM = " << CMTSOM << endl ;
488 //
489 if ( readHeader( CMTNODE, "CMTNODE", sameEndian ) != 0 ) {
490 os << LogIO::WARN << "Error while reading data CMTNODE." << LogIO::POST ;
491 return -1 ;
492 }
493 // DEBUG
494 //cout << "CMTNODE = " << CMTNODE << endl ;
495 //
496 if ( readHeader( CMTI, "CMTI", sameEndian ) != 0 ) {
497 os << LogIO::WARN << "Error while reading data CMTI." << LogIO::POST ;
498 return -1 ;
499 }
500 // DEBUG
501 //cout << "CMTI = " << CMTI << endl ;
502 //
503 if ( readHeader( CMTTM, "CMTTM" ) != 0 ) {
504 os << LogIO::WARN << "Error while reading data CMTTM." << LogIO::POST ;
505 return -1 ;
506 }
507 // DEBUG
508 //cout << "CMTTM = \'" << CMTTM << "\'" << endl ;
509 //
510 if ( readHeader( SBDX, "SDBX", sameEndian ) != 0 ) {
511 os << LogIO::WARN << "Error while reading data SBDX." << LogIO::POST ;
512 return -1 ;
513 }
514 // DEBUG
515 //cout << "SBDX = " << SBDX << endl ;
516 //
517 if ( readHeader( SBDY, "SDBY", sameEndian ) != 0 ) {
518 os << LogIO::WARN << "Error while reading data SBDY." << LogIO::POST ;
519 return -1 ;
520 }
521 // DEBUG
522 //cout << "SBDY = " << SBDY << endl ;
523 //
524 if ( readHeader( SBDZ1, "SDBZ1", sameEndian ) != 0 ) {
525 os << LogIO::WARN << "Error while reading data SBDZ1." << LogIO::POST ;
526 return -1 ;
527 }
528 // DEBUG
529 //cout << "SBDZ1 = " << SBDZ1 << endl ;
530 //
531 if ( readHeader( SBDZ2, "SDBZ2", sameEndian ) != 0 ) {
532 os << LogIO::WARN << "Error while reading data SBDZ2." << LogIO::POST ;
533 return -1 ;
534 }
535 // DEBUG
536 //cout << "SBDZ2 = " << SBDZ2 << endl ;
537 //
538 if ( readHeader( DAZP, "DAZP", sameEndian ) != 0 ) {
539 os << LogIO::NORMAL << "DAZP set to 0." << LogIO::POST ;
540 DAZP = 0.0 ;
541 }
542 // DEBUG
543 //cout << "DAZP = " << DAZP << endl ;
544 //
545 if ( readHeader( DELP, "DELP", sameEndian ) != 0 ) {
546 os << LogIO::NORMAL << "DELP set to 0." << LogIO::POST ;
547 DELP = 0.0 ;
548 }
549 // DEBUG
550 //cout << "DELP = " << DELP << endl ;
551 //
552 if ( readHeader( CHBIND, "CHBIND", sameEndian ) != 0 ) {
553 os << LogIO::NORMAL << "CHBIND set to 1." << LogIO::POST ;
554 CHBIND = 1 ;
555 }
556 // DEBUG
557 //cout << "CHBIND = " << CHBIND << endl ;
558 //
559 if ( readHeader( NUMCH, "NCH", sameEndian ) != 0 ) {
560 if ( readTable( NUMCH, "NCH", sameEndian ) != 0 ) {
561 os << LogIO::NORMAL << "NUMCH set to " << chmax_ << "." << LogIO::POST ;
562 NUMCH = chmax_ ;
563 }
564 }
565 // DEBUG
566 //cout << "NUMCH = " << NUMCH << endl ;
567 //
568 if ( readHeader( CHMIN, "CHMIN", sameEndian ) != 0 ) {
569 os << LogIO::NORMAL << "CHMIN set to 1." << LogIO::POST ;
570 CHMIN = 1 ;
571 }
572 // DEBUG
573 //cout << "CHMIN = " << CHMIN << endl ;
574 //
575 if ( readHeader( CHMAX, "CHMAX", sameEndian ) != 0 ) {
576 os << LogIO::NORMAL << "CHMAX set to " << chmax_ << "." << LogIO::POST ;
577 CHMAX = chmax_ ;
578 }
579 // DEBUG
580 //cout << "CHMAX = " << CHMAX << endl ;
581 //
582 if ( readHeader( ALCTM, "ALCTM", sameEndian ) != 0 ) {
583 if ( readTable( ALCTM, "ALCTM", sameEndian ) != 0 ) {
584 os << LogIO::WARN << "Error while reading data ALCTM." << LogIO::POST ;
585 return -1 ;
586 }
587 }
588 // DEBUG
589 //cout << "ALCTM = " << ALCTM << endl ;
590 //
591 int itmp ;
592 if ( readHeader( itmp, "INTEG", sameEndian ) != 0 ) {
593 if ( readTable( itmp, "INTEG", sameEndian ) != 0 ) {
594 os << LogIO::WARN << "Error while reading data IPTIM." << LogIO::POST ;
595 return -1 ;
596 }
597 }
598 IPTIM = (double)itmp ;
599 // DEBUG
600 //cout << "IPTIM = " << IPTIM << endl ;
601 //
602 if ( readHeader( PA, "PA", sameEndian ) != 0 ) {
603 if ( readTable( PA, "PA", sameEndian ) != 0 ) {
604 os << LogIO::WARN << "Error while reading data PA." << LogIO::POST ;
605 return -1 ;
606 }
607 }
608 // DEBUG
609 //cout << "PA = " << PA << endl ;
610 //
611
612 // find data index for each ARYTP
613 findData() ;
614 vector<char *> v( ARYNM ) ;
615 if ( readColumn( RX, "RX" ) != 0 ) {
616 os << LogIO::WARN << "Error while reading data RX." << LogIO::POST ;
617 return -1 ;
618 }
619 // DEBUG
620// for ( int i = 0 ; i < ARYNM ; i++ ) {
621// if ( i == 0 ) {
622// //cout << "RX " ;
623// }
624// else if ( ( i % 5 ) == 0 ) {
625// //cout << endl << " " ;
626// }
627// //cout << "\'" << RX[i] << "\' " ;
628// }
629// //cout << endl ;
630 //
631 if ( readColumn( HPBW, "HPBW", sameEndian ) != 0 ) {
632 os << LogIO::WARN << "Error while reading data HPBW." << LogIO::POST ;
633 return -1 ;
634 }
635 // DEBUG
636// for ( int i = 0 ; i < ARYNM ; i++ ) {
637// if ( i == 0 ) {
638// //cout << "HPBW " ;
639// }
640// else if ( ( i % 5 ) == 0 ) {
641// //cout << endl << " " ;
642// }
643// //cout << HPBW[i] << " " ;
644// }
645// //cout << endl ;
646 //
647 if ( readColumn( EFFA, "EFFA", sameEndian ) != 0 ) {
648 os << LogIO::WARN << "Error while reading data EFFA." << LogIO::POST ;
649 return -1 ;
650 }
651 // DEBUG
652// for ( int i = 0 ; i < ARYNM ; i++ ) {
653// if ( i == 0 ) {
654// //cout << "EFFA " ;
655// }
656// else if ( ( i % 5 ) == 0 ) {
657// //cout << endl << " " ;
658// }
659// //cout << EFFA[i] << " " ;
660// }
661// //cout << endl ;
662 //
663 if ( readColumn( EFFB, "EFFB", sameEndian ) != 0 ) {
664 os << LogIO::WARN << "Error while reading data EFFB." << LogIO::POST ;
665 return -1 ;
666 }
667 // DEBUG
668// for ( int i = 0 ; i < ARYNM ; i++ ) {
669// if ( i == 0 ) {
670// //cout << "EFFB " ;
671// }
672// else if ( ( i % 5 ) == 0 ) {
673// //cout << endl << " " ;
674// }
675// //cout << EFFB[i] << " " ;
676// }
677// //cout << endl ;
678 //
679 if ( readColumn( EFFL, "EFFL", sameEndian ) != 0 ) {
680 os << LogIO::WARN << "Error while reading data EFFL." << LogIO::POST ;
681 return -1 ;
682 }
683 // DEBUG
684// for ( int i = 0 ; i < ARYNM ; i++ ) {
685// if ( i == 0 ) {
686// //cout << "EFFL " ;
687// }
688// else if ( ( i % 5 ) == 0 ) {
689// //cout << endl << " " ;
690// }
691// //cout << EFFL[i] << " " ;
692// }
693// //cout << endl ;
694 //
695 if ( readColumn( EFSS, "EFSS", sameEndian ) != 0 ) {
696 os << LogIO::WARN << "Error while reading data EFSS." << LogIO::POST ;
697 return -1 ;
698 }
699 // DEBUG
700// for ( int i = 0 ; i < ARYNM ; i++ ) {
701// if ( i == 0 ) {
702// //cout << "EFSS " ;
703// }
704// else if ( ( i % 5 ) == 0 ) {
705// //cout << endl << " " ;
706// }
707// //cout << EFSS[i] << " " ;
708// }
709// //cout << endl ;
710 //
711 if ( readColumn( GAIN, "GAIN", sameEndian ) != 0 ) {
712 os << LogIO::WARN << "Error while reading data GAIN." << LogIO::POST ;
713 return -1 ;
714 }
715 // DEBUG
716// for ( int i = 0 ; i < ARYNM ; i++ ) {
717// if ( i == 0 ) {
718// //cout << "GAIN " ;
719// }
720// else if ( ( i % 5 ) == 0 ) {
721// //cout << endl << " " ;
722// }
723// //cout << GAIN[i] << " " ;
724// }
725// //cout << endl ;
726 //
727 if ( readColumn( HORN, "HORN" ) != 0 ) {
728 os << LogIO::WARN << "Error while reading data HORN." << LogIO::POST ;
729 return -1 ;
730 }
731 // DEBUG
732// for ( int i = 0 ; i < ARYNM ; i++ ) {
733// if ( i == 0 ) {
734// //cout << "HORN " ;
735// }
736// else if ( ( i % 5 ) == 0 ) {
737// //cout << endl << " " ;
738// }
739// //cout << "\'" << HORN[i] << "\' " ;
740// }
741// //cout << endl ;
742 //
743 if ( readColumn( POLTP, "POLTP" ) != 0 ) {
744 os << LogIO::WARN << "Error while reading data POLTP." << LogIO::POST ;
745 return -1 ;
746 }
747 // DEBUG
748// for ( int i = 0 ; i < ARYNM ; i++ ) {
749// if ( i == 0 ) {
750// //cout << "POLTP " ;
751// }
752// else if ( ( i % 5 ) == 0 ) {
753// //cout << endl << " " ;
754// }
755// //cout << "\'" << POLTP[i] << "\' " ;
756// }
757// //cout << endl ;
758 //
759 vector<int> ipoldr( ARYNM, 0 ) ;
760 if ( readColumn( ipoldr, "POLDR", sameEndian ) != 0 ) {
761 os << LogIO::WARN << "Error while reading data POLDR." << LogIO::POST ;
762 return -1 ;
763 }
764 for ( int i = 0 ; i < ARYNM ; i++ )
765 POLDR[i] = (double)ipoldr[i] ;
766 // DEBUG
767// for ( int i = 0 ; i < ARYNM ; i++ ) {
768// if ( i == 0 ) {
769// //cout << "POLDR " ;
770// }
771// else if ( ( i % 5 ) == 0 ) {
772// //cout << endl << " " ;
773// }
774// //cout << POLDR[i] << " " ;
775// }
776// //cout << endl ;
777 //
778 if ( readColumn( POLAN, "POLAN", sameEndian ) != 0 ) {
779 os << LogIO::WARN << "Error while reading data POLAN." << LogIO::POST ;
780 return -1 ;
781 }
782 // DEBUG
783// for ( int i = 0 ; i < ARYNM ; i++ ) {
784// if ( i == 0 ) {
785// //cout << "POLAN " ;
786// }
787// else if ( ( i % 5 ) == 0 ) {
788// //cout << endl << " " ;
789// }
790// //cout << POLAN[i] << " " ;
791// }
792// //cout << endl ;
793 //
794 if ( readColumn( DFRQ, "DFRQ", sameEndian ) != 0 ) {
795 os << LogIO::WARN << "Error while reading data DFRQ." << LogIO::POST ;
796 return -1 ;
797 }
798 // DEBUG
799// for ( int i = 0 ; i < ARYNM ; i++ ) {
800// if ( i == 0 ) {
801// //cout << "DFRQ " ;
802// }
803// else if ( ( i % 5 ) == 0 ) {
804// //cout << endl << " " ;
805// }
806// //cout << DFRQ[i] << " " ;
807// }
808// //cout << endl ;
809 //
810 if ( readColumn( SIDBD, "SIDBD" ) != 0 ) {
811 os << LogIO::WARN << "Error while reading data SIDBD." << LogIO::POST ;
812 return -1 ;
813 }
814 // DEBUG
815// for ( int i = 0 ; i < ARYNM ; i++ ) {
816// if ( i == 0 ) {
817// //cout << "SIDBD " ;
818// }
819// else if ( ( i % 5 ) == 0 ) {
820// //cout << endl << " " ;
821// }
822// //cout << "\'" << SIDBD[i] << "\' " ;
823// }
824// //cout << endl ;
825 //
826 if ( readColumn( REFN, "REFN", sameEndian ) != 0 ) {
827 os << LogIO::WARN << "Error while reading data REFN." << LogIO::POST ;
828 return -1 ;
829 }
830 // DEBUG
831// for ( int i = 0 ; i < ARYNM ; i++ ) {
832// if ( i == 0 ) {
833// //cout << "REFN " ;
834// }
835// else if ( ( i % 5 ) == 0 ) {
836// //cout << endl << " " ;
837// }
838// //cout << REFN[i] << " " ;
839// }
840// //cout << endl ;
841 //
842 if ( readColumn( IPINT, "IPINT", sameEndian ) != 0 ) {
843 os << LogIO::WARN << "Error while reading data IPINT." << LogIO::POST ;
844 return -1 ;
845 }
846 // DEBUG
847// for ( int i = 0 ; i < ARYNM ; i++ ) {
848// if ( i == 0 ) {
849// //cout << "IPINT " ;
850// }
851// else if ( ( i % 5 ) == 0 ) {
852// //cout << endl << " " ;
853// }
854// //cout << IPINT[i] << " " ;
855// }
856// //cout << endl ;
857 //
858 if ( readColumn( MULTN, "MULTN", sameEndian ) != 0 ) {
859 os << LogIO::WARN << "Error while reading data MULTN." << LogIO::POST ;
860 return -1 ;
861 }
862 // DEBUG
863// for ( int i = 0 ; i < ARYNM ; i++ ) {
864// if ( i == 0 ) {
865// //cout << "MULTN " ;
866// }
867// else if ( ( i % 5 ) == 0 ) {
868// //cout << endl << " " ;
869// }
870// //cout << MULTN[i] << " " ;
871// }
872// //cout << endl ;
873 //
874 if ( readColumn( MLTSCF, "MLTSCF", sameEndian ) != 0 ) {
875 os << LogIO::WARN << "Error while reading data MLTSCF." << LogIO::POST ;
876 return -1 ;
877 }
878 // DEBUG
879// for ( int i = 0 ; i < ARYNM ; i++ ) {
880// if ( i == 0 ) {
881// //cout << "MLTSCF " ;
882// }
883// else if ( ( i % 5 ) == 0 ) {
884// //cout << endl << " " ;
885// }
886// //cout << MLTSCF[i] << " " ;
887// }
888// //cout << endl ;
889 //
890 if ( readColumn( LAGWIND, "LAGWIN" ) != 0 ) {
891 os << LogIO::WARN << "Error while reading data LAGWIND." << LogIO::POST ;
892 return -1 ;
893 }
894 // DEBUG
895// for ( int i = 0 ; i < ARYNM ; i++ ) {
896// if ( i == 0 ) {
897// //cout << "LAGWIND " ;
898// }
899// else if ( ( i % 5 ) == 0 ) {
900// //cout << endl << " " ;
901// }
902// //cout << "\'" << LAGWIND[i] << "\' " ;
903// }
904// //cout << endl ;
905 //
906 if ( readColumn( BEBW, "BEBW", sameEndian ) != 0 ) {
907 os << LogIO::WARN << "Error while reading data BEBW." << LogIO::POST ;
908 return -1 ;
909 }
910 // DEBUG
911// for ( int i = 0 ; i < ARYNM ; i++ ) {
912// if ( i == 0 ) {
913// //cout << "BEBW " ;
914// }
915// else if ( ( i % 5 ) == 0 ) {
916// //cout << endl << " " ;
917// }
918// //cout << BEBW[i] << " " ;
919// }
920// //cout << endl ;
921 //
922 if ( readColumn( BERES, "BERES", sameEndian ) != 0 ) {
923 os << LogIO::WARN << "Error while reading data BERES." << LogIO::POST ;
924 return -1 ;
925 }
926 // DEBUG
927// for ( int i = 0 ; i < ARYNM ; i++ ) {
928// if ( i == 0 ) {
929// //cout << "BERES " ;
930// }
931// else if ( ( i % 5 ) == 0 ) {
932// //cout << endl << " " ;
933// }
934// //cout << BERES[i] << " " ;
935// }
936// //cout << endl ;
937 //
938 if ( readColumn( CHWID, "CHWID", sameEndian ) != 0 ) {
939 os << LogIO::WARN << "Error while reading data CHWID." << LogIO::POST ;
940 return -1 ;
941 }
942 // DEBUG
943// for ( int i = 0 ; i < ARYNM ; i++ ) {
944// if ( i == 0 ) {
945// //cout << "CHWID " ;
946// }
947// else if ( ( i % 5 ) == 0 ) {
948// //cout << endl << " " ;
949// }
950// //cout << CHWID[i] << " " ;
951// }
952// //cout << endl ;
953 //
954 if ( readARRY() != 0 ) {
955 os << LogIO::WARN << "Error while reading data ARRY." << LogIO::POST ;
956 return -1 ;
957 }
958 // DEBUG
959// for ( int i = 0 ; i < NRO_FITS_ARYMAX ; i++ ) {
960// if ( i == 0 ) {
961// //cout << "ARRY " ;
962// }
963// else if ( ( i % 20 ) == 0 ) {
964// //cout << endl << " " ;
965// }
966// //cout << ARRY[i] << " " ;
967// }
968// //cout << endl ;
969 //
970 if ( readColumn( NFCAL, "NFCAL", sameEndian ) != 0 ) {
971 os << LogIO::WARN << "Error while reading data NFCAL." << LogIO::POST ;
972 return -1 ;
973 }
974 // DEBUG
975// for ( int i = 0 ; i < ARYNM ; i++ ) {
976// if ( i == 0 ) {
977// //cout << "NFCAL " ;
978// }
979// else if ( ( i % 5 ) == 0 ) {
980// //cout << endl << " " ;
981// }
982// //cout << NFCAL[i] << " " ;
983// }
984// //cout << endl ;
985 //
986 if ( readColumn( F0CAL, "F0CAL", sameEndian ) != 0 ) {
987 os << LogIO::WARN << "Error while reading data F0CAL." << LogIO::POST ;
988 return -1 ;
989 }
990 // DEBUG
991// for ( int i = 0 ; i < ARYNM ; i++ ) {
992// if ( i == 0 ) {
993// //cout << "F0CAL " ;
994// }
995// else if ( ( i % 5 ) == 0 ) {
996// //cout << endl << " " ;
997// }
998// //cout << F0CAL[i] << " " ;
999// }
1000// //cout << endl ;
1001 //
1002 for ( int i= 0 ; i < 10 ; i++) {
1003 vector<double> vv( ARYNM, 0 ) ;
1004 if ( readColumn( vv, "FQCAL", sameEndian, i ) != 0 ) {
1005 os << LogIO::WARN << "Error while reading data FQCAL." << LogIO::POST ;
1006 return -1 ;
1007 }
1008 for ( int j = 0 ; j < ARYNM ; j++ ) {
1009 FQCAL[j][i] = vv[j] ;
1010 }
1011 }
1012 // DEBUG
1013// for ( int i = 0 ; i < ARYNM ; i++ ) {
1014// for ( int j = 0 ; j < 10 ; j++ ) {
1015// if ( j == 0 ) {
1016// if ( i < 10 )
1017// //cout << "FQCAL0" << i << " " ;
1018// else
1019// //cout << "FQCAL" << i << " " ;
1020// }
1021// else if ( ( j % 5 ) == 0 ) {
1022// //cout << endl << " " ;
1023// }
1024// //cout << FQCAL[i][j] << " " ;
1025// }
1026// //cout << endl ;
1027// }
1028 //
1029 for ( int i= 0 ; i < 10 ; i++) {
1030 vector<double> vv( ARYNM, 0 ) ;
1031 if ( readColumn( vv, "CHCAL", sameEndian, i ) != 0 ) {
1032 os << LogIO::WARN << "Error while reading data CHCAL." << LogIO::POST ;
1033 return -1 ;
1034 }
1035 for ( int j = 0 ; j < ARYNM ; j++ ) {
1036 CHCAL[j][i] = vv[j] ;
1037 }
1038 }
1039 // DEBUG
1040// for ( int i = 0 ; i < ARYNM ; i++ ) {
1041// for ( int j = 0 ; j < 10 ; j++ ) {
1042// if ( j == 0 ) {
1043// if ( i < 10 )
1044// //cout << "CHCAL0" << i << " " ;
1045// else
1046// //cout << "CHCAL" << i << " " ;
1047// }
1048// else if ( ( j % 5 ) == 0 ) {
1049// //cout << endl << " " ;
1050// }
1051// //cout << CHCAL[i][j] << " " ;
1052// }
1053// //cout << endl ;
1054// }
1055 //
1056 for ( int i= 0 ; i < 10 ; i++) {
1057 vector<double> vv( ARYNM, 0 ) ;
1058 if ( readColumn( vv, "CWCAL", sameEndian, i ) != 0 ) {
1059 os << LogIO::WARN << "Error while reading data CWCAL." << LogIO::POST ;
1060 return -1 ;
1061 }
1062 for ( int j = 0 ; j < ARYNM ; j++ ) {
1063 CWCAL[j][i] = vv[j] ;
1064 }
1065 }
1066 // DEBUG
1067// for ( int i = 0 ; i < ARYNM ; i++ ) {
1068// for ( int j = 0 ; j < 10 ; j++ ) {
1069// if ( j == 0 ) {
1070// if ( i < 10 )
1071// //cout << "CWCAL0" << i << " " ;
1072// else
1073// //cout << "CWCAL" << i << " " ;
1074// }
1075// else if ( ( j % 5 ) == 0 ) {
1076// //cout << endl << " " ;
1077// }
1078// //cout << CWCAL[i][j] << " " ;
1079// }
1080// //cout << endl ;
1081// }
1082 //
1083 if ( readHeader( SCNLEN, "NAXIS1", sameEndian ) != 0 ) {
1084 os << LogIO::WARN << "Error while reading data SCNLEN." << LogIO::POST ;
1085 return -1 ;
1086 }
1087 // DEBUG
1088 //cout << "SCNLEN = " << SCNLEN << endl ;
1089 //
1090 if ( readHeader( SBIND, "SBIND", sameEndian ) != 0 ) {
1091 os << LogIO::NORMAL << "SBIND set to 0." << LogIO::POST ;
1092 SBIND = 0 ;
1093 }
1094 // DEBUG
1095 //cout << "SBIND = " << SBIND << endl ;
1096 //
1097 if ( readHeader( IBIT, "IBIT", sameEndian ) != 0 ) {
1098 os << LogIO::NORMAL << "IBIT set to 8." << LogIO::POST ;
1099 IBIT = 8 ; // 8 bit? 12 bit?
1100 }
1101 // DEBUG
1102 //cout << "IBIT = " << IBIT << endl ;
1103 //
1104 if ( readHeader( SITE, "TELESCOP" ) != 0 ) {
1105 os << LogIO::WARN << "Error while reading data SITE." << LogIO::POST ;
1106 return -1 ;
1107 }
1108 // DEBUG
1109 //cout << "SITE = \'" << SITE << "\'" << endl ;
1110 //
1111 if ( readColumn( DSBFC, "DSBFC", sameEndian ) != 0 ) {
1112 os << LogIO::NORMAL << "All DSBFC elements set to 1." << LogIO::POST ;
1113 for ( int i = 0 ; i < ARYNM ; i++ )
1114 DSBFC[i] = 1.0 ;
1115 }
1116 // DEBUG
1117// for ( int i = 0 ; i < ARYNM ; i++ ) {
1118// if ( i == 0 ) {
1119// //cout << "DSBFC " ;
1120// }
1121// else if ( ( i % 5 ) == 0 ) {
1122// //cout << endl << " " ;
1123// }
1124// //cout << DSBFC[i] << " " ;
1125// }
1126// //cout << endl ;
1127 //
1128
1129 show() ;
1130
1131 return 0 ;
1132}
1133
1134int NROFITSDataset::fillRecord( int i )
1135{
1136 LogIO os( LogOrigin( "NROFITSDataset", "fillRecord()", WHERE ) ) ;
1137
1138 int status = 0 ;
1139 string str4( 4, ' ' ) ;
1140 string str8( 8, ' ' ) ;
1141 string str24( 24, ' ' ) ;
1142
1143 strcpy( record_->LSFIL, str4.c_str() ) ;
1144 status = readTable( record_->LSFIL, "LSFIL", 4, i ) ;
1145 if ( status ) {
1146 os << LogIO::WARN << "Error while reading LSFIL." << LogIO::POST ;
1147 return status ;
1148 }
1149 // DEBUG
1150 //cout << "LSFIL(" << i << ") = " << record_->LSFIL << endl ;
1151 //
1152 status = readTable( record_->ISCAN, "ISCN", same_, i ) ;
1153 if ( status ) {
1154 os << LogIO::WARN << "Error while reading ISCAN." << LogIO::POST ;
1155 return status ;
1156 }
1157 // DEBUG
1158 //cout << "ISCAN(" << i << ") = " << record_->ISCAN << endl ;
1159 //
1160 vector<int> itmp( 6, 0 ) ;
1161 status = readTable( itmp, "LAVST", same_, i ) ;
1162 if ( status ) {
1163 os << LogIO::WARN << "Error while reading LAVST." << LogIO::POST ;
1164 return status ;
1165 }
1166 else {
1167 sprintf( record_->LAVST, "%d%d%d%d%d%d.000", itmp[0], itmp[1], itmp[2], itmp[3], itmp[4], itmp[5] ) ;
1168 }
1169 // DEBUG
1170 //cout << "LAVST(" << i << ") = " << record_->LAVST << endl ;
1171 //
1172 strcpy( record_->SCANTP, str8.c_str() ) ;
1173 status = readTable( record_->SCANTP, "SCNTP", strlen(record_->SCANTP), i ) ;
1174 if ( status ) {
1175 os << LogIO::WARN << "Error while reading SCANTP." << LogIO::POST ;
1176 return status ;
1177 }
1178 // DEBUG
1179 //cout << "SCANTP(" << i << ") = " << record_->SCANTP << endl ;
1180 //
1181 char *name1 = "" ;
1182 char *name2 = "" ;
1183 if ( SCNCD == 0 ) {
1184 name1 = "DRA" ;
1185 name2 = "DDEC" ;
1186 }
1187 else if ( SCNCD == 1 ) {
1188 name1 = "DGL" ;
1189 name2 = "DGB" ;
1190 }
1191 else {
1192 name1 = "DAZ" ;
1193 name2 = "DEL" ;
1194 }
1195 status = readTable( record_->DSCX, name1, same_, i ) ;
1196 if ( status ) {
1197 os << LogIO::WARN << "Error while reading DSCX." << LogIO::POST ;
1198 return status ;
1199 }
1200 // DEBUG
1201 //cout << "DSCX(" << i << ") = " << record_->DSCX << endl ;
1202 //
1203 status = readTable( record_->DSCY, name2, same_, i ) ;
1204 if ( status ) {
1205 os << LogIO::WARN << "Error while reading DSCY." << LogIO::POST ;
1206 return status ;
1207 }
1208 // DEBUG
1209 //cout << "DSCY(" << i << ") = " << record_->DSCY << endl ;
1210 //
1211 if ( SCNCD == 0 ) {
1212 name1 = "RA" ;
1213 name2 = "DEC" ;
1214 }
1215 else if ( SCNCD == 1 ) {
1216 name1 = "GL" ;
1217 name2 = "GB" ;
1218 }
1219 else {
1220 name1 = "AZ" ;
1221 name2 = "EL" ;
1222 }
1223 status = readTable( record_->SCX, name1, same_, i ) ;
1224 if ( status ) {
1225 os << LogIO::WARN << "Error while reading SCX." << LogIO::POST ;
1226 return status ;
1227 }
1228 // DEBUG
1229 //cout << "SCX(" << i << ") = " << record_->SCX << endl ;
1230 //
1231 status = readTable( record_->SCY, name2, same_, i ) ;
1232 if ( status ) {
1233 os << LogIO::WARN << "Error while reading SCY." << LogIO::POST ;
1234 return status ;
1235 }
1236 // DEBUG
1237 //cout << "SCY(" << i << ") = " << record_->SCY << endl ;
1238 //
1239 status = readTable( record_->PAZ, "PAZ", same_, i ) ;
1240 if ( status ) {
1241 os << LogIO::WARN << "Error while reading PAZ." << LogIO::POST ;
1242 return status ;
1243 }
1244 // DEBUG
1245 //cout << "PAZ(" << i << ") = " << record_->PAZ << endl ;
1246 //
1247 status = readTable( record_->PEL, "PEL", same_, i ) ;
1248 if ( status ) {
1249 os << LogIO::WARN << "Error while reading PEL." << LogIO::POST ;
1250 return status ;
1251 }
1252 // DEBUG
1253 //cout << "PEL(" << i << ") = " << record_->PEL << endl ;
1254 //
1255 status = readTable( record_->RAZ, "RAZ", same_, i ) ;
1256 if ( status ) {
1257 os << LogIO::WARN << "Error while reading RAZ." << LogIO::POST ;
1258 return status ;
1259 }
1260 // DEBUG
1261 //cout << "RAZ(" << i << ") = " << record_->RAZ << endl ;
1262 //
1263 status = readTable( record_->REL, "REL", same_, i ) ;
1264 if ( status ) {
1265 os << LogIO::WARN << "Error while reading REL." << LogIO::POST ;
1266 return status ;
1267 }
1268 // DEBUG
1269 //cout << "REL(" << i << ") = " << record_->REL << endl ;
1270 //
1271 status = readTable( record_->XX, "XX", same_, i ) ;
1272 if ( status ) {
1273 os << LogIO::WARN << "Error while reading XX." << LogIO::POST ;
1274 return status ;
1275 }
1276 // DEBUG
1277 //cout << "XX(" << i << ") = " << record_->XX << endl ;
1278 //
1279 status = readTable( record_->YY, "YY", same_, i ) ;
1280 if ( status ) {
1281 os << LogIO::WARN << "Error while reading YY." << LogIO::POST ;
1282 return status ;
1283 }
1284 // DEBUG
1285 //cout << "YY(" << i << ") = " << record_->YY << endl ;
1286 //
1287 strcpy( record_->ARRYT, str4.c_str() ) ;
1288 status = readTable( record_->ARRYT, "ARRYT", strlen(record_->ARRYT), i ) ;
1289 record_->ARRYT[3] = '\0' ;
1290 if ( status ) {
1291 os << LogIO::WARN << "Error while reading ARRYT." << LogIO::POST ;
1292 return status ;
1293 }
1294 // DEBUG
1295 //cout << "ARRYT(" << i << ") = " << record_->ARRYT << endl ;
1296 //
1297 double dtmp ;
1298 status = readTable( dtmp, "TEMP", same_, i ) ;
1299 if ( status ) {
1300 os << LogIO::WARN << "Error while reading TEMP." << LogIO::POST ;
1301 return status ;
1302 }
1303 else {
1304 record_->TEMP = dtmp ;
1305 }
1306 // DEBUG
1307 //cout << "TEMP(" << i << ") = " << record_->TEMP << endl ;
1308 //
1309 status = readTable( dtmp, "PATM", same_, i ) ;
1310 if ( status ) {
1311 os << LogIO::WARN << "Error while reading PATM." << LogIO::POST ;
1312 return status ;
1313 }
1314 else {
1315 record_->PATM = dtmp ;
1316 }
1317 // DEBUG
1318 //cout << "PATM(" << i << ") = " << record_->PATM << endl ;
1319 //
1320 status = readTable( dtmp, "PH2O", same_, i ) ;
1321 if ( status ) {
1322 os << LogIO::WARN << "Error while reading PH2O." << LogIO::POST ;
1323 return status ;
1324 }
1325 else {
1326 record_->PH2O = dtmp ;
1327 }
1328 // DEBUG
1329 //cout << "PH2O(" << i << ") = " << record_->PH2O << endl ;
1330 //
1331 status = readTable( dtmp, "VWIND", same_, i ) ;
1332 if ( status ) {
1333 os << LogIO::WARN << "Error while reading VWIND." << LogIO::POST ;
1334 return status ;
1335 }
1336 else {
1337 record_->VWIND = dtmp ;
1338 }
1339 // DEBUG
1340 //cout << "VWIND(" << i << ") = " << record_->VWIND << endl ;
1341 //
1342 status = readTable( dtmp, "DWIND", same_, i ) ;
1343 if ( status ) {
1344 os << LogIO::WARN << "Error while reading DWIND." << LogIO::POST ;
1345 return status ;
1346 }
1347 else {
1348 record_->DWIND = dtmp ;
1349 }
1350 // DEBUG
1351 //cout << "DWIND(" << i << ") = " << record_->DWIND << endl ;
1352 //
1353 status = readTable( dtmp, "TAU", same_, i ) ;
1354 if ( status ) {
1355 os << LogIO::WARN << "Error while reading TAU." << LogIO::POST ;
1356 return status ;
1357 }
1358 else {
1359 record_->TAU = dtmp ;
1360 }
1361 // DEBUG
1362 //cout << "TAU(" << i << ") = " << record_->TAU << endl ;
1363 //
1364 status = readTable( dtmp, "TSYS", same_, i ) ;
1365 if ( status ) {
1366 os << LogIO::WARN << "Error while reading TSYS." << LogIO::POST ;
1367 return status ;
1368 }
1369 else {
1370 record_->TSYS = dtmp ;
1371 }
1372 // DEBUG
1373 //cout << "TSYS(" << i << ") = " << record_->TSYS << endl ;
1374 //
1375 status = readTable( dtmp, "BATM", same_, i ) ;
1376 if ( status ) {
1377 os << LogIO::WARN << "Error while reading BATM." << LogIO::POST ;
1378 return status ;
1379 }
1380 else {
1381 record_->BATM = dtmp ;
1382 }
1383 // DEBUG
1384 //cout << "BATM(" << i << ") = " << record_->BATM << endl ;
1385 //
1386 status = readTable( record_->VRAD, "VRAD", same_, i ) ;
1387 if ( status ) {
1388 os << LogIO::WARN << "Error while reading TEMP." << LogIO::POST ;
1389 return status ;
1390 }
1391 // DEBUG
1392 //cout << "VRAD(" << i << ") = " << record_->VRAD << endl ;
1393 //
1394 status = readTable( record_->FREQ0, "FRQ0", same_, i ) ;
1395 if ( status ) {
1396 os << LogIO::WARN << "Error while reading FREQ0." << LogIO::POST ;
1397 return status ;
1398 }
1399 // DEBUG
1400 //cout << "FREQ0(" << i << ") = " << record_->FREQ0 << endl ;
1401 //
1402 status = readTable( record_->FQTRK, "FQTRK", same_, i ) ;
1403 if ( status ) {
1404 os << LogIO::WARN << "Error while reading FQTRK." << LogIO::POST ;
1405 return status ;
1406 }
1407 // DEBUG
1408 //cout << "FQTRK(" << i << ") = " << record_->FQTRK << endl ;
1409 //
1410 status = readTable( record_->FQIF1, "FQIF1", same_, i ) ;
1411 if ( status ) {
1412 os << LogIO::WARN << "Error while reading FQIF1." << LogIO::POST ;
1413 return status ;
1414 }
1415 // DEBUG
1416 //cout << "FQIF1(" << i << ") = " << record_->FQIF1 << endl ;
1417 //
1418 status = readTable( record_->ALCV, "ALCV", same_, i ) ;
1419 if ( status ) {
1420 os << LogIO::WARN << "Error while reading ALCV." << LogIO::POST ;
1421 return status ;
1422 }
1423 // DEBUG
1424 //cout << "ALCV(" << i << ") = " << record_->ALCV << endl ;
1425 //
1426 record_->IDMY0 = 0 ;
1427 status = readTable( record_->DPFRQ, "DPFRQ", same_, i ) ;
1428 if ( status ) {
1429 //os << LogIO::WARN << "Error DPFRQ set to 0." << LogIO::POST ;
1430 record_->DPFRQ = 0.0 ;
1431 }
1432 // DEBUG
1433 //cout << "DPFRQ(" << i << ") = " << record_->DPFRQ << endl ;
1434 //
1435 status = readTable( record_->SFCTR, "SFCTR", same_, i ) ;
1436 if ( status ) {
1437 os << LogIO::WARN << "Error while reading SFCTR." << LogIO::POST ;
1438 return status ;
1439 }
1440 // DEBUG
1441 //cout << "SFCTR(" << i << ") = " << record_->SFCTR << endl ;
1442 //
1443 status = readTable( record_->ADOFF, "ADOFF", same_, i ) ;
1444 if ( status ) {
1445 os << LogIO::WARN << "Error while reading ADOFF." << LogIO::POST ;
1446 return status ;
1447 }
1448 // DEBUG
1449 //cout << "ADOFF(" << i << ") = " << record_->ADOFF << endl ;
1450 //
1451 //status = readTable( record_->JDATA, "LDATA", same_, i ) ;
1452 status = readTable( JDATA, "LDATA", same_, i ) ;
1453 if ( status ) {
1454 os << LogIO::WARN << "Error while reading JDATA." << LogIO::POST ;
1455 return status ;
1456 }
1457 // DEBUG
1458// for ( int i = 0 ; i < chmax_ ; i++ )
1459// //cout << "JDATA[" << i << "] = " << JDATA[i] << " " ;
1460// //cout << endl ;
1461 //
1462 return status ;
1463}
1464
1465vector< vector<double> > NROFITSDataset::getSpectrum()
1466{
1467 vector< vector<double> > spec;
1468
1469 for ( int i = 0 ; i < rowNum_ ; i++ ) {
1470 spec.push_back( getSpectrum( i ) ) ;
1471 }
1472
1473 return spec ;
1474}
1475
1476vector<double> NROFITSDataset::getSpectrum( int i )
1477{
1478 vector<double> spec( chmax_, 0.0 ) ;
1479 vector<double> specout( chmax_, 0.0 ) ;
1480 NRODataRecord *record = getRecord( i ) ;
1481 double scale = record->SFCTR ;
1482 double offset = record->ADOFF ;
1483 double dscale = MLTSCF[getIndex( i )] ;
1484 //vector<int> ispec = record->JDATA ;
1485 vector<int> ispec = JDATA ;
1486 for ( int ii = 0 ; ii < chmax_ ; ii++ ) {
1487 spec[ii] = (double)( ispec[ii] * scale + offset ) * dscale ;
1488 }
1489
1490 // for AOS, re-gridding is needed
1491 if ( strncmp( record->ARRYT, "H", 1 ) == 0
1492 || strncmp( record->ARRYT, "W", 1 ) == 0
1493 || strncmp( record->ARRYT, "U", 1 ) == 0 ) {
1494
1495 string arryt = string( record->ARRYT ) ;
1496 uInt ib = getArrayId( arryt ) ;
1497 vector<double> fqcal = getFQCAL()[ib] ;
1498 vector<double> chcal = getCHCAL()[ib] ;
1499 int ncal = getNFCAL()[ib] ;
1500
1501// //cout << "NRODataset::getFrequencies() ncal = " << ncal << endl ;
1502 while ( ncal < (int)fqcal.size() ) {
1503 fqcal.pop_back() ;
1504 chcal.pop_back() ;
1505 }
1506 Vector<Double> xin( chcal ) ;
1507 Vector<Double> yin( fqcal ) ;
1508 int nchan = getNUMCH() ;
1509 Vector<Double> xout( nchan ) ;
1510 indgen( xout ) ;
1511 Vector<Double> yout ;
1512 InterpolateArray1D<Double, Double>::interpolate( yout, xout, xin, yin, InterpolateArray1D<Double,Double>::cubic ) ;
1513 // debug
1514 //cout << "i=" << i << endl ;
1515 if ( i == 16 ) {
1516 ofstream ofs0( "spgrid0.dat" ) ;
1517 for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
1518 ofs0 << xout[ii] << "," ;
1519 ofs0 << endl ;
1520 for ( int ii = 0 ; ii < getNUMCH() ; ii++ )
1521 ofs0 << setprecision(16) << record->FREQ0+yout[ii] << "," ;
1522 ofs0 << endl ;
1523 ofs0.close() ;
1524 }
1525 //
1526 Vector<Double> z( nchan ) ;
1527 Double bw = abs( yout[nchan-1] - yout[0] ) ;
1528 bw += 0.5 * abs( yout[nchan-1] - yout[nchan-2] + yout[1] - yout[0] ) ;
1529 Double dz = bw / (Double)nchan ;
1530 if ( yout[0] > yout[nchan-1] )
1531 dz = - dz ;
1532 z[0] = yout[0] - 0.5 * ( yout[1] - yout[0] - dz ) ;
1533 for ( int ii = 1 ; ii < nchan ; ii++ )
1534 z[ii] = z[ii-1] + dz ;
1535 Vector<Double> zi( nchan+1 ) ;
1536 Vector<Double> yi( nchan+1 ) ;
1537 zi[0] = z[0] - 0.5 * dz ;
1538 zi[1] = z[0] + 0.5 * dz ;
1539 yi[0] = yout[0] - 0.5 * ( yout[1] - yout[0] ) ;
1540 yi[1] = yout[0] + 0.5 * ( yout[1] - yout[0] ) ;
1541 for ( int ii = 2 ; ii < nchan ; ii++ ) {
1542 zi[ii] = zi[ii-1] + dz ;
1543 yi[ii] = yi[ii-1] + 0.5 * ( yout[ii] - yout[ii-2] ) ;
1544 }
1545 zi[nchan] = z[nchan-1] + 0.5 * dz ;
1546 yi[nchan] = yout[nchan-1] + 0.5 * ( yout[nchan-1] - yout[nchan-2] ) ;
1547// // debug
1548// //cout << "nchan=" << nchan << ", bw=" << bw << ", dz=" << dz
1549// << ", y[1]-y[0]=" << yout[1]-yout[0] << endl ;
1550// //cout << "z: " << z[0] << " - " << z[nchan-1]
1551// << ", zi: " << zi[0] << " - " << zi[nchan] << endl ;
1552// //cout << "y: " << yout[0] << " - " << yout[nchan-1]
1553// << ", yi: " << yi[0] << " - " << yi[nchan] << endl ;
1554// ofstream ofs1( "spgrid1.dat", ios::out | ios::app ) ;
1555// ofs1 << "spid=" << i << ", ARRYT=" << record->ARRYT << endl ;
1556// ofs1 << "z[0]=" << z[0] << ", yout[0]=" << yout[0] << endl ;
1557// for ( int ii = 1; ii < nchan ; ii++ ) {
1558// ofs1 << " dz=" << z[ii]-z[ii-1] << ", dy=" << yout[ii]-yout[ii-1] << endl ;
1559// ofs1 << "z[" << ii << "]=" << z[ii] << ", yout[" << ii << "]=" << yout[ii] << endl ;
1560// }
1561// ofs1.close() ;
1562// ofstream ofs2( "spgrid2.dat", ios::out | ios::app ) ;
1563// ofs2 << "spid=" << i << ", ARRYT=" << record->ARRYT << endl ;
1564// for ( int ii = 0 ; ii < nchan+1 ; ii++ )
1565// ofs2 << "zi[" << ii << "]=" << zi[ii] << ", yi[" << ii << "]=" << yi[ii] << endl ;
1566// ofs2.close() ;
1567// //
1568 int ichan = 0 ;
1569 double wsum = 0.0 ;
1570 // debug
1571 //ofstream ofs3( "spgrid3.dat", ios::out | ios::app ) ;
1572 if ( dz > 0.0 ) {
1573 for ( int ii = 0 ; ii < nchan ; ii++ ) {
1574 double zl = zi[ii] ;
1575 double zr = zi[ii+1] ;
1576 for ( int j = ichan ; j < nchan ; j++ ) {
1577 double yl = yi[j] ;
1578 double yr = yi[j+1] ;
1579 if ( yl <= zl ) {
1580 if ( yr <= zl ) {
1581 continue ;
1582 }
1583 else if ( yr <= zr ) {
1584 specout[ii] += spec[j] * ( yr - zl ) ;
1585 wsum += ( yr - zl ) ;
1586 }
1587 else {
1588 specout[ii] += spec[j] * dz ;
1589 wsum += dz ;
1590 ichan = j ;
1591 break ;
1592 }
1593 }
1594 else if ( yl < zr ) {
1595 if ( yr <= zr ) {
1596 specout[ii] += spec[j] * ( yr - yl ) ;
1597 wsum += ( yr - yl ) ;
1598 }
1599 else {
1600 specout[ii] += spec[j] * ( zr - yl ) ;
1601 wsum += ( zr - yl ) ;
1602 ichan = j ;
1603 break ;
1604 }
1605 }
1606 else {
1607 ichan = j - 1 ;
1608 break ;
1609 }
1610 }
1611 specout[ii] /= wsum ;
1612 wsum = 0.0 ;
1613 }
1614 }
1615 else if ( dz < 0.0 ) {
1616 for ( int ii = 0 ; ii < nchan ; ii++ ) {
1617 double zl = zi[ii] ;
1618 double zr = zi[ii+1] ;
1619 for ( int j = ichan ; j < nchan ; j++ ) {
1620 double yl = yi[j] ;
1621 double yr = yi[j+1] ;
1622 if ( yl >= zl ) {
1623 if ( yr >= zl ) {
1624 continue ;
1625 }
1626 else if ( yr >= zr ) {
1627 specout[ii] += spec[j] * abs( yr - zl ) ;
1628 wsum += abs( yr - zl ) ;
1629 }
1630 else {
1631 specout[ii] += spec[j] * abs( dz ) ;
1632 wsum += abs( dz ) ;
1633 ichan = j ;
1634 break ;
1635 }
1636 }
1637 else if ( yl > zr ) {
1638 if ( yr >= zr ) {
1639 specout[ii] += spec[j] * abs( yr - yl ) ;
1640 wsum += abs( yr - yl ) ;
1641 }
1642 else {
1643 specout[ii] += spec[j] * abs( zr - yl ) ;
1644 wsum += abs( zr - yl ) ;
1645 ichan = j ;
1646 break ;
1647 }
1648 }
1649 else {
1650 ichan = j - 1 ;
1651 break ;
1652 }
1653 }
1654 specout[ii] /= wsum ;
1655 wsum = 0.0 ;
1656 }
1657 }
1658 //specout = spec ;
1659 //ofs3.close() ;
1660 }
1661 else {
1662 specout = spec ;
1663 }
1664
1665 return specout ;
1666}
1667
1668int NROFITSDataset::getIndex( int irow )
1669{
1670 NRODataRecord *record = getRecord( irow ) ;
1671 string str = record->ARRYT ;
1672 string::size_type pos = str.find( " " ) ;
1673 if ( pos != string::npos )
1674 str = str.substr( 0, pos ) ;
1675 int index = -1 ;
1676 for ( int i = 0 ; i < ARYNM ; i++ ) {
1677 if ( str.compare( 0, 3, ARYTP[i] ) == 0 ) {
1678 index = i ;
1679 break ;
1680 }
1681 }
1682 return index ;
1683}
1684
1685double NROFITSDataset::radRA( string ra )
1686{
1687 int pos1 = ra.find( ':' ) ;
1688 int pos2 ;
1689 string ch = ra.substr( 0, pos1 ) ;
1690 //cout << "ch = \'" << ch << "\'" << endl ;
1691 pos2 = pos1 + 1 ;
1692 pos1 = ra.find( ':', pos2 ) ;
1693 string cm = ra.substr( pos2, pos1 - pos2 ) ;
1694 //cout << "cm = \'" << cm << "\'" << endl ;
1695 pos2 = pos1 + 1 ;
1696 pos1 = ra.size() ;
1697 string cs = ra.substr( pos2, pos1 - pos2 ) ;
1698 //cout << "cs = \'" << cs << "\'" << endl ;
1699 double h ;
1700 if ( ra[0] != '-' )
1701 h = atof( ch.c_str() ) + atof( cm.c_str() ) / 60.0 + atof( cs.c_str() ) / 3600.0 ;
1702 else
1703 h = atof( ch.c_str() ) - atof( cm.c_str() ) / 60.0 - atof( cs.c_str() ) / 3600.0 ;
1704 double rra = h * M_PI / 12.0 ;
1705 return rra ;
1706}
1707
1708double NROFITSDataset::radDEC( string dec )
1709{
1710 int pos1 = dec.find( ':' ) ;
1711 int pos2 ;
1712 string cd = dec.substr( 0, pos1 ) ;
1713 //cout << "cd = \'" << cd << "\'" << endl ;
1714 pos2 = pos1 + 1 ;
1715 pos1 = dec.find( ':', pos2 ) ;
1716 string cm = dec.substr( pos2, pos1 - pos2 ) ;
1717 //cout << "cm = \'" << cm << "\'" << endl ;
1718 pos2 = pos1 + 1 ;
1719 pos1 = dec.size() ;
1720 string cs = dec.substr( pos2, pos1 - pos2 ) ;
1721 //cout << "cs = \'" << cs << "\'" << endl ;
1722 double h ;
1723 if ( dec[0] != '-' )
1724 h = atof( cd.c_str() ) + atof( cm.c_str() ) / 60.0 + atof( cs.c_str() ) / 3600.0 ;
1725 else
1726 h = atof( cd.c_str() ) - atof( cm.c_str() ) / 60.0 - atof( cs.c_str() ) / 3600.0 ;
1727 double rdec = h * M_PI / 180.0 ;
1728 return rdec ;
1729}
1730
1731void NROFITSDataset::getField()
1732{
1733 for ( int i = 0 ; i < numField_ ; i++ ) {
1734 char key1[9] ;
1735 char key2[9] ;
1736 char key3[9] ;
1737 if ( i < 9 ) {
1738 sprintf( key1, "TFORM%d ", i+1 ) ;
1739 sprintf( key2, "TTYPE%d ", i+1 ) ;
1740 sprintf( key3, "TUNIT%d ", i+1 ) ;
1741 //cout << "key1 = " << key1 << ", key2 = " << key2 << ", key3 = " << key3 << endl ;
1742 }
1743 else if ( i < 99 ) {
1744 sprintf( key1, "TFORM%2d ", i+1 ) ;
1745 sprintf( key2, "TTYPE%2d ", i+1 ) ;
1746 sprintf( key3, "TUNIT%2d ", i+1 ) ;
1747 //cout << "key1 = " << key1 << ", key2 = " << key2 << ", key3 = " << key3 << endl ;
1748 }
1749 else {
1750 sprintf( key1, "TFORM%3d", i+1 ) ;
1751 sprintf( key2, "TTYPE%3d", i+1 ) ;
1752 sprintf( key3, "TUNIT%3d", i+1 ) ;
1753 //cout << "key1 = " << key1 << ", key2 = " << key2 << ", key3 = " << key3 << endl ;
1754 }
1755 //char tmp[9] ;
1756 string tmp ;
1757 //strcpy( tmp, " " ) ;
1758 if ( readHeader( tmp, key1 ) != 0 ) {
1759 cerr << "Error while reading field keyword for scan header." << endl ;
1760 return ;
1761 }
1762 //forms_[i] = string( tmp ) ;
1763 forms_[i] = tmp ;
1764 string::size_type spos = forms_[i].find( " " ) ;
1765 if ( spos != string::npos )
1766 forms_[i] = forms_[i].substr( 0, spos ) ;
1767 //strcpy( tmp, " " ) ;
1768 if ( readHeader( tmp, key2 ) != 0 ) {
1769 cerr << "Error while reading field type for scan header." << endl ;
1770 return ;
1771 }
1772 //names_[i] = string( tmp ) ;
1773 names_[i] = tmp ;
1774 spos = names_[i].find( " " ) ;
1775 if ( spos != string::npos )
1776 names_[i] = names_[i].substr( 0, spos ) ;
1777 //strcpy( tmp, " " ) ;
1778 if ( forms_[i].find( "A" ) != string::npos ) {
1779 //cout << "skip to get unit: name = " << forms_[i] << endl ;
1780 //strcpy( tmp, "none " ) ;
1781 tmp = "none" ;
1782 }
1783 else {
1784 //cout << "get unit: name = " << forms_[i] << endl ;
1785 if ( readHeader( tmp, key3 ) != 0 ) {
1786 //strcpy( tmp, "none " ) ;
1787 tmp = "none" ;
1788 }
1789 }
1790 //units_[i] = string( tmp ) ;
1791 units_[i] = tmp ;
1792 spos = units_[i].find( " " ) ;
1793 if ( spos != string::npos )
1794 units_[i] = units_[i].substr( 0, spos ) ;
1795 //cout << "i = " << i << ": name=" << forms_[i] << " type=" << names_[i] << " unit=" << units_[i] << endl ;
1796 }
1797}
1798
1799void NROFITSDataset::fillARYTP()
1800{
1801 string arry ;
1802 int count = 0 ;
1803 string arry1 ;
1804 string arry2 ;
1805 string arry3 ;
1806 string arry4 ;
1807 if ( readHeader( arry, "ARRY1" ) == 0 )
1808 arry1 = arry ;
1809 else
1810 arry1 = "00000000000000000000" ;
1811 for ( int i = 0 ; i < 20 ; i++ ) {
1812 if ( arry1[i] == '1' ) {
1813 char arytp[4] ;
1814 sprintf( arytp, "H%d", i+1 ) ;
1815 ARYTP[count++] = string( arytp ) ;
1816 //cout << "ARYTP[" << count-1 << "] = " << ARYTP[count-1] << endl ;
1817 }
1818 }
1819 if ( readHeader( arry, "ARRY2" ) == 0 )
1820 arry2 = arry ;
1821 else
1822 arry2 = "00000000000000000000" ;
1823 for ( int i = 0 ; i < 20 ; i++ ) {
1824 if ( arry2[i] == '1' ) {
1825 if ( i < 10 ) {
1826 char arytp[4] ;
1827 sprintf( arytp, "W%d", i+1 ) ;
1828 ARYTP[count++] = string( arytp ) ;
1829 //cout << "ARYTP[" << count-1 << "] = " << ARYTP[count-1] << endl ;
1830 }
1831 else if ( i < 15 ) {
1832 char arytp[4] ;
1833 sprintf( arytp, "U%d", i-9 ) ;
1834 ARYTP[count++] = string( arytp ) ;
1835 //cout << "ARYTP[" << count-1 << "] = " << ARYTP[count-1] << endl ;
1836 }
1837 else {
1838 char arytp[4] ;
1839 sprintf( arytp, "X%d", i-14 ) ;
1840 ARYTP[count++] = string( arytp ) ;
1841 //cout << "ARYTP[" << count-1 << "] = " << ARYTP[count-1] << endl ;
1842 }
1843 }
1844 }
1845 if ( readHeader( arry, "ARRY3" ) == 0 )
1846 arry3 = arry ;
1847 else
1848 arry3 = "00000000000000000000" ;
1849 for ( int i = 0 ; i < 20 ; i++ ) {
1850 if ( arry3[i] == '1' ) {
1851 char arytp[4] ;
1852 sprintf( arytp, "A%d", i+1 ) ;
1853 ARYTP[count++] = string( arytp ) ;
1854 //cout << "ARYTP[" << count-1 << "] = " << ARYTP[count-1] << endl ;
1855 }
1856 }
1857 if ( readHeader( arry, "ARRY4" ) == 0 )
1858 arry4 = arry ;
1859 else
1860 arry4 = "00000000000000000000" ;
1861 for ( int i = 0 ; i < 20 ; i++ ) {
1862 if ( arry4[i] == '1' ) {
1863 char arytp[4] ;
1864 sprintf( arytp, "A%d", i+21 ) ;
1865 ARYTP[count++] = string( arytp ) ;
1866 //cout << "ARYTP[" << count-1 << "] = " << ARYTP[count-1] << endl ;
1867 }
1868 }
1869}
1870
1871int NROFITSDataset::readARRY()
1872{
1873 LogIO os( LogOrigin( "NROFITSDataset", "readARRY()", WHERE ) ) ;
1874
1875 string arry1 ;
1876 string arry2 ;
1877 string arry3 ;
1878 string arry4 ;
1879 int status = readHeader( arry1, "ARRY1" ) ;
1880 if ( status ) {
1881 os << LogIO::SEVERE << "Error while reading ARRY1" << LogIO::POST ;
1882 return status ;
1883 }
1884 status = readHeader( arry2, "ARRY2" ) ;
1885 if ( status ) {
1886 os << LogIO::SEVERE << "Error while reading ARRY2" << LogIO::POST ;
1887 return status ;
1888 }
1889 status = readHeader( arry3, "ARRY3" ) ;
1890 if ( status ) {
1891 os << LogIO::SEVERE << "Error while reading ARRY3" << LogIO::POST ;
1892 return status ;
1893 }
1894 status = readHeader( arry4, "ARRY4" ) ;
1895 if ( status ) {
1896 os << LogIO::SEVERE << "Error while reading ARRY4" << LogIO::POST ;
1897 return status ;
1898 }
1899 int index = 0 ;
1900 for ( int i = 0 ; i < 20 ; i++ ) {
1901 // AOSH
1902 if ( arry1[i] == '1' )
1903 ARRY[index] = 1 ;
1904 else
1905 ARRY[index] = 0 ;
1906 // AOSW, AOSU, FX
1907 if ( arry2[i] == '1' )
1908 ARRY[index+20] = 1 ;
1909 else
1910 ARRY[index+20] = 0 ;
1911 // AC45 (1-20)
1912 if ( arry3[i] == '1' )
1913 ARRY[index+40] = 1 ;
1914 else
1915 ARRY[index+40] = 0 ;
1916 // AC45 (21-35)
1917 if ( i < 15 ) {
1918 if ( arry4[i] == '1' )
1919 ARRY[index+60] = 1 ;
1920 else
1921 ARRY[index+60] = 0 ;
1922 }
1923 index++ ;
1924 }
1925 return status ;
1926}
1927
1928void NROFITSDataset::findData()
1929{
1930 LogIO os( LogOrigin( "NROFITSDataset", "findData()", WHERE ) ) ;
1931
1932 // skip header
1933 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
1934
1935 // get offset
1936 int offset = getOffset( "ARRYT" ) ;
1937 if ( offset == -1 ) {
1938 //cerr << "Error, ARRYT is not found in the name list." << endl ;
1939 return ;
1940 }
1941 //cout << "offset for ARRYT is " << offset << " bytes." << endl ;
1942 fseek( fp_, offset, SEEK_CUR ) ;
1943 int count = 0 ;
1944 int index = 0 ;
1945 while ( count < ARYNM && index < rowNum_ ) {
1946 char ctmp[5] ;
1947 fread( ctmp, 1, 4, fp_ ) ;
1948 ctmp[4] = '\0' ;
1949 //cout << "ctmp = " << ctmp << endl ;
1950 for ( int i = 0 ; i < ARYNM ; i++ ) {
1951 if ( arrayid_[i] != -1 )
1952 continue ;
1953 else if ( strncmp( ctmp, ARYTP[i].c_str(), ARYTP[i].size() ) == 0 ) {
1954 //cout << "matched: i = " << i << ", ARYTP = " << ARYTP[i] << ", ctmp = " << ctmp << endl ;
1955 arrayid_[i] = index ;
1956 count++ ;
1957 }
1958 }
1959 fseek( fp_, scanLen_-4, SEEK_CUR ) ;
1960 index++ ;
1961 }
1962
1963 if ( count != ARYNM ) {
1964 os << LogIO::WARN << "NROFITSDataset::findData() failed to find rows for " ;
1965 for ( int i = 0 ; i < ARYNM ; i++ ) {
1966 if ( arrayid_[i] == -1 ) {
1967 os << LogIO::WARN << ARYTP[i] << " " ;
1968 }
1969 }
1970 os.post() ;
1971 }
1972
1973// for ( int i = 0 ; i < ARYNM ; i++ )
1974// //cout << "arrayid_[" << i << "] = " << arrayid_[i] << endl ;
1975 }
1976
1977int NROFITSDataset::getOffset( char *name )
1978{
1979 int offset = 0 ;
1980 string sname( name ) ;
1981 bool found = false ;
1982 for ( int i = 0 ; i < numField_ ; i++ ) {
1983 // escape if name is found
1984 //cout << "names_[" << i << "] = " << names_[i] << " sname = " << sname << endl ;
1985 if ( names_[i] == sname ) {
1986 found = true ;
1987 break ;
1988 }
1989
1990 // form analysis
1991 string substr1 = forms_[i].substr( 0, forms_[i].size()-1 ) ;
1992 string substr2 = forms_[i].substr( forms_[i].size()-1, 1 ) ;
1993 //cout << "substr1 = " << substr1 << ", substr2 = " << substr2 << endl ;
1994 int o1 = atoi( substr1.c_str() ) ;
1995 int o2 = 0 ;
1996 if ( substr2 == "A" )
1997 o2 = sizeof(char) ;
1998 else if ( substr2 == "J" )
1999 o2 = sizeof(int) ;
2000 else if ( substr2 == "F" )
2001 o2 = sizeof(float) ;
2002 else if ( substr2 == "D" )
2003 o2 = sizeof(double) ;
2004 //cout << "o1 = " << o1 << ", o2 = " << o2 << endl ;
2005 offset += o1 * o2 ;
2006 }
2007
2008 if ( !found )
2009 offset = -1 ;
2010
2011 return offset ;
2012}
2013
2014int NROFITSDataset::getPolarizationNum()
2015{
2016 int npol = 0 ;
2017
2018 vector<char> type( 2 ) ;
2019 type[0] = 'C' ;
2020 type[1] = 'L' ;
2021 vector<double> crot ;
2022 vector<double> lagl ;
2023
2024 for ( int i = 0 ; i < ARYNM ; i++ ) {
2025 if ( POLTP[i][0] == type[0] ) {
2026 // circular polarization
2027 if( count( crot.begin(), crot.end(), POLDR[i] ) != 0 ) {
2028 crot.push_back( POLDR[i] ) ;
2029 npol++ ;
2030 }
2031 }
2032 else if ( POLTP[i][0] == type[1] ) {
2033 // linear polarization
2034 if ( count( lagl.begin(), lagl.end(), POLAN[i] ) != 0 ) {
2035 lagl.push_back( POLAN[i] ) ;
2036 npol++ ;
2037 }
2038 }
2039 }
2040
2041 if ( npol == 0 )
2042 npol = 1 ;
2043
2044
2045 return npol ;
2046}
2047
2048int NROFITSDataset::readHeader( string &v, char *name )
2049{
2050 //
2051 // Read 'name' attribute defined as char from the FITS Header
2052 //
2053 int status = 0 ;
2054
2055 char buf[81] ;
2056 strcpy( buf, " " ) ;
2057 fseek( fp_, 0, SEEK_SET ) ;
2058 int count = 0 ;
2059 while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
2060 fread( buf, 1, 80, fp_ ) ;
2061 buf[80] = '\0' ;
2062 count++ ;
2063 }
2064 if ( strncmp( buf, "END", 3 ) == 0 ) {
2065 //cerr << "NROFITSDataset::readHeader() keyword " << name << " not found." << endl ;
2066 //cerr << "count = " << count << endl ;
2067 status = -1 ;
2068 return status ;
2069 }
2070 string str( buf ) ;
2071 int pos1 = str.find( '\'' ) + 1 ;
2072 int pos2 = str.find( '\'', pos1 ) ;
2073 unsigned int clen = pos2 - pos1 ;
2074 //cout << "string: " << str << endl ;
2075 //cout << "value: " << str.substr( pos1, clen ).c_str() << endl ;
2076 //cout << "clen = " << clen << endl ;
2077 v = str.substr( pos1, clen ) ;
2078 //cout << "v = \'" << v << "\'" << endl ;
2079
2080 return status ;
2081}
2082
2083int NROFITSDataset::readHeader( int &v, char *name, int b )
2084{
2085 //
2086 // Read 'name' attribute defined as int from the FITS Header
2087 //
2088 int status = 0 ;
2089
2090 char buf[81] ;
2091 strcpy( buf, " " ) ;
2092 fseek( fp_, 0, SEEK_SET ) ;
2093 while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
2094 fread( buf, 1, 80, fp_ ) ;
2095 buf[80] = '\0' ;
2096 //char bufo[9] ;
2097 //strncpy( bufo, buf, 8 ) ;
2098 //bufo[8] = '\0' ;
2099 //cout << "header: " << bufo << endl ;
2100 }
2101 if ( strncmp( buf, "END", 3 ) == 0 ) {
2102 //cerr << "NROFITSDataset::readHeader() keyword " << name << " not found." << endl ;
2103 status = -1 ;
2104 return status ;
2105 }
2106 string str( buf ) ;
2107 int pos1 = str.find( '=' ) + 1 ;
2108 int pos2 = str.find( '/' ) ;
2109 //cout << "string: " << str << endl ;
2110 //cout << "value: " << str.substr( pos1, pos2 - pos1 ).c_str() << endl ;
2111 v = atoi( str.substr( pos1, pos2 - pos1 ).c_str() ) ;
2112 //cout << "v = " << v << endl ;
2113
2114 //cout << "NROFITSDataset::readHeader() end to read" << endl ;
2115 return status ;
2116}
2117
2118
2119int NROFITSDataset::readHeader( float &v, char *name, int b )
2120{
2121 //
2122 // Read 'name' attribute defined as float from the FITS Header
2123 //
2124 int status = 0 ;
2125
2126 char buf[81] ;
2127 strcpy( buf, " " ) ;
2128 fseek( fp_, 0, SEEK_SET ) ;
2129 while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
2130 fread( buf, 1, 80, fp_ ) ;
2131 buf[80] = '\0' ;
2132 //char bufo[9] ;
2133 //strncpy( bufo, buf, 8 ) ;
2134 //bufo[8] = '\0' ;
2135 //cout << "header: " << bufo << endl ;
2136 }
2137 if ( strncmp( buf, "END", 3 ) == 0 ) {
2138 //cerr << "NROFITSDataset::readHeader() keyword " << name << " not found." << endl ;
2139 status = -1 ;
2140 return status ;
2141 }
2142 string str( buf ) ;
2143 int pos1 = str.find( '=' ) + 1 ;
2144 int pos2 = str.find( '/' ) ;
2145 //cout << "string: " << str << endl ;
2146 //cout << "value: " << str.substr( pos1, pos2 - pos1 ).c_str() << endl ;
2147 v = atof( str.substr( pos1, pos2 - pos1 ).c_str() ) ;
2148 //cout << "v = " << v << endl ;
2149
2150 return status ;
2151}
2152
2153int NROFITSDataset::readHeader( double &v, char *name, int b )
2154{
2155 //
2156 // Read 'name' attribute defined as double from the FITS Header
2157 //
2158 int status = 0 ;
2159
2160 char buf[81] ;
2161 strcpy( buf, " " ) ;
2162 fseek( fp_, 0, SEEK_SET ) ;
2163 while ( strncmp( buf, name, strlen(name) ) != 0 && strncmp( buf, "END", 3 ) != 0 ) {
2164 fread( buf, 1, 80, fp_ ) ;
2165 buf[80] = '\0' ;
2166 char bufo[9] ;
2167 strncpy( bufo, buf, 8 ) ;
2168 bufo[8] = '\0' ;
2169 //cout << "header: \'" << bufo << "\' bufo = \'" << bufo << "\' ";
2170 //cout << strncmp( buf, name, strlen(name) ) << endl ;
2171 }
2172 if ( strncmp( buf, "END", 3 ) == 0 ) {
2173 //cerr << "NROFITSDataset::readHeader() keyword " << name << " not found." << endl ;
2174 status = -1 ;
2175 return status ;
2176 }
2177 string str( buf ) ;
2178 int pos1 = str.find( '=' ) + 1 ;
2179 int pos2 = str.find( '/' ) ;
2180 //cout << "string: " << str << endl ;
2181 //cout << "value: " << str.substr( pos1, pos2 - pos1 ).c_str() << endl ;
2182 v = atof( str.substr( pos1, pos2 - pos1 ).c_str() ) ;
2183 //cout << "v = " << v << endl ;
2184
2185 return status ;
2186}
2187
2188int NROFITSDataset::readTable( char *v, char *name )
2189{
2190 //
2191 // Read 'name' attribute defined as char from the 0-th row
2192 // of the FITS Scan Record
2193 //
2194 int status = readTable( v, name, (int)strlen( v ), 0 ) ;
2195
2196 return status ;
2197}
2198
2199int NROFITSDataset::readTable( char *v, char *name, int clen, int idx )
2200{
2201 //
2202 // Read 'name' attribute defined as char from the idx-th row
2203 // of the FITS Scan Record
2204 //
2205 int status = 0 ;
2206
2207 // skip header
2208 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2209
2210 // find offset
2211 int offset = getOffset( name ) ;
2212 if ( offset == -1 ) {
2213 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2214 return -1 ;
2215 }
2216 offset += idx * scanLen_ ;
2217 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2218 fseek( fp_, offset, SEEK_CUR ) ;
2219
2220 // get length of char
2221 int index = -1 ;
2222 for ( int i = 0 ; i < numField_ ; i++ ) {
2223 if ( names_[i] == name ) {
2224 index = i ;
2225 break ;
2226 }
2227 }
2228 string substr = forms_[index].substr( 0, forms_[index].size()-1 ) ;
2229 int xsize = atoi( substr.c_str() ) ;
2230 //cout << "xsize = " << xsize << endl ;
2231
2232 // read data
2233 if ( xsize < clen ) {
2234 fread( v, 1, xsize, fp_ ) ;
2235 //v[xsize] = '\0' ;
2236 }
2237 else {
2238 fread( v, 1, clen-1, fp_ ) ;
2239 //v[clen-1] = '\0' ;
2240 }
2241
2242 return status ;
2243}
2244
2245int NROFITSDataset::readTable( int &v, char *name, int b )
2246{
2247 //
2248 // Read 'name' attribute defined as int from the 0-th row
2249 // of the FITS Scan Record
2250 //
2251 int status = readTable( v, name, b, 0 ) ;
2252
2253 return status ;
2254}
2255
2256int NROFITSDataset::readTable( int &v, char *name, int b, int idx )
2257{
2258 //
2259 // Read 'name' attribute defined as int from the idx-th row
2260 // of the FITS Scan Record
2261 //
2262 int status = 0 ;
2263 // skip header
2264 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2265
2266 // find offset
2267 int offset = getOffset( name ) ;
2268 if ( offset == -1 ) {
2269 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2270 return -1 ;
2271 }
2272 offset += idx * scanLen_ ;
2273 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2274 fseek( fp_, offset, SEEK_CUR ) ;
2275
2276 // read data
2277 fread( &v, sizeof(int), 1, fp_ ) ;
2278 if ( b == 0 )
2279 convertEndian( v ) ;
2280
2281 return status ;
2282}
2283
2284int NROFITSDataset::readTable( float &v, char *name, int b )
2285{
2286 //
2287 // Read 'name' attribute defined as float from the 0-th row
2288 // of the FITS Scan Record
2289 //
2290 int status = readTable( v, name, b, 0 ) ;
2291
2292 return status ;
2293}
2294
2295int NROFITSDataset::readTable( float &v, char *name, int b, int idx )
2296{
2297 //
2298 // Read 'name' attribute defined as float from the idx-th row
2299 // of the FITS Scan Record
2300 //
2301 int status = 0 ;
2302
2303 // skip header
2304 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2305
2306 // find offset
2307 int offset = getOffset( name ) ;
2308 if ( offset == -1 ) {
2309 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2310 return -1 ;
2311 }
2312 offset += idx * scanLen_ ;
2313 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2314 fseek( fp_, offset, SEEK_CUR ) ;
2315
2316 // read data
2317 fread( &v, sizeof(float), 1, fp_ ) ;
2318 if ( b == 0 )
2319 convertEndian( v ) ;
2320
2321 return status ;
2322}
2323
2324int NROFITSDataset::readTable( double &v, char *name, int b )
2325{
2326 //
2327 // Read 'name' attribute defined as double from the 0-th row
2328 // of the FITS Scan Record
2329 //
2330 int status = readTable( v, name, b, 0 ) ;
2331
2332 return status ;
2333}
2334
2335int NROFITSDataset::readTable( double &v, char *name, int b, int idx )
2336{
2337 //
2338 // Read 'name' attribute defined as double from the idx-th row
2339 // of the FITS Scan Record
2340 //
2341 int status = 0 ;
2342
2343 // skip header
2344 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2345
2346 // find offset
2347 int offset = getOffset( name ) ;
2348 if ( offset == -1 ) {
2349 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2350 return -1 ;
2351 }
2352 offset += idx * scanLen_ ;
2353 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2354 fseek( fp_, offset, SEEK_CUR ) ;
2355
2356 // read data
2357 fread( &v, sizeof(double), 1, fp_ ) ;
2358 if ( b == 0 )
2359 convertEndian( v ) ;
2360
2361 return status ;
2362}
2363
2364int NROFITSDataset::readTable( vector<char *> &v, char *name, int idx )
2365{
2366 //
2367 // Read 'name' attribute defined as char array from the FITS Scan Record
2368 //
2369 int status = 0 ;
2370
2371 // skip header
2372 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2373
2374 // find offset
2375 int offset = getOffset( name ) ;
2376 if ( offset == -1 ) {
2377 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2378 return -1 ;
2379 }
2380 offset += idx * scanLen_ ;
2381 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2382 fseek( fp_, offset, SEEK_CUR ) ;
2383
2384 // get length of char
2385 int index = -1 ;
2386 for ( int i = 0 ; i < numField_ ; i++ ) {
2387 if ( names_[i] == name ) {
2388 index = i ;
2389 break ;
2390 }
2391 }
2392 string substr = forms_[index].substr( 0, forms_[index].size()-1 ) ;
2393 int xsize = atoi( substr.c_str() ) ;
2394 //cout << "xsize = " << xsize << endl ;
2395
2396 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2397 int clen = strlen( v[i] ) ;
2398 if ( clen > xsize ) {
2399 fread( v[i], 1, xsize, fp_ ) ;
2400 //v[i][xsize] = '\0' ;
2401 }
2402 else {
2403 fread( v[i], 1, clen, fp_ ) ;
2404 //v[i][clen-1] = '\0' ;
2405 }
2406 //cout << "v[" << i << "] = " << v[i] << endl ;
2407 }
2408
2409 return status ;
2410}
2411
2412int NROFITSDataset::readTable( vector<int> &v, char *name, int b, int idx )
2413{
2414 //
2415 // Read 'name' attribute defined as int array from the FITS Scan Record
2416 //
2417 int status = 0 ;
2418
2419 // skip header
2420 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2421
2422 // find offset
2423 int offset = getOffset( name ) ;
2424 if ( offset == -1 ) {
2425 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2426 return -1 ;
2427 }
2428 offset += idx * scanLen_ ;
2429 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2430 fseek( fp_, offset, SEEK_CUR ) ;
2431
2432 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2433 fread( &v[i], 1, sizeof(int), fp_ ) ;
2434 if ( b == 0 )
2435 convertEndian( v[i] ) ;
2436 //cout << "v[" << i << "] = " << v[i] << endl ;
2437 }
2438
2439 return status ;
2440}
2441
2442int NROFITSDataset::readTable( vector<float> &v, char *name, int b, int idx )
2443{
2444 //
2445 // Read 'name' attribute defined as float array from the FITS Scan Record
2446 //
2447 int status = 0 ;
2448
2449 // skip header
2450 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2451
2452 // find offset
2453 int offset = getOffset( name ) ;
2454 if ( offset == -1 ) {
2455 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2456 return -1 ;
2457 }
2458 offset += idx * scanLen_ ;
2459 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2460 fseek( fp_, offset, SEEK_CUR ) ;
2461
2462 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2463 fread( &v[i], 1, sizeof(float), fp_ ) ;
2464 if ( b == 0 )
2465 convertEndian( v[i] ) ;
2466 //cout << "v[" << i << "] = " << v[i] << endl ;
2467 }
2468
2469 return status ;
2470}
2471
2472int NROFITSDataset::readTable( vector<double> &v, char *name, int b, int idx )
2473{
2474 //
2475 // Read 'name' attribute defined as double array from the FITS Scan Record
2476 //
2477 int status = 0 ;
2478
2479 // skip header
2480 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2481
2482 // find offset
2483 int offset = getOffset( name ) ;
2484 if ( offset == -1 ) {
2485 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2486 return -1 ;
2487 }
2488 offset += idx * scanLen_ ;
2489 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2490 fseek( fp_, offset, SEEK_CUR ) ;
2491
2492 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2493 fread( &v[i], 1, sizeof(double), fp_ ) ;
2494 if ( b == 0 )
2495 convertEndian( v[i] ) ;
2496 //cout << "v[" << i << "] = " << v[i] << endl ;
2497 }
2498
2499 return status ;
2500}
2501
2502int NROFITSDataset::readColumn( vector<string> &v, char *name )
2503{
2504 //
2505 // Read 0-th column of ARRYTP-dependent 'name' attributes
2506 // defined as char array from the FITS Scan Record
2507 //
2508 int status = readColumn( v, name, 0 ) ;
2509
2510 return status ;
2511}
2512
2513int NROFITSDataset::readColumn( vector<string> &v, char *name, int idx )
2514{
2515 //
2516 // Read idx-th column of ARRYTP-dependent 'name' attributes
2517 // defined as char array from the FITS Scan Record
2518 //
2519 int status = 0 ;
2520
2521 // skip header
2522 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2523
2524 // find offset
2525 int offset = getOffset( name ) ;
2526 if ( offset == -1 ) {
2527 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2528 return -1 ;
2529 }
2530 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2531 fseek( fp_, offset, SEEK_CUR ) ;
2532
2533 // get length of char
2534 int index = -1 ;
2535 for ( int i = 0 ; i < numField_ ; i++ ) {
2536 if ( names_[i] == name ) {
2537 index = i ;
2538 break ;
2539 }
2540 }
2541 string substr = forms_[index].substr( 0, forms_[index].size()-1 ) ;
2542 int xsize = atoi( substr.c_str() ) ;
2543 //cout << "xsize = " << xsize << endl ;
2544
2545 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2546 offset = scanLen_ * arrayid_[i] + xsize * idx ;
2547 fseek( fp_, offset, SEEK_CUR ) ;
2548// int clen = (int)strlen( v[i] ) ;
2549// if ( clen > xsize ) {
2550// fread( v[i], 1, xsize, fp_ ) ;
2551// //v[i][xsize] = '\0' ;
2552// }
2553// else {
2554// fread( v[i], 1, clen-1, fp_ ) ;
2555// //v[i][clen-1] = '\0' ;
2556// }
2557 char c[xsize+1] ;
2558 fread( c, 1, xsize, fp_ ) ;
2559 c[xsize] = '\0' ;
2560 v[i] = string( c ) ;
2561 //cout << "v[" << i << "] = \'" << v[i] << "\'" << endl ;
2562 fseek( fp_, -xsize-offset, SEEK_CUR ) ;
2563 }
2564
2565 return status ;
2566}
2567
2568int NROFITSDataset::readColumn( vector<int> &v, char *name, int b )
2569{
2570 //
2571 // Read 0-th column of ARRYTP-dependent 'name' attributes
2572 // defined as int array from the FITS Scan Record
2573 //
2574 int status = readColumn( v, name, b, 0 ) ;
2575
2576 return status ;
2577}
2578
2579int NROFITSDataset::readColumn( vector<int> &v, char *name, int b, int idx )
2580{
2581 //
2582 // Read idx-th column of ARRYTP-dependent 'name' attributes
2583 // defined as int array from the FITS Scan Record
2584 //
2585 int status = 0 ;
2586
2587 // skip header
2588 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2589
2590 // find offset
2591 int offset = getOffset( name ) ;
2592 if ( offset == -1 ) {
2593 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2594 return -1 ;
2595 }
2596 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2597 fseek( fp_, offset, SEEK_CUR ) ;
2598
2599 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2600 offset = scanLen_ * arrayid_[i] + sizeof(int) * idx ;
2601 fseek( fp_, offset, SEEK_CUR ) ;
2602 fread( &v[i], 1, sizeof(int), fp_ ) ;
2603 if ( b == 0 )
2604 convertEndian( v[i] ) ;
2605 //cout << "v[" << i << "] = " << v[i] << endl ;
2606 fseek( fp_, -sizeof(int)-offset, SEEK_CUR ) ;
2607 }
2608
2609 return status ;
2610}
2611
2612int NROFITSDataset::readColumn( vector<float> &v, char *name, int b )
2613{
2614 //
2615 // Read 0-th column of ARRYTP-dependent 'name' attributes
2616 // defined as float array from the FITS Scan Record
2617 //
2618 int status = readColumn( v, name, b, 0 ) ;
2619
2620 return status ;
2621}
2622
2623int NROFITSDataset::readColumn( vector<float> &v, char *name, int b, int idx )
2624{
2625 //
2626 // Read idx-th column of ARRYTP-dependent 'name' attributes
2627 // defined as float array from the FITS Scan Record
2628 //
2629 int status = 0 ;
2630
2631 // skip header
2632 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2633
2634 // find offset
2635 int offset = getOffset( name ) ;
2636 if ( offset == -1 ) {
2637 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2638 return -1 ;
2639 }
2640 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2641 fseek( fp_, offset, SEEK_CUR ) ;
2642
2643 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2644 offset = scanLen_ * arrayid_[i] + sizeof(float) * idx ;
2645 fseek( fp_, offset, SEEK_CUR ) ;
2646 fread( &v[i], 1, sizeof(float), fp_ ) ;
2647 if ( b == 0 )
2648 convertEndian( v[i] ) ;
2649 //cout << "v[" << i << "] = " << v[i] << endl ;
2650 fseek( fp_, -sizeof(float)-offset, SEEK_CUR ) ;
2651 }
2652
2653 return status ;
2654}
2655
2656int NROFITSDataset::readColumn( vector<double> &v, char *name, int b )
2657{
2658 //
2659 // Read 0-th column of ARRYTP-dependent 'name' attributes
2660 // defined as double array from the FITS Scan Record
2661 //
2662 int status = readColumn( v, name, b, 0 ) ;
2663
2664 return status ;
2665}
2666
2667int NROFITSDataset::readColumn( vector<double> &v, char *name, int b, int idx )
2668{
2669 //
2670 // Read idx-th column of ARRYTP-dependent 'name' attributes
2671 // defined as double array from the FITS Scan Record
2672 //
2673 int status = 0 ;
2674
2675 // skip header
2676 fseek( fp_, FITS_HEADER_SIZE, SEEK_SET ) ;
2677
2678 // find offset
2679 int offset = getOffset( name ) ;
2680 if ( offset == -1 ) {
2681 //cerr << "Error, " << name << " is not found in the name list." << endl ;
2682 return -1 ;
2683 }
2684 //cout << "offset for " << name << " is " << offset << " bytes." << endl ;
2685 fseek( fp_, offset, SEEK_CUR ) ;
2686
2687 for ( unsigned int i = 0 ; i < v.size() ; i++ ) {
2688 offset = scanLen_ * arrayid_[i] + sizeof(double) * idx ;
2689 fseek( fp_, offset, SEEK_CUR ) ;
2690 fread( &v[i], 1, sizeof(double), fp_ ) ;
2691 if ( b == 0 )
2692 convertEndian( v[i] ) ;
2693 //cout << "offset = " << offset << ", v[" << i << "] = " << v[i] << endl ;
2694 fseek( fp_, -sizeof(double)-offset, SEEK_CUR ) ;
2695 }
2696
2697// //cout << "v: " << endl ;
2698// for ( vector<double>::iterator i = v.begin() ; i != v.end() ; i++ )
2699// //cout << *i << " " ;
2700// //cout << endl ;
2701
2702 return status ;
2703}
2704
2705uInt NROFITSDataset::getArrayId( string type )
2706{
2707 uInt ib = 99 ;
2708 for ( uInt i = 0 ; i < arrayid_.size() ; i++ ) {
2709 uInt len = ARYTP[i].size() ;
2710 if ( type.compare( 0, len, ARYTP[i], 0, len ) == 0 ) {
2711 ib = i ;
2712 break ;
2713 }
2714 }
2715 return ib ;
2716}
Note: See TracBrowser for help on using the repository browser.