files = list(folder.glob("Tara_HyperPro*.txt"))
data = table.vstack([read(file, data_start=35) for file in files])
header = read(files[0], data_start=32, data_end=33)
header["col1"][0] = "year"
header = header[0].as_void()
for key, new_key in zip(data.keys(), header):
data.rename_column(key, new_key)
data.remove_columns(get_keys_with_label(data, "LU"))
data.rename_column("lat", "Latitude")
data.rename_column("lon", "Longitude")
rename_columns(data, "ES", "Ed_", exclude="None")
rename_columns(data, "Rrs", "R_rs_", exclude="None")
convert_to_unit(data, "Ed", u.microwatt / (u.centimeter**2 * u.nanometer),
u.watt / (u.meter**2 * u.nanometer))
convert_to_unit(data, "R_rs", 1 / u.steradian)
data = add_Lw_from_Ed_Rrs(data)
map_data(data, data_label="TaraO", lon_0=0, resolution="i")
plot_spectra(data, data_label="TaraO", alpha=0.1)
write_data(data, label="TaraO")
from sba.data_processing import split_spectrum, get_keys_with_label, remove_negative_R_rs, convert_to_unit, rename_columns
Ed = read("data/SOP4/SO-P4_irrad.tab", data_start=142, header_start=141)
Lu = read("data/SOP4/SO-P4_rad_up_40deg.tab", data_start=142, header_start=141)
Ls = read("data/SOP4/SO-P4_sky_rad_40deg.tab",
data_start=142,
header_start=141)
data = table.join(Ed, Lu, keys=["Date/Time"])
data = table.join(data, Ls, keys=["Date/Time"])
rename_columns(data, "Ed", "Ed", strip=True)
rename_columns(data, "Lu", "Lu", strip=True)
rename_columns(data, "Ls", "Ls", strip=True)
convert_to_unit(data, "Ed", u.watt / (u.meter**2 * u.nanometer))
convert_to_unit(data, "Lu",
u.microwatt / (u.centimeter**2 * u.nanometer * u.steradian),
u.watt / (u.meter**2 * u.nanometer * u.steradian))
convert_to_unit(data, "Ls", u.watt / (u.meter**2 * u.nanometer * u.steradian))
Ed_keys, Lu_keys, Ls_keys = get_keys_with_label(data, "Ed", "Lu", "Ls")
for Ed_k, Lu_k, Ls_k in zip(Ed_keys, Lu_keys, Ls_keys):
Lw_k = Lu_k.replace("Lu", "Lw")
Lw = data[Lu_k] - 0.028 * data[Ls_k]
Lw.name = Lw_k
data.add_column(Lw)
R_rs = data[Lw_k] / data[Ed_k]
R_rs.name = Lw_k.replace("Lw", "R_rs")
R_rs.unit = 1 / u.steradian