Ignore:
Timestamp:
Jan 11, 2007, 7:07:52 PM (19 years ago)
Author:
cmv
Message:

gestion gros fichiers FITS dans fabtcolread et fabtwriter cmv 11/01/2007

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/SophyaExt/FitsIOServer/fabtcolread.cc

    r3114 r3128  
    209209 int sta=0; long val=0;
    210210 if(fits_read_key(fitsptr,TLONG,keyname,&val,NULL,&sta))
     211   printerror(sta);
     212 return val;
     213}
     214
     215/*!
     216  Read a fitsheader key into long
     217  \param fitsptr : cfitio pointer to Fits file
     218  \param keyname : name of the key
     219  \return value into long long
     220*/
     221LONGLONG FitsOpenFile::ReadKeyLL(fitsfile *fitsptr,char *keyname)
     222{
     223 if(keyname==NULL || fitsptr==NULL) return 0;
     224 int sta=0; LONGLONG val=0;
     225 if(fits_read_key(fitsptr,TLONGLONG,keyname,&val,NULL,&sta))
    211226   printerror(sta);
    212227 return val;
    … …  
    267282  fbt.Print(3);
    268283  // Read element by element
    269   for(long i=0;i<fbt.GetNbLine();i++) {
     284  for(LONGLONG i=0;i<fbt.GetNbLine();i++) {
    270285    double x = fbt.Read(i);
    271286    double y = fbt2.Read(i);
    … …  
    274289  // Read into a vector
    275290  TVector<double> data;
    276   long n = fbt.Read(32,50,data);
     291  LONGLONG n = fbt.Read(32,50,data);
    277292    cout<<"Number of values read: "<<n<<endl;
    278293  data.ReSize(100);
    … …  
    402417
    403418 // Get number of rows
    404  if(fits_get_num_rows(GetFitsPtr(),&NBline,&sta)) {
     419 if(fits_get_num_rowsll(GetFitsPtr(),&NBline,&sta)) {
    405420   FitsOpenFile::printerror(sta);
    406421   throw NotAvailableOperation("FitsABTColRd::Init: Error getting number of rows\n");
    … …  
    425440
    426441 // Get column type
    427  if(fits_get_coltype(GetFitsPtr(),ColNum+1,&ColTypeCode,NULL,NULL,&sta)) {
     442 if(fits_get_coltypell(GetFitsPtr(),ColNum+1,&ColTypeCode,NULL,NULL,&sta)) {
    428443   FitsOpenFile::printerror(sta);
    429444   throw ParmError("FitsABTColRd::Init: Error getting column type\n");
    … …  
    436451 // Get column name back, tunit, tform
    437452 char tunit[64], tform[64], tdisp[64];
    438  long repeat=0; double tscale=1., tzero=0.;
     453 LONGLONG repeat=0; double tscale=1., tzero=0.;
    439454 int rc=0;
    440455 if(HDUType()==BINARY_TBL) {
    441    fits_get_bcolparms(GetFitsPtr(),ColNum+1,labelcol,tunit,tform
     456   fits_get_bcolparmsll(GetFitsPtr(),ColNum+1,labelcol,tunit,tform
    442457                     ,&repeat,&tscale,&tzero,NULL,tdisp,&sta);
    443458 } else {
    444    fits_get_acolparms(GetFitsPtr(),ColNum+1,labelcol,&repeat,tunit,tform
     459   long repeatlng;
     460   fits_get_acolparms(GetFitsPtr(),ColNum+1,labelcol,&repeatlng,tunit,tform
    445461                     ,&tscale,&tzero,NULL,tdisp,&sta);
     462   repeat = repeatlng;
    446463 }
    447464 if(rc) {
    … …  
    525542{
    526543 return FitsOpenFile::ReadKeyL(GetFitsPtr(),keyname);
     544}
     545
     546/*!
     547  Read a fitsheader key into long long
     548  \param keyname : name of the key
     549  \return value into long long
     550*/
     551LONGLONG FitsABTColRd::ReadKeyLL(char *keyname)
     552{
     553 return FitsOpenFile::ReadKeyLL(GetFitsPtr(),keyname);
    527554}
    528555
    … …  
    549576  \endverbatim
    550577*/
    551 double FitsABTColRd::Read(long n,bool usebuffer)
     578double FitsABTColRd::Read(LONGLONG n,bool usebuffer)
    552579// Attention: n=nline [0,NBline[, cfistio veut [1,NBline]
    553580// Attention: colnum  [0,NBcol[ , cfistio veut [1,NBcol]
    … …  
    581608 if(n<LineDeb || n>LineFin) {
    582609   NFitsRead++;
    583    long row1,row2,nrow;
     610   LONGLONG row1,row2,nrow;
    584611   if(BuffSens>0) { // Cas remplissage forward
    585612     row1 = n+1;
    … …  
    640667  \endverbatim
    641668*/
    642 long FitsABTColRd::Read(long n1,long n2,TVector<double>& data)
     669LONGLONG FitsABTColRd::Read(LONGLONG n1,LONGLONG n2,TVector<double>& data)
    643670{
    644671 if(n1<0 || n1>=NBline)
    … …  
    649676 if(n2>=NBline) n2 = NBline-1;
    650677
    651  sa_size_t nread = n2-n1+1;
     678 LONGLONG nread = n2-n1+1;
    652679 if(data.Size()<nread) data.SetSize(nread);
    653680
    654  //for(long i=n1;i<=n2;i++) data(i-n1) = Read(i);
     681 //for(LONGLONG i=n1;i<=n2;i++) data(i-n1) = Read(i);
    655682 int sta=0;
    656683 fits_read_col(GetFitsPtr(),TDOUBLE,ColNum+1,n1+1,1,nread,NULL,data.Data(),NULL,&sta);
    … …  
    664691
    665692/*! idem before but for TVector of float */
    666 long FitsABTColRd::Read(long n1,long n2,TVector<float>& data)
     693LONGLONG FitsABTColRd::Read(LONGLONG n1,LONGLONG n2,TVector<float>& data)
    667694{
    668695 if(n1<0 || n1>=NBline)
    … …  
    673700 if(n2>=NBline) n2 = NBline-1;
    674701
    675  sa_size_t nread = n2-n1+1;
     702 LONGLONG nread = n2-n1+1;
    676703 if(data.Size()<nread) data.SetSize(nread);
    677704
    678  //for(long i=n1;i<=n2;i++) data(i-n1) = Read(i);
     705 //for(LONGLONG i=n1;i<=n2;i++) data(i-n1) = Read(i);
    679706 int sta=0;
    680707 fits_read_col(GetFitsPtr(),TFLOAT,ColNum+1,n1+1,1,nread,NULL,data.Data(),NULL,&sta);
    … …  
    688715
    689716/*! idem before but for TVector of unsigned short */
    690 long FitsABTColRd::Read(long n1,long n2,TVector<uint_2>& data)
     717LONGLONG FitsABTColRd::Read(LONGLONG n1,LONGLONG n2,TVector<uint_2>& data)
    691718{
    692719 if(n1<0 || n1>=NBline)
    … …  
    697724 if(n2>=NBline) n2 = NBline-1;
    698725
    699  sa_size_t nread = n2-n1+1;
     726 LONGLONG nread = n2-n1+1;
    700727 if(data.Size()<nread) data.SetSize(nread);
    701728
    … …  
    711738
    712739/*! idem before but for TVector of int_4 */
    713 long FitsABTColRd::Read(long n1,long n2,TVector<int_4>& data)
     740LONGLONG FitsABTColRd::Read(LONGLONG n1,LONGLONG n2,TVector<int_4>& data)
    714741{
    715742 if(n1<0 || n1>=NBline)
    … …  
    720747 if(n2>=NBline) n2 = NBline-1;
    721748
    722  sa_size_t nread = n2-n1+1;
     749 LONGLONG nread = n2-n1+1;
    723750 if(data.Size()<nread) data.SetSize(nread);
    724751
    725  //for(long i=n1;i<=n2;i++) data(i-n1) = Read(i);
     752 //for(LONGLONG i=n1;i<=n2;i++) data(i-n1) = Read(i);
    726753 int sta=0;
    727754 int T = (sizeof(long)==4) ? TLONG: TINT;
    … …  
    736763
    737764/*! idem before but for TVector of int_8 */
    738 long FitsABTColRd::Read(long n1,long n2,TVector<int_8>& data)
     765LONGLONG FitsABTColRd::Read(LONGLONG n1,LONGLONG n2,TVector<int_8>& data)
    739766{
    740767#ifdef TLONGLONG
    … …  
    746773 if(n2>=NBline) n2 = NBline-1;
    747774
    748  sa_size_t nread = n2-n1+1;
     775 LONGLONG nread = n2-n1+1;
    749776 if(data.Size()<nread) data.SetSize(nread);
    750777
    … …  
    773800  \return <0 means not found
    774801*/
    775 long FitsABTColRd::FirstRow(double val1,double val2,long rowstart)
    776 {
    777  long row = -1;
     802LONGLONG FitsABTColRd::FirstRow(double val1,double val2,LONGLONG rowstart)
     803{
     804 LONGLONG row = -1;
    778805 if(NBline==0) return row;
    779806 // Change buffer for efficiency
    … …  
    782809 if(rowstart<0) rowstart = 0;
    783810 if(rowstart>=NBline) rowstart = NBline-1;
    784  for(long i=rowstart;i<NBline;i++) {
     811 for(LONGLONG i=rowstart;i<NBline;i++) {
    785812   double val = Read(i);
    786813   if(val<val1 || val>val2) continue;
    … …  
    802829  \endverbatim
    803830*/
    804 long FitsABTColRd::LastRow(double val1,double val2,long rowstart)
    805 {
    806  long row = -1;
     831LONGLONG FitsABTColRd::LastRow(double val1,double val2,LONGLONG rowstart)
     832{
     833 LONGLONG row = -1;
    807834 if(NBline==0) return row;
    808835 // Change buffer for efficiency
    … …  
    810837 if(bsens>=0) {ChangeBuffer(BuffLen,-1); bchange=true;}
    811838 if(rowstart<0 || rowstart>=NBline) rowstart = NBline-1;
    812  for(long i=rowstart;i>=0;i--) {
     839 for(LONGLONG i=rowstart;i>=0;i--) {
    813840   double val = Read(i);
    814841   if(val<val1 || val>val2) continue;
    … …  
    853880  fbt.Print(3);
    854881  // Read element by element
    855   for(long i=0;i<fbt.GetNbLine();i++) {
     882  for(LONGLONG i=0;i<fbt.GetNbLine();i++) {
    856883    double x = fbt.Read(i);
    857884    if(i%lpmod==0) cout<<i<<": "<<x<<endl;
    … …  
    859886  // Read into a vector
    860887  TVector<double> data;
    861   long n = fbt.Read(32,50,data);
     888  LONGLONG n = fbt.Read(32,50,data);
    862889    cout<<"Number of values read: "<<n<<endl;
    863890  data.ReSize(100);
    … …  
    10271054 // Get NAXIS 1 et 2
    10281055 int nfound=0;
    1029  if(fits_read_keys_lng(GetFitsPtr(),"NAXIS",1,2,Naxis,&nfound,&sta)) {
     1056 // car fits_read_keys_lnglng n'est pas prototype dans longnam.h
     1057 if(ffgknjj(GetFitsPtr(),"NAXIS",1,2,Naxis,&nfound,&sta)) {
    10301058   FitsOpenFile::printerror(sta);
    10311059   throw RangeCheckError("FitsImg2DRd::Init: Error reading NAXIS cards\n");
    … …  
    10651093{
    10661094 return FitsOpenFile::ReadKeyL(GetFitsPtr(),keyname);
     1095}
     1096
     1097/*!
     1098  Read a fitsheader key into long long
     1099  \param keyname : name of the key
     1100  \return value into long long
     1101*/
     1102LONGLONG FitsImg2DRd::ReadKeyLL(char *keyname)
     1103{
     1104 return FitsOpenFile::ReadKeyLL(GetFitsPtr(),keyname);
    10671105}
    10681106
    … …  
    10881126\warning TMatrix data(Naxis2,Naxis1)
    10891127*/
    1090 size_t FitsImg2DRd::Read(TMatrix<uint_2>& data)
     1128LONGLONG FitsImg2DRd::Read(TMatrix<uint_2>& data)
    10911129{
    10921130 int sta=0;
    … …  
    10941132 data.ReSize(Naxis[1],Naxis[0]);
    10951133
    1096  for(int j=0;j<Naxis[1];j++) {
    1097    long deb = j*Naxis[0]+1, nel = Naxis[0];
     1134 for(LONGLONG j=0;j<Naxis[1];j++) {
     1135   LONGLONG deb = j*Naxis[0]+1, nel = Naxis[0];
    10981136   fits_read_img(GetFitsPtr(),TUSHORT,deb,nel,&NulVal,arr,NULL,&sta);
    10991137   if(sta) {
    … …  
    11021140       NotAvailableOperation("FitsImg2DRd::Read(TMatrix<uint_2>): Error Reading Fits file\n");
    11031141   }
    1104    for(int i=0;i<Naxis[0];i++) data(j,i) = arr[i];
     1142   for(LONGLONG i=0;i<Naxis[0];i++) data(j,i) = arr[i];
    11051143 }
    11061144
    … …  
    11101148
    11111149/*! Read image into a TMatrix<int_4> */
    1112 size_t FitsImg2DRd::Read(TMatrix<int_4>& data)
     1150LONGLONG FitsImg2DRd::Read(TMatrix<int_4>& data)
    11131151{
    11141152 int sta=0;
    … …  
    11171155 int T = (sizeof(long)==4) ? TLONG: TINT;
    11181156
    1119  for(int j=0;j<Naxis[1];j++) {
    1120    long deb = j*Naxis[0]+1, nel = Naxis[0];
     1157 for(LONGLONG j=0;j<Naxis[1];j++) {
     1158   LONGLONG deb = j*Naxis[0]+1, nel = Naxis[0];
    11211159   fits_read_img(GetFitsPtr(),T,deb,nel,&NulVal,arr,NULL,&sta);
    11221160   if(sta) {
    … …  
    11251163       NotAvailableOperation("FitsImg2DRd::Read(TMatrix<int_4>): Error Reading Fits file\n");
    11261164   }
    1127    for(int i=0;i<Naxis[0];i++) data(j,i) = arr[i];
     1165   for(LONGLONG i=0;i<Naxis[0];i++) data(j,i) = arr[i];
    11281166 }
    11291167
    … …  
    11331171
    11341172/*! Read image into a TMatrix<int_8> */
    1135 size_t FitsImg2DRd::Read(TMatrix<int_8>& data)
     1173LONGLONG FitsImg2DRd::Read(TMatrix<int_8>& data)
    11361174{
    11371175 int sta=0;
    … …  
    11391177 data.ReSize(Naxis[1],Naxis[0]);
    11401178
    1141  for(int j=0;j<Naxis[1];j++) {
    1142    long deb = j*Naxis[0]+1, nel = Naxis[0];
     1179 for(LONGLONG j=0;j<Naxis[1];j++) {
     1180   LONGLONG deb = j*Naxis[0]+1, nel = Naxis[0];
    11431181   fits_read_img(GetFitsPtr(),TLONGLONG,deb,nel,&NulVal,arr,NULL,&sta);
    11441182   if(sta) {
    … …  
    11471185       NotAvailableOperation("FitsImg2DRd::Read(TMatrix<int_8>): Error Reading Fits file\n");
    11481186   }
    1149    for(int i=0;i<Naxis[0];i++) data(j,i) = arr[i];
     1187   for(LONGLONG i=0;i<Naxis[0];i++) data(j,i) = arr[i];
    11501188 }
    11511189
    … …  
    11551193
    11561194/*! Read image into a TMatrix<float> */
    1157 size_t FitsImg2DRd::Read(TMatrix<float>& data)
     1195LONGLONG FitsImg2DRd::Read(TMatrix<float>& data)
    11581196{
    11591197 int sta=0;
    … …  
    11611199 data.ReSize(Naxis[1],Naxis[0]);
    11621200
    1163  for(int j=0;j<Naxis[1];j++) {
    1164    long deb = j*Naxis[0]+1, nel = Naxis[0];
     1201 for(LONGLONG j=0;j<Naxis[1];j++) {
     1202   LONGLONG deb = j*Naxis[0]+1, nel = Naxis[0];
    11651203   fits_read_img(GetFitsPtr(),TFLOAT,deb,nel,&NulVal,arr,NULL,&sta);
    11661204   if(sta) {
    … …  
    11691207       NotAvailableOperation("FitsImg2DRd::Read(TMatrix<float>): Error Reading Fits file\n");
    11701208   }
    1171    for(int i=0;i<Naxis[0];i++) data(j,i) = arr[i];
     1209   for(LONGLONG i=0;i<Naxis[0];i++) data(j,i) = arr[i];
    11721210 }
    11731211
    … …  
    11771215
    11781216/*! Read image into a TMatrix<double> */
    1179 size_t FitsImg2DRd::Read(TMatrix<double>& data)
     1217LONGLONG FitsImg2DRd::Read(TMatrix<double>& data)
    11801218{
    11811219 int sta=0;
    … …  
    11831221 data.ReSize(Naxis[1],Naxis[0]);
    11841222
    1185  for(int j=0;j<Naxis[1];j++) {
    1186    long deb = j*Naxis[0]+1, nel = Naxis[0];
     1223 for(LONGLONG j=0;j<Naxis[1];j++) {
     1224   LONGLONG deb = j*Naxis[0]+1, nel = Naxis[0];
    11871225   fits_read_img(GetFitsPtr(),TDOUBLE,deb,nel,&NulVal,arr,NULL,&sta);
    11881226   if(sta) {
    … …  
    11911229       NotAvailableOperation("FitsImg2DRd::Read(TMatrix<double>): Error Reading Fits file\n");
    11921230   }
    1193    for(int i=0;i<Naxis[0];i++) data(j,i) = arr[i];
     1231   for(LONGLONG i=0;i<Naxis[0];i++) data(j,i) = arr[i];
    11941232 }
    11951233
    … …  
    13431381 // Get NAXIS 1, 2 et 3
    13441382 int nfound=0;
    1345  if(fits_read_keys_lng(GetFitsPtr(),"NAXIS",1,3,Naxis,&nfound,&sta)) {
     1383 // car fits_read_keys_lnglng n'est pas prototype dans longnam.h
     1384 if(ffgknjj(GetFitsPtr(),"NAXIS",1,3,Naxis,&nfound,&sta)) {
    13461385   FitsOpenFile::printerror(sta);
    13471386   throw RangeCheckError("FitsImg3DRd::Init: Error reading NAXIS cards\n");
    … …  
    13841423
    13851424/*!
     1425  Read a fitsheader key into long long
     1426  \param keyname : name of the key
     1427  \return value into long long
     1428*/
     1429LONGLONG FitsImg3DRd::ReadKeyLL(char *keyname)
     1430{
     1431 return FitsOpenFile::ReadKeyLL(GetFitsPtr(),keyname);
     1432}
     1433
     1434/*!
    13861435  Read a fitsheader key into string
    13871436  \param keyname : name of the key
    … …  
    14011450Read 3D image into a TArray<uint_2>
    14021451*/
    1403 size_t FitsImg3DRd::Read(TArray<uint_2>& data)
     1452LONGLONG FitsImg3DRd::Read(TArray<uint_2>& data)
    14041453{
    14051454 int sta=0;
    … …  
    14071456 sa_size_t ndim[3] = {Naxis[0],Naxis[1],Naxis[2]}; data.ReSize(3,ndim);
    14081457
    1409  for(int k=0;k<Naxis[2];k++) for(int j=0;j<Naxis[1];j++) {
    1410    long deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
     1458 for(LONGLONG k=0;k<Naxis[2];k++) for(LONGLONG j=0;j<Naxis[1];j++) {
     1459   LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
    14111460   fits_read_img(GetFitsPtr(),TUSHORT,deb,nel,&NulVal,arr,NULL,&sta);
    14121461   if(sta) {
    … …  
    14151464       NotAvailableOperation("FitsImg3DRd::Read(TArray<uint_2>): Error Reading Fits file\n");
    14161465   }
    1417    for(int i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
     1466   for(LONGLONG i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
    14181467 }
    14191468
    … …  
    14231472
    14241473/*! Read 3D image into a TArray<int_4> */
    1425 size_t FitsImg3DRd::Read(TArray<int_4>& data)
     1474LONGLONG FitsImg3DRd::Read(TArray<int_4>& data)
    14261475{
    14271476 int sta=0;
    … …  
    14301479 int T = (sizeof(long)==4) ? TLONG: TINT;
    14311480
    1432  for(int k=0;k<Naxis[2];k++) for(int j=0;j<Naxis[1];j++) {
    1433    long deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
     1481 for(LONGLONG k=0;k<Naxis[2];k++) for(LONGLONG j=0;j<Naxis[1];j++) {
     1482   LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
    14341483   fits_read_img(GetFitsPtr(),T,deb,nel,&NulVal,arr,NULL,&sta);
    14351484   if(sta) {
    … …  
    14381487       NotAvailableOperation("FitsImg3DRd::Read(TArray<int_4>): Error Reading Fits file\n");
    14391488   }
    1440    for(int i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
     1489   for(LONGLONG i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
    14411490 }
    14421491
    … …  
    14461495
    14471496/*! Read 3D image into a TArray<int_8> */
    1448 size_t FitsImg3DRd::Read(TArray<int_8>& data)
     1497LONGLONG FitsImg3DRd::Read(TArray<int_8>& data)
    14491498{
    14501499 int sta=0;
    … …  
    14521501 sa_size_t ndim[3] = {Naxis[0],Naxis[1],Naxis[2]}; data.ReSize(3,ndim);
    14531502
    1454  for(int k=0;k<Naxis[2];k++) for(int j=0;j<Naxis[1];j++) {
    1455    long deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
     1503 for(LONGLONG k=0;k<Naxis[2];k++) for(LONGLONG j=0;j<Naxis[1];j++) {
     1504   LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
    14561505   fits_read_img(GetFitsPtr(),TLONGLONG,deb,nel,&NulVal,arr,NULL,&sta);
    14571506   if(sta) {
    … …  
    14601509       NotAvailableOperation("FitsImg3DRd::Read(TArray<int_8>): Error Reading Fits file\n");
    14611510   }
    1462    for(int i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
     1511   for(LONGLONG i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
    14631512 }
    14641513
    … …  
    14681517
    14691518/*! Read 3D image into a TArray<float> */
    1470 size_t FitsImg3DRd::Read(TArray<float>& data)
     1519LONGLONG FitsImg3DRd::Read(TArray<float>& data)
    14711520{
    14721521 int sta=0;
    … …  
    14741523 sa_size_t ndim[3] = {Naxis[0],Naxis[1],Naxis[2]}; data.ReSize(3,ndim);
    14751524
    1476  for(int k=0;k<Naxis[2];k++) for(int j=0;j<Naxis[1];j++) {
    1477    long deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
     1525 for(LONGLONG k=0;k<Naxis[2];k++) for(LONGLONG j=0;j<Naxis[1];j++) {
     1526   LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
    14781527   fits_read_img(GetFitsPtr(),TFLOAT,deb,nel,&NulVal,arr,NULL,&sta);
    14791528   if(sta) {
    … …  
    14821531       NotAvailableOperation("FitsImg3DRd::Read(TArray<float>): Error Reading Fits file\n");
    14831532   }
    1484    for(int i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
     1533   for(LONGLONG i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
    14851534 }
    14861535
    … …  
    14901539
    14911540/*! Read 3D image into a TArray<double> */
    1492 size_t FitsImg3DRd::Read(TArray<double>& data)
     1541LONGLONG FitsImg3DRd::Read(TArray<double>& data)
    14931542{
    14941543 int sta=0;
    … …  
    14961545 sa_size_t ndim[3] = {Naxis[0],Naxis[1],Naxis[2]}; data.ReSize(3,ndim);
    14971546
    1498  for(int k=0;k<Naxis[2];k++) for(int j=0;j<Naxis[1];j++) {
    1499    long deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
     1547 for(LONGLONG k=0;k<Naxis[2];k++) for(LONGLONG j=0;j<Naxis[1];j++) {
     1548   LONGLONG deb = Naxis[0]*(j+Naxis[1]*k)+1, nel = Naxis[0];
    15001549   fits_read_img(GetFitsPtr(),TDOUBLE,deb,nel,&NulVal,arr,NULL,&sta);
    15011550   if(sta) {
    … …  
    15041553       NotAvailableOperation("FitsImg3DRd::Read(TArray<double>): Error Reading Fits file\n");
    15051554   }
    1506    for(int i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
     1555   for(LONGLONG i=0;i<Naxis[0];i++) data(i,j,k) = arr[i];
    15071556 }
    15081557
Note: See TracChangeset for help on using the changeset viewer.