Scribd
Upload
119 views

Tutorial Openpiv - Version 01

The document provides an overview of OpenPIV, an open source Particle Image Velocimetry (PIV) analysis software. It describes the graphical user interface, how to load PIV data, the available analysis parameters including interrogation window size, signal-to-noise ratio filters, and more. It also outlines the output files generated, how to convert the velocity units from pixels/time to physical units like m/s, and references related materials.

Uploaded by

Abhijit Kushwaha
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
119 views

Tutorial Openpiv - Version 01

The document provides an overview of OpenPIV, an open source Particle Image Velocimetry (PIV) analysis software. It describes the graphical user interface, how to load PIV data, the available analysis parameters including interrogation window size, signal-to-noise ratio filters, and more. It also outlines the output files generated, how to convert the velocity units from pixels/time to physical units like m/s, and references related materials.

Uploaded by

Abhijit Kushwaha
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 7

Tutorial OpenPIV - version 01

ISTA - Institut für Strömungsmechanik und Technische Akustik


TU - Berlin, Germany

September 6, 2012

Contents
1 Inroduccion 2

2 Graphical User Interface GUI 2

3 Load Data 2

4 Analysis parameters 3
4.1 Interrogation window size . . . . . . . . . . . . . . . . . . . . . . 3
4.2 Spacing/Overlap . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
4.3 Signal-to-Noise ratio and filters . . . . . . . . . . . . . . . . . . . 3
4.4 Scale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
4.5 Outlier filter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
4.6 Jump . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
4.7 Select and reset Region Of Interest(ROI) . . . . . . . . . . . . . . 4
4.8 Preprocessing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5

5 Start and stop 5

6 Output files 5

7 Units convert 6

1
1 Inroduccion
OpenPIV is an open source Particle Image Velocimetry (PIV) analysis software
written in Python, MATLAB and C++. The aim of this document is to de-
liver a fast overview of OpenPIV. The next sections will refer to the selectable
parameters, the ouput files and other features of OpenPIV. Some basics of the
PIV technique will be mentioned as well. All the distributions of OpenPIV and
other related material can be found at [1, 2].

This tutorial runs OpenPIV on MATLAB, nevertheless the same applies for
other distributions of OpenPIV. The PIV-Data used as example in this tutorial
was extracted from [3].

2 Graphical User Interface GUI


Figure 1 depicts the Graphical User Interface (GUI) of OpenPIV after loading
some example data.

Figure 1: GUI OpenPIV

3 Load Data
The first step before to run the PIV-Analysis is to load the PIV-Data. On the
top left the menu bar is located. Under the drop-down menu “File” the option
“Load” is available, either this option or its shortcut Ctrl+L opens a window
where to select the data to be loaded. Once the PIV-Images are loaded the but-
tons  and  at the bottom jump between the loaded images, the indicator at
the right-bottom displays the number of the currently shown image.

Note that at least two PIV-Images are needed for a PIV-Analysis.

2
4 Analysis parameters
In this section the selectable parameters on the right side of the GUI are de-
scribed.

4.1 Interrogation window size


The PIV-Image will be divided in interrogation windows. The width and the
height control the size in pixels of the interrogation windows. Since the analysis
is based on statistical a large interrogation windows are more robust against
background noise and outlier vectors.

4.2 Spacing/Overlap
The horizontal and vertical parameters control the size of the correlation win-
dows and set the resolution of the two components of the velocity field ux uy .
In figure 2 is depicted an example.

Figure 2: Image of 256x256 pixels divided by 32x32 interrogation window.


Space/Overlap: red 16x16, green 32x32, blue 64x64 and yellow 128x128

4.3 Signal-to-Noise ratio and filters


In OpenPIV are implemented two types of Signal-to-Noise ratio (S/N): peak to
the mean and peak to the second peak. In the drop-down menu 1 corresponds
to peak to the mean and 2 to peak to the second peak, option 1 is faster.

The parameter S/N value is the threshold between the correct/wrong vec-
tors according to the S/N type. There are no default values, a try-out has to
be perform in order to find out the best value.

