Changeset 2613 in Sophya for trunk/SophyaPI/PIext/nomskymapadapter.cc

Timestamp:

Sep 7, 2004, 3:41:54 PM (21 years ago)

Author:

ansari

Message:

MAJ Makefile et piaversion.h, modification NomSkyMapAdapter apres ajout SphereECP<T> - Reza , 7 Septembre 2004

File:

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

trunk/SophyaPI/PIext/nomskymapadapter.cc

@@ -99 +99 @@
 SphereHEALPix<T>* sg = dynamic_cast< SphereHEALPix<T> * >(mMap); 
 if (sg != NULL) type = "SphereHEALPix< ";
+SphereECP<T>* se = dynamic_cast< SphereECP<T> * >(mMap); 
+if (se != NULL) type = "SphereECP< ";
 
 // type +=  DecodeTypeIdName(typeid(T).name());
@@ -116 +118 @@
 SphereHEALPix<T>* sg = dynamic_cast< SphereHEALPix<T> * >(mMap); 
 if (sg != NULL) return( new SphereHEALPix<T>(*sg, share) );
+SphereECP<T>* se = dynamic_cast< SphereECP<T> * >(mMap); 
+if (se != NULL) return( new SphereECP<T>(*se, share) );
 return(NULL);
 }
@@ -210 +214 @@
 if (sg != NULL) {
   FIO_SphereHEALPix<T> fio(sg);
+  fio.Write(pos, nom);
+  return;
+  }
+SphereECP<T>* se = dynamic_cast< SphereECP<T> * >(mMap); 
+if (se != NULL) {
+  FIO_SphereECP<T> fio(se);
   fio.Write(pos, nom);
   return;
@@ -233 +243 @@
 /* --Methode-- */
 template <class T> 
-P2DArrayAdapter* NOMAdapter_PixelMap<T>::Get2DArray(string &)
+P2DArrayAdapter* NOMAdapter_PixelMap<T>::Get2DArray(string & opt)
 {
 LocalMap<T>* lm = dynamic_cast< LocalMap<T> * >(mMap); 
@@ -239 +249 @@
 int nr = 250;
 int nc = 500;
+SphereECP<T>* se = dynamic_cast< SphereECP<T> * >(mMap); 
+size_t olen = opt.length();
+if ((se != NULL) && (opt.find("ecparray") < olen)) {
+  // On veut la carte bi-dim complete de la sphere ECP
+  TMatrix<T> * mtxecp = new TMatrix<T>(se->GetPixelArray(), true);  // on partage les donnees
+  POTMatrixAdapter<T> * potmtxecp = new POTMatrixAdapter<T>(mtxecp, true);
+  if (opt.find("ecpscaledeg") < olen)  // Echelle angulaire en degre
+    potmtxecp->DefineXYCoordinates(se->MinPhi()*180./M_PI, se->MinTheta()*180./M_PI, 
+                                   se->DeltaPhi()*180./M_PI, se->DeltaTheta()*180./M_PI);
+  else if (opt.find("ecpscalerad") < olen)  // Echelle angulaire en radian
+    potmtxecp->DefineXYCoordinates(se->MinPhi(), se->MinTheta(), se->DeltaPhi(), se->DeltaTheta());
+
+  return ( potmtxecp );
+}
 SphericalMap<T>* sm = dynamic_cast< SphericalMap<T> *>(mMap);
 if (sm != NULL) { nr = sqrt(0.75*mMap->NbPixels());  nc = 2*nr; }