Changeset 518 in Sophya for trunk/SophyaLib/Samba/pixelmap.h

Timestamp:

Oct 26, 1999, 4:31:44 PM (26 years ago)

Author:

ansari

Message:

Introduction de SpherePosition and SphereCoordSys, and Initiator for module Samba - Reza+I. Grivell 26/10/99

File:

• trunk/SophyaLib/Samba/pixelmap.h (modified) (6 diffs)

trunk/SophyaLib/Samba/pixelmap.h

@@ +1 @@
+// This may look like C code, but it is really -*- C++ -*-
 #ifndef PIXELMAP_SEEN
 #define PIXELMAP_SEEN
@@ -4 +5 @@
 #include "ppersist.h" 
 #include "dvlist.h"
+#include "spherepos.h"
 #include <iostream.h>
 
@@ -21 +23 @@
 public:
   
-PixelMap():mInfo_(NULL) {};
-virtual ~PixelMap() {};
-  
+              PixelMap(SphereCoordSys* cs = NULL) : mInfo_(NULL) 
+                    { if (cs) cs_ = cs; else cs = new SphereCoordSys; }
+virtual       ~PixelMap() 
+                    { if (mInfo_) delete mInfo_; if (cs_) delete cs_; }
+
+// Set/Change/Get the coordinate system
+virtual void  SetCoordSys(SphereCoordSys* cs)
+                    { if (cs) { delete cs_;  cs_ = cs; } }
+inline SphereCoordSys* GetCoordSys() const { return(cs_); }
+
 // Number of pixels
 virtual int NbPixels() const=0;
@@ -30 +39 @@
 virtual T& PixVal(int k)=0;
 virtual T const& PixVal(int k) const=0;
-   
+
+// Map s coverage
+virtual bool ContainsSph(double theta, double phi) const=0;
+virtual bool Contains(const SpherePosition& spos) const;
+
 // Index of pixel at (theta,phi)
 virtual int PixIndexSph(double theta, double phi) const=0;
+// Index of pixel at a sky-position 
+virtual int PixIndex(const SpherePosition& spos);
+
 
 // Value of pixel number at (theta,phi)
@@ -55 +71 @@
 inline T const& operator()(int k) const {return(PixVal(k));}
 
+// Overloading of () to access pixel at a sky position .
+inline T& operator()(const SpherePosition& spos) 
+                     { return(PixVal(PixIndex(spos))); }
+inline T const& operator()(const SpherePosition& spos) const 
+                     { return(PixVal(PixIndex(spos))); }
+
+
 // Note : no overloading of (double,double) to access pixel (theta,phi).
 // overloading of (double,double) in SphericalMap
@@ -75 +98 @@
 
 protected :
-
+    SphereCoordSys *cs_;   // Coordinate system used in the map
     DVList* mInfo_;        // Infos (variables) attachees 
 };
+
+
+template <class T>
+int PixelMap<T>::PixIndex(const SpherePosition& spos)
+{
+UnitVector v = spos.Transform(*cs_);
+return(PixIndexSph(v.Theta(), v.Phi()));
+}
+
+template <class T>
+bool PixelMap<T>::Contains(const SpherePosition& spos) const
+{
+UnitVector v = spos.Transform(*cs_);
+return(ContainsSph(v.Theta(), v.Phi()));
+}
+
 #endif