Commit c22b192d authored by MAURICE's avatar MAURICE
Browse files

correction indentation

parent d2f3cb56
......@@ -153,13 +153,13 @@ Les variables suivantes sont à initialiser:
"""
# Chemins windows
path_evt='H:/WBAT/CasinoEchosondeBis_Leg2.txt'
path_config ='E:/WBAT/configMovies3Dter_chgt_offsets'
path_evt='H:/EchoSonde2021-1/CasinoManuelEchoSonde2021-1.txt'
path_config ='H:/Config_M3D_bis/configsAnalyseM3D_Echosonde/EIlay/'
thresholds = [-80]
path_config = path_config + "/" + str(thresholds[0])
path_hac_survey='E:/WBAT/HERMES/LEG2/SBES_MBES/ECHOSONDE_BIS/RUN127/'
path_hac_survey='D:/Hac_Echosonde2021/SBES_MBES/ECHOSONDE2021-1/RUN067/'
type_campagne='echosonde_bis'
type_campagne='echosonde'
if type_campagne=='echosonde':
analyse_couple_sondeur='non'
analyse_par_saison='oui'
......@@ -184,11 +184,11 @@ clustering = True
methode_EI = 'netcdf'
#Chemin de sauvegarde
fpath='H:/WBAT/Result_leg2_bis/Leg2/'
fpath='H:/Rearchivage/ECHOSONDE2021/LEG1/'
#fpath='H:/Rearchivage/ECHOSONDE2019/LEG2/120/Echosonde2019_netcdf/Result/'
# Sélection des radiales, transducteurs et fréquences d'interêt
ind_decoup=list(range(0,721,1))
ind_decoup=list(range(0,4,1))
indice_freq_max=342
if analyse_couple_sondeur=='oui':
......@@ -213,15 +213,15 @@ if analyse_couple_sondeur=='oui':
freq_nominale=120000
freq_custom =[freq_nominale]
else:
nbfreqs = 1
indices_transduc_choisis = np.array([0])
nb_transduc=1
freq_custom =[70000]
nbfreqs = 2
indices_transduc_choisis = np.array([0,1])
nb_transduc=2
freq_custom =[120000,200000]
# Classification de la donnée
nombre_de_categories_sv=3 #classif par portion de radiale
nombre_de_categories_sv_campEntiere=5 #classif par recombinaison à l'échelle de la campagne entière
pas=200
pas=1000
methode_clustering="LLE_KMeans"
if methode_clustering=="LLE_KMeans":
nombre_composantes=4
......@@ -248,25 +248,25 @@ offset=2
# dir_vent_vrai,
# ) = ic.infos_radiales(C)
# C, Entete, index_events, date_heure = ic.import_casino(path_evt)
# (
# datedebut,
# datefin,
# heuredebut,
# heurefin,
# nameTransect,
# vit_vent_vrai,
# dir_vent_vrai,
# ) = ic.infos_radiales_echosonde(C)
C, Entete, index_events, date_heure = ic.import_casino_echosonde_bis(path_evt)
C, Entete, index_events, date_heure = ic.import_casino(path_evt)
(
datedebut,
datefin,
heuredebut,
heurefin,
nameTransect,
) = ic.infos_radiales_echosonde_bis(C)
vit_vent_vrai,
dir_vent_vrai,
) = ic.infos_radiales_echosonde(C)
# C, Entete, index_events, date_heure = ic.import_casino_echosonde_bis(path_evt)
# (
# datedebut,
# datefin,
# heuredebut,
# heurefin,
# nameTransect,
# ) = ic.infos_radiales_echosonde_bis(C)
# supression des variables qui ne sont plus utiles
del C, Entete, date_heure
......@@ -321,7 +321,7 @@ for indexd in ind_decoup:
#path windows
survey_name = 'ECHOSONDE_BIS'+'_'+str(premiere_date)+'_'+str(derniere_date)
survey_name = 'ECHOSONDE2021'+'_'+str(premiere_date)+'_'+str(derniere_date)
path_save = fpath + survey_name + "/"
......@@ -440,6 +440,21 @@ for indexd in ind_decoup:
save=save_data,
filename=filename_data_EK80,
)
with open(filename_data_EK80, "rb") as f:
(
freq_MFR,
freqs_moy_transducteur,
Sv,
Sv_moy_transducteur_x,
Sa,
Depth,
Lat,
Lon,
Time,
nb_transduc,
tableau_radiales,
nom_numero_radiales,
) = pickle.load(f)
del (
time_EK80_db,
......@@ -470,6 +485,17 @@ for indexd in ind_decoup:
premiere_date,
nameTransect,
)
file=path_netcdf+nameTransect[premiere_date:][0]+'_1.xsf.nc'
dataset=nc.Dataset(file)
Sv=dataset.groups['Sonar'].groups['Grid_group_1'].variables['integrated_backscatter'][:]
freq_MFR=dataset.groups['Sonar'].groups['Grid_group_1'].variables['frequency'][:]
freq_MFR=np.sort(freq_MFR)
Time=dataset.groups['Sonar'].groups['Grid_group_1'].variables['cell_ping_time'][:]
Depth=dataset.groups['Sonar'].groups['Grid_group_1'].variables['cell_depth'][:]
Depth=Depth.T
Lat=dataset.groups['Sonar'].groups['Grid_group_1'].variables['cell_latitude'][:]
Lon=dataset.groups['Sonar'].groups['Grid_group_1'].variables['cell_longitude'][:]
nom_numero_radiales=([nameTransect[premiere_date]],[premiere_date])
else:
if methode_EI=='pickle':
......@@ -503,30 +529,30 @@ for indexd in ind_decoup:
Lon=dataset.groups['Sonar'].groups['Grid_group_1'].variables['cell_longitude'][:]
nom_numero_radiales=([nameTransect[premiere_date]],[premiere_date])
if run_once<1:
tr=[]
for i in range(len(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:])):
if int(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:][i][9:]) in freq_custom:
tr.append(i)
else:
pass
if len(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:])!=nb_transduc:
run=input('Le nombre de transducteurs rentré par l\'utilisateur est différent de celui contenu dans les fichiers d\'EI, voulez vous poursuivre quand même ?')
if run=='n':
raise ValueError("Modifiez les paramètres suivants: nb_transduc, nbfreqs, indices_transduc_choisis et freq_custom")
else:
pass
elif len(tr)!=len(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:]):
run=input('Certaines fréquences rentrées par l\'utilisateur sont différentes de celles contenues dans les fichiers d\'EI ou il manque des fréquences, voulez vous poursuivre quand même ?')
if run=='n':
raise ValueError("Modifiez les paramètres suivants: nb_transduc, nbfreqs, indices_transduc_choisis et freq_custom")
else:
pass
if run_once<1:
tr=[]
for i in range(len(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:])):
if int(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:][i][9:]) in freq_custom:
tr.append(i)
else:
pass
run_once+=1
if len(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:])!=nb_transduc:
run=input('Le nombre de transducteurs rentré par l\'utilisateur est différent de celui contenu dans les fichiers d\'EI, voulez vous poursuivre quand même ?')
if run=='n':
raise ValueError("Modifiez les paramètres suivants: nb_transduc, nbfreqs, indices_transduc_choisis et freq_custom")
else:
pass
elif len(tr)!=len(dataset.groups['Sonar'].groups['Grid_group_1'].variables['beam'][:]):
run=input('Certaines fréquences rentrées par l\'utilisateur sont différentes de celles contenues dans les fichiers d\'EI ou il manque des fréquences, voulez vous poursuivre quand même ?')
if run=='n':
raise ValueError("Modifiez les paramètres suivants: nb_transduc, nbfreqs, indices_transduc_choisis et freq_custom")
else:
pass
else:
pass
run_once+=1
else:
pass
len_Sv1.append(len(Sv[1]))
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment