Changes between Version 22 and Version 23 of Tutorial/ParticleImageVelocimetry


Ignore:
Timestamp:
Jan 15, 2015, 4:20:25 PM (10 years ago)
Author:
vaillant1p
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • Tutorial/ParticleImageVelocimetry

    v22 v23  
    103103To improve the results, come back to the GUI '''CIV''', and follow these successive steps, corresponding to a sequence of operations.
    104104
    105 ''' * Optimizing the Civ1 parameters''': Select the '''[CIV1]''' check box so the corresponding parameters show up. Improve the spatial resolution by selecting smaller correlation boxes in the '''civ1''' menu, for instance [Bx],[By] to 19 and 13 (image pixels). This is possible because of the good image quality and high particle density. The use of a smaller box in y allows to optimize the resolution in this direction, to deal with the strong vertical shear. It is now possible to adjust the search range, using knowledge on extremal velocities, see histograms displayed by '''uvmat''' in phys coordinates. We introduce estimated bounds on each velocity component, [min] to -2 and [max] to 6 for ’u’ and (-3, 3) for ’v’, and press the button '''[Search Range]'''. The optimum search ranges and shifts are now displayed. In consistency with higher resolution, we set the grid mesh to [Dx] and [Dy] to 10. Finally select the '''Mask option''' : the mask name should be displayed in the corresponding edit box, else a browser leads to the appropriate mask file.
     105'''* Optimizing the Civ1 parameters''': Select the '''[CIV1]''' check box so the corresponding parameters show up. Improve the spatial resolution by selecting smaller correlation boxes in the '''civ1''' menu, for instance [Bx],[By] to 19 and 13 (image pixels). This is possible because of the good image quality and high particle density. The use of a smaller box in y allows to optimize the resolution in this direction, to deal with the strong vertical shear. It is now possible to adjust the search range, using knowledge on extremal velocities, see histograms displayed by '''uvmat''' in phys coordinates. We introduce estimated bounds on each velocity component, [min] to -2 and [max] to 6 for ’u’ and (-3, 3) for ’v’, and press the button '''[Search Range]'''. The optimum search ranges and shifts are now displayed. In consistency with higher resolution, we set the grid mesh to [Dx] and [Dy] to 10. Finally select the '''Mask option''' : the mask name should be displayed in the corresponding edit box, else a browser leads to the appropriate mask file.
    106106
    107 ''' * FIX1''': Select the ’'''FIX1'''’ operation, which eliminates some false vectors using several criteria. Use the default parameters.
     107'''* FIX1''': Select the ’'''FIX1'''’ operation, which eliminates some false vectors using several criteria. Use the default parameters.
    108108
    109 ''' * PATCH1''': Select the ’'''PATCH1'''’ operation, to interpolate the vectors on a regular grid and calculate spatial derivatives. Select a high resolution [nx]=100, [ny]=30, number of grid points in which the velocity will be interpolated. These values correspond to dx=dy=10 pixels for an image 1000x300. Choose the default value 10 for the smoothing parameter ’rho’. You can later try different values, the smoothing effect increasing with rho. Keep the default values for the other parameters.
     109'''* PATCH1''': Select the ’'''PATCH1'''’ operation, to interpolate the vectors on a regular grid and calculate spatial derivatives. Select a high resolution [nx]=100, [ny]=30, number of grid points in which the velocity will be interpolated. These values correspond to dx=dy=10 pixels for an image 1000x300. Choose the default value 10 for the smoothing parameter ’rho’. You can later try different values, the smoothing effect increasing with rho. Keep the default values for the other parameters.
    110110
    111 ''' * CIV2''': Select the ’'''CIV2'''’ operation to improve the correlation results, using the information on local image deformation, provided by the previous knowledge on velocity spatial derivatives (calculated in patch1). Use a finer grid dx= dy=5 than for civ1. The spatial resolution can be slightly improved by decreasing the correlation box, using for instance Bx,By=(15,11). The search range is determined automatically using the prior knowledge on velocity, obtained with the civ1 and patch1 operations. Note that ’civ2’ corresponds to a new measurement from the images, the previous civ1 and patch operations being used only as an initial guess for the search of optimal correlations.
     111'''* CIV2''': Select the ’'''CIV2'''’ operation to improve the correlation results, using the information on local image deformation, provided by the previous knowledge on velocity spatial derivatives (calculated in patch1). Use a finer grid dx= dy=5 than for civ1. The spatial resolution can be slightly improved by decreasing the correlation box, using for instance Bx,By=(15,11). The search range is determined automatically using the prior knowledge on velocity, obtained with the civ1 and patch1 operations. Note that ’civ2’ corresponds to a new measurement from the images, the previous civ1 and patch operations being used only as an initial guess for the search of optimal correlations.
    112112
    113113Then select ’'''FIX2'''’ and ’'''PATCH2'''’ with the same parameters as ’'''FIX1'''’ and ’'''PATCH1'''’.
    114114
    115 ''' * Running the calculation''': Press '''[RUN]''' to run the calculation. The results are stored in a new subdirectory, CIV_1, so the previous results are not erased (you can also set the name of the subdirectory in the GUI, by the edit boxes '''[SubDirCiv1]''' and '''[SubDirCiv2]'''. The existing subdirectories are listed above.
     115'''* Running the calculation''': Press '''[RUN]''' to run the calculation. The results are stored in a new subdirectory, CIV_1, so the previous results are not erased (you can also set the name of the subdirectory in the GUI, by the edit boxes '''[SubDirCiv1]''' and '''[SubDirCiv2]'''. The existing subdirectories are listed above.
    116116
    117117The status of the calculation is displayed in a new window which refreshes automatically. Close it to avoid blocking of new Matlab operations. This status window can be opened again by pushing the button '''[STATUS]''' in the GUI '''civ'''. Do not close civ until the calculation is finished.