OpenPIV applies global and local filtering. By global filtering the obvious
outlier vectors will be supposedly removed, i.e. the vectors which length is larger
than the mean of the flow field plus 3 times its standard deviation. These are
global outliers in the statistical sense.

In order to find local vector outliers, i.e. vectors that are dissimilar from the
close neighbors, a local filtering is performed on small neighborhoods of vectors,

3
e.g. 3x3 or 5x5. Typically there are about 5 % of erroneous vectors, these are
removed and later the missing values are interpolated from the neighbor vector
values.

4.4 Scale
The scale value is a factor that will multiply the resulting vectors. It helps with
the visualization of the vectors but it has to be taken into account when the
units are transform to physical units.

4.5 Outlier filter


The outlier filter value is the threshold of the global outlier filter and indicates
how many times the standard deviation of the whole vector field is exceeded
before the vector is considered as outlier.

4.6 Jump
The jump value indicates how many PIV-Images are evaluated in the anaylisis.
The count starts at the next image from the current one, i.e. if jump is equal to
n then OpenPIV takes the current image plus the next n PIV-Images for one
analysis. Note that at least two PIV-Images are needed for the analysis, what
makes jump values less than one pointless.

4.7 Select and reset Region Of Interest(ROI)


The option Select ROI provides the possibility of evaluate only part of the whole
PIV-Image. The option Reset ROI release the selected region. In figure 3 one
example is shown.

Figure 3: Results for a selected region

4
4.8 Preprocessing
The preprocessing option applies transformation-image function to the PIV-
Images. The current MATLAB version uses the function adapthisteq which
enhances the contrast of images. For more details see MATLAB help.

5 Start and stop


The start and stop buttons either start or stop the analyze process.

6 Output files
The first vectors results are depicted on the PIV-Images in the GUI, as figure
3 and 4 illustrate. In the folder containing the PIV-Data OpenPIV saves three
lists of files: no filtered results dataName noflt.txt, filtered results after global
and local filters dataName flt.txt and final results after filtering and interpola-
tion dataName.txt.

The output file has 5 columns of data, i.e. x y ux uy S/R. The first and
second correspond to the pixel coordinates, the third and fourth to the two ve-
locity vector components and the fifth column is the S/N ration. Note that the
value is different if the user chose peak to second peak or peak to mean ratio as
S/N type. The value of peak to second peak or peak to mean ratio is stored for
further processing.

The number of output files will be the number of loaded PIV-Images minus
the selected jump value.

Figure 4: Output results

5
7 Units convert
After the program is executed, the velocity vector units of ux and uy are not
in physical units [m/s]. OpenPIV results are delivered in [scale × (pixels/∆t)],
with scale as the selected scale value and ∆t as the time interval between two
images. In order to convert the units to [m/s] is needed:

• the vector components ux , uy in [scale × (pixels/∆t)]

• the ∆t value in s
• the dimensionless scale value
• the relation between meters and pixels in that image [m/pixels], i.e. the
image size and the number of pixels.

The following example demonstrates a particular case of unit conversion.

• ux = 5, uy = 0 in [scale × (pixels/∆t)]
• the ∆t = 0.5 in s

• scale value = 2
• image size 0.25 x 0.25 m and 256 x 256 pixels, i.e. the relation between
meters and pixels is 0.25/256 [m/pixels].
• the conversion of ux from [scale × (pixels/∆t)] to [m/s]:
ux m 5 0.25
× rel. = × = 0.0049m/s
scale × ∆t pixels 2 × 0.5 256

6
References
[1] R. Gurka A. Liberzon and Z. Taylor. OpenPIV Home Page. http://www.
openpiv.net, 2009.
[2] G.A. Kopp A. Liberzon Z.J. Taylor, R. Gurka. Long-Duration
Time-Resolved PIV to Study Unsteady Aerodynamics. Instrumenta-
tion and Measurement, IEEE Transactions on, 59(12):3262–3269, Dec.
2010. http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=
5464317&isnumber=5609237.
[3] J. Carlier. Second set of fluid mechanics image sequences. European Project
’Fluid image analisys and description’ (FLUID) – http://www.fluid.irisa.fr/,
2005.

You might also like