1
! @(#)fripple.prg 19.1 (ES0-DMD) 02/25/03 14:29:18
2
!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
6
! MIDAS COMMAND PROCEDURE : FRIPPLE.PRG
8
! D. Baade, ST-ECF, Garching version 1.0 291085
9
! version 1.1 P. Ballester 220291
14
! execute the command :
16
! FILTER/RIPPLE inimage outimage period [start,end]
18
! inimage: input frame
19
! outimage: output frame
20
! period: ripple period
21
! start,end: first and last pixel of range over which analysis is to be performe
23
!-------------------------------------------------------------------------------
25
DEFINE/PARAM P1 ? IMA "Enter input image:"
26
DEFINE/PARAM P2 ? IMA "Enter output image:"
27
DEFINE/PARAM P3 ? NUMBER "Enter period:"
28
DEFINE/LOCAL DELTA/R/1/1 0.0
30
WRITE/KEYW HISTORY "FILTER/RIPPLE "
33
WRITE/KEYW INPUTR/R/1/1 {P3}
35
IF P4(1:1) .EQ. "?" THEN
37
COPY/DK {P1} NPIX/I/1/1 INPUTI/I/2/1
39
WRITE/KEYW INPUTI/I/1/2 {P4}
42
! Check consistency of input data 910222
44
IF INPUTI(2) .LE. INPUTI(1) THEN
45
WRITE/OUT "Condition End > Start not fulfilled"
47
ELSEIF INPUTR(1) .LT. 1. THEN
48
WRITE/OUT "Period must be larger than 1. pixel"
52
DELTA = INPUTI(2) - INPUTI(1)
54
IF DELTA .LE. INPUTR(1) THEN
55
WRITE/OUT "Condition (End-Start)>Period not fulfilled"