希尔伯特谱分析matlab代码-PIV-PTV-oilcoatedBubbleTracking:结合PTV-PTV气泡轨迹跟踪和

  • V9_569691
    了解作者
  • 45.7MB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • VIP专享
    资源类型
  • 0
    下载次数
  • 2022-05-11 11:57
    上传日期
希尔伯特谱分析matlab代码PIV-PTV-油膜气泡跟踪 步骤进度 生成ptv piv分析文件夹(数据文件夹) 运行python create文件夹,其中包括所有情况。 在所有情况下运行piv处理。(09/19/2020) 检测上升周期。 消除背景噪声以进行气泡轨迹跟踪。 调用图像序列:“文件”-“导入”-“图像序列” 计算背景:“图像”-“堆栈”-“ z工程”-“投影类型:最小强度” 删除背景:“处理”-“图像计算器”-“操作:减” 锐化对比度:“处理-增强对比度”-“饱和像素:0.3%” 保存过滤后的序列:“文件”-“另存为”-“图像序列”-“格式:jpeg” PTV点检测,高斯核斑点检测。 面对某些帧的激光太亮的问题,在某些情况下很难检测到斑点形状,这会导致轨迹发散和混合。 通过使用预测器方法解决了09/24/2020问题,该方法使用最近附近的速度来预测下一帧气泡的位置。 对于每种情况,都需要通过应用imageJ局部最大值检测来重新检查气泡大小。 轨迹计算。 泡泡盒3失败。 在中的轨迹的俯视图,请使用我发送给您的密码访问数据。 如果您无法打开matlab图文件,请在matlab
PIV-PTV-oilcoatedBubbleTracking-master.zip
  • PIV-PTV-oilcoatedBubbleTracking-master
  • .vscode
  • settings.json
    39B
  • smoothspline.m
    391B
  • rt2traj.m
    1.5KB
  • sevenMoreCaseVelocity.fig
    274.1KB
  • VerticalV100Oil.fig
    26.1KB
  • trajp2traj.m
    1.3KB
  • MatlabImportData.mlx
    4.2KB
  • TransverseSpeedAmplitudeFiveMoreCase.fig
    290.2KB
  • Velocity10oil.fig
    460.5KB
  • TransversePositionAmplitude10oil.fig
    273.2KB
  • .gitattributes
    66B
  • TransversePositionAmplitudeAir.fig
    71.6KB
  • ImportData.m
    1.2KB
  • nearestneighborlinker.m
    5.4KB
  • freqSevenNewCase.fig
    6.6MB
  • freq100oil.fig
    11.2MB
  • freqAir.fig
    3.1MB
  • Velocity100oil.fig
    488.6KB
  • TransverseSpeed100oil.fig
    361KB
  • 100122233X-1Project.csv
    7.8KB
  • TransversePositionAmplitude100oil.fig
    267.6KB
  • TransverseSpeedAir.fig
    100.1KB
  • 10010226X-2Project.csv
    6.5KB
  • TransverseSpeedAmplitude10oilTime.fig
    273.6KB
  • VerticalV10Oil.fig
    26.1KB
  • PIV_PTV_combine_LiuHong.mlx
    105.2KB
  • 10010226X-4Project.csv
    8.9KB
  • PIVProcessing.mlx
    88.4KB
  • 1082230X-1Project.csv
    6.2KB
  • 10142230X-2Project.csv
    7.6KB
  • experiment_infoMatLab.csv
    1.8KB
  • Air_TopViewTraj.fig
    92.4KB
  • 1082230X-4Project.csv
    6.2KB
  • 1082210X-3Project.csv
    6.1KB
  • 10142220X-5Project.csv
    8.4KB
  • LiuHongProcessExtraAirBubble.py
    4.8KB
  • munkres.m
    5.4KB
  • VelocityAir.fig
    133.6KB
  • AirBubble-2Project.csv
    8.8KB
  • sevenMoreCaseTransverseSpeed.fig
    208.6KB
  • 1082220X-3Project.csv
    9.9KB
  • 1082230X-3Project.csv
    6.4KB
  • BestPlaneView10Oil.fig
    325.6KB
  • 10142220X-3Project.csv
    6.8KB
  • 100102230X-3Project.csv
    6.8KB
  • freqCompare.fig
    82.4KB
  • trajp2traj_nopositionfilter.m
    1.4KB
  • TransverseSpeedAmplitude10oil.fig
    496.6KB
  • simpletracker.m
    13.2KB
  • 100102226X-3Project.csv
    6.8KB
  • 100142230X-3Project.csv
    7.2KB
  • BestPlaneViewAir.fig
    92.7KB
  • TransverseSpeedAmplitudeSevenMoreCaseTime.fig
    185.2KB
  • Air_New_TopViewTraj.fig
    143.3KB
  • 100102226X-1Project.csv
    6.6KB
  • 10142220X-2Project.csv
    6.8KB
  • AirBubble-4Project.csv
    9.5KB
  • plot_line.m
    440B
  • VelocityAir.png
    67.7KB
  • SevenMoreCaseBestPlaneView.fig
    195.5KB
  • 100122233X-3Project.csv
    7.3KB
  • TransverseSpeedAmplitude100oil.fig
    471.8KB
  • LiuHongPTVProcess.mlx
    77.3KB
  • 1082220X-1Project.csv
    6.1KB
  • 10142230X-4Project.csv
    8.1KB
  • freq10oil.fig
    11.1MB
  • BubbleDetect.mlx
    131.2KB
  • 1082210X-2Project.csv
    5.7KB
  • 100102230X-2Project.csv
    9.7KB
  • BestPlaneView100Oil.fig
    338.4KB
  • 10_oil_TopViewTraj.fig
    330.8KB
  • freqAirNew.fig
    4.8MB
  • README.md
    13.5KB
  • TransverseSpeedAmplitudeSevenMoreCase.fig
    188.8KB
  • CreateDatafolder.py
    703B
  • FiveMoreCase.fig
    143.4KB
  • experiment_infoNew.csv
    1.7KB
  • trajProcessing.mlx
    482.4KB
  • 100142230X-2Project.csv
    11KB
  • WaveLet10oil.fig
    169.6KB
  • 10122230X-3Project.csv
    6.9KB
  • BestPlaneViewAirNew.fig
    142KB
  • experiment_info.csv
    1.7KB
  • 100102213X-3Project.csv
    6.4KB
  • 1082220X-2Project.csv
    6.5KB
  • AirBubble-1Project.csv
    6.3KB
  • 100102213X-2Project.csv
    6.5KB
  • TransverseSpeed10oil.fig
    345.7KB
  • Velocity100oil.png
    102.6KB
  • 10122230X-4Project.csv
    6.7KB
  • sevenMoreCaseVelocity.png
    72.6KB
  • .gitignore
    5B
  • TransverseSpeedAmplitude100oilTime.fig
    266.2KB
  • LiuHongPTVProcess.m
    2.3KB
  • 100102230X-1Project.csv
    7.7KB
  • preprocess.py
    2.1KB
  • LiuHongProcess.py
    5KB
  • hungarianlinker.m
    4.5KB
