# Author: G. Hennenfent
#         Seismic Laboratory for Imaging and Modeling
#         Department of Earch & Ocean Sciences
#         The University of British Columbia
#         
# Date  : January, 07

from rsf.proj import *
import os

Flow('model',None,
     'sfsigmoid n1=300 n2=400')
Result('model','sfgrey title="Model" ')

Flow('data','model','sfnoise var=1e-4 seed=12010')
Result('data','sfgrey title="Data" ')

Flow('res','data',
     '''sffdct nbs=6 nba=32 ac=1 adj=n |
     sfthr thr=4e-3 |
     sffdct nbs=6 nba=32 ac=1 adj=y'''%vars() )
Result('res','sfgrey title="Denoised image" ')

End()