source: trunk/external-alma/atnf/PKSIO/NROFITSDataset.cc@ 2436

Last change on this file since 2436 was 2436, checked in by Takeshi Nakazato, 13 years ago

New Development: No

JIRA Issue: No

Ready for Test: Yes

Interface Changes: No

What Interface Changed: Please list interface changes

Test Programs: List test programs

Put in Release Notes: Yes/No

Module(s): Module Names change impacts.

Description: Describe your changes here...

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