args.bbw = {
'file_name': r'.\00_Dados\Dados_Victor\bbw_zhang_s2a.xlsx',
'sheet_name': 'Sheet1'
}
args.Kd_ZEU = {
'file_name': r'.\00_Dados\Dados_Victor\Kd_ZEU\Kd_ZEU_TRM.xlsx',
'sheet_name': 'Sheet1'
}
args.Kd_3m = {
'file_name': r'.\00_Dados\Dados_Victor\Kd_3m\Kd_3m_TRM.xlsx',
'sheet_name': 'Sheet1'
}
args.Kd_6m = {
'file_name': r'.\00_Dados\Dados_Victor\Kd_6m\Kd_6m_TRM.xlsx',
'sheet_name': 'Sheet1'
}
data_obj = read_files(args)
lb0 = 560
Opt = '_OptY'
qaa_v5_s2a_560 = Qaa_v5_model_RUN(data_obj, 'TRM_2019', lb0)
qaa_v5_s2a_560 = load_data('.\Results_S2A\QAAv5_TRM_2019_S2a_' + str(lb0) +
Opt + '.pickle')
args19.kd_6m = {
'file_name': r'.\00_Dados\Dados_Victor\Kd_6m\Kd_6m_TRM.xlsx',
'sheet_name': 'Sheet1'
}
data_obj19 = read_files(args19)
qaa_v6_fit19_665 = Qaa_v6_model_RUN_original(data_obj19, 'original_TRM_2019',
665)
###############################################################################
######################### QAA V6 - Campanhas 2013/19 ########################
###############################################################################
## Abrindo arquivos contendo dicionários de dados resultados da aplicação do QAA/Kd;
qaa_v6_fit13_665 = load_data(
r'.\Results_Espectral\QAAv6_original_TRM_2013_665_Ajustado.pickle')
qaa_v6_fit19_665 = load_data(
r'.\Results_Espectral\QAAv6_original_TRM_2019_665_Ajustado.pickle')
kd_ref = pd.concat(
[qaa_v6_fit13_665['Kd_ZEU'].copy(), qaa_v6_fit19_665['Kd_ZEU'].copy()],
axis=1)
kd_pred = pd.concat([
qaa_v6_fit13_665['kd_qaa_665'].copy(),
qaa_v6_fit19_665['kd_qaa_665'].copy()
],
axis=1)
kd_ref13 = qaa_v6_fit13_665['Kd_ZEU'].copy().drop(['P25', 'P26'], axis=1)
kd_pred13 = qaa_v6_fit13_665['kd_qaa_665'].copy().drop(['P25', 'P26'], axis=1)
args.coleta = {
'file_name': r'.\00_Dados\Dados_Renato\Coleta\TRM_2013_Datas-2.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_zeu = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_ZEU\Kd_TRM_2013_ZEU.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_3m = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_3m\Kd_TRM_2013_3m.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_6m = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_6m\Kd_TRM_2013_6m.xlsx',
'sheet_name': 'Sheet1'
}
data_obj = read_files(args)
lb0 = 555
Opt = '_OptY'
qaa_v5_555_2 = Qaa_v5_model_RUN(data_obj, 'TRM_2013', 555)
qaa_v5_555_2 = load_data('.\Results_Espectral\QAAv5_TRM_2013_' + str(lb0) +
Opt + '.pickle')
'file_name': r'.\00_Dados\Dados_Victor\Kd_3m\Kd_3m_TRM.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_6m = {
'file_name': r'.\00_Dados\Dados_Victor\Kd_6m\Kd_6m_TRM.xlsx',
'sheet_name': 'Sheet1'
}
data_obj = read_files(args)
## Rodando o modelo QAA/Kd com base nos arquivos chamados pelos argumentos.
qaa_v6_fit19_560 = Qaa_v6_model_RUN(data_obj, 'TRM_2019', 560)
## Abrindo umarquivo contendo dicionário de dados resultados da aplicação do QAA/Kd.
qaa_v6_fit19_560 = load_data(r'.\Results_Espectral\QAAv6_TRM_2019_560_Ajustado.pickle')
plt.figure(figsize = (10, 5))
plt.subplot(121)
plt.plot(qaa_v6_fit13_560['Kd_ZEU']['P15'], 'r-')
plt.plot(qaa_v6_fit13_560['kd_qaa_560_acs_492']['P15'], 'b-')
plt.xlabel('Comprimento de Onda [nm]', fontsize = 14)
plt.ylabel('Kd [m ^ -1]', fontsize = 14)
plt.xlim([400, 700])
plt.ylim([0, 2])
plt.rc('xtick', labelsize = 14)
plt.rc('ytick', labelsize = 14)
plt.title('2013 - P15', fontsize = 14)
plt.show()
'sheet_name': 'Sheet1'
}
args.kd_zeu = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_ZEU\Kd_TRM_2013_ZEU.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_3m = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_3m\Kd_TRM_2013_3m.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_6m = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_6m\Kd_TRM_2013_6m.xlsx',
'sheet_name': 'Sheet1'
}
data_obj = read_files(args)
lb0 = 704
Opt = '_OptY'
OptAw = '_OptAw'
qaa_v5_704_2 = Qaa_v5_model_RUN(data_obj, 'TRM_2013', 704)
qaa_v5_704_2 = load_data('.\Results_Espectral\QAAv5_TRM_2013_' + str(lb0) +
Opt + '.pickle')
qaa_v5_704_2aw = load_data('.\Results_Espectral\QAAv5_TRM_2013_' + str(lb0) +
OptAw + '.pickle')
'file_name': r'.\00_Dados\Dados_Renato\bbw_zhang_PAR.xlsx',
'sheet_name': 'Sheet1'
}
args.coleta = {
'file_name': r'.\00_Dados\Dados_Renato\Coleta\TRM_2013_Datas-2.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_zeu = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_ZEU\Kd_TRM_2013_ZEU.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_3m = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_3m\Kd_TRM_2013_3m.xlsx',
'sheet_name': 'Sheet1'
}
args.kd_6m = {
'file_name': r'.\00_Dados\Dados_Renato\Kd_6m\Kd_TRM_2013_6m.xlsx',
'sheet_name': 'Sheet1'
}
data_obj = read_files(args)
qaa_v6_670_2 = Qaa_v6_model_RUN(data_obj, 'TRM_2013', 670)
qaa_v6_670_2 = load_data('.\Results_Espectral\QAAv6_TRM_2013_670.pickle')