Up | Example 01 - Calculate Area | Example 02 - Calculate mean value | Example 03 - Manipulate the results from two LineFinders | Example 04 - Dynamic threshold | Example 05 - Auto Exposure | Example 06 - DrawLine | Example 07 - Overlays | Example 08 - Python Objects | Example 09 - Timing | Example 10 - Image Averaging | Example 11 - Image Resampling | Example 12 - Constant Contrast | Example 13 - Serial Communication | Example 14 - Python results | Example 15 - Making a result string | Example 16 - Running tools from a script | Example 17 - Image Manipulation using Python Scripting | Example 18 - Calculating the median angle | Example 19 - Iterating objects located by a blob | Example 20 - Resampling using non-linear calibration | Example 21 - Custom Scorpion Python extension | Example 22 - Accessing Image Pixels | Example 23 - Implementing a tcp/ip Socket Server | Example 24 - Setting ExternalReference from calculated four points | Example 25 - Rotating a reference around in a circle | Example 26 - Grabbing an image from an MOXA Video IP Server | Example 27 - Toolbox Switch | Example 28 - Color Matcher Iteration | Example 29 - Audio Notification | Example 30 - Windows MessageBox | Example 31 -  Client to tcp Socket Server | Example 32 -  Read / Write External Data from / to file | Example 33 - Changing a tool's ROI | Example 34 - Histogram Equalization | Example 35 - Robust Adam 6060 scripts | Example 36 - Bubblesort | Example 37 - Element Statistics | Example 38 - Saving 3D Image | Example 39 - Disabling Zoom in Image Windows | Example 40 - Filtering timeseries | Example 41 - Scorpion Watchdog keep system running | Example 42 - Binary Search | Example 43 - Creating an ordered pointcloud | Example 44 - UDP Socket Communication

 

  

 
Example 40 - Filtering timeseries

This example contains two use filters for arrays:

def MedianFilter(x,n):
  # median filter of array x, length n
  # filter behaves pretty well at the edges (shortens the length as needed)
  def med(x):
    x.sort()
    return x[len(x)/2]
  # force odd filter order (better filters then)
  if n%2==0:
    n += 1
  nh = n/2
  l = len(x)
  p = list(x)
  for i in range(nh):
    v = x[0:n-nh+i]
    p[i] = med(v)
  for i in range(nh,l-nh):
    v = x[i-nh:i-nh+n]
    p[i] = med(v)
  for i in range(l-nh,l):
    v = x[i-nh:l]
    p[i] = med(v)
  return p
def LowpassFilter(p,N,C):
  # N-tap running average, run C times
  l = len(p)
  if (N%2): N += 1 # N odd: zero phase
  NH = N / 2
  r = list(p)
  for c in range(C):
    lpad = [r[0] for i in range(NH)]
    rpad = [r[l-1] for i in range(NH)]
    rr = lpad + r + rpad
    r = [sum(rr[i:i+N])/N for i in range(l)]
  return r

 

 

Scorpion Vision Version XII : Build 646 - Date: 20170225
Scorpion Vision Software® is a registered trademark of Tordivel AS.
Copyright © 2000 - 2017 Tordivel AS.