Changeset 1176 in Sophya

Timestamp:

Sep 4, 2000, 11:17:21 AM (25 years ago)

Author:

ansari

Message:

Mollweide et sinus projections

Location:

trunk/SophyaExt/FitsIOServer

Files:

• fitsspherehealpix.cc (modified) (2 diffs)
• fitsspherehealpix.h (modified) (1 diff)

trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc

@@ -1 +1 @@
 #include "pexceptions.h"
 #include "fitsspherehealpix.h"
-
+#include "tarray.h"
+#include "fitstarray.h" 
 
 
@@ -198 +199 @@
 }
 
+template <class T>
+void FITS_SphereHEALPix<T>::Mollweide_picture_projection(char filename[])
+{
+  int ni = 300;
+  int nj = 600;
+  TMatrix<float> map(ni, nj);
+
+   int i,j;
+   for(i=0; i < ni; i++) 
+     {
+       double yd = (i+0.5)/ni-0.5;
+       double facteur=2.*Pi/sin(acos(yd*2));
+       double theta = (0.5-yd)*Pi;
+       for(j=0; j < nj; j++)  
+         {
+           double xa = (j+0.5)/nj-0.5;
+           double phi =  xa*facteur+Pi;
+           float th=float(theta);
+           float ff=float(phi);
+           if (phi<2*Pi && phi>= 0) 
+             {
+               map(i,j) = (float)dobj_->PixValSph(th, ff);
+             }
+         }
+     }
+   FitsOutFile  fout(filename);
+   fout.firstImageOnPrimaryHeader();
+   FITS_TArray<r_4> fta(map);
+   fta.Write(fout);
+   
+}
+template <class T>
+void FITS_SphereHEALPix<T>::sinus_picture_projection(char filename[])
+{
+  int ni = 300;
+  int nj = 600;
+  TMatrix<float> map(ni, nj);
+
+  int i,j;
+   for(i=0; i < ni; i++) 
+     {
+       double yd = (i+0.5)/ni-0.5;
+       double theta = (0.5-yd)*Pi;
+        double facteur=1./sin(theta);
+       for(j=0; j < nj; j++)  
+         {
+           double xa = (j+0.5)/nj-0.5;
+           double phi = 2.*Pi*xa*facteur+Pi;
+           float th=float(theta);
+           float ff=float(phi);
+           if (phi<2*Pi && phi>= 0) 
+             {
+               map(i,j) = (float)dobj_->PixValSph(th, ff);
+             }
+         }
+     }
+   FitsOutFile  fout(filename);
+   fout.firstImageOnPrimaryHeader();
+   FITS_TArray<r_4> fta(map);
+   fta.Write(fout);
+
+   
+}
+
 
 #ifdef __CXX_PRAGMA_TEMPLATES__

trunk/SophyaExt/FitsIOServer/fitsspherehealpix.h

@@ -29 +29 @@
 virtual void        SetDataObj(AnyDataObj & o);
 
+void Mollweide_picture_projection(char flnm[]);
+void sinus_picture_projection(char flnm[]);
+
+
 inline operator SphereHEALPix<T>() { return(*dobj_); }