← Back to branch summary

~ubuntu-branches/ubuntu/wily/eso-midas/wily-proposed

« back to all changes in this revision

Viewing changes to contrib/wavelet/src/wa_entr.c

  • Committer: Package Import Robot
  • Author(s): Ole Streicher
  • Date: 2014-04-22 14:44:58 UTC
  • Revision ID: package-import@ubuntu.com-20140422144458-okiwi1assxkkiz39
Tags: upstream-13.09pl1.2+dfsg
Import upstream version 13.09pl1.2+dfsg

Show diffs side-by-side

added added

removed removed

Lines of Context:
contrib/wavelet/src/wa_entr.c
 
1
/*===========================================================================
 
2
  Copyright (C) 1993-2009 European Southern Observatory (ESO)
 
3
 
 
4
  This program is free software; you can redistribute it and/or 
 
5
  modify it under the terms of the GNU General Public License as 
 
6
  published by the Free Software Foundation; either version 2 of 
 
7
  the License, or (at your option) any later version.
 
8
 
 
9
  This program is distributed in the hope that it will be useful,
 
10
  but WITHOUT ANY WARRANTY; without even the implied warranty of
 
11
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
12
  GNU General Public License for more details.
 
13
 
 
14
  You should have received a copy of the GNU General Public 
 
15
  License along with this program; if not, write to the Free 
 
16
  Software Foundation, Inc., 675 Massachusetts Ave, Cambridge, 
 
17
  MA 02139, USA.
 
18
 
 
19
  Correspondence concerning ESO-MIDAS should be addressed as follows:
 
20
        Internet e-mail: midas@eso.org
 
21
        Postal address: European Southern Observatory
 
22
                        Data Management Division 
 
23
                        Karl-Schwarzschild-Strasse 2
 
24
                        D 85748 Garching bei Muenchen 
 
25
                        GERMANY
 
26
===========================================================================*/
 
27
 
 
28
/******************************************************************************
 
29
**
 
30
**    UNIT
 
31
**
 
32
**    Version: 19.1
 
33
**
 
34
**    Author: Jean-Luc Starck
 
35
**
 
36
**    Date:  03/02/25
 
37
**    
 
38
**    File:  wa_entr.c
 
39
**
 
40
.VERSION
 
41
 090810         last modif
 
42
 
 
43
*******************************************************************************
 
44
**
 
45
**    DESCRIPTION  Extraction of a plane from a wavelet transform
 
46
**    ----------- 
 
47
**                 
 
48
******************************************************************************
 
49
**
 
50
**    DESCRIPTION   program which allows to replace a plane of a wavelet 
 
51
**    -----------   transform by an image
 
52
**
 
53
**    PARAMETER
 
54
**    ---------
 
55
**   
 
56
**     File_Name_Transform (keyword: IN_A) = 
 
57
**          input file name of the wavelet transform
 
58
**     File_Name_Imag (keyword: IN_B) = input file name of the image
 
59
**     File_Name_Transform (keyword: OUT_A) = 
 
60
**          output file name of the wavelet transform
 
61
**     Num_Plan (keyword: INPUTI) =
 
62
**          plane number to extract (Num_Plan = 1 .. Number_of_Planes)
 
63
**
 
64
******************************************************************************/
 
65
 
 
66
 
 
67
#include <stdio.h>
 
68
#include <math.h>
 
69
#include <string.h>
 
70
 
 
71
#include <midas_def.h>
 
72
 
 
73
#include "Def_Math.h"
 
74
#include "Def_Mem.h"
 
75
#include "Def_Wavelet.h"
 
76
 
 
77
#define VISU_PARAM FALSE
 
78
 
 
79
extern void wavelet_enter_plan_file();
 
80
 
 
81
 
 
82
/*****************************************************************************/
 
83
 
 
84
int main()
 
85
{
 
86
    string File_Name_Transform, File_Name_Imag, File_Trans_Out;
 
87
    int Unit; 
 
88
    int Num_Plan;
 
89
    int Stat, Null, Actvals, Buffer_Int, Maxvals, Felem;
 
90
 
 
91
    /* Initialisation */  
 
92
    SCSPRO("enter_plan");
 
93
 
 
94
    /* read the image file name */
 
95
    Felem = 1;
 
96
    Maxvals = 60;
 
97
    Stat = SCKGETC("IN_A", Felem, Maxvals, &Actvals, File_Name_Transform);
 
98
 
 
99
    /* read the  output image file name */
 
100
    Stat = SCKGETC("IN_B", Felem, Maxvals, &Actvals, File_Name_Imag);
 
101
 
 
102
    Felem = 1;
 
103
    Maxvals = 1;
 
104
    Stat= SCKRDI("INPUTI", Felem, Maxvals, &Actvals, &Buffer_Int, &Unit, &Null);
 
105
    Num_Plan = Buffer_Int;
 
106
 
 
107
    /* read the  output image file name */
 
108
    Felem = 1;
 
109
    Maxvals = 60;
 
110
    Stat = SCKGETC("OUT_A", Felem, Maxvals, &Actvals, File_Trans_Out);
 
111
 
 
112
#if VISU_PARAM
 
113
    { char Send[100];
 
114
    sprintf (Send,"Nom du fichier de la transformee   =%s\n",File_Name_Transform);
 
115
    SCTPUT(Send);
 
116
    sprintf (Send,"Nom de fichier image               =%s\n",File_Name_Imag);
 
117
    SCTPUT(Send);
 
118
    sprintf (Send,"Nom du fichier ondelette en sortie =%s\n",File_Trans_Out);
 
119
    SCTPUT(Send);
 
120
    sprintf(Send,"Numero du plan=%d\n",Num_Plan); 
 
121
    SCTPUT(Send);
 
122
    }
 
123
#endif
 
124
 
 
125
    wavelet_enter_plan_file (File_Name_Imag, File_Name_Transform, 
 
126
                             Num_Plan, File_Trans_Out);
 
127
 
 
128
 
 
129
   /* End */
 
130
   return SCSEPI();
 
131
}