Changeset 1553 in Sophya for trunk/SophyaPI/PIext

Timestamp:

Jun 27, 2001, 12:22:19 PM (24 years ago)

Author:

ansari

Message:

Info position angulaire ds affichage SphericalMap - Reza 27/6/2001

File:

• trunk/SophyaPI/PIext/nomskymapadapter.cc (modified) (4 diffs)

trunk/SophyaPI/PIext/nomskymapadapter.cc

@@ -215 +215 @@
 }
 
+// -- Fonction pour convertir un X,Y en Theta-Phi sur projection MolleWeide --
+char * _SphMollweideProj_XY2ThetaPhi(P2DArrayAdapter * aa, int ix, int iy);
 
 /* --Methode-- */
@@ -228 +230 @@
 TMatrix<T> * mtx = new TMatrix<T>(nr, nc);
 Project_Mol(*mtx);
-return (new POTMatrixAdapter<T>(mtx, true) );
+POTMatrixAdapter<T> * potma = new POTMatrixAdapter<T>(mtx, true);
+potma->SetInfoStringFunction(_SphMollweideProj_XY2ThetaPhi);
+return ( potma );
 }
 
@@ -271 +275 @@
     //    teta = (0.5-yd)*M_PI;
     for(c=0; c<nc; c++)  {
-      xa = (r_8)(c+0.5)/(r_8)nc-0.5;
+      xa = (r_8)(c+0.5)/(r_8)nc-0.5;  
       phi = xa*facteur+M_PI;
       if ( (phi <= 2*M_PI) && (phi >= 0.) ) {
@@ -279 +283 @@
     }
   }
+}
+
+static char m_buff[128];
+char * _SphMollweideProj_XY2ThetaPhi(P2DArrayAdapter * aa, int ix, int iy)
+{
+  r_8 xa, yd, teta,phi, facteur;
+  int_4 nl = aa->YSize();
+  int_4 nc = aa->XSize();
+  int l, c;
+  aa->IJ2xy(ix, iy, c, l);
+  yd = (r_8)(l+0.5)/(r_8)nl-0.5;
+  facteur=2.*M_PI/sin(acos((double)yd*2));
+  teta = (yd+0.5)*Pi;
+  xa = (r_8)(c+0.5)/(r_8)nc-0.5;  
+  phi = xa*facteur+M_PI;
+  if ( (phi <= 2*M_PI) && (phi >= 0.) ) 
+    sprintf(m_buff,"(l,b=%5.1f,%5.1f)", 90.-teta*180./M_PI, phi*180./M_PI);
+  else 
+    sprintf(m_buff,"(l,b=?,?)", teta*180./M_PI, phi*180./M_PI);
+
+  return (m_buff);
 }