parouby
/
acoused
mirror of https://forgemia.inra.fr/theophile.terraz/acoused
1
0
Fork 0
Code Issues Packages Projects Releases Activity

ubsediflow raw extract is cleaned from useless test

dev-brahim
brahim 2025-03-05 13:36:37 +01:00
parent 77f28462e6
commit 9fcbe0e306
1 changed files with 2 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
Model/AquascatDataLoader.py
View File

@ -721,14 +721,14 @@ class RawAquascatData:
elif pktType == 41: elif pktType == 41:
chan = st.unpack("H", f.read(2))[0] + 1 chan = st.unpack("H", f.read(2))[0] + 1
# Increase the Index # Increase the Index
AbsIndex[chan-1] = AbsIndex[chan-1] + 1 AbsIndex[chan-1] = AbsIndex[chan-1] + 1
Tmp = [st.unpack( Tmp = [st.unpack(
"h", f.read(2))[0] for ui in np.arange( "h", f.read(2))[0] for ui in np.arange(
1, 2 * AbsNumBins[chan-1]+1)] 1, 2 * AbsNumBins[chan-1]+1)]
sss = [x / 32768 for x in Tmp] sss = [x / 32768 for x in Tmp]
AbsData[:, int( AbsData[:, int(
AbsIndex[chan-1]), int(chan-1)] = np.abs(np.array(sss[0::2])+complex(0,1)*np.array(sss[1::2])) AbsIndex[chan-1]), int(chan-1)] = np.abs(np.array(sss[0::2])+complex(0,1)*np.array(sss[1::2]))
@ -1097,19 +1097,3 @@ class MeanAquascatProfile:
file.close() file.close()
# ------------------------- Test --------------------------------------#
start_time = time.time()
if __name__ == "__main__":
# path1 = r'C:\Users\vergne\Documents\Donnees_aquascat\2017_juillet - experience cuve sediments fins\2017_07_19 - mercredi eau claire\20170719114700.aqa'
path1 = r'//home/brahim.moudjed/Documents/3 Software_Project/river_inversion_project/Data/Aquascat data test/20171213135800.aqa'
data1 = RawAquascatData(path1)
# path2 = r'C:\Users\vergne\Documents\Donnees_aquascat\2017_juillet - experience cuve sediments fins\2017_07_19 - mercredi eau claire\20170719114700.aqa.txt'
path2 = r'//home/brahim.moudjed/Documents/3 Software_Project/river_inversion_project/Data/Aquascat data test/20171213135800.txt'
data2 = RawAquascatData(path2)
print(data1.PingRate)
print(data2.PingRate)
print("Computational time: %.2f min" %((time.time() - start_time)/60) )