Changeset 3772 in Sophya for trunk/SophyaExt/FitsIOServer/fabtcolread.cc

Timestamp:

May 9, 2010, 12:07:49 AM (15 years ago)

Author:

cmv

Message:

lecture fits3D axe 1 dans TVector, cmv 09/05/2010

File:

• trunk/SophyaExt/FitsIOServer/fabtcolread.cc (modified) (1 diff)

trunk/SophyaExt/FitsIOServer/fabtcolread.cc

@@ -1925 +1925 @@
 }
 
+/*! Read 3D image into a TVector<int_4> */
+LONGLONG FitsImg3DRd::Read(LONGLONG j, LONGLONG k, TVector<int_4>& data)
+{
+ if(j<0 || k<0 || j>=Naxis[1] || k>=Naxis[2])
+   throw NotAvailableOperation("FitsImg3DRd::Read(TVector<int_4>): bad j/k number\n");
+ int sta=0;
+ if(data.Size() != Naxis[0]) data.ReSize(Naxis[0]);
+ int T = (sizeof(long)==4) ? TLONG: TINT;
+ LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
+ fits_read_img(GetFitsPtr(),T,deb,nel,&NulVal,data.Data(),NULL,&sta);
+ if(sta) {
+   FitsOpenFile::printerror(sta);
+   throw
+     NotAvailableOperation("FitsImg3DRd::Read(TVector<int_4>): Error Reading Fits file\n");
+  }
+ return Naxis[0];
+}
+
+/*! Read 3D image into a TVector<float> */
+LONGLONG FitsImg3DRd::Read(LONGLONG j, LONGLONG k, TVector<float>& data)
+{
+ if(j<0 || k<0 || j>=Naxis[1] || k>=Naxis[2])
+   throw NotAvailableOperation("FitsImg3DRd::Read(TVector<float>): bad j/k number\n");
+ int sta=0;
+ if(data.Size() != Naxis[0]) data.ReSize(Naxis[0]);
+ LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
+ fits_read_img(GetFitsPtr(),TFLOAT,deb,nel,&NulVal,data.Data(),NULL,&sta);
+ if(sta) {
+   FitsOpenFile::printerror(sta);
+   throw
+     NotAvailableOperation("FitsImg3DRd::Read(TVector<float>): Error Reading Fits file\n");
+ }
+ return Naxis[0];
+}
+
+/*! Read 3D image into a TVector<double> */
+LONGLONG FitsImg3DRd::Read(LONGLONG j, LONGLONG k, TVector<double>& data)
+{
+ if(j<0 || k<0 || j>=Naxis[1] || k>=Naxis[2])
+   throw NotAvailableOperation("FitsImg3DRd::Read(TVector<double>): bad j/k number\n");
+ int sta=0;
+ if(data.Size() != Naxis[0]) data.ReSize(Naxis[0]);
+ LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
+ fits_read_img(GetFitsPtr(),TDOUBLE,deb,nel,&NulVal,data.Data(),NULL,&sta);
+ if(sta) {
+   FitsOpenFile::printerror(sta);
+   throw
+     NotAvailableOperation("FitsImg3DRd::Read(TVector<double>): Error Reading Fits file\n");
+ }
+ return Naxis[0];
+}
+
 /*! Read 3D image pixel i,j,k with i=[0,Naxis1[ , j=[0,Naxis2[ , k=[0,Naxis3[ */
 double FitsImg3DRd::Read(LONGLONG i, LONGLONG j, LONGLONG k)