Changeset 1047 in Sophya for trunk/SophyaExt/FitsIOServer/fitstarray.cc
- Timestamp:
- Jun 15, 2000, 12:22:23 PM (25 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaExt/FitsIOServer/fitstarray.cc
r972 r1047 60 60 61 61 template <class T> 62 void FITS_TArray<T>::ReadFromFits( FitsFile& fn)62 void FITS_TArray<T>::ReadFromFits() 63 63 { 64 64 65 if (!fn.IsFitsImage()) 65 // if (!fn.IsFitsImage()) 66 if (!IsFitsImage()) 66 67 { 67 68 throw PException("ReadFromFits: the fits file seems not to be an image"); 68 69 } 69 int dimension = fn.nbDimOfImage(); 70 // int dimension = fn.nbDimOfImage(); 71 int dimension = nbDimOfImage(); 70 72 cout << " dimension de l'image a lire: " << dimension << endl; 71 73 72 74 uint_4* siz = new uint_4[dimension]; 73 for (int k=0; k< dimension; k++) siz[k] = fn.dimOfImageAxes()[k]; 75 // for (int k=0; k< dimension; k++) siz[k] = fn.dimOfImageAxes()[k]; 76 for (int k=0; k< dimension; k++) siz[k] = dimOfImageAxes()[k]; 74 77 if(dobj_ == NULL) 75 78 dobj_ = new TArray<T>(dimension,siz); … … 78 81 79 82 delete [] siz; 80 if (dobj_->Size() != fn.nbOfImageData() ) 83 // if (dobj_->Size() != fn.nbOfImageData() ) 84 if (dobj_->Size() != nbOfImageData() ) 81 85 { 82 86 cout << " total size of TArray: " << dobj_->Size() << endl; 83 cout << " total size from fits file: " << fn.nbOfImageData() << endl;87 cout << " total size from fits file: " << nbOfImageData() << endl; 84 88 throw PException("ReadFromFits: size conflict"); 85 89 } 86 90 // On lit le tableau 87 fn.GetSingleColumn( dobj_->Data(),dobj_->Size()); 91 // fn.GetSingleColumn( dobj_->Data(),dobj_->Size()); 92 GetSingleColumn( dobj_->Data(),dobj_->Size()); 88 93 89 94 } … … 91 96 92 97 template <class T> 93 void FITS_TArray<T>::WriteToFits( FitsFile& fn)98 void FITS_TArray<T>::WriteToFits() 94 99 { 95 100 if(dobj_ == NULL) return; … … 125 130 } 126 131 cout << " nombre total d'elements a copier " << nbels << endl; 127 fn.makeHeaderImageOnFits(type, nbdim, naxisn); 132 // fn.makeHeaderImageOnFits(type, nbdim, naxisn); 133 makeHeaderImageOnFits(type, nbdim, naxisn); 128 134 if (!dobj_->IsPacked()) 129 135 { … … 138 144 139 145 } 140 fn.putImageToFits(nbels, dobj_->Data()); 146 // fn.putImageToFits(nbels, dobj_->Data()); 147 putImageToFits(nbels, dobj_->Data()); 141 148 142 149 delete [] naxisn;
Note:
See TracChangeset
for help on using the changeset viewer.