Changeset 2621 in Sophya for trunk/SophyaLib/SkyMap/sphereecp.cc

Timestamp:

Sep 23, 2004, 10:38:41 AM (21 years ago)

Author:

ansari

Message:

Correction pour codage pos theta,phi ds index pour pixels hors couverture de la carte SphereECP - Reza 23/09/2004

File:

• trunk/SophyaLib/SkyMap/sphereecp.cc (modified) (19 diffs)

trunk/SophyaLib/SkyMap/sphereecp.cc

@@ -15 +15 @@
   _phi2 = 2*M_PI;
   _outofmapidx = -99;
+  _outofmapnphi = 0;
+  _outofmapntet = 0;
   _outofmappix = 0;
   _outofmapval = 0;
@@ -31 +33 @@
   _phi2 = 2*M_PI;
   _outofmapidx = -99;
+  _outofmapnphi = 0;
+  _outofmapntet = 0;
   _outofmappix = 0;
   _outofmapval = 0;
@@ -46 +50 @@
   _phi2 = 2*M_PI; 
   _outofmapidx = -99;
+  _outofmapnphi = 0;
+  _outofmapntet = 0;
   _outofmappix = 0;
   _outofmapval = 0;
@@ -67 +73 @@
   _phi2 = phi2; 
   _outofmapidx = _pixels.Size()+659;
+  // Reza, 23 sep 2004 
+  // pour des cartes partielles, et des pixels sur la sphere, mais 
+  // en dehors de la couverture, on code le numero de tranche en theta, 
+  // compte a partir de theta=0,phi=0 +  _outofmapidx
+  // theta -> idx =  _outofmapidx + theta/_dtheta*_outofmapnphi +  _outofmapntet
+  _outofmapntet = 0;
   _outofmappix = 0;
   _outofmapval = 0;
   _dtheta = (_theta2-_theta1)/ntet;
   _dphi = (_phi2-_phi1)/nphi;
+
+  _outofmapntet = M_PI/_dtheta;
+  _outofmapnphi = 2*M_PI/_dphi;
 }
 
@@ -83 +98 @@
   _phi2 = a._phi2;   
   _outofmapidx = a._outofmapidx;
+  _outofmapnphi = a._outofmapnphi;
+  _outofmapntet = a._outofmapntet;
   _outofmappix = a._outofmappix;
   _outofmapval = a._outofmapval;
@@ -99 +116 @@
   _phi2 = a._phi2;   
   _outofmapidx = a._outofmapidx;
+  _outofmapnphi = a._outofmapnphi;
+  _outofmapntet = a._outofmapntet;
   _outofmappix = a._outofmappix;
   _outofmapval = a._outofmapval;
