[Thread Prev][Thread Next][Index]
Re: Locate good data point
Hi Jens-Olaf,
Try this demo
! define a 3-d variable
def axis/x=1:10:1 xax ; def axis/y=1:10:1 yax ; def axis/depth/z=1:5:1 zax
def grid/x=xax/y=yax/z=zax grd ; let v=x[g=grd]+y[g=grd]-z[g=grd]
! and suppose only values greater than 10 are good
let vg=if(v gt 10)then v
! look at the first layer or two to see what its like
shade/k=1 vg
...
! define a mask function for "good" values ...
let mask=if(vg)then 1
! ... and an x-edge function that locates the first good point
let xedge=mask[x=@loc:1]
! see if it is what you want
shade/k=1 vg ; list/k=1 xedge
! you could also define a y-edge function
let yedge=mask[x=@loc:1]
list/k=1 yedge
Hope this helps,
Mick
|--****--****-*---*---***--***--|____spillane@pmel.noaa.gov____|
|-*__---*-----*--*-*--*--*-*--*-|_SCIENCE APPLICATIONS SUPPORT_|
|--***--*-----*-*---*-***--***--|____EPIC/Ferret/PlotPlus______|
|-----*-*-----*-*****-*----*----|__Room 2070 Bldg#3 NOAA/PMEL__|
|-****---****-*-*---*-*----*----|____Phone_:_(206)526-6780_____|
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement