datetick 0.9.0__py3-none-any.whl

datetick/__init__.py

@@ -0,0 +1,7 @@
+from datetick.datetick import datetick
+from importlib.metadata import version, PackageNotFoundError
+
+try:
+    __version__ = version('datetick')
+except PackageNotFoundError:
+    __version__ = 'unknown'

datetick/datetick.py

@@ -0,0 +1,694 @@
+def datetick(*args,
+             axes=None,
+             set_cb=True,
+             adjust_last_xlabel=False,
+             adjust_first_xlabel=False,
+             adjust_xrange=False,
+             adjust_yrange=False,
+             debug=False):
+  """
+  datetick('x') or datetick('y') formats the major and minor tick labels
+  of the current figure.
+
+  datetick('x', axes=ax) or datetick('y', axes=ax) formats the given axes `ax`.
+
+  Example:
+  --------
+    import datetime as dt
+    import matplotlib.pyplot as plt
+    from datetick import datetick
+    d1 = dt.datetime(1900, 1, 2)
+    d2 = dt.datetime.fromordinal(10 + dt.datetime.toordinal(d1))
+    x = [d1, d2]
+    y = [0.0,1.0]
+    plt.clf()
+    plt.plot(x, y)
+    datetick('x')
+  """
+
+  # Based on spacepy/plot/utils.py on 07/10/2017, but many additions.
+  # See also https://github.com/JouleCai/geospacelab/blob/master/geospacelab/visualization/mpl/axis_ticks.py
+
+  # TODO: Use _numsize() to determine if figure width and height
+  #       will cause overlap when default number major tick labels is used.
+  # TODO: If time[0].day > 28, need to make first tick at time[0].day = 28
+  #       as needed.
+  # TODO: If first data point has fractional seconds, the plot won't have
+  #       a major x-label right below it. This is due to the fact that
+  #       MicrosecondLocator() does not take a keyword argument of
+  #       "bymicroseconds".
+  # TODO: Adjust lower and upper limits as in 366*8 span
+
+  # Get all kwargs passed using locals
+  kwargs = {k: v for k, v in locals().items() if k != 'args'}
+
+  import warnings
+  from datetime import datetime
+
+  import matplotlib.dates as mpld
+
+  if len(args) == 0:
+    dir = 'x'
+  else:
+    dir = args[0]
+
+  plt = _get_plt(debug=debug)
+
+  if 'axes' in kwargs:
+    axes = kwargs['axes']
+    fig = axes.figure
+  else:
+    axes = plt.gca()
+    fig = plt.gcf()
+
+  fig.canvas.draw()
+  bbox = axes.dataLim
+
+  if dir == 'x':
+    datamin = bbox.x0
+    datamax = bbox.x1
+    lim = axes.get_xlim()
+    ticks = axes.get_xticks()
+  else:
+    datamin = bbox.y0
+    datamax = bbox.y1
+    lim = axes.get_ylim()
+    ticks = axes.get_yticks()
+
+  try:
+    mpld.num2date(lim[0])
+  except:
+    msg = 'Lower axis limit of %f is not a valid Matplotlib datenum' % lim[0]
+    raise ValueError(msg)
+  try:
+    mpld.num2date(lim[1])
+  except:
+    msg = 'Upper axis limit of %f is not a valid Matplotlib datenum' % lim[1]
+    raise ValueError(msg)
+
+  # If all values are NaN, datamin = np.inf and datamax = -np.inf.
+  msg = "is not a valid Matplotlib datenum. Are all data values NaN? Cannot use datetick()"
+  try:
+    mpld.num2date(datamin)
+  except:
+    warnings.warn(f'Minimum value of %f {msg}' % datamin)
+    return
+  try:
+    mpld.num2date(datamax)
+  except:
+    warnings.warn(f'Maximum value of %f {msg}' % datamax)
+    return
+
+  if datamin == datamax:
+    axes.set_xticks([mpld.num2date(datamin)])
+    xticklabel = datetime.strftime(mpld.num2date(datamin),'%Y-%m-%dT%H:%M:%S')
+    axes.set_xticklabels([xticklabel])
+    return
+
+  # Need to document why this was used. It creates
+  # problems if labels extend beyond the axis limits.
+  #tmin = np.min((lim[0], datamin))
+  #tmax = np.max((lim[1], datamax))
+  tmin = lim[0]
+  tmax = lim[1]
+
+  tspan = mpld.num2date((tmin, tmax))
+
+  deltaT = tspan[-1] - tspan[0]
+  if debug:
+    print("Total seconds: %s" % deltaT.total_seconds())
+
+  """
+  fmt1 is format of the tick labels
+
+  fmt2 contains additional information that is used for the first tick label
+  or when there is a major change. For example, if
+    fmt1 = %M:%S and fmt2 = %H,
+  the labels will have only minute and hour and the first tick will have a
+  label of %M:%S\n%H. If there is a change in hour somewhere on the axis,
+  that label will include the new hour.
+
+  Note that interval=... is specified even when it would seem to be redundant.
+  It is needed to workaround the bug discussed at stackoverflow.com/q/31072589
+  """
+
+  Mtick, mtick, fmt1, fmt2 = _locator(axes, deltaT, debug=debug)
+
+  if debug:
+    print(f'{dir} data min:         {mpld.num2date(datamin)}')
+    print(f'Default {dir}lim[0]:    {mpld.num2date(lim[0])}')
+    print(f'Default {dir}ticks[0]:  {mpld.num2date(ticks[0])}')
+    print(f'{dir} data max:         {mpld.num2date(datamax)}')
+    print(f'Default {dir}lim[-1]:   {mpld.num2date(lim[-1])}')
+    print(f'Default {dir}ticks[-1]: {mpld.num2date(ticks[-1])}')
+    print(f'Default {dir}ticks:')
+    for i in range(0,len(ticks)):
+      print(f' {mpld.num2date(ticks[i])}')
+
+  if dir == 'x':
+    axes.xaxis.set_major_locator(Mtick)
+    axes.xaxis.set_minor_locator(mtick)
+    axes.xaxis.set_major_formatter(fmt1)
+    if adjust_xrange:
+      _adjust_range(dir, fig, axes, datamin, datamax, debug=debug)
+    fig.canvas.draw() # Render new labels so updated for next line
+    labels = [item.get_text() for item in axes.get_xticklabels()]
+    ticks = axes.get_xticks()
+  else:
+    axes.yaxis.set_major_locator(Mtick)
+    axes.yaxis.set_minor_locator(mtick)
+    axes.yaxis.set_major_formatter(fmt1)
+    if adjust_xrange:
+      _adjust_range(dir, fig, axes, datamin, datamax, debug=debug)
+    fig.canvas.draw() # Render new labels so updated for next line
+    labels = [item.get_text() for item in axes.get_yticklabels()]
+    ticks = axes.get_yticks()
+
+  if debug:
+    xl = axes.get_xlim()
+    print(f'New {dir}ticks:')
+    for i in range(0,len(ticks)):
+      note = ''
+      if ticks[i] < xl[0] or ticks[i] > xl[1]:
+        note = ' (will be clipped by mpl b/c outside of axis limits)'
+      print(f' {mpld.num2date(ticks[i])} {note}')
+
+  if len(labels) == 0:
+    if debug:
+      print('No labels.')
+    return
+
+  if debug:
+    print(f'fmt1 {dir}labels:')
+    for i in range(0,len(labels)):
+      print(f' {labels[i]}')
+
+  if fmt2 != '':
+    labels = _add_fmt2(fmt2, lim, ticks, deltaT, labels, dir, debug=debug)
+
+    if debug:
+      print(f'fmt1 + fmt2 {dir}labels:')
+      for i in range(0,len(labels)):
+        print(f' {labels[i].replace("\n", "\\n")}')
+
+    if dir == 'x':
+      # Without the set_xticks(), warning is generated:
+      #   UserWarning: set_ticklabels() should only be used.
+      # Additional discussion: https://github.com/matplotlib/matplotlib/issues/18848
+      # The correct way to avoid the warning: https://stackoverflow.com/a/69126185
+      # Better: Create custom class:
+      #   https://matplotlib.org/stable/gallery/ticks/date_index_formatter.html
+      axes.set_xticks(axes.get_xticks())
+      axes.set_xticklabels(labels)
+      _adjust_xlabels(axes,
+                      adjust_first_xlabel=adjust_first_xlabel,
+                      adjust_last_xlabel=adjust_last_xlabel,
+                      debug=debug)
+
+    if dir == 'y':
+      axes.set_yticks(axes.get_yticks())
+      axes.set_yticks(axes.get_yticks())
+      axes.set_yticklabels(labels)
+
+  # Trigger update of ticks when limits change due to user interaction.
+  if 'set_cb':
+    def on_xlims_change(ax):
+      datetick('x', **{**kwargs, 'set_cb': False})
+
+    def on_ylims_change(ax):
+      datetick('y', **{**kwargs, 'set_cb': False})
+
+    if dir == 'x':
+      axes.callbacks.connect('xlim_changed', on_xlims_change)
+    else:
+      axes.callbacks.connect('ylim_changed', on_ylims_change)
+
+
+def _adjust_range(dir, fig, axes, datamin, datamax, debug=False):
+  if dir == 'x':
+    ticks = axes.get_xticks()
+  else:
+    ticks = axes.get_yticks()
+  if len(ticks) >= 2:
+    fig.canvas.draw()
+    dt = ticks[1] - ticks[0]
+    pad = 0.05 * dt
+    first_candidates  = ticks[ticks <= datamin]
+    last_candidates = ticks[ticks >= datamax]
+    first  = first_candidates[-1]  if len(first_candidates)  > 0 else ticks[0] - dt
+    last = last_candidates[0]  if len(last_candidates) > 0 else ticks[-1] + dt
+    if dir == 'x':
+      axes.set_xlim(first - pad, last + pad)
+    else:
+      axes.set_ylim(first - pad, last + pad)
+
+
+def _add_fmt2(fmt2, lim, ticks, deltaT, labels, dir, debug=False):
+  from datetime import datetime
+  import matplotlib.dates as mpld
+
+  time = mpld.num2date(ticks)
+
+  nDays  = deltaT.days
+  nSecs  = deltaT.total_seconds()
+
+  first = 0
+  if ticks[0] < lim[0]:
+    if debug:
+      msg = 'First tick is less than lower axis limit. Applying fmt2 to second tick label.'
+      print(msg)
+    # Work-around for bug in Matplotlib where left-most tick is less than
+    # lower x-limit. Could more than one tick be less than lower x-limit?
+    first = 1
+
+  # Always apply fmt2 to first tick label
+  if debug:
+    print(f'Applying fmt2 to first tick label at {mpld.num2date(ticks[first])}.')
+  labels[first] = '%s\n%s' % (labels[first], datetime.strftime(time[first], fmt2))
+
+  for i in range(first+1, len(time)):
+    # First label will always have fmt1 applied.
+    # Modify labels after first under certain conditions.
+    modify = False
+
+    if time[i].year > time[i-1].year:
+      modify = True
+    if nDays < 60 and time[i].month > time[i-1].month:
+      modify = True
+    if nDays < 4 and time[i].day > time[i-1].day:
+      modify = True
+    if nSecs < 60*30 and time[i].hour > time[i-1].hour:
+      modify = True
+    if nSecs < 1 and time[i].minute > time[i-1].minute:
+      modify = True
+    if nSecs < 1 and time[i].second > time[i-1].second:
+      modify = True
+
+    if not modify:
+      continue
+
+    if i == first + 1 and dir == 'x':
+      # If first two major tick labels have fmt2 applied, the will
+      # likely run together. This keeps fmt2 label for second major
+      # tick.
+      if debug:
+        print(f'Removing fmt2 to first tick label at {mpld.num2date(ticks[first])} to avoid overlap with second label.')
+        print(f'Applying fmt2 to tick label at           {mpld.num2date(ticks[i])}.')
+      labels[first] = labels[first].split('\n')[0]
+      labels[i] = '%s\n%s' % (labels[i], datetime.strftime(mpld.num2date(ticks[i]), fmt2))
+    else:
+      if debug:
+        print(f'Applying fmt2 to tick label at       {mpld.num2date(ticks[i])}.')
+      labels[i] = '%s\n%s' % (labels[i], datetime.strftime(mpld.num2date(ticks[i]), fmt2))
+
+  return labels
+
+
+def _get_plt(debug=False):
+  import warnings
+  import matplotlib
+  import matplotlib.rcsetup
+
+  # TODO: Determine earliest version of Matplotlib where the problems
+  #       addressed here are fixed.
+  try:    # Available in Matplotlib >= 3.9
+    if matplotlib.__version__.split(".") >= (3, 9):
+      import matplotlib.backends.backend_registry
+      backends = matplotlib.backends.backend_registry.list_builtin()
+    else:
+      backends = matplotlib.rcsetup.all_backends
+  except Exception:
+    # In earlier versions of Matplotlib (at least 3.3), the list of backends
+    # was not available and get_backend() did not return lower-case names.
+    backends = ['Qt5Agg', 'QT4Agg', 'GTKAgg', 'TKAgg', 'WXAgg']
+
+  if debug:
+    print(f'Matplotlib version: {matplotlib.__version__}')
+    print(f'Matplotlib backend: {matplotlib.get_backend()}')
+    print(f'Available backends: {backends}')
+
+  if matplotlib.get_backend() == 'MacOSX':
+    """
+    With MacOSX backend, draw() does not update the ticks. See warning at
+    https://matplotlib.org/3.3.0/tutorials/advanced/blitting.html
+    """
+    import sys
+    if debug:
+      print('Matplotlib backend is MacOSX. Switching backend globally to workaround draw() not updating ticks.')
+    if sys.version_info[0:2] < (3, 6):
+      # warnings.filterwarnings("ignore", '.*backend.*', category=UserWarning)
+      # the above should work and is better because more specific.
+      warnings.simplefilter("ignore", category=UserWarning)
+    for backend in backends:
+      cmd = f"matplotlib.use('{backend}', force=True)"
+      try:
+        if debug:
+          print(f"Trying {cmd}")
+        matplotlib.use(backend, force=True)
+        import matplotlib.pyplot as plt
+        if debug:
+          print("  Success.")
+        break
+      except:
+        if debug:
+          print(" Failure.")
+        continue
+  else:
+    try:
+      import matplotlib.pyplot as plt
+    except:
+      print('Failed: "import matplotlib.pyplot as plt". Switching backend globally.')
+      cmd = f"matplotlib.use('{backends}', force=True)"
+      for backend in backends:
+        try:
+          if debug:
+            print(f"Trying {cmd}")
+          matplotlib.use(backend, force=True)
+          import matplotlib.pyplot as plt
+          if debug:
+            print("  Success.")
+          break
+        except:
+          if debug:
+            print(" Failure.")
+          continue
+
+  return plt
+
+
+def _locator(axes, deltaT, debug=False):
+  import matplotlib
+  import matplotlib.dates as mpld
+
+  def _millis(x, pos):
+    x = matplotlib.dates.num2date(x)
+    label = x.strftime('.%f')
+    label = label[0:3]
+    #label = label.rstrip(".")
+    return label
+
+  nHours = deltaT.days * 24.0 + deltaT.seconds/3600.0
+  if deltaT.total_seconds() < 0.1:
+    # < 0.1 second
+    Mtick = mpld.MicrosecondLocator(interval=10000)
+    mtick = mpld.MicrosecondLocator(interval=2000)
+    fmt1 = matplotlib.ticker.FuncFormatter(_millis)
+    fmt2  = '%H:%M:%S\n%Y-%m-%d'
+  if deltaT.total_seconds() < 0.5:
+    # < 0.5 seconds
+    # Locators don't locate at this resolution.
+    # Need to do this manually. See comment above.
+    Mtick = mpld.MicrosecondLocator(interval=50000)
+    mtick = mpld.MicrosecondLocator(interval=10000)
+    fmt1 = matplotlib.ticker.FuncFormatter(_millis)
+    fmt2  = '%H:%M:%S\n%Y-%m-%d'
+  if deltaT.total_seconds() < 1:
+    # < 1 second
+    # https://matplotlib.org/api/dates_api.html#matplotlib.dates.MicrosecondLocator
+    # MircosecondLocator() does not have a "bymicrosecond" option. If
+    # first point is not at zero microseconds, it won't be labeled.
+    Mtick = mpld.MicrosecondLocator(interval=100000)
+    mtick = mpld.MicrosecondLocator(interval=20000)
+    fmt1 = matplotlib.ticker.FuncFormatter(_millis)
+    #fmt1  = mpld.DateFormatter('%M:%S.%f')
+    fmt2  = '%H:%M:%S\n%Y-%m-%d'
+  elif deltaT.total_seconds() < 5:
+    # < 5 seconds
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 1)) )
+    mtick = mpld.MicrosecondLocator(interval=200000)
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 10:
+    # < 10 seconds
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 1)) )
+    mtick = mpld.MicrosecondLocator(interval=500000)
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 20:
+    # < 20 seconds
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 2)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 1)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 30:
+    # < 30 seconds
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 5)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 1)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60:
+    # < 1 minute
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 10)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 2)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60*2:
+    # < 2 minutes
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 20)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 5)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60*3:
+    # < 3 minutes
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 20)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 5)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60*5:
+    # < 5 minutes
+    Mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 30)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 10)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60*10:
+    # < 10 minutes
+    Mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 1)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 15)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60*20:
+    # < 20 minutes
+    Mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 2)) )
+    mtick = mpld.SecondLocator(bysecond=list(range(0, 60, 30)) )
+    fmt1  = mpld.DateFormatter('%M:%S')
+    fmt2  = '%Y-%m-%dT%H'
+  elif deltaT.total_seconds() < 60*30:
+    # < 30 minutes
+    Mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 5)) )
+    mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 1)) )
+    fmt1  = mpld.DateFormatter('%H:%M')
+    fmt2  = '%Y-%m-%d'
+  elif deltaT.total_seconds() < 60*60:
+    # < 60 minutes
+    Mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 10)) )
+    mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 2)) )
+    fmt1  = mpld.DateFormatter('%H:%M')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 2:
+    Mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 15)) )
+    mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 5)) )
+    fmt1  = mpld.DateFormatter('%H:%M')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 4:
+    Mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 20)) )
+    mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 5)) )
+    fmt1  = mpld.DateFormatter('%H:%M')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 6:
+    Mtick = mpld.HourLocator(byhour=list(range(0,24,1)) )
+    mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 10)) )
+    fmt1  = mpld.DateFormatter('%H:%M')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 12:
+    Mtick = mpld.HourLocator(byhour=list(range(0,24,2)) )
+    mtick = mpld.MinuteLocator(byminute=list(range(0, 60, 30)) )
+    fmt1  = mpld.DateFormatter('%H:%M')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 24:
+    # < 1 day
+    Mtick = mpld.HourLocator(byhour=list(range(0, 24, 3)) )
+    mtick = mpld.HourLocator(byhour=list(range(0, 24, 1)) )
+    fmt1  = mpld.DateFormatter('%H')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 48:
+    # < 2 days
+    Mtick = mpld.HourLocator(byhour=list(range(0, 24, 4)) )
+    mtick = mpld.HourLocator(byhour=list(range(0, 24, 2)) )
+    fmt1  = mpld.DateFormatter('%H')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 72:
+    # < 3 days
+    Mtick = mpld.HourLocator(byhour = list(range(0, 24, 6)))
+    mtick = mpld.HourLocator(byhour = list(range(0, 24, 3)))
+    fmt1  = mpld.DateFormatter('%H')
+    fmt2  = '%Y-%m-%d'
+  elif nHours < 96:
+    # < 4 days
+    Mtick = mpld.HourLocator(byhour = list(range(0, 24, 12)))
+    mtick = mpld.HourLocator(byhour = list(range(0, 24, 3)))
+    fmt1  = mpld.DateFormatter('%H')
+    fmt2  = '%Y-%m-%d'
+  elif deltaT.days < 8:
+    Mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 1)))
+    mtick = mpld.HourLocator(byhour=list(range(0, 24, 4)))
+    fmt1  = mpld.DateFormatter('%d')
+    fmt2  = '%Y-%m'
+  elif deltaT.days < 16:
+    Mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 1)))
+    mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 1)))
+    fmt1  = mpld.DateFormatter('%d')
+    fmt2  = '%Y-%m'
+  elif deltaT.days < 32:
+    Mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 4)))
+    mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 1)))
+    fmt1  = mpld.DateFormatter('%d')
+    fmt2  = '%Y-%m'
+  elif deltaT.days < 60:
+    Mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 7)))
+    mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 1)))
+    fmt1  = mpld.DateFormatter('%d')
+    fmt2  = '%Y-%m'
+  elif deltaT.days < 183:
+    Mtick = mpld.MonthLocator(bymonth=list(range(1, 13, 1)))
+    mtick = mpld.DayLocator(bymonthday=list(range(1, 32, 7)))
+    fmt1  = mpld.DateFormatter('%m')
+    fmt2  = '%Y'
+  elif deltaT.days < 367:
+    Mtick = mpld.MonthLocator(bymonth=list(range(1, 13, 1)))
+    mtick = mpld.MonthLocator(bymonth=list(range(1, 13, 1)))
+    fmt1  = mpld.DateFormatter('%m')
+    fmt2  = '%Y'
+  elif deltaT.days < 366*2:
+    Mtick = mpld.MonthLocator(bymonth=list(range(1, 13, 2)))
+    mtick = mpld.MonthLocator(bymonth=list(range(1, 13, 1)))
+    fmt1  = mpld.DateFormatter('%m')
+    fmt2  = '%Y'
+  elif deltaT.days < 366*8:
+    Mtick = mpld.YearLocator(1)
+    mtick = mpld.MonthLocator(bymonth=list(range(1, 13, 4)))
+    fmt1  = mpld.DateFormatter('%Y')
+    fmt2  = ''
+  elif deltaT.days < 366*15:
+    to = axes.lines[0].get_xdata()[0]
+    tf = axes.lines[0].get_xdata()[-1]
+    # Ideally would set byyear=list(range(to.year, tf.year,2)) but
+    # byyear is not a kwarg. Would need to something like
+    # https://stackoverflow.com/questions/48428729/matplotlib-dates-yearlocator-with-odd-intervals
+    Mtick = mpld.YearLocator(1)
+    mtick = mpld.YearLocator(1)
+    fmt1  = mpld.DateFormatter('%Y')
+    fmt2  = ''
+    if False:
+      xl = axes.get_xlim()
+      a = mpld.num2date(xl[0])
+      print(a)
+      import pdb;pdb.set_trace()
+      a = mpld.date2num(a.replace(month=1, day=1, hour=0, minute=0, second=0, microsecond=0))
+      b = mpld.num2date(xl[1])
+      b = mpld.date2num(b.replace(year=(b.year+1), month=1, day=1, hour=0, minute=0, second=0, microsecond=0))
+      axes.set_xlim([a, b])
+  elif deltaT.days < 366*40:
+    Mtick = mpld.YearLocator(5)
+    mtick = mpld.YearLocator(1)
+    fmt1  = mpld.DateFormatter('%Y')
+    fmt2  = ''
+  elif deltaT.days < 366*100:
+    Mtick = mpld.YearLocator(10)
+    mtick = mpld.YearLocator(2)
+    fmt1  = mpld.DateFormatter('%Y')
+    fmt2  = ''
+  elif deltaT.days < 366*200:
+    Mtick = mpld.YearLocator(20)
+    mtick = mpld.YearLocator(5)
+    fmt1  = mpld.DateFormatter('%Y')
+    fmt2  = ''
+  else:
+    Mtick = mpld.YearLocator(50)
+    mtick = mpld.YearLocator(10)
+    fmt1  = mpld.DateFormatter('%Y')
+    fmt2  = ''
+
+  return Mtick, mtick, fmt1, fmt2
+
+
+def _adjust_xlabels(axes, adjust_first_xlabel=False, adjust_last_xlabel=False, debug=False):
+  xticklabels = axes.get_xticklabels()
+  lastlabel_text = xticklabels[-1].get_text()
+  firstlabel_text = xticklabels[0].get_text()
+  adjusted = False
+
+  if adjust_first_xlabel and '\n' in firstlabel_text:
+    if len(firstlabel_text.split('\n')[-1]) > 7:
+      if debug:
+        print(f'Adjusting first x-label: "{firstlabel_text.replace("\n", "\\n")}"')
+      adjusted = True
+      # If fmt1 in first label longer than YYYY-MM, set justification to left.
+      xticklabels[0].set_ha('left')
+      # Shift to right by 1/2 width of fmt1 in first label.
+      first_fmt1 = firstlabel_text.split('\n')[0]
+      offset = _numsize(axes, first_fmt1, +1, debug=debug)
+      xticklabels[0].set_transform(xticklabels[0].get_transform() - offset)
+
+  if adjust_last_xlabel and '\n' in lastlabel_text:
+    if len(lastlabel_text.split('\n')[-1]) > 7:
+      if debug:
+        print(f'Adjusting last x-label: "{lastlabel_text.replace("\n", "\\n")}"')
+      adjusted = True
+      # If fmt1 iin last label longer than YYYY-MM, set justification to right.
+      xticklabels[-1].set_ha('right')
+      # Shift to left by 1/2 width of fmt1 in last label.
+      last_fmt1 = lastlabel_text.split('\n')[0]
+      offset = _numsize(axes, last_fmt1, -1, debug=debug)
+      xticklabels[-1].set_transform(xticklabels[-1].get_transform() - offset)
+
+  if adjusted:
+    # Make all labels without newline slightly smaller than default fontsize 
+    # so it is clearer that fmt2 applies to larger number.
+    if debug:
+      print('Adjusting font size of x-labels without newline')
+    for label in xticklabels:
+      if '\n' not in label.get_text():
+        label.set_fontsize(label.get_fontsize()*0.85)
+
+
+def _numsize(ax, num, sign, debug=False):
+  '''Returns (width, height) of str(num) in pixels.
+
+  If ax is given, measures against that axes' renderer and DPI (correct).
+  Otherwise creates a temporary figure using rcParams figure.dpi.
+  '''
+  import re
+  import matplotlib
+  import matplotlib.figure
+  import matplotlib.backends.backend_agg
+
+  num = str(num)
+
+  #dpi = ax.figure.get_dpi() if ax is not None else matplotlib.rcParams['figure.dpi']
+  dpi = 72 # Why not use above dpi? On OS-X when dpi = 200 is returned, offset is wrong.
+  fig = matplotlib.figure.Figure(dpi=dpi)
+  canvas = matplotlib.backends.backend_agg.FigureCanvasAgg(fig)
+  ax_tmp = fig.add_subplot(111)
+  renderer = canvas.get_renderer()
+  fontsize = matplotlib.rcParams['xtick.labelsize']
+  t = ax_tmp.text(0.5, 0.5, num, fontsize=fontsize)
+
+  w, h, d = renderer.get_text_width_height_descent(num, t.get_fontproperties(), ismath=False)
+  dpi = ax.figure.get_dpi() if ax is not None else matplotlib.rcParams['figure.dpi']
+  delta = w/len(num)
+  # It seems like sign should not be needed, but does not work if
+  # + offset instead of - offset is used code that uses offset.
+  offset =  matplotlib.transforms.ScaledTranslation(sign*delta/dpi, 0, ax.figure.dpi_scale_trans)
+
+  if debug:
+    print('_numsize():')
+    print(f'  num      = "{num}"')
+    print(f'  fontsize = {fontsize}')
+    print(f'  dpi      = {dpi}')
+    print(f'  width    = {w}')
+    print(f'  height   = {h}')
+    print(f'  descent  = {d}')
+    print(f'  delta    = {delta}')
+    print(f'  offset   = {re.sub(r"\n\s+", "", str(offset))}')
+
+  return offset

