Ejemplo n.º 1
0
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)
Ejemplo n.º 3
0
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()
Ejemplo n.º 5
0
    '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')
Ejemplo n.º 6
0
                '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')