Ignore:
Timestamp:
Aug 24, 2000, 12:18:02 PM (25 years ago)
Author:
ansari
Message:

structuration a la ppersist+ convention sur hdu

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc

    r1047 r1136  
    1515{
    1616  dobj_= new SphereHEALPix<T>;
    17   ownobj= true;
     17  ownobj_= true;
    1818}
    1919 
     
    2222{
    2323  dobj_= new SphereHEALPix<T>;
    24   ownobj= true;
    25 
    26   ReadF(inputfile,hdunum);
    27   dobj_->SetTemp(true);
     24  ownobj_= true;
     25
     26  Read(inputfile,hdunum);
     27  //  dobj_->SetTemp(true);
    2828}
    2929
     
    3333{
    3434  dobj_= new SphereHEALPix<T>(obj, true);
    35   dobj_->SetTemp(true);
    36   ownobj= true;
     35  //  dobj_->SetTemp(true);
     36  ownobj_= true;
    3737}
    3838
     
    4141{
    4242  dobj_= obj;
    43   ownobj= false;
     43  ownobj_= false;
    4444}
    4545
     
    4747FITS_SphereHEALPix<T>::~FITS_SphereHEALPix()
    4848{
    49   if (ownobj && dobj_) delete dobj_;
     49  if (ownobj_ && dobj_) delete dobj_;
    5050}
    5151
     
    6161  SphereHEALPix<T> * po = dynamic_cast< SphereHEALPix<T> * >(&o);
    6262  if (po == NULL) return;
    63   if (ownobj && dobj_) delete dobj_;
     63  if (ownobj_ && dobj_) delete dobj_;
    6464  dobj_ = po;
    65   ownobj = false;
     65  ownobj_ = false;
    6666}
    6767
    6868
    6969
    70 template <class T>
    71 void FITS_SphereHEALPix<T>::Write(char outputfile[], bool OldFile)
    72 {
    73   WriteF(outputfile, OldFile);
    74 }
    75 template <class T>
    76 void FITS_SphereHEALPix<T>::Read(char inputfile[],int hdunum)
    77 {
    78   ReadF(inputfile,hdunum);
    79 }
    80 
    81 template <class T>
    82 void FITS_SphereHEALPix<T>::WriteToFits()
     70
     71template <class T>
     72void FITS_SphereHEALPix<T>::WriteToFits(FitsOutFile& os)
    8373{
    8474  if(dobj_ == NULL)
     
    130120  Type[1]='\0';
    131121  vector<int> dummy;
    132   //  fn.makeHeaderBntblOnFits(Type, Noms, nPix, 1, dvl, extname, dummy);
    133   makeHeaderBntblOnFits(Type, Noms, nPix, 1, dvl, extname, dummy);
     122  os.makeHeaderBntblOnFits(Type, Noms, nPix, 1, dvl, extname, dummy);
    134123  delete [] Noms[0];
    135124  delete [] Noms;
    136   putColToFits(0,  nPix, dobj_->pixels_.Data());
     125  os.putColToFits(0,  nPix, dobj_->pixels_.Data());
    137126}
    138127
    139128template <class T>
    140 void FITS_SphereHEALPix<T>::ReadFromFits()
     129void FITS_SphereHEALPix<T>::ReadFromFits(FitsInFile& is)
    141130{
    142131  if(dobj_ == NULL)
    143132    {
    144133      dobj_= new SphereHEALPix<T>;
    145       dobj_->SetTemp(true);
    146       ownobj= true;     
     134      ownobj_= true;     
    147135    }
    148136
    149137 
    150138  int nbcols, nbentries;
    151   //  nbcols = fn.NbColsFromFits();
    152   nbcols = NbColsFromFits();
     139  nbcols = is.NbColsFromFits();
    153140  if (nbcols != 1)
    154141    {
    155142      throw IOExc("le fichier fits n'est pas une sphere Healpix");
    156143    }
    157   //   const DVList* dvl = &fn.DVListFromFits();
    158   //  DVList dvl=fn.DVListFromFits();
    159   DVList dvl=DVListFromFits();
    160   //  dvl.Print();
    161   //   nbentries = fn.NentriesFromFits(0);
    162    nbentries = NentriesFromFits(0);
     144  DVList dvl=is.DVListFromFits();
     145   nbentries = is.NentriesFromFits(0);
    163146  int lastpix=dvl.GetI("LASTPIX");
    164147  if (lastpix>0)
     
    211194  // On lit les DataBlocks;
    212195  dobj_->pixels_.ReSize(nPix);
    213   //  fn.GetSingleColumn(dobj_->pixels_.Data(),nPix);
    214   GetSingleColumn(dobj_->pixels_.Data(),nPix);
     196  is.GetSingleColumn(dobj_->pixels_.Data(),nPix);
    215197
    216198  // on effectue le decoupage en tranches
    217199  dobj_->SetThetaSlices();
    218   //  dobj_->Info()=fn.DVListFromFits();
    219   dobj_->Info()=DVListFromFits();
     200  dobj_->Info()=is.DVListFromFits();
    220201
    221202
Note: See TracChangeset for help on using the changeset viewer.