datetick-0.9.0.dist-info/METADATA

@@ -0,0 +1,216 @@
+Metadata-Version: 2.4
+Name: datetick
+Version: 0.9.0
+Summary: Sensible date tick locator for matplotlib
+Author: Brendan Gallagher
+Author-email: Bob Weigel <rweigel@gmu.edu>
+License-Expression: MIT
+Project-URL: Homepage, http://pypi.python.org/pypi/datetick/
+Requires-Python: >=3.7
+Description-Content-Type: text/markdown
+License-File: LICENSE.txt
+Requires-Dist: matplotlib
+Requires-Dist: numpy
+Requires-Dist: python-dateutil
+Dynamic: license-file
+
+# datetick
+
+Sensible numeric time and date tick labels for Matplotlib
+
+# Motivation
+
+Matplotlib's default time tick labels are often poor, and adjusting them requires using [locators and formatters](https://matplotlib.org/stable/api/ticker_api.html) on an ad-hoc basis. In addition, the interfaces for locators and formatters complex and non-intuitive and require study and experimentation.
+
+`datetick()` contains logic for locators and formatters that apply to plots with arbitrary time ranges. One only needs to add the command `datetick()` after the usual Matplotlib `plt.plot(...)` command to have sensible and useable time tick labels.
+
+# Usage
+
+```
+import datetime as dt
+import matplotlib.pyplot as plt
+from datetick import datetick
+
+dt1 = dt.datetime(2011, 1, 2)
+dt2 = dt1 + dt.timedelta(days=1, hours=1, minutes=1)
+
+plt.plot([dt1, dt2], [0.0,1.0])
+datetick()
+plt.show()
+# or
+# datetick('x') (use 'y' if y variable is a time)
+# or
+# datetick('x', axes=plt.gca())
+# or
+# fig, axes = plt.subplots(2)
+# plt.plot([dt1, dt2], [0.0, 1.0])
+# datetick('x', axes=axes[0])
+```
+
+# Comparison to default Matplotlib
+
+![test/mpl-3.10.9/20010101000000.0-20010101000000.1.svg](test/mpl-3.10.9/20010101000000.0-20010101000000.1.svg)
+
+![test/mpl-3.10.9/20010101000000.0-20010101000000.2.svg](test/mpl-3.10.9/20010101000000.0-20010101000000.2.svg)
+
+![test/mpl-3.10.9/20010101000000.0-20010101000000.5.svg](test/mpl-3.10.9/20010101000000.0-20010101000000.5.svg)
+
+![test/mpl-3.10.9/20010101000000.01-20010101000000.99.svg](test/mpl-3.10.9/20010101000000.01-20010101000000.99.svg)
+
+![test/mpl-3.10.9/20010101000059.8-20010101000100.3.svg](test/mpl-3.10.9/20010101000059.8-20010101000100.3.svg)
+
+![test/mpl-3.10.9/20010101005959.8-20010101010003.0.svg](test/mpl-3.10.9/20010101005959.8-20010101010003.0.svg)
+
+![test/mpl-3.10.9/20010101235959.8-20010102000003.0.svg](test/mpl-3.10.9/20010101235959.8-20010102000003.0.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101000001.svg](test/mpl-3.10.9/20010101000000-20010101000001.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101000005.svg](test/mpl-3.10.9/20010101000000-20010101000005.svg)
+
+![test/mpl-3.10.9/20010101000059-20010101000103.svg](test/mpl-3.10.9/20010101000059-20010101000103.svg)
+
+![test/mpl-3.10.9/20010101005958-20010101010003.svg](test/mpl-3.10.9/20010101005958-20010101010003.svg)
+
+![test/mpl-3.10.9/20010101005959-20010101010004.svg](test/mpl-3.10.9/20010101005959-20010101010004.svg)
+
+![test/mpl-3.10.9/20010101235958-20010102000003.svg](test/mpl-3.10.9/20010101235958-20010102000003.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101000007.svg](test/mpl-3.10.9/20010101000000-20010101000007.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101000009.svg](test/mpl-3.10.9/20010101000000-20010101000009.svg)
+
+![test/mpl-3.10.9/20010101000001-20010101000006.svg](test/mpl-3.10.9/20010101000001-20010101000006.svg)
+
+![test/mpl-3.10.9/20010101235958-20010102000005.svg](test/mpl-3.10.9/20010101235958-20010102000005.svg)
+
+![test/mpl-3.10.9/20010101235959-20010102000005.svg](test/mpl-3.10.9/20010101235959-20010102000005.svg)
+
+![test/mpl-3.10.9/20010101005958-20010101010005.svg](test/mpl-3.10.9/20010101005958-20010101010005.svg)
+
+![test/mpl-3.10.9/20010101005958-20010101010007.svg](test/mpl-3.10.9/20010101005958-20010101010007.svg)
+
+![test/mpl-3.10.9/20010101005958-20010101010003_2.svg](test/mpl-3.10.9/20010101005958-20010101010003_2.svg)
+
+![test/mpl-3.10.9/20010101235956-20010102000010.svg](test/mpl-3.10.9/20010101235956-20010102000010.svg)
+
+![test/mpl-3.10.9/20010101235958-20010102000010.svg](test/mpl-3.10.9/20010101235958-20010102000010.svg)
+
+![test/mpl-3.10.9/20010101005956-20010101010010.svg](test/mpl-3.10.9/20010101005956-20010101010010.svg)
+
+![test/mpl-3.10.9/20010101005956-20010101010015.svg](test/mpl-3.10.9/20010101005956-20010101010015.svg)
+
+![test/mpl-3.10.9/20010101005956-20010101010006.svg](test/mpl-3.10.9/20010101005956-20010101010006.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101000021.svg](test/mpl-3.10.9/20010101000000-20010101000021.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101060000.svg](test/mpl-3.10.9/20010101000000-20010101060000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101090000.svg](test/mpl-3.10.9/20010101000000-20010101090000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101110000.svg](test/mpl-3.10.9/20010101000000-20010101110000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101110000_2.svg](test/mpl-3.10.9/20010101000000-20010101110000_2.svg)
+
+![test/mpl-3.10.9/20010101005958-20010101010028.svg](test/mpl-3.10.9/20010101005958-20010101010028.svg)
+
+![test/mpl-3.10.9/20010101000058-20010101000118.svg](test/mpl-3.10.9/20010101000058-20010101000118.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101180000.svg](test/mpl-3.10.9/20010101000000-20010101180000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101120000.svg](test/mpl-3.10.9/20010101000000-20010101120000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010101230000.svg](test/mpl-3.10.9/20010101000000-20010101230000.svg)
+
+![test/mpl-3.10.9/20010101020000-20010102010000.svg](test/mpl-3.10.9/20010101020000-20010102010000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010102010000.svg](test/mpl-3.10.9/20010101000000-20010102010000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010102230000.svg](test/mpl-3.10.9/20010101000000-20010102230000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010102230000_2.svg](test/mpl-3.10.9/20010101000000-20010102230000_2.svg)
+
+![test/mpl-3.10.9/20010101000000-20010102230000_3.svg](test/mpl-3.10.9/20010101000000-20010102230000_3.svg)
+
+![test/mpl-3.10.9/20010101000000-20010102230000_4.svg](test/mpl-3.10.9/20010101000000-20010102230000_4.svg)
+
+![test/mpl-3.10.9/20010101060000-20010102070000.svg](test/mpl-3.10.9/20010101060000-20010102070000.svg)
+
+![test/mpl-3.10.9/20010101003000-20010102010000.svg](test/mpl-3.10.9/20010101003000-20010102010000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010103000000.svg](test/mpl-3.10.9/20010101000000-20010103000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010103123000.svg](test/mpl-3.10.9/20010101000000-20010103123000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010103235959.svg](test/mpl-3.10.9/20010101000000-20010103235959.svg)
+
+![test/mpl-3.10.9/20010101000000-20010105000000.svg](test/mpl-3.10.9/20010101000000-20010105000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010105000000_2.svg](test/mpl-3.10.9/20010101000000-20010105000000_2.svg)
+
+![test/mpl-3.10.9/20010130000000-20010201230000.svg](test/mpl-3.10.9/20010130000000-20010201230000.svg)
+
+![test/mpl-3.10.9/20011230000000-20020101230000.svg](test/mpl-3.10.9/20011230000000-20020101230000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010109000000.svg](test/mpl-3.10.9/20010101000000-20010109000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010116230000.svg](test/mpl-3.10.9/20010101000000-20010116230000.svg)
+
+![test/mpl-3.10.9/20010130000000-20010204230000.svg](test/mpl-3.10.9/20010130000000-20010204230000.svg)
+
+![test/mpl-3.10.9/20011230000000-20020104230000.svg](test/mpl-3.10.9/20011230000000-20020104230000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010131000000.svg](test/mpl-3.10.9/20010101000000-20010131000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010116230000_2.svg](test/mpl-3.10.9/20010101000000-20010116230000_2.svg)
+
+![test/mpl-3.10.9/20010130000000-20010215230000.svg](test/mpl-3.10.9/20010130000000-20010215230000.svg)
+
+![test/mpl-3.10.9/20011230000000-20010115230000.svg](test/mpl-3.10.9/20011230000000-20010115230000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010202000000.svg](test/mpl-3.10.9/20010101000000-20010202000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010227230000.svg](test/mpl-3.10.9/20010101000000-20010227230000.svg)
+
+![test/mpl-3.10.9/20010115000000-20010216230000.svg](test/mpl-3.10.9/20010115000000-20010216230000.svg)
+
+![test/mpl-3.10.9/20011231000000-20020226230000.svg](test/mpl-3.10.9/20011231000000-20020226230000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010502000000.svg](test/mpl-3.10.9/20010101000000-20010502000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20010227230000_2.svg](test/mpl-3.10.9/20010101000000-20010227230000_2.svg)
+
+![test/mpl-3.10.9/20011231000000-20020226230000_2.svg](test/mpl-3.10.9/20011231000000-20020226230000_2.svg)
+
+![test/mpl-3.10.9/20010101000000-20010702000000.svg](test/mpl-3.10.9/20010101000000-20010702000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20011231000000.svg](test/mpl-3.10.9/20010101000000-20011231000000.svg)
+
+![test/mpl-3.10.9/20010212000000-20020131000000.svg](test/mpl-3.10.9/20010212000000-20020131000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20020103000000.svg](test/mpl-3.10.9/20010101000000-20020103000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20021231000000.svg](test/mpl-3.10.9/20010101000000-20021231000000.svg)
+
+![test/mpl-3.10.9/20010401000000-20020430000000.svg](test/mpl-3.10.9/20010401000000-20020430000000.svg)
+
+![test/mpl-3.10.9/20011001000000-20031004000000.svg](test/mpl-3.10.9/20011001000000-20031004000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20081231000000.svg](test/mpl-3.10.9/20010101000000-20081231000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20030104000000.svg](test/mpl-3.10.9/20010101000000-20030104000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20090104000000.svg](test/mpl-3.10.9/20010101000000-20090104000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20120104000000.svg](test/mpl-3.10.9/20010101000000-20120104000000.svg)
+
+![test/mpl-3.10.9/20000101000000-20170104000000.svg](test/mpl-3.10.9/20000101000000-20170104000000.svg)
+
+![test/mpl-3.10.9/20010101000000-20180104000000.svg](test/mpl-3.10.9/20010101000000-20180104000000.svg)
+
+![test/mpl-3.10.9/20020101000000-20190104000000.svg](test/mpl-3.10.9/20020101000000-20190104000000.svg)
+
+![test/mpl-3.10.9/20030101000000-20200104000000.svg](test/mpl-3.10.9/20030101000000-20200104000000.svg)
+
+![test/mpl-3.10.9/20040101000000-20300104000000.svg](test/mpl-3.10.9/20040101000000-20300104000000.svg)
+
+![test/mpl-3.10.9/19500101000000-20120104000000.svg](test/mpl-3.10.9/19500101000000-20120104000000.svg)

