ratt-ru / shadeMS

Rapid Measurement Set plotting with dask-ms and datashader
20 stars 6 forks source link

shadems crashes when an antenna is fully flagged (marked faulty) and you are making one plot. #103

Open bennahugo opened 1 year ago

bennahugo commented 1 year ago

@SharmilaGoedhart reports that this only works when iterating.

Backlog to be posted

SharmilaGoedhart commented 1 year ago

(venv-shadems) sharmila@com06:/scratch2/sharmila/nb-artefacts$ cd sci-obs/ (venv-shadems) sharmila@com06:/scratch2/sharmila/nb-artefacts/sci-obs$ shadems --xaxis CHAN --yaxis DATA:amp bpcal.ms 2022-10-05 12:37:16 - shadems - INFO - using colourmap colorcet.bkr 2022-10-05 12:37:16 - shadems - INFO - using colourmap cmasher.pride 2022-10-05 12:37:16 - shadems - INFO - using colourmap colorcet.glasbey_dark 2022-10-05 12:37:16 - shadems - INFO - /scratch2/sharmila/venv-shadems/bin/shadems --xaxis CHAN --yaxis DATA:amp bpcal.ms 2022-10-05 12:37:16 - shadems - INFO - ------------------------------------------------------ 2022-10-05 12:37:16 - shadems - INFO - : MS bpcal.ms contains 135420 rows 2022-10-05 12:37:16 - shadems - INFO - : (1, 32768) spectral windows and channels 2022-10-05 12:37:16 - shadems - INFO - : 1 fields: J1939-6342 2022-10-05 12:37:16 - shadems - INFO - : 1 scans: 1 2022-10-05 12:37:16 - shadems - INFO - : 60/60 antennas: 0:m000 1:m001 2:m002 3:m003 4:m004 5:m005 6:m006 7:m007 8:m008 9:m009 10:m010 11:m011 12:m012 13:m013 14:m014 15:m015 16:m016 17:m017 18:m018 19:m019 20:m020 21:m021 22:m022 23:m023 24:m024 25:m025 26:m026 27:m030 28:m031 29:m032 30:m033 31:m034 32:m035 33:m036 34:m037 35:m038 36:m039 37:m040 38:m041 39:m042 40:m043 41:m044 42:m045 43:m046 44:m047 45:m048 46:m049 47:m050 48:m051 49:m052 50:m053 51:m054 52:m055 53:m057 54:m058 55:m059 56:m060 57:m061 58:m062 59:m063 2022-10-05 12:37:16 - shadems - INFO - : 1830/1830 baselines present 2022-10-05 12:37:16 - shadems - INFO - : corrs/Stokes XX YY I Q 2022-10-05 12:37:16 - shadems - INFO - ------------------------------------------------------ 2022-10-05 12:37:16 - shadems - INFO - : Data selected for plotting: 2022-10-05 12:37:16 - shadems - INFO - Antenna(s) : all 2022-10-05 12:37:16 - shadems - INFO - Baseline(s) : all except autocorrelations 2022-10-05 12:37:16 - shadems - INFO - Field(s) : all 2022-10-05 12:37:16 - shadems - INFO - SPW(s) : all 2022-10-05 12:37:16 - shadems - INFO - Scan(s) : all 2022-10-05 12:37:16 - shadems - INFO - Channels : all 2022-10-05 12:37:16 - shadems - INFO - Corr/Stokes : XX YY 2022-10-05 12:37:16 - shadems - INFO - ------------------------------------------------------ 2022-10-05 12:37:16 - shadems - INFO - loading minmax cache from bpcal-minmax-cache.json 2022-10-05 12:37:16 - shadems - INFO - axis: CHAN, range (None, None), discretization None 2022-10-05 12:37:16 - shadems - INFO - axis: amp(DATA), corr None, range (None, None), discretization None 2022-10-05 12:37:16 - shadems - INFO - : you have asked for 1 plots employing 2 unique datums 2022-10-05 12:37:16 - shadems - INFO - : Indexing MS and building dataframes (130980 rows, chunk size is 5000) 2022-10-05 12:37:16 - shadems - INFO - : complete 2022-10-05 12:37:16 - shadems - INFO - : rendering 1 dataframes with 8.58e+09 points into 1 plot types 2022-10-05 12:37:16 - shadems - INFO - : rendering plot-bpcal-DATA-XX-YY-amp-CHAN.png 2022-10-05 12:37:16 - shadems - INFO - : scanning axis min/max for DATA_amp_None /scratch2/sharmila/venv-shadems/lib/python3.8/site-packages/shade_ms/data_plots.py:254: RuntimeWarning: All-NaN slice encountered np.array([[(np.nanmin(df[axis].values).item() if bounds[axis][0] is None else bounds[axis][0]) for axis in unknowns]+ /scratch2/sharmila/venv-shadems/lib/python3.8/site-packages/shade_ms/data_plots.py:255: RuntimeWarning: All-NaN slice encountered [(np.nanmax(df[axis].values).item() if bounds[axis][1] is None else bounds[axis][1]) for axis in unknowns]]), /scratch2/sharmila/venv-shadems/lib/python3.8/site-packages/shade_ms/data_plots.py:254: RuntimeWarning: All-NaN slice encountered np.array([[(np.nanmin(df[axis].values).item() if bounds[axis][0] is None else bounds[axis][0]) for axis in unknowns]+ /scratch2/sharmila/venv-shadems/lib/python3.8/site-packages/shade_ms/data_plots.py:255: RuntimeWarning: All-NaN slice encountered [(np.nanmax(df[axis].values).item() if bounds[axis][1] is None else bounds[axis][1]) for axis in unknowns]]), 2022-10-05 12:41:51 - shadems - INFO - : wrote plot-bpcal-DATA-XX-YY-amp-CHAN.png 2022-10-05 12:41:51 - shadems - INFO - Total time : 275.65 seconds 2022-10-05 12:41:51 - shadems - INFO - Saved minmax cache to bpcal-minmax-cache.json (disable with --no-lim-save) 2022-10-05 12:41:51 - shadems - INFO - Finished 2022-10-05 12:41:51 - shadems - INFO - ------------------------------------------------------

SharmilaGoedhart commented 1 year ago

Apologies - it seems like it did work.. I missed seeing the output file last night. Shouldn't be working so late :-/