Another try at using median (actually 45%-ile) for base level in waterfall.

git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/map65@577 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
This commit is contained in:
Joe Taylor 2007-11-12 14:55:30 +00:00
parent 30c1a0427f
commit 7c8adacb71
2 changed files with 12 additions and 4 deletions

View File

@ -1,4 +1,4 @@
#---------------------------------------------------------------------- MAP65 #--------------------------------------------------------------------- MAP65
# $Date$ $Revision$ # $Date$ $Revision$
# #
from Tkinter import * from Tkinter import *

View File

@ -17,7 +17,7 @@ subroutine spec(brightness,contrast,ngain,nspeed,a,a2)
! Could save memory by doing the averaging-by-7 (or 10?) of ss5 in symspec. ! Could save memory by doing the averaging-by-7 (or 10?) of ss5 in symspec.
include 'spcom.f90' include 'spcom.f90'
real s(NFFT,NY) real s(NFFT,NY),savg2(NFFT),tmp(NFFT)
include 'gcom1.f90' include 'gcom1.f90'
include 'gcom2.f90' include 'gcom2.f90'
include 'gcom3.f90' include 'gcom3.f90'
@ -45,6 +45,15 @@ subroutine spec(brightness,contrast,ngain,nspeed,a,a2)
enddo enddo
enddo enddo
enddo enddo
call zero(savg2,NFFT)
do j=1,nlines
call add(savg2,s(1,j),savg2,NFFT)
enddo
ia=0.08*NFFT
nn=0.84*NFFT
call pctile(savg2(ia),tmp,nn,45,base)
base=base/nlines
newpts=NX*nlines newpts=NX*nlines
do i=newpts+1,NX*NY do i=newpts+1,NX*NY
@ -55,7 +64,7 @@ subroutine spec(brightness,contrast,ngain,nspeed,a,a2)
gain=40*sqrt(nstep(nspeed)/5.0) * 5.0**(0.01*contrast) gain=40*sqrt(nstep(nspeed)/5.0) * 5.0**(0.01*contrast)
gamma=1.3 + 0.01*contrast gamma=1.3 + 0.01*contrast
offset=(brightness+64.0)/2 offset=(brightness+64.0)/2
fac=20.0/nadd fac=3.2/base
foffset=0.001*(1270+nfcal) foffset=0.001*(1270+nfcal)
nbpp=(nfb-nfa)*NFFT/(96.0*NX) !Bins per pixel in wideband (upper) waterfall nbpp=(nfb-nfa)*NFFT/(96.0*NX) !Bins per pixel in wideband (upper) waterfall
fselect=mousefqso + foffset fselect=mousefqso + foffset
@ -73,7 +82,6 @@ subroutine spec(brightness,contrast,ngain,nspeed,a,a2)
do j=nlines,1,-1 !Reverse order so last will be on top do j=nlines,1,-1 !Reverse order so last will be on top
do i=1,NX do i=1,NX
k=k+1 k=k+1
n=0 n=0
x=0. x=0.
iia=(i-1)*nbpp + ii0 + 1 iia=(i-1)*nbpp + ii0 + 1