datetick-0.9.0.dist-info/RECORD

@@ -0,0 +1,9 @@
+datetick/__init__.py,sha256=Q_1eN-e8ewVpwdzExrVDYCIVwFJ4Le_AueDsdjoLZBA,200
+datetick/datetick.py,sha256=g_ejtah9mvsctC9n27-jOczeCwvSClNhNRv8pIoJims,24860
+datetick-0.9.0.dist-info/licenses/LICENSE.txt,sha256=UyVxGUH2kA6ZPDrjS3dtWGk5cUnsLnY7OlNSabyeCYA,1087
+etc/release.py,sha256=-ApOOmegDxt4QUBc_ATr6YgHdZufM5ecXyj_eOn41U8,2737
+test/datetick_test.py,sha256=GHHeMA7SUNCkBupOfI56dSE-l3BfoybMTokBQzf39LU,3412
+datetick-0.9.0.dist-info/METADATA,sha256=DQXL_ajkgncha8NfMfGdAYfIHFe-P1NKOnvy2MJEdzA,10368
+datetick-0.9.0.dist-info/WHEEL,sha256=aeYiig01lYGDzBgS8HxWXOg3uV61G9ijOsup-k9o1sk,91
+datetick-0.9.0.dist-info/top_level.txt,sha256=1ZRIT2tCJ-2z6LU-rrxHoJe2JmzN2uIN4WCZYolEAxg,18
+datetick-0.9.0.dist-info/RECORD,,

