Changeset 972 in Sophya for trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc

Timestamp:

Apr 27, 2000, 6:12:10 PM (25 years ago)

Author:

ansari

Message:

mise a jour 27/04/00 GLM

File:

• trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc (modified) (6 diffs)

trunk/SophyaExt/FitsIOServer/fitsspherehealpix.cc

@@ -69 +69 @@
 
 template <class T>
-void FITS_SphereHEALPix<T>::Write(char outputfile[],int hdunum)
-{
-  WriteF(outputfile, hdunum);
+void FITS_SphereHEALPix<T>::Write(char outputfile[], bool OldFile)
+{
+  WriteF(outputfile, OldFile);
 }
 template <class T>
@@ -80 +80 @@
 
 template <class T>
-void FITS_SphereHEALPix<T>::WriteToFits(const FitsFile& fn)
+void FITS_SphereHEALPix<T>::WriteToFits(FitsFile& fn)
 {
   if(dobj_ == NULL) 
@@ -94 +94 @@
   dvl["COORDSYS"]= description;
 
+  dvl["PIXTYPE"] = "HEALPIX";
+  dvl.SetComment("PIXTYPE", "HEALPIX Pixelization");
   dvl["ORDERING"]= dobj_->TypeOfMap();
+  dvl.SetComment("ORDERING", "Pixel ordering scheme, either RING or NESTED");
   
   int nSide= dobj_->SizeIndex();
   dvl["NSIDE"]= nSide;
+  dvl.SetComment("NSIDE","Resolution parameter for HEALPIX" );
 
   int nPix= dobj_->NbPixels();
   dvl["FIRSTPIX"]= 0;
+  dvl.SetComment("FIRSTPIX", "First pixel # (0 based)");
   dvl["LASTPIX"]= nPix-1;
+  dvl.SetComment("LASTPIX", "Last pixel # (0 based)");
   dvl["Content"]= "SphereHEALPix";
   dvl["EXTNAME"]= "SIMULATION";
@@ -131 +137 @@
 
 template <class T> 
-void FITS_SphereHEALPix<T>::ReadFromFits(const FitsFile& fn)
+void FITS_SphereHEALPix<T>::ReadFromFits(FitsFile& fn)
 {
   if(dobj_ == NULL) 
@@ -164 +170 @@
   int id= SphereCoordSys_NEUTRAL;
   string description= "NEUTRAL SphereCoordSystem";
-  const char* coordsys= dvl.GetS("COORDSYS").c_str();    
-  if(strncmp(coordsys,"ROTATION",8) == 0)
+  string coordsys= dvl.GetS("COORDSYS");    
+  if(coordsys.compare(0,7,"ROTATION",0,7) == 0)
     {
       id= SphereCoordSys_ROTATION;
       description= "ROTATION SphereCoordSystem";
     } 
-  else if(strncmp(coordsys,"OTHER",5) == 0)
+  else if(coordsys.compare(0,4,"OTHER",0,4) == 0)
     {
       id= SphereCoordSys_OTHER;
@@ -177 +183 @@
   dobj_->SetCoordSys(new SphereCoordSys(id,description));   
 
-  const char* ordering= dvl.GetS("ORDERING").c_str();
-  if(strncmp(ordering,"RING",4) != 0)
+  string ordering= dvl.GetS("ORDERING");
+  if(ordering.compare(0,3,"RING",0,3) != 0)
     { 
       cout << " ORDERING= " << ordering << endl;