diff --git a/bstplot.py b/bstplot.py index 0a14772b42e123e62e75687a8d4792e2397e5fb8..b37eb7613e97b2e39fc785633acf043942f4e12a 100644 --- a/bstplot.py +++ b/bstplot.py @@ -33,8 +33,6 @@ if __name__ == "__main__": print(data.shape) print(metadata) - t = Time("2024-10-10T09:00:00", format="isot", scale="utc") + np.arange(86305) * u.s - metadata["timestamps"] = t.isot plot_bst_dynamic_spectrum(metadata, data, args.pol) diff --git a/lofty/io.py b/lofty/io.py index e65139a439c2bf781980f37915e63c1fde5f33a6..de81411e9b527a7054c8b6615b10fd6b84a94857 100644 --- a/lofty/io.py +++ b/lofty/io.py @@ -481,12 +481,15 @@ def read_single_bst(bst_filename): metadata = {} # for key in ["antenna_names", "antenna_quality", "antenna_type", "antenna_usage_mask", "station_name"]: # metadata[key] = h5["/"].attrs[key] - for key in ["antenna_names", "antenna_type", "station_name"]: + for key in ["antenna_names", "antenna_type", "station_name", "antenna_usage_mask"]: metadata[key] = h5["/"].attrs[key] # Antennafield - metadata["antennafield"] = h5["/"].attrs["antennafield_device"].replace("STAT/AFH/", "").replace("STAT/AFL/", "") - + if "antennafield_device" in h5["/"].attrs: + metadata["antennafield"] = h5["/"].attrs["antennafield_device"].replace("STAT/AFH/", "").replace("STAT/AFL/", "") + elif "antenna_field" in h5["/"].attrs: + metadata["antennafield"] = h5["/"].attrs["antenna_field"] + # Number of subintegrations nsub = len(h5.keys()) metadata["nsub"] = nsub diff --git a/lofty/plotting.py b/lofty/plotting.py index a2152be3a99ed5d5f396ff1b0b5d2ede73fe7123..e6652e8f2722368c89fcc61f4c8b9ff9bd48b50b 100644 --- a/lofty/plotting.py +++ b/lofty/plotting.py @@ -354,8 +354,8 @@ def plot_bst_dynamic_spectrum(metadata, data, pol): data_mean = np.mean(data, axis=0) - #data /= data_mean[np.newaxis, :, :] - data = decibel(data) + data /= data_mean[np.newaxis, :, :] + #data = decibel(data) if pol.lower() == "x": p = data[:, 0, :].squeeze() diff --git a/xstimager.py b/xstimager.py index 61e403ab8535c9ae606d07ab76f44fb41970596f..6372a685a95d93715ddeeb276c407129f25a1f48 100644 --- a/xstimager.py +++ b/xstimager.py @@ -28,6 +28,8 @@ if __name__ == "__main__": else: metadata, data = read_xst(sorted(args.filenames), fill_empty=True) + metadata["timestamps"] = [timestamp.replace("+0000", "") for timestamp in metadata["timestamps"]] + print(metadata["timestamps"]) # Get information station_type = get_station_type(metadata["station_name"].upper()) antenna_set = metadata["antennafield"]