@@ -131 +150 @@
 T& SphereECP<T>::PixVal(int_4 k)
 {
-  if ((_partial) && (k==_outofmapidx))   return _outofmappix;
+  // On pourrait etre plus strict en demandant aussi k<_outofmapidx+_outofmapnphi*_outofmapntet (Reza 23/09/04)
+  if ((_partial) && (k>=_outofmapidx))   return _outofmappix;
   if((k < 0) || (k >= _pixels.Size()) )  
      throw RangeCheckError("SphereECP::PIxVal() Pixel index out of range");
@@ -140 +160 @@
 T const& SphereECP<T>::PixVal(int_4 k)  const
 {
-  if ((_partial) && (k==_outofmapidx))   return _outofmapval;
+  // On pourrait etre plus strict en demandant aussi k<_outofmapidx+_outofmapnphi*_outofmapntet (Reza 23/09/04)
+  if ((_partial) && (k>=_outofmapidx))   return _outofmapval;
   if((k < 0) || (k >= _pixels.Size()) )  
      throw RangeCheckError("SphereECP::PIxVal() Pixel index out of range");
@@ -165 +186 @@
   if (_partial) {
     if ( (theta < _theta1) || (theta > _theta2) ||
-         (phi < _phi1) || (phi > _phi2) )  return _outofmapidx;
+         (phi < _phi1) || (phi > _phi2) )  return _outofmapidx+theta/_dtheta*_outofmapnphi+phi/_dphi;
   }
   int_4 it = (theta-_theta1)/_dtheta;
@@ -175 +196 @@
 void  SphereECP<T>::PixThetaPhi(int_4 k, double& theta, double& phi) const
 {
-  if ((_partial) && (k==_outofmapidx)) {
-    theta = M_PI;
-    phi = 2*M_PI;
+  if ((_partial) && (k>=_outofmapidx)) {
+    theta = (((k-_outofmapidx) / _outofmapnphi)+0.5)*_dtheta;
+    phi = (((k-_outofmapidx) % _outofmapnphi)+0.5)*_dphi;;
     return;
   }
@@ -185 +206 @@
   int_4 it = k / _pixels.SizeX();
   int_4 jp = k % _pixels.SizeX();
-  theta = _theta1+it*_dtheta;
-  phi = _phi1*jp*_dphi;
+  theta = _theta1+(it+0.5)*_dtheta;
+  phi = _phi1*(jp+0.5)*_dphi;
   return;
 }
@@ -200 +221 @@
 double  SphereECP<T>::PixSolAngle(int_4 k) const
 {
-  if ((k < 0) || (k >= _pixels.Size()))  return (_dtheta*_dphi);
+  if ((_partial) && (k>=_outofmapidx)) {
+    double theta = (((k-_outofmapidx) / _outofmapnphi)+0.5)*_dtheta;
+    return (_dtheta*_dphi*sin(theta));
+  }
+  if((k < 0) || (k >= _pixels.Size()) )  
+     throw RangeCheckError("SphereECP::PixSolAngle(int_4 k) Pixel index out of range");
+
   int_4 it = k / _pixels.SizeX();
   double theta = _theta1+it*_dtheta;
@@ -229 +256 @@
   _dtheta = (_theta2-_theta1)/m;
   _dphi = (_phi2-_phi1)/mphi;
+
+  _outofmapntet = M_PI/_dtheta;
+  _outofmapnphi = 2*M_PI/_dphi;
   return;
 }
@@ -245 +275 @@
      throw RangeCheckError("SphereECP::GetThetaSlice() theta index out of range");
 
-  theta = _theta1 + index*_dtheta;
+  theta = _theta1 + (index+0.5)*_dtheta;
   int_4 nphit = 2.*M_PI/_dphi;
   phi.ReSize(nphit); 
@@ -252 +282 @@
   int_4 i;
   for(i=0; i<ioff; i++) {
-    phi(i) = _phi1 + (i-ioff)*_dphi;
+    phi(i) = _phi1 + (i-ioff+0.5)*_dphi;
     value(i) = _outofmapval;
   }
   for(i=ioff; i<ioff+_pixels.SizeX(); i++) {
-    phi(i) = _phi1 + (i-ioff)*_dphi;
+    phi(i) = _phi1 + (i-ioff+0.5)*_dphi;
     value(i) = _pixels(i-ioff,index,0);
   }
   for(i=ioff+_pixels.SizeX(); i<nphit; i++) {
-    phi(i) = _phi1 + (i-ioff)*_dphi;
+    phi(i) = _phi1 + (i-ioff+0.5)*_dphi;
     value(i) = _outofmapval;
   }
@@ -273 +303 @@
     throw RangeCheckError("SphereECP::GetThetaSlice() theta index out of range");
   
-  theta = _theta1 + index*_dtheta;
+  theta = _theta1 + (index+0.5)*_dtheta;
   int_4 ioff = _phi1/_dphi;
-  phi0 = _phi1 - ioff*_dphi;
+  phi0 = _phi1 - (ioff-0.5)*_dphi;
 
   int_4 nphit = 2.*M_PI/_dphi;
@@ -283 +313 @@
   int_4 i;
   for(i=0; i<ioff; i++) {
-    pixidx(i) =  _outofmapidx;
+    pixidx(i) = _outofmapidx+theta/_dtheta*_outofmapnphi+phi0/_dphi+i;
     value(i) = _outofmapval;
   }
@@ -291 +321 @@
   }
   for(i=ioff+_pixels.SizeX(); i<nphit; i++) {
-    pixidx(i) =  _outofmapidx;
+    pixidx(i) = _outofmapidx+theta/_dtheta*_outofmapnphi+phi0/_dphi+i;
     value(i) = _outofmapval;
   }
@@ -322 +352 @@
   _phi2 = a._phi2;   
   _outofmapidx = a._outofmapidx;
+  _outofmapnphi = a._outofmapnphi;
+  _outofmapntet = a._outofmapntet;
   _outofmappix = a._outofmappix;
   _outofmapval = a._outofmapval;