datetick-0.9.0.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (82.0.1)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

datetick-0.9.0.dist-info/licenses/LICENSE.txt

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026- Bob Weigel, Brendan Gallagher
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

datetick-0.9.0.dist-info/top_level.txt

@@ -0,0 +1,3 @@
+datetick
+etc
+test

etc/release.py

@@ -0,0 +1,93 @@
+"""
+Usage:
+  python release.py
+    [path/to/pyproject.toml]
+    [--pypi-config-file path/to/pypirc]
+
+Creates GitHub and PyPi releases with the same version number in pyproject.toml.
+
+If path/to/pyproject.toml is not specified, it will look for it in the parent
+directory of this script.
+
+Requires:
+  A GitHub token file at ~/git/admin/etc/github_token (plain text, one line).
+  A PyPI config file at ~/git/admin/etc/pypirc
+"""
+
+import os
+import sys
+import subprocess
+
+try:
+  import tomllib
+except ImportError:
+  try:
+    import tomli as tomllib
+  except ImportError:
+    cmd_list = [sys.executable, '-m', 'pip', 'install', 'tomli']
+    print("Executing:", " ".join(cmd_list))
+    subprocess.run(cmd_list, check=True)
+    import tomli as tomllib
+
+def _install_deps():
+  try:
+    import build  # noqa: F401
+  except ImportError:
+    cmd_list = [sys.executable, '-m', 'pip', 'install', 'build']
+    print("Executing:", " ".join(cmd_list))
+    subprocess.run(cmd_list, check=True)
+
+  try:
+    import twine  # noqa: F401
+  except ImportError:
+    cmd_list = [sys.executable, '-m', 'pip', 'install', 'twine']
+    print("Executing:", " ".join(cmd_list))
+    subprocess.run(cmd_list, check=True)
+
+def main(pypi_config_file):
+  if len(sys.argv) > 1:
+    toml_path = sys.argv[1]
+  else:
+    toml_path = os.path.join(os.path.dirname(__file__), '..', 'pyproject.toml')
+  toml_path = os.path.abspath(toml_path)
+
+  if not os.path.exists(toml_path):
+    print(f"Could not find pyproject.toml at {toml_path}")
+    sys.exit(1)
+
+  with open(toml_path, 'rb') as f:
+    data = tomllib.load(f)
+
+  version = data.get('project', {}).get('version')
+  if not version:
+    print(f"Could not find [project] version in {toml_path}")
+    sys.exit(1)
+
+  print(f"Creating release {version} on GitHub and PyPi")
+
+  # Check if gh CLI is installed
+  try:
+    subprocess.run(['gh', '--version'], check=True, stdout=subprocess.DEVNULL)
+  except (subprocess.CalledProcessError, FileNotFoundError):
+    print("Error: gh CLI is not installed. Please install it from https://cli.github.com/")
+    sys.exit(1)
+
+  cmd_list = ['gh', 'release', 'create', version, '-t', version, '-n', f'Release {version}']
+  #print("Executing:", " ".join(cmd_list))
+  #subprocess.run(cmd_list, check=True)
+
+  # Build package
+  cmd_list = [sys.executable, '-m', 'build']
+  print("Executing:", " ".join(cmd_list))
+  subprocess.run(cmd_list, check=True)
+  exit()
+
+  # Upload package to PyPi
+  cmd_list = ['twine', 'upload', '--config-file', pypi_config_file, 'dist/*']
+  subprocess.run(cmd_list, check=True)
+
+if __name__ == '__main__':
+  _install_deps()
+  pypi_config_path = os.path.expanduser('~/git/admin/etc')
+  pypi_config_file = os.path.join(pypi_config_path, 'pypirc')
+  main(pypi_config_file)