内容介绍
# PIV-PTV-oilcoatedBubbleTracking # Step progress - [x] build ptv piv analysis folder (data folder) Run python [script](https://github.com/lipilian/PIV-PTV-oilcoatedBubbleTracking/blob/master/CreateDatafolder.py) create [data](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) folder, which includes all the cases. - [x] run piv processing for all case.(09/19/2020) - [x] Detect the Raising period. - [x] Remove the background noise for bubble trajectory tracking. 1) Call image sequence: ‘file’-‘import’-‘image sequence’ 2) Compute background: ‘image’-‘stacks’-‘z project’-‘projection type:min intensity’ 3) Remove background: ‘process’-‘image calculator’-‘operatioin:substract’ 4) Sharpen the contrast: ‘process-‘enhance contrast’-‘saturated pixels:0.3%’ 5) Save filtered sequence: ‘file’-‘save as’-‘image sequence’-‘format:jpeg’ - [x] PTV point detection, Gaussian kernel blob detection. 1) Face the problem about the laser light too light for some frames that it hard to detect the blob shape for certain case, which cause trajectory divergent and mixing. <ul><li>- [x] </li> solved 09/24/2020 by using predictor method, which use the recently nearby velocity to predict the next frame bubble location. For each case, the bubble size need to be rechecked by applying imageJ local maximum detection. - [x] Trajectory computing. Bubble case 3 failed. 1) Top view of trajectory in [TopViewTraj folder](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn), please use the password I send you to access the data. If you can't open the matlab figure file Try these two lines code in matlab: eg. >>fig = openfig('Graph.fig'); >>set(fig, 'Visible', 'on') - [ ] PIV-PTV Combine. TODO-------------------------------------------------------- Using POD to calculate the freq of vortex, test adding mass effect, sloshing effect, for the trajectory frequency change. (10/28/2020) - [x] Trajectory analysis. **Start** (10/07/2020) 1) trajectory analysis - [x] PCA best fit plane finding. Update [TopViewTraj folder](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn), the projection 2D plane function are shown in x label for each figure, and plotted as red line in the figure.(10/07/2020). The red arrow represents the view point of the following figures for 2D fit plane. - [x] Side 2d view of trajectory on best fit plane. Update [SideViewTraj folder](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn), the side view is 3D trajectory which are viewed from norm vector of best fitting 2D Plane. (10/07/2020) 2) velocity analysis. - [x] Transverse velocity, vertical speed, and total speed. Velocity profile [VelocityProfile](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn), the u, v hortizontal velocity of camera view plane; w is the vertical velocity. Speed is the overall velocity. (10/07/2020) - [x] Transerse speed with vertical speed [TransverseSpeed](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn). **I found a very interesting behavior about damping. For higher oil fraction, the damping motion start disspiating, and higher viscousity oil will trigger it much earlier. My guess this motion of oil with in the bubble can balance the vortex induced vibration. Higher oil fraction have higher inertia and momantum. And lower viscousity oil can be easier to move inside the bubble, which can balance this vibration motion.** Based on this, we can study **actual damping motion, wavelet analysis, Hilbert transform.** - [x] vertical average speed with oil fraction for each case. Vertical velocity plot with error bar for two cases average. [VerticalVelocity](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) (10/09/2020) 3) frequency analysis - [x] Calualte instantaneous frequency based on first conditional spectral moment. Overall frequency seems to decrease when oil fraction increase. But pure air bubble also has low frequency. **Hard to conclude** Need discuss. (10/09/2020) [freqAnalysis](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) Update the freq vs oil fraction plot. (10/14/2020) [freqAnalysis](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) (shot time Fourier transform) - Compute spectrogram power spectrum P(t,f). - Estimate the joint time-frequency moment. 4) Amplitude analysis - [x] Position amplitude. (10/14/2020) FFT only work for periodic case, for dissipation case. May be we can use underdamping function fitting method to calculate the trend of amplitude. For some case the projection to 2D Plane method didn't work vary well to detection the main amplitude. Maybe due to the small disturbance from 3D, which is consistant with top view. (For example, some case the amplitude peak from frequency domain is not that obvious. I chekc their top view, it is much more like 3D trajectory.)[TransversePositionAmplitude](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) - [x] Fit underdamped harmonic oscillator to the model. Red line indicate the fitting curve, blue line is the original data. [DampingAmplitude](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) (10/19/2020) - [x] remake 5 airbubble tracking experiment. This time, use regionprops for tracking. (10/22/2022) [NewAirBubbleExperimentByBing](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) (10/23/2022) Excel sheet for vertical speed, transverse speed, freq information, Bing will add amplitude information. (10/28/2020) [ExcelSheetForInformation](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) 5) Five more cases, - [x] process 6 more cases for trajectory tracking. [FiveNewCase](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) 6) Excel sheet for 2D projected trajectory [Project2DExcel](https://uofi.box.com/s/sw9vvcwwh0pjfxlhspk6wisby87mabtn) # Overall process progress report | Outer needle | Inner needle | Initial oil volume (?L) | deq (mm) | foil | ID | deq | Processing Progress | StartFrame-EndFrame | | ------------ | ------------ | ----------------------- | -------- | ---- | ------------ | ----------- | -------------------- | ------------------- | | 8G | 22G | 10 | 4.04 | 0.08 | 1082210X-2 | 3.92869399 | <ul><li>- [x] finished </li> | 711-899 (The last frame is around 890 for bubble hitting the top edge, which may caused piv data missing about induced vortex) | | 8G | 22G | 10 | 4.04 | 0.08 | 1082210X-3 | 3.92755005 | <ul><li>- [x] Left seperate trajectory, try Link(09/25/2020) **Solved** (09/28/2020) </li> |1024-1349 | | 8G | 22G | 20 | 4.01 | 0.22 | 1082220X-1 | 3.8863055 | <ul><li>- [x] double check right camera 2 traj (09/25/2020) **Solved** (09/28/2020) </li> |129-449 | | 8G | 22G | 20 | 4.01 | 0.22 | 1082220X-3 | 3.88859096 | <ul><li>- [x] Left Camera 2 trajs (09/25/2020) **Solved** (09/28/2020) </li> |1056-1349 | | 8G | 22G | 30 | 4.08 | 0.31 | 1082230X-1 | 3.96551573 | <ul><li>- [x] finished </li> |190-449 | | 8G | 22G | 30 | 4.08 | 0.31 | 1082230X-4 | 3.97278205 | <ul><li>- [x] finished </li> |1558- 1799 | | 12G | 22G | 30 | 4.02 | 0.38 | 10122230X-3 | 4.17866011 | <ul><li>- [x] Right Camera has two trajectories. need further process **Solved** (09/28/2020) </li> |998-1349 | | 12G | 22G | 30 | 4.02 | 0.38 | 10122230X-4 | 4.28105654 | <ul><li>- [x] Right Cmaera has two trajectories
评论
    相关推荐
    • Matlab-Python-Pictures-master.zip
      利用matlab绘制节点,以及他们之间的关系,利用Excel表格 adjoin.xlsx,描述点之间的关系,1表示连接,0表示不连接,pos.xlsx表示要绘制的屏幕上的位置
    • MATLABPython接口
      介绍MATLABPython接口的多种方法,包括MATLABPython的相互调用以及将MATLAB代码打包成组件然后在Python中集成。
    • matlab导入mat代码-MATLAB-PYthon-Combine-Code:MATLAB-PYthon-Combine-C
      MATLAB的官方网站,有关于如何为Python安装MATLAB API的 2. 相关代码 2.1 Python当中的 import matlab.engine eng = matlab.engine.start_matlab() res = eng.GetRowColOfImage('test_img.jpg') print res 2.2 ...
    • Matlab Python UDP Socket:Matlab中使用Python Socket实现UDP通信的基本示例-matl
      这个基本示例提供了一个使用 python 套接字实现的 UDP 通信接口。 我体验过这种方法在时间关键应用程序中运行比 matlab/java UDP 套接字更... 据我所知,所有 python 依赖项都应该由 Matlab 附带的 python 版本解决。
    • matlab2python
      git clone https://gitee.com/zongfang/matlab2python.git python matlab2python.py file.m -o file.py
    • matlabfunction代码生成-python:Python
      matlab function 代码生成 install python 在window下下载python安装包并安装(会自动安装pip和setuptools,如没有可以然后安装:python get-pip.py),然后添加脚本(如C:\Python27\Scripts或C:\Python36\Scripts)到...
    • matlab特征点代码-matlab-python:MatlabPython
      matlab特征点代码 matlab-python #1.将matlab代码转成python #2.特征点匹配图
    • TripInference:用于行程推断的MATLABPython代码
      该存储库包含用于Trip Trip的MATLAB / Python代码。 TripParseTest.m:演示TripParse.m功能的MATLAB脚本 TripParse.m:解析行程的主要功能 getActivities.py:用于活动匹配和Google查询的Python库 mainActivities....
    • Codes:MatlabPython脚本
      MatlabPython的脚本
    • matlab-with-python:结合使用MATLABPython的示例
      结合使用MATLABPython的示例 该演示演示了如何从MATLAB调用Python,反之亦然,以及如何从PythonMATLAB Compiler SDK)调用打包的MATLAB代码。 主要示例调用了开放天气API: : 没有提供访问密钥(ha),但是您...