Changes between Version 1 and Version 2 of Tutorial/ComparaisonNumSolution


Ignore:
Timestamp:
Feb 19, 2015, 3:13:10 PM (6 years ago)
Author:
vaillant1p
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Tutorial/ComparaisonNumSolution

    v1 v2  
    33= [wiki:Tutorial] / Comparaison with a Numerical Solution =
    44
     5== Introduction ==
     6Using synthetic images allows us to compare our PIV solution with the numerical one used to create the images. It is this way easy to determine the errors made during the PIC calculations.
     7
     8The first part in this tutorial shows how to analyse perfect noise-free images and compare them with the synthetic solution.
     9The next part will be focusing on creating noise in these same synthetic images and analysing them.
     10
     11== Synthetic Images ==
     12In the folder 'AnalyticFlows', select one of the images from the 'a3b' series. It shows a diagonal uniform flow running from the bottom left corner to the top right one. '''[RUN]''' a PIV on the images series only using '''[CIV1]''', '''[FIX1]''' and '''[PATCH1]''' with the default parameters.
     13Then '''[Open]''' one of the vector files (.civ) created. We need to substract the numerical solution to this image in order to keep only the errors between our calculation and the theorical one. Select the mall red square'''[SubField]''' under the '''[Input]''' table to open a second file in '''uvmat'''. Open the netcdf file of the synthetic solution 'a3b_sol_num.nc'. The GUI '''get_field''' pops up, select 'vectors' in the '''[!FieldOptions]''' table and 'U' and 'V' in the right and left tables under '''[Ordinate]'''. The '''[Coordinates]''' should still remain as 'X' and 'Y'. Click on '''[OK]'''. An 'ERROR' window pops up, select '''[OK]''' and close the GUI '''set_object]''' which automatically popped up.
     14
     15On the left hand side of '''uvmat''' select the '''[Tranform]''' function 'diff_vel'. If this one does not appear on the list, click on 'more...' and select 'diff_vel.m' from the folder 'UVMAT/transform_field'. The difference between the two vectors should run automatically and the result appears in '''uvmat'''. Selecting '''[<-]''' and '''[+>]''' we can go through all the series observing the difference result on screen. A single file will be easier to observe the errors between the analytic and synthetic solutions. As we do not want to focus on the borders of the images we are going to create a '''[Projection object/rectangle]''' of '''[Max]''' (100, 100) with a center at x=128 and y=128. You can change this parameters once you drew the rectangle in '''set_object'''.
     16
     17Then in order to obtain a single result file, go to '''series''' ('''[RUN/field series]'''). The two files should appear in the '''[Input Files]''' table. Select the '''[Action]''' function 'aver_stat'. A 'WARNING' window pops up ask you if you want to use a projection object, select '''[OK]''', the rectangle newly created should appear on the left hand corner under '''[Projection Object]'''. Another 'INPUT_TXT' window pops up to set the bin size for the histograms, enter '0.005' instead of 'auto' and select '''[OK]'''. Make sure that the i indices go from 1 to 41, that the '''[Field Transform]''' on the left hand side is 'diff_vel' and that the '''[!VelType]''' in '''[Fields]''' (right hand side) is 'civ1' instead of '*'. It is possible to change the name of the file under '''[output !SubDir]''', otherwise '''[RUN]''' the calculation.
     18
     19
     20
     21
     22
     23Once it is done, the GUI '''get_field''' pops up again. Keep '1D plot' in the '''[!FieldOptions]''' table and select 'UHisto' or 'Vhisto' in the '''[Ordinate]''' one. The '''[Coordinates]''' should change according to the input in the '''[Ordinate]''' ('U' if 'UHisto' is selected, same for 'V' and 'VHisto'). Select '''[OK]'''.
     24
     25A curve, close to a Gaussian, (gathering all the vectors from the difference between all the analytic and synthetic files of the series) appears in '''uvmat'''. On the right hand corner a table shows all the parameters of this curve. In the scroll bar '''[Field/get_field]''' it is possible to get back to the GUI '''get_field''' and change the '''[Ordinate]'''. You can also '''[Export]''' the curve '''[/in new figure]''' or '''[/on existing axis]'''.
     26
     27This is useful to determine and optimize some PIV parameters. The figure below shows the RMS, !FirstCentile, !FirstDecile, !LastDecile and !LastCentile of PIV calculations with different '''[CorrBoxSize]'''.
     28
     29
     30
     31
     32
     33== Adding Noise ==
     34However as synthetic images do not exactly relate reality it is necessary to add noise to them in order to obtain a more realystic result. '''[Open]''' once again one of the 'a3b' images, then select the '''[Tranform]''' function 'ima_noise' (again if it does not appear on the list, click on 'more...' and select 'ima_noise.m' from the folder 'UVMAT/transform_field'). Some multiplicative noise is then added to the image.
     35
     36As it is not possible to keep the '''[Transform]''' function 'ima_noise' while running a PIV we need to create another image serie with noise. To do that '''[RUN/field series]''' and select the '''[Action]''' function 'merge_proj'. Do not add any projection object, make sure that the '''[Field Transform]''' function is set as 'ima_noise' and that the calculation will run for all the images (from indice i=1 to i=41) and '''[RUN]'''. A new folder '.mproj' is then created with all the ''noisy'' images.
     37
     38
     39
     40It is now possible to '''[Run]''' a '''[/PIV]''' on this serie and by following all the steps from part1 (Synthetic Images) to obtain the difference between the analytical and noise solution and the synthetic one. The curve is not as smooth as in part1, you can indeed notice the effect of noise in a PIV calculation.
     41
     42
     43
     44