test/datetick_test.py

@@ -0,0 +1,130 @@
+# Create plots with varying time ranges.
+
+import os
+import yaml
+import dateutil.parser
+import matplotlib.pyplot as plt
+
+from datetick import datetick
+
+debug = False
+
+def plot(ds1, ds2, **kwargs):
+  _, axes = plt.subplots(2, figsize=(8, 2))
+  plt.subplots_adjust(hspace=1.0)
+
+  dt1 = dateutil.parser.parse(ds1)
+  dt2 = dateutil.parser.parse(ds2)
+  x = [dt1, dt2]
+  xt = x[0] + (x[1] - x[0])/2
+  y = [0.0,0.0]
+
+  bbox = dict(boxstyle='round,pad=0.3', facecolor='white', edgecolor='gray', alpha=0.8)
+
+  axes[0].set_title(ds1 + ' - ' + ds2, fontfamily='monospace')
+  axes[0].plot(x, y, '*')
+  axes[0].text(xt, 0.00, 'matplotlib', ha='center', bbox=bbox)
+
+  text = 'datetick'
+  if kwargs:
+    bbox['facecolor'] = 'lightblue'
+    for key, value in kwargs.items():
+      text += f'\n{key}={value}'
+  axes[1].plot(x, y, '*')
+  axes[1].text(xt, 0.00, text, ha='center', bbox=bbox)
+  datetick('x', axes=axes[1], debug=debug, **kwargs)
+
+  for axis in axes:
+    axis.grid()
+    axis.spines[['top', 'right', 'left']].set_visible(False)
+    axis.yaxis.set_visible(False)
+
+  file = _savefig(ds1, ds2)
+
+  return file
+
+
+def _savefig(ds1, ds2):
+  ds1 = ds1.replace(":","").replace("-","").replace("T","").replace("Z","")
+  ds2 = ds2.replace(":","").replace("-","").replace("T","").replace("Z","")
+
+  file = f'{_script_dir()}/{_mpl_version()}/{ds1}-{ds2}.svg'
+  if file in files:
+    v = 2
+    while file in files:
+      file = f'{_script_dir()}/{_mpl_version()}/{ds1}-{ds2}_{v}.svg'
+      v += 1
+
+  print("Writing", file)
+  dirname = os.path.dirname(file)
+  if not os.path.exists(dirname):
+    os.makedirs(dirname, exist_ok=True)
+
+  plt.savefig(file, bbox_inches='tight')
+  plt.close()
+
+  return file
+
+
+def _append_to_readme(files):
+
+  readme = 'README.md' # Repo README
+  readme = os.path.join(_script_dir(), "..", readme)
+
+  with open(readme, 'r+') as file:
+    lines = file.readlines()
+
+  index = next(i for i, line in enumerate(lines) if "Comparison to default Matplotlib" in line)
+  del lines[index+1:]
+
+  image_links = []
+  for file in files:
+    # Make path relative to README
+    file = os.path.relpath(file, os.path.dirname(readme))
+    image_links.append(f'![{file}]({file})')
+
+  lines.append("\n" + "\n\n".join(image_links))
+
+  print(f"Updating {readme} with {len(files)} images")
+  with open(readme, 'w') as file:
+    file.writelines(lines)
+
+def _create_subdir_readme(files):
+  # Create README in mpl subdir
+  image_links = []
+  for file in files:
+    # Make path relative to README
+    file = os.path.basename(file)
+    image_links.append(f'![{file}]({file})')
+
+  readme = os.path.join(_script_dir(), _mpl_version(), 'README.md')
+  print(f"Writing {readme} with {len(files)} images")
+  with open(readme, 'w') as file:
+    file.writelines("\n" + "\n\n".join(image_links))
+
+
+def _mpl_version():
+  return f"mpl-{plt.matplotlib.__version__}"
+
+
+def _script_dir():
+  return os.path.dirname(os.path.realpath(__file__))
+
+
+if False:
+  plot('2001-01-01T00:00:58Z','2001-01-01T00:01:18Z', adjust_xrange=True)
+  #plot('2001-01-01T00:00:00Z','2001-01-02T23:00:00Z', adjust_first_xlabel=True)
+  exit()
+
+test_file = os.path.join(_script_dir(), 'datetick_test.yaml')
+with open(test_file, 'r') as file:
+  tests = yaml.safe_load(file)
+
+files = []
+for test in tests['main']:
+  kwargs = test[2] if len(test) > 2 else {}
+  file = plot(test[0], test[1], **kwargs)
+  files.append(file)
+
+_append_to_readme(files)
+_create_subdir_readme(files)