Changeset 1136 in Sophya for trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc
- Timestamp:
- Aug 24, 2000, 12:18:02 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc
r1047 r1136 15 15 { 16 16 dobj_= new SphereHEALPix<T>; 17 ownobj = true;17 ownobj_= true; 18 18 } 19 19 … … 22 22 { 23 23 dobj_= new SphereHEALPix<T>; 24 ownobj = true;25 26 Read F(inputfile,hdunum);27 dobj_->SetTemp(true);24 ownobj_= true; 25 26 Read(inputfile,hdunum); 27 // dobj_->SetTemp(true); 28 28 } 29 29 … … 33 33 { 34 34 dobj_= new SphereHEALPix<T>(obj, true); 35 dobj_->SetTemp(true);36 ownobj = true;35 // dobj_->SetTemp(true); 36 ownobj_= true; 37 37 } 38 38 … … 41 41 { 42 42 dobj_= obj; 43 ownobj = false;43 ownobj_= false; 44 44 } 45 45 … … 47 47 FITS_SphereHEALPix<T>::~FITS_SphereHEALPix() 48 48 { 49 if (ownobj && dobj_) delete dobj_;49 if (ownobj_ && dobj_) delete dobj_; 50 50 } 51 51 … … 61 61 SphereHEALPix<T> * po = dynamic_cast< SphereHEALPix<T> * >(&o); 62 62 if (po == NULL) return; 63 if (ownobj && dobj_) delete dobj_;63 if (ownobj_ && dobj_) delete dobj_; 64 64 dobj_ = po; 65 ownobj = false;65 ownobj_ = false; 66 66 } 67 67 68 68 69 69 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 71 template <class T> 72 void FITS_SphereHEALPix<T>::WriteToFits(FitsOutFile& os) 83 73 { 84 74 if(dobj_ == NULL) … … 130 120 Type[1]='\0'; 131 121 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); 134 123 delete [] Noms[0]; 135 124 delete [] Noms; 136 putColToFits(0, nPix, dobj_->pixels_.Data());125 os.putColToFits(0, nPix, dobj_->pixels_.Data()); 137 126 } 138 127 139 128 template <class T> 140 void FITS_SphereHEALPix<T>::ReadFromFits( )129 void FITS_SphereHEALPix<T>::ReadFromFits(FitsInFile& is) 141 130 { 142 131 if(dobj_ == NULL) 143 132 { 144 133 dobj_= new SphereHEALPix<T>; 145 dobj_->SetTemp(true); 146 ownobj= true; 134 ownobj_= true; 147 135 } 148 136 149 137 150 138 int nbcols, nbentries; 151 // nbcols = fn.NbColsFromFits(); 152 nbcols = NbColsFromFits(); 139 nbcols = is.NbColsFromFits(); 153 140 if (nbcols != 1) 154 141 { 155 142 throw IOExc("le fichier fits n'est pas une sphere Healpix"); 156 143 } 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); 163 146 int lastpix=dvl.GetI("LASTPIX"); 164 147 if (lastpix>0) … … 211 194 // On lit les DataBlocks; 212 195 dobj_->pixels_.ReSize(nPix); 213 // fn.GetSingleColumn(dobj_->pixels_.Data(),nPix); 214 GetSingleColumn(dobj_->pixels_.Data(),nPix); 196 is.GetSingleColumn(dobj_->pixels_.Data(),nPix); 215 197 216 198 // on effectue le decoupage en tranches 217 199 dobj_->SetThetaSlices(); 218 // dobj_->Info()=fn.DVListFromFits(); 219 dobj_->Info()=DVListFromFits(); 200 dobj_->Info()=is.DVListFromFits(); 220 201 221 202
Note:
See TracChangeset
for help on using the changeset viewer.