GuyRon1 / FFT-Hydrodynamic-sensing

This code is a fft algorithm used to analyze the vortex shedding freq
0 stars 0 forks source link

Frequency from PIV #8

Open GuyRon1 opened 1 month ago

GuyRon1 commented 1 month ago

https://github.com/GuyRon1/FFT-Hydrodynamic-sensing/blame/e2facaf5f6fed375f44eeb6f2dd795c17269c9e9/PIV_Folder/PIV%20FREQ%20DETECTION.ipynb#L1

Hi Alex, I have finished the code for calculating the frequency for the PIV measurements. I went over a lot of the files manually to select the thresholds and area of interest. In the end, I got 54 cycles in a 250 sec run which correlates to a frequency of 0.21 [Hz]. I will add more explanation to the ReaserchNotes on the overleaf file.

alexlib commented 1 month ago

cannot run without the data.

GuyRon1 commented 1 month ago

https://drive.google.com/drive/folders/1bRSO8nZyVJikAWYJFiQbYdb3sHZRdVIK?usp=drive_link Here is a link to a google drive with all the data. Would you like me to upload some of the data here?

alexlib commented 1 month ago

no, drive is enough

alexlib commented 1 month ago

I suggest to use .nc file instead of reading every time the large directory of .vec files. It's so much faster. please upload the .nc file to the Google drive and use it for the vortex detection.

GuyRon1 commented 1 month ago

When saving the .nc file and working with it, I get a Permission Error while loading vfField. Because I need to load vfField multiple times in the "for loop" " I decided not to work with the .nc file but to load all the files once and then only run the for loop for different time slices as I see fit. In issue #7 I understood from you that I need to delete the .nc file to load a new vfField manually.

alexlib commented 1 month ago

I'm talking only about the first reading

raw_data_nc_file = path / 'data.nc'
if raw_data_nc_file.exists():
    d = xr.load_dataset(raw_data_nc_file) #Loading the .netcdf file
else:
    d = io.load_directory(path) #Loading all files in the directory on a single xarray called "d"
    print(type(d)) 
    d.to_netcdf(path / 'data.nc') #Saves the xarray data from the files to a .netcdf file