Changeset 1626 in Sophya for trunk/SophyaProg/PrgMap/extrapsph.cc

Timestamp:

Aug 7, 2001, 6:20:32 PM (24 years ago)

Author:

cmv

Message:

Uniformisation des programmes + doc cmv 7/8/01

File:

• trunk/SophyaProg/PrgMap/extrapsph.cc (modified) (8 diffs)

trunk/SophyaProg/PrgMap/extrapsph.cc

@@ -1 @@
-// Pour boucher les trous d'une sphere HEALPIX en clusterisant dans une sphere
-// plus petite             cmv 13/6/01
+// Pour boucher les trous d'une sphere HEALPIX en la clusterisant
+// dans une sphere plus petite puis en la re-extrapolant.
+//                    cmv 13/6/01
 // extrapsph -r 2 sph143k05.fits sph143k05_e.fits
 //                sphout.fits sphoutw.fits sphred.fits sphredw.fits
@@ -13 +14 @@
 #include "fitsspherehealpix.h" 
 
+/*!
+  \ingroup PrgMap
+  \file extrapsph.cc
+  \brief \b extrapsph: Fill holes in a sphere of datas by reducing
+  it size and re-extrapolating.
+  \verbatim
+csh> extrapsph -h
+extrapsph [-r reduc_factor] sphin.fits sphinw.fits
+          sphout.fits [sphout_w.fits sphred.fits sphred_w.fits]
+ -r reduc : reduction factor for clustering (must be 2^n, def=2)
+ sphin.fits    : input sphere
+ sphin_w.fits  : input sphere filling weight
+ sphout.fits   : output sphere
+ sphout_w.fits : output sphere filling weight
+ sphred.fits   : output reducted sphere
+ sphred_w.fits : output reducted sphere filling weight
+  \endverbatim
+*/
+
 void usage();
 void usage()
 {
-cout<<"extrapsph [-r reduc] sphin.fits sphinw.fits"<<endl
-    <<"                     sphout.fits [sphout_w.fits"
+cout<<"extrapsph [-r reduc_factor] sphin.fits sphinw.fits"<<endl
+    <<"          sphout.fits [sphout_w.fits"
     <<" sphred.fits sphred_w.fits]"<<endl
     <<" -r reduc : reduction factor for clustering (must be 2^n, def=2)"<<endl
@@ -31 +51 @@
 {
 int red=2;
+string dum;
 int c;
 while((c = getopt(narg,arg,"hr:")) != -1) {
@@ -122 +143 @@
   if(sphredw(ir)<=0.) continue;
   sphin(i)  = sphred(ir);
+  // On passe en negatif les pixels qui ont ete extrapoles
   sphinw(i) = -sphredw(ir);
   n++;
@@ -129 +151 @@
 // Ecriture des spheres
 {
-string dum = "rm -f "; dum += fsphout; system(dum.c_str());
+dum = "rm -f "; dum += fsphout; system(dum.c_str());
 FitsOutFile sfits(fsphout);
 cout<<"Writing Output Sphere Fits file"<<endl;
@@ -135 +157 @@
 }
 if(fsphoutw) {
-string dum = "rm -f "; dum += fsphoutw; system(dum.c_str());
+dum = "rm -f "; dum += fsphoutw; system(dum.c_str());
 FitsOutFile sfits(fsphoutw);
 cout<<"Writing Output Sphere Weight Fits file"<<endl;
@@ -141 +163 @@
 }
 if(fsphred) {
-string dum = "rm -f "; dum += fsphred; system(dum.c_str());
+dum = "rm -f "; dum += fsphred; system(dum.c_str());
 FitsOutFile sfits(fsphred);
 cout<<"Writing Reducted Sphere Fits file"<<endl;
@@ -147 +169 @@
 }
 if(fsphredw) {
-string dum = "rm -f "; dum += fsphredw; system(dum.c_str());
+dum = "rm -f "; dum += fsphredw; system(dum.c_str());
 FitsOutFile sfits(fsphredw);
 cout<<"Writing Reducted Sphere Weight Fits file"<<endl;