use hms/dms format for plot axes

Author: PaulHancock

robbie_viewer_server/main.py

@@ -3,6 +3,8 @@
 from astropy.io import fits
 from astropy.io.votable import parse
 from astropy.wcs import WCS
+from astropy.coordinates import Angle
+import astropy.units as u
 import numpy as np
 import glob
 import pandas as pd
@@ -15,6 +17,7 @@
 from bokeh.layouts import layout
 from bokeh.io import curdoc
 from bokeh.models import ColumnDataSource, Circle, CustomJS
+from bokeh.models.formatters import FuncTickFormatter
 from utils import get_joined_table_source, get_scatter_plots, get_mean_image_plot,\
     get_light_curve_plot, get_light_curve_plot, get_epoch_image_plots, load_mean_image, mean_image_data, get_imdata
 
@@ -69,6 +72,42 @@
 sky.x_range = mean_image.x_range = epochs.x_range
 sky.y_range = mean_image.y_range = epochs.y_range
 
+# make the plot coords be in hms/dms format
+dms = """
+var s = Math.sign(tick);
+var t = Math.abs(tick);
+var deg = Math.floor(t);
+t = (t -deg)*60;
+var min = Math.floor(t);
+t = (t-min)*60;
+var sec = Math.floor(t);
+var tk = ""
+if (s < 0) {
+  tk = String(s)[0]
+}
+tk = tk + String(deg).padStart(2,'0') + ":"+String(min).padStart(2,'0')+":"+String(sec).padStart(2,'0');
+return tk;
+"""
+hms = """
+var s = Math.sign(tick);
+var t = Math.abs(tick/15);
+var deg = Math.floor(t);
+t = (t -deg)*60;
+var min = Math.floor(t);
+t = (t-min)*60;
+var sec = Math.floor(t);
+var tk = ""
+if (s < 0) {
+  tk = String(s)[0]
+}
+tk = tk + String(deg).padStart(2,'0') + ":"+String(min).padStart(2,'0')+":"+String(sec).padStart(2,'0');
+return tk;
+"""
+
+sky.yaxis.formatter = mean_image.yaxis.formatter = epochs.yaxis.formatter = FuncTickFormatter(code=dms)
+sky.xaxis.formatter = mean_image.xaxis.formatter = epochs.xaxis.formatter = FuncTickFormatter(code=hms)
+
+
 # # Callback to update epochs
 # def update_epochs():
 

robbie_viewer_server/utils.py

@@ -139,7 +139,7 @@ def get_scatter_plots(source):
     left = figure(
         tools=TOOLS,
         title='Sky Plot',
-        x_axis_label='RA (deg)', y_axis_label='DEC (deg)',
+        x_axis_label='RA (hour)', y_axis_label='DEC (deg)',
         x_range=(0,1),
         y_range=(0,1),
     )
@@ -186,8 +186,7 @@ def get_mean_image_plot(source, result_dir, ra_ref=None, dec_ref=None, degrees_a
             ("RA", "$x{0.00}°"),
             ("DEC", "$y{0.00}°")
         ],
-        x_axis_label='RA',
-        y_axis_label='DEC',
+        x_axis_label='RA (hour)', y_axis_label='DEC (deg)',
         x_range=(min(source.data['ref_ra']),  max(source.data['ref_ra'])),
         y_range=(min(source.data['ref_dec']), max(source.data['ref_dec'])),
     )
@@ -218,8 +217,7 @@ def get_epoch_image_plots(epoch_files, mean_source, ra_ref=None, dec_ref=None, d
             ("RA", "$x{0.00}°"),
             ("DEC", "$y{0.00}°")
         ],
-        x_axis_label='RA',
-        y_axis_label='DEC',
+        x_axis_label='RA (hour)', y_axis_label='DEC (deg)',
         x_range=(min(mean_source.data['ref_ra']),  max(mean_source.data['ref_ra'])),
         y_range=(min(mean_source.data['ref_dec']), max(mean_source.data['ref_dec'])),
     )