Changeset 1119 in Sophya for trunk/SophyaProg

Timestamp:

Jul 31, 2000, 3:18:35 PM (25 years ago)

Author:

ansari

Message:

Remplacement des FitsIoServer par des FITS_XXX<T> ds tspm2.cc , Reza 31/7/2000

File:

• trunk/SophyaProg/Tests/tspm2.cc (modified) (7 diffs)

trunk/SophyaProg/Tests/tspm2.cc

@@ -5 +5 @@
 
 #include "tod.h"
-#include "fitsioserver.h"
+#include "fitstarray.h"
+#include "fitsspherehealpix.h"
 #include "nbrandom.h"
 
@@ -50 +51 @@
 
 }
-template <class T>
-void  Project_Mol_Old(PixelMap<T> const & map, ImageR4 & img, float defval=-999.)
-{
-  r_8 xa, yd, teta,phi, facteur;
-  int_4 i,j,k,n;
-  printf("Xsize= %d   Ysize= %d  NPix= %d\n",img.XSize(),img.YSize(),img.XSize()*img.YSize() );
-  n = 0;
-  img.Zero();
-  for(j=0; j<img.YSize(); j++) {
-    yd = (r_4)(j+0.5)/(r_4)img.YSize()-0.5;
-    facteur=2.*M_PI/sin(acos((double)yd*2));
-    teta = (yd+0.5)*M_PI;
-    //    teta = (0.5-yd)*M_PI;
-    for(i=0; i<img.XSize(); i++)  {
-      xa = (r_4)(i+0.5)/(r_4)img.XSize()-0.5;
-      phi = xa*facteur+M_PI;
-      if ( (phi <= 2*M_PI) && (phi >= 0.) ) {
-        k = map.PixIndexSph(teta, phi);
-        img(i,j) = map(k);
-        //        if (n<20) {
-        //          printf("i,j= %d %d -> %g %g -> k=%d -> val= %g \n", 
-        //               i,j,teta,phi,k,map(k)); n++; }
-      }
-    }
-  }
-
-}
 
 int main(int narg, char* arg[]) 
@@ -94 +68 @@
   cout <<  " ===== HEALPix Spherical Map Test M= " << m << endl;
 
+  try {
+
   POutPersist s("spheres.ppf"); 
   string nomobj;
@@ -115 +91 @@
    }
 
-  FitsIoServer fios;
-  // On projete dans un fichier FITS
-  fios.sinus_picture_projection(sph, "gsin1.fits");
-  fios.save(sph, "sph1.fits");
 
   TMatrix<double> mtx(3*m, 6*m);
+  cout << " Project_Mol(sph, mtx) " << endl;
   Project_Mol(sph, mtx);
-  fios.save(mtx, "mtx1.fits");
-
-  cout << " Project_Mol_Old(sph, img) " << endl;
-  ImageR4 img(6*m, 3*m);
-  Project_Mol_Old(sph, img);
-  img.CheckDyn();
-  cout << img; 
-  fios.save(img, "img1.fits");
-  cout << " Apres writeFits " << endl;
-  img.CheckDyn();
-  cout << img; 
-  /*
-  ImageR4 imgr;
-  fios.load(imgr, "img1.fits");
-  cout << " Apres readFits " << endl;
-  imgr.CheckDyn();
-  cout << imgr; 
-  */
+  {
+  cout << " Writing to FITS prjmol1.fits " << endl;  
+  FITS_TArray<double> fios(mtx);
+  fios.Write("prjmol1.fits");
+  }
 
   // Computing mean and sigma on the sphere 
@@ -173 +133 @@
 
   // On projete dans un fichier FITS
-  FitsIoServer fios;
    {
    cout << "Test of Write/Read SphereHEALPix<float> to FITS (sphg_r4.fits) " << endl;
-   fios.save(sph, "sphg_r4.fits");
- 
-   SphereHEALPix<float> sphr(4);
-   fios.load(sphr, "sphg_r4.fits", 2);
+   FITS_SphereHEALPix<r_4> fiosw(sph);
+   fiosw.Write("sphg_r4.fits");
+   }
+
+   
+   SphereHEALPix<float> sphr;
+   {
+   FITS_SphereHEALPix<r_4> fiosr("sphg_r4.fits");
+   sphr = fiosr;
    cout << " Read from file - SphereHEALPix<float> NPixels= " << sphr.NbPixels() << endl;
    int ndiff = 0;
@@ -193 +157 @@
 
 
-  fios.sinus_picture_projection(sph, "gsin2.fits");
-  fios.save(sph, "sph2.fits");
-
   TMatrix<float> mtx(3*m, 6*m);
   Project_Mol(sph, mtx);
-  fios.save(mtx, "mtx2.fits");
-
-  
-  cout << " Project_Mol_Old(sph, img) " << endl;
-  ImageR4 img(6*m, 3*m);
-  Project_Mol_Old(sph, img);
-  img.CheckDyn();
-  cout << img;
-  fios.save(img, "img2.fits");
-  cout << " Apres writeFits " << endl;
-  img.CheckDyn();
-  cout << img; 
+  {
+  cout << " Writing to FITS prjmol2.fits " << endl;  
+  FITS_TArray<double> fios(mtx);
+  fios.Write("prjmol2.fits");
+  }
 
   // Computing mean and sigma on the sphere 
@@ -241 +195 @@
    }
   }
+  }
+  catch (PThrowable & exc) {
+    cerr << "tspm2: Catched Exception " << (string)typeid(exc).name()
+         << " - Msg= " << exc.Msg() << endl;
+  }
+  catch (...) {
+    cerr << "tspm2:  some other exception was caught ! " << endl;
+  }
+
   cout <<  " ===== Fin de TSPM2_Test ======== " << endl;
   return 0;