OPA_App_v2.3.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Observability of Primordial Asteroids Application\n",
    "## (Primo Pole Project)\n",
    "In the frame of <a href=\"https://mp3c.oca.eu/catalogue/projects.htm\" target=\"_blank\">Primo Pole Project </a> leaded by Dr. Marco Delbo from Observatoire de la Côte d'Azur (OCA), asteroids that belong to the Primordial family need to be observed, thus their spin pole to be estimated. This family was discovered by Delbo et al. in 2017 and is published as \"*Delbo, M., Walsh, K., Bolin, B., Avdellidou, C., & Morbidelli, A. (2017). Identification of a primordial asteroid family constrains the original planetesimal population. Science, 357(6355).*\". This family could be as old as the Solar System, i.e. ~4.5 Gyears old.\n",
    "The working hypothesis is that according to theories of asteroid orbital evolution, under the influence of the Yarkovsky effect, members of the inward side of the Primordial family should have a statistical predominance of retrograde objects.\n",
    "\n",
    "**This code finds the observable asteroids, members of the Primordial family, based on the location, the observation period and instrument limitations.** In this way, a short list of candidate targets will be exported automatically providing us the capability to do a smart sampling.   "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "code_folding": [
     0,
     36
    ]
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/srv/conda/envs/notebook/lib/python3.7/site-packages/torch/cuda/__init__.py:52: UserWarning: CUDA initialization: Found no NVIDIA driver on your system. Please check that you have an NVIDIA GPU and installed a driver from http://www.nvidia.com/Download/index.aspx (Triggered internally at  /pytorch/c10/cuda/CUDAFunctions.cpp:100.)\n",
      "  return torch._C._cuda_getDeviceCount() > 0\n"
     ]
    }
   ],
   "source": [
    "#Useful Packages and Functions\n",
    "from fastai.vision.all import *\n",
    "from fastai.vision.widgets import *\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import astropy\n",
    "from astroquery.jplhorizons import Horizons\n",
    "import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.ticker import MultipleLocator\n",
    "import datetime as dt #Package for the dates\n",
    "from scipy.stats import norm #Package for statistics\n",
    "#import seaborn as sns #Package Styling for Statistics\n",
    "from bisect import bisect_left\n",
    "\n",
    "#To export plots as images\n",
    "import io\n",
    "import urllib, base64\n",
    "import time\n",
    "\n",
    "#Packages for interaction\n",
    "import ipywidgets as ipw\n",
    "from ipywidgets import interact, interactive, fixed, interact_manual\n",
    "from IPython.core.display import HTML, display, Javascript, clear_output\n",
    "\n",
    "\n",
    "progress=ipw.IntProgress(\n",
    "    value=0,\n",
    "    min=0,\n",
    "    max=10,\n",
    "    step=1,\n",
    "    description='Searching:',\n",
    "    bar_style='', # 'success', 'info', 'warning', 'danger' or ''\n",
    "    orientation='horizontal'\n",
    ")\n",
    "\n",
    "\n",
    "def take_closest(myList, myNumber):\n",
    "    \"\"\"\n",
    "    Assumes myList is sorted. Returns closest value to myNumber.\n",
    "\n",
    "    If two numbers are equally close, return the smallest number.\n",
    "    \"\"\"\n",
    "    pos = bisect_left(myList, myNumber)\n",
    "    if pos == 0:\n",
    "        return myList[0]\n",
    "    if pos == len(myList):\n",
    "        return myList[-1]\n",
    "    before = pos - 1\n",
    "    after = pos\n",
    "    if myList[after] - myNumber < myNumber - myList[before]:\n",
    "       return after\n",
    "    else:\n",
    "       return before\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "code_folding": [
     98
    ]
   },
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "fcd6d3269ce945ce98ac9b86ab88d21c",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "VBox(children=(HTML(value='<p style=\"text-align: center;\"><span style=\"font-size: 40px; font-family: Helvetica…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "2e83b3fa0f9e45ac85f3542f8192533b",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "VBox(children=(HTML(value='', layout=Layout(height='20px')), VBox(children=(HTML(value='<p style=\"text-align: …"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "bc3f1845cd694ed388d3b63bab2f36ae",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "VBox(children=(HTML(value='', layout=Layout(height='20px')), HTML(value='<p style=\"text-align: left;\"><strong>…"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "6e4aae63f77a4ae3b73cc3aed41c6633",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "Output()"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "################################################ MAIN CODE ########################################\n",
    "path = Path()\n",
    "output = widgets.Output()\n",
    "\n",
    "#TITLE\n",
    "Title=ipw.HTML('<p style=\"text-align: center;\"><span style=\"font-size: 40px; font-family: Helvetica;\"><span style=\"color: rgb(85, 57, 130);\"><strong>O</strong>bservability of <strong>P</strong>rimordial <strong>A</strong>steroids</span> <strong>App</strong>lication</span></p>')\n",
    "#Title.layout.height = '80px'\n",
    "#Title.layout.align_content = 'center'\n",
    "\n",
    "#COPYRIGHTS\n",
    "Copyright=ipw.HTML('<footer style=\"text-align: right;\">&copy; 2021 Dimitrios Athanasopoulos </footer>')\n",
    "Copyright.layout.height = '15px'\n",
    "#CONTACT\n",
    "Contact=ipw.HTML('<footer style=\"text-align: right;\">Contact: dimathanaso@phys.uoa.gr </footer>')\n",
    "Contact.layout.height = '30px'\n",
    "\n",
    "################################################ A PART ####################################################\n",
    "\n",
    "#TAB_1\n",
    "Location_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 18px; font-family: Helvetica;\">Insert your Location and Timezone</span></strong></p>')\n",
    "LON_ft=ipw.FloatText(description=\"Longitude (deg):\", value=23.783368, style={'description_width': '100px'})\n",
    "LAT_ft=ipw.FloatText(description=\"Latitude (deg):\", value=37.968561, style={'description_width': '100px'})\n",
    "ALT_ft=ipw.FloatText(description=\"Altidute (km):\", value=0.250, style={'description_width': '100px'})\n",
    "TimeZone_dd=ipw.Dropdown(description=\"TimeZone:\", options=['Etc/GMT+12', 'Etc/GMT+11', 'Etc/GMT+10', 'Etc/GMT+9', 'Etc/GMT+8', 'Etc/GMT+7', 'Etc/GMT+6',\n",
    "                               'Etc/GMT+5', 'Etc/GMT+4', 'Etc/GMT+3', 'Etc/GMT+2', 'Etc/GMT+1', 'UTC',\n",
    "                               'Etc/GMT-1', 'Etc/GMT-2', 'Etc/GMT-3', 'Etc/GMT-4', 'Etc/GMT-5', 'Etc/GMT-6', 'Etc/GMT-7', \n",
    "                               'Etc/GMT-8', 'Etc/GMT-9', 'Etc/GMT-10', 'Etc/GMT-11', 'Etc/GMT-12', 'Etc/GMT-13', 'Etc/GMT-14'],\n",
    "                      value='UTC', style={'description_width': '100px'})\n",
    "OutputMessage1a=ipw.HTML(value=\"\")\n",
    "OutputMessage1b=ipw.HTML(value=\"\")\n",
    "\n",
    "def location_input(b=None):\n",
    "    global LON, LAT, ALT, TimeZone, OutputMessage1a, InsertLocation\n",
    "    LON=LON_ft.value; LAT=LAT_ft.value; ALT=ALT_ft.value; TimeZone=TimeZone_dd.value\n",
    "    clear_output(wait=True)\n",
    "    OutputMessage1a.value=\"Your Location is: Longitude: \" + str(LON) + \" deg --- Latitude: \" + str(LAT) + \" deg --- Altitude: \" + str(ALT) + \" km\"\n",
    "    OutputMessage1b.value=\"Your TimeZone is: \" + TimeZone\n",
    "    display(DashBoard1)\n",
    "\n",
    "Insert_btn1=ipw.Button(description=\"Insert\", button_style='Success')\n",
    "\n",
    "Insert_btn1.on_click(location_input)\n",
    "\n",
    "InsertLocation=ipw.VBox([Location_text, LON_ft, LAT_ft, ALT_ft, TimeZone_dd, Insert_btn1, OutputMessage1a, OutputMessage1b])\n",
    "\n",
    "#TAB_2\n",
    "ObservationPeriod_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 18px; font-family: Helvetica;\">Insert your Observation period </span></strong></p>')\n",
    "START_t=ipw.Text(description=\"Start date (YYYY-MM-DD):\", style={'description_width': '180px'})\n",
    "STOP_t=ipw.Text(description=\"Stop date (YYYY-MM-DD):\", style={'description_width': '180px'})\n",
    "OutputMessage2=ipw.HTML(value=\"\")\n",
    "\n",
    "def observationPeriod_input(b=None):\n",
    "    global START, STOP\n",
    "    START=START_t.value; STOP=STOP_t.value\n",
    "    clear_output(wait=True)\n",
    "    OutputMessage2.value=\"Your Observation Period is from \" + START + \" at 12:00 (UTC) to \" + STOP + \" at 12:00 (UTC).\"\n",
    "    display(DashBoard1)\n",
    "\n",
    "\n",
    "Insert_btn2=ipw.Button(description=\"Insert\", button_style='Success')\n",
    "\n",
    "Insert_btn2.on_click(observationPeriod_input)\n",
    "\n",
    "InsertObservationPeriod=ipw.VBox([ObservationPeriod_text, START_t, STOP_t, Insert_btn2, OutputMessage2])\n",
    "\n",
    "#TAB_3\n",
    "Preferences_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 18px; font-family: Helvetica;\">Insert your instrument limitations and preferences</span></strong></p>')\n",
    "APmag_min_ft=ipw.FloatText(description=\"The target max. ap. magnitude (mag):\", value=17, style={'description_width': '230px'})\n",
    "MinObsTime_ft=ipw.FloatText(description=\"The minimum observation time (hours):\", value=4, style={'description_width': '230px'})\n",
    "OutputMessage3=ipw.HTML(\"\")\n",
    "\n",
    "def preferences_input(b=None):\n",
    "    global APmag_min, MinObsTime\n",
    "    \n",
    "    APmag_min=APmag_min_ft.value; MinObsTime=MinObsTime_ft.value\n",
    "    clear_output(wait=True)\n",
    "    OutputMessage3.value=\"You want to observe targets with ≤ \" + str(APmag_min) + \" mag for at least \" + str(MinObsTime) + \" hours.\"\n",
    "    display(DashBoard1)\n",
    "    #display(InsertPreferences)\n",
    "    #print(\"You want to observe targets with ≤ \" + str(APmag_min) + \" mag for at least \" + str(MinObsTime) + \" hours.\")\n",
    "\n",
    "\n",
    "Insert_btn3=ipw.Button(description=\"Insert\", button_style='Success')\n",
    "Insert_btn3.on_click(preferences_input)\n",
    "\n",
    "InsertPreferences=ipw.VBox([Preferences_text, APmag_min_ft, MinObsTime_ft, Insert_btn3, OutputMessage3])\n",
    "\n",
    "\n",
    "#ALL_TABS\n",
    "InsertTab=ipw.Tab([InsertLocation, InsertObservationPeriod, InsertPreferences])\n",
    "InsertTab.set_title(0, 'Your Location')\n",
    "InsertTab.set_title(1, 'Observation Period')\n",
    "InsertTab.set_title(2, 'Preferences')\n",
    "\n",
    "IDtargetList=[]\n",
    "\n",
    "#FIND_OBSERVABLE_ASTEROIDS\n",
    "@output.capture(clear_output=False,wait=True)\n",
    "def find_asteroids(b=None):\n",
    "    time.sleep(0.1)\n",
    "    DashBoard2.layout.display = 'None'\n",
    "    DashBoard3.layout.display = 'None'\n",
    "    \n",
    "    global Family, KnownFamily, TargetList, IDtargetList, IDtargetListPR, Observatory, ObsDays, OptionTargets, FIG1\n",
    "    ct1=0 #This is for the final button\n",
    "    \n",
    "    if (OutputMessage1a.value==\"\"):\n",
    "        OutputMessage4.value='<p><strong><span style=\"color: red;\">Please insert your Location and Timezone!</span></strong></p>'\n",
    "    elif (OutputMessage2.value==\"\") or (START_t.value==\"\") or (STOP_t.value==\"\"):\n",
    "        OutputMessage4.value='<p><strong><span style=\"color: red;\">Please insert your Observation Period!</span></strong></p>'\n",
    "    elif (OutputMessage3.value==\"\"):\n",
    "        OutputMessage4.value='<p><strong><span style=\"color: red;\">Please insert your Preferences!</span></strong></p>'\n",
    "    else:\n",
    "        Observatory={'lon': LON, 'lat': LAT, 'elevation': ALT}\n",
    "        step=1 #Time Step\n",
    "        units=\"h\"#Time units, type h: for hours, m: for minutes, d: for days\n",
    "        ObsDays={'start':START+' 12:00', 'stop':STOP+' 12:00','step':(str(step)+units)}\n",
    "        \n",
    "        if units=='h':\n",
    "            step=step\n",
    "        elif units=='m':\n",
    "            step=step/60\n",
    "        elif units=='d':\n",
    "            step=step*24\n",
    "        else:\n",
    "            print('WARNING!: Please insert again a time-step with units in minutes(\"m\") or hours (\"h\") or days (\"d\")')\n",
    "        \n",
    "        #READ FILE\n",
    "        Family = pd.read_excel('PrimordialFamily.xls', sheet_name='All_Members')\n",
    "        \n",
    "        OutputMessage4.value='<p><strong><span style=\"color: green;\">Initialisation completed!</span></strong></p>'\n",
    "        \n",
    "        #DASHBOARD_1 (with progress bar)\n",
    "        ObservableAsteroids_text.layout.display = 'None'\n",
    "        TableList.value=\"\"\n",
    "        TableFootNote.layout.display = 'None'\n",
    "        FIG1.value=\"\"\n",
    "        FIG2.value=\"\"\n",
    "        Message1.value=\"\"\n",
    "        Message2.value=\"\"\n",
    "        TableOP.value=\"\"\n",
    "        FIG3.value=\"\"\n",
    "    \n",
    "        #clear_output(wait=True)\n",
    "        progress.layout.visibility = 'visible'\n",
    "        #DashBoard2=ipw.VBox([Title, Copyright, Contact, InsertTab, FindAsteroids, progress])\n",
    "        #display(DashBoard1)\n",
    "        \n",
    "        #Search for Observable Asteroids\n",
    "        TargetList=[]\n",
    "        IDtargetList=[]\n",
    "        IDtargetListPR=[]\n",
    "        \n",
    "        it=0\n",
    "        progress.max=(len(Family['Number'])-1)\n",
    "        for IDtarget in Family['Number']:\n",
    "            progress.value=it\n",
    "            progress.value=it\n",
    "            \n",
    "            #Find epherides from Horizon - JPL Small Bodies Solar System database\n",
    "            obj = Horizons(id=int(IDtarget), location=Observatory, epochs=ObsDays)\n",
    "            eph = obj.ephemerides()\n",
    "            eph=eph.to_pandas()\n",
    "            \n",
    "            #This is the full name of the target\n",
    "            TargetName=eph['targetname'][0]\n",
    "            \n",
    "            #Datetime formatting\n",
    "            Date=eph['datetime_str']\n",
    "            DateT=pd.to_datetime(Date)\n",
    "            if (TimeZone=='UTC') | (TimeZone=='Europe/London'):\n",
    "                LocalDate=DateT\n",
    "            else:\n",
    "                LocalDate = DateT.dt.tz_localize('UTC').dt.tz_convert(TimeZone)\n",
    "            \n",
    "            #Determine the nightime\n",
    "            Sun=eph['solar_presence']\n",
    "            NightTime=[]\n",
    "            n=0\n",
    "            for x in range(len(Sun)):\n",
    "                if (Sun[x]==\"\"):\n",
    "                    NightTime.append(LocalDate[x])\n",
    "                    n=+1\n",
    "            \n",
    "            #The target at nightime\n",
    "            Target=eph[['datetime_str', 'solar_presence','V', 'airmass','RA','DEC','alpha_true','PABLon','PABLat']].copy(deep=True) #Suppress the SettingWithCopyWarning\n",
    "            Target['datetime_str']=pd.to_datetime(Target['datetime_str'])\n",
    "            pd.options.mode.chained_assignment = None  # default='warn'\n",
    "            Target.rename(columns={'datetime_str':'datetime'}, inplace=True)\n",
    "            Target=Target[(Target['airmass']<=3) & (Target['solar_presence']==\"\")]\n",
    "            Target.reset_index(drop=True, inplace=True)\n",
    "            \n",
    "            #To release memory\n",
    "            del eph\n",
    "            \n",
    "            #Duration of observable target\n",
    "            DurationT=[]\n",
    "            if len(Target)>1 and abs(Target[\"datetime\"][0]-Target[\"datetime\"][len(Target[\"datetime\"])-1])<dt.timedelta(hours=24):\n",
    "                check=1\n",
    "                c1=0\n",
    "                c2=len(Target[\"datetime\"])-1\n",
    "                D0=(pd.Timestamp.date(Target[\"datetime\"][c2]))\n",
    "                dd=Target[\"datetime\"][c2]-Target[\"datetime\"][c1]\n",
    "                D1=(dd.total_seconds()/3600)\n",
    "                D2=(Target[\"V\"][int((c2+c1)/2)])\n",
    "                if (TimeZone=='UTC') | (TimeZone=='Europe/London'):\n",
    "                    D3=(Target[\"datetime\"][c1])\n",
    "                    D4=(Target[\"datetime\"][c2])\n",
    "                else:\n",
    "                    D3=(Target[\"datetime\"][c1].tz_localize(TimeZone))\n",
    "                    D4=(Target[\"datetime\"][c2].tz_localize(TimeZone))\n",
    "                D5=([Target[\"RA\"][c1],Target[\"DEC\"][c1]])\n",
    "                D6=([Target[\"RA\"][c2],Target[\"DEC\"][c2]])\n",
    "                D7=(Target[\"alpha_true\"][int((c2+c1)/2)])\n",
    "                D8=(Target[\"PABLon\"][int((c2+c1)/2)])\n",
    "                D9=(Target[\"PABLat\"][int((c2+c1)/2)])\n",
    "                DurationT.append([D0,D1,D2,D3,D4,D5,D6,D7,D8,D9])\n",
    "            else:\n",
    "                c1=0\n",
    "                c2=0\n",
    "                check=0\n",
    "                for i in range(len(Target[\"datetime\"])-1):\n",
    "                    if  abs(Target[\"datetime\"][i]-Target[\"datetime\"][i+1])>dt.timedelta(hours=step):\n",
    "                        check=1\n",
    "                        c2=i\n",
    "                        D0=(pd.Timestamp.date(Target[\"datetime\"][c2]))\n",
    "                        dd=Target[\"datetime\"][c2]-Target[\"datetime\"][c1]\n",
    "                        D1=(dd.total_seconds()/3600)\n",
    "                        D2=(Target[\"V\"][int((c2+c1)/2)])\n",
    "                        if (TimeZone=='UTC') | (TimeZone=='Europe/London'):\n",
    "                            D3=(Target[\"datetime\"][c1])\n",
    "                            D4=(Target[\"datetime\"][c2])\n",
    "                        else:\n",
    "                            D3=(Target[\"datetime\"][c1].tz_localize(TimeZone))\n",
    "                            D4=(Target[\"datetime\"][c2].tz_localize(TimeZone))\n",
    "                        D5=([Target[\"RA\"][c1],Target[\"DEC\"][c1]])\n",
    "                        D6=([Target[\"RA\"][c2],Target[\"DEC\"][c2]])\n",
    "                        D7=(Target[\"alpha_true\"][int((c2+c1)/2)])\n",
    "                        D8=(Target[\"PABLon\"][int((c2+c1)/2)])\n",
    "                        D9=(Target[\"PABLat\"][int((c2+c1)/2)])\n",
    "                        DurationT.append([D0,D1,D2,D3,D4,D5,D6,D7,D8,D9])\n",
    "                        c1=c2+1\n",
    "                    \n",
    "            if (check==1):\n",
    "                colD=['Mid-Date', 'Duration (h)', 'V (mag)', 'Observation Start (UTC + Timezone)', 'Observation End (UTC + Timezone)', 'Start Point [RA,DEC]', 'End Point [RA,DEC]', 'Solar Phase (\\u00b0)','PABLon (\\u00b0)','PABLat (\\u00b0)']\n",
    "                DurationT=pd.DataFrame(DurationT, columns=colD)\n",
    "                if(max(DurationT[\"Duration (h)\"])>=MinObsTime) & (min(DurationT[\"V (mag)\"])<=APmag_min):\n",
    "                    TargetList.append([TargetName, max(DurationT[\"Duration (h)\"]), min(DurationT[\"V (mag)\"]), Family[\"P (h)\"][it], Family[\"lambda 1 (deg)\"][it], Family[\"beta 1 (deg)\"][it], Family[\"lambda 2 (deg)\"][it], Family[\"beta 2 (deg)\"][it]])\n",
    "                    IDtargetList.append(IDtarget)\n",
    "        \n",
    "            it=it+1\n",
    "                    \n",
    "        \n",
    "        colTL=['Target', 'Max. Duration (h)', 'Min. V (mag)', 'P (h)', '$\\lambda_1$ (\\u00b0)', '$\\\\beta_1$ (\\u00b0)', '$\\lambda_2$ (\\u00b0)', '$\\\\beta_2$ (\\u00b0)' ]\n",
    "        #colTL=['Target', 'Maximum Duration (h)', 'Minimum V (mag)', 'P (h)', \n",
    "        #       '<p>&lambda;<sub>1</sub> (deg)</p>', '<p>&beta;<sub>1</sub></p> (deg)', '&lambda;<sub>2</sub> (deg)', '&beta;<sub>2</sub> (deg)' ]\n",
    "        TargetList=pd.DataFrame(TargetList, columns=colTL)\n",
    "        #OptionTargets=list(zip(TargetList['Target'].tolist(),IDtargetList))\n",
    "        \n",
    "        #PRIORITY_RATE\n",
    "        ApAll=Family[\"a_p (au)\"]\n",
    "        KnownFamily=Family[Family[\"lambda 1 (deg)\"]!='-'] #Creating a new dataframe with the family members with known spin\n",
    "        ApK=KnownFamily[\"a_p (au)\"]\n",
    "        \n",
    "        #Density distribution of all family members\n",
    "        Ap = np.linspace(min(ApAll), max(ApAll), 1000)\n",
    "        densityA = sum(norm.pdf(Ap, xi, 0.02) for xi in ApAll)\n",
    "        \n",
    "        #Density distribution of family members with known spin \n",
    "        densityK = sum(norm.pdf(Ap,xi, 0.02) for xi in ApK)\n",
    "        \n",
    "        #Definition of priority rate for an homogenous observing sample\n",
    "        PR=np.floor(100*(max(densityK/densityA)-densityK/densityA))\n",
    "        TL=[]\n",
    "        for ta in range(len(IDtargetList)):\n",
    "            cin=Family[Family[\"Number\"]==IDtargetList[ta]].index\n",
    "            xa=take_closest(Ap, Family[\"a_p (au)\"][cin].tolist())\n",
    "            if Family[\"lambda 1 (deg)\"][cin].tolist()==['-']:\n",
    "                TL.append(PR[xa])\n",
    "            else:\n",
    "                TL.append(0)\n",
    "        \n",
    "        TargetList[\"Priority Rate (%)\"]=TL\n",
    "        \n",
    "        TargetList=TargetList.sort_values([\"Priority Rate (%)\", \"Max. Duration (h)\", \"Min. V (mag)\"], ascending = (False, False, True))\n",
    "        #pd.set_option('display.max_columns', None)  \n",
    "        #print('The targets are display here with a priority row:')\n",
    "        #TargetListBest.head(len(TargetList))\n",
    "        for row in TargetList.index:\n",
    "            IDtargetListPR.append(IDtargetList[row])\n",
    "        \n",
    "        OptionTargets=list(zip(TargetList['Target'].tolist(), IDtargetListPR))\n",
    "        \n",
    "        # Turn interactive plotting off\n",
    "        plt.ioff()\n",
    "        \n",
    "        CenterOfFamily=2.366\n",
    "        fig1 = plt.figure(figsize=(15,3))\n",
    "        \n",
    "        #Make room for x label\n",
    "        #plt.tight_layout()\n",
    "        plt.gcf().subplots_adjust(bottom=0.25)\n",
    "        \n",
    "        plt.axvline(0, color='grey', linestyle='--')\n",
    "        plt.text(-0.01,2*max(100*densityK/densityA)/3,\"$a_c=\"+str(CenterOfFamily)+\"\\ au$\", bbox=dict(facecolor='white', alpha=0.95))\n",
    "        plt.fill_between(Ap-CenterOfFamily, 100*(max(densityK/densityA)-densityK/densityA), alpha=0.5)\n",
    "        plt.plot(ApAll-CenterOfFamily, np.full_like(ApAll, -1), '+k', markersize=10 , markeredgewidth=1, label='All members')\n",
    "        plt.plot(ApK-CenterOfFamily, np.full_like(ApK, -1), '+r', markersize=10 , markeredgewidth=3, label='Members with known spin')\n",
    "        plt.title(\"Priority Rate for the Primordial Family members for an homogenous observing sample\")\n",
    "        plt.ylabel(\"Priority Rate (%)\")\n",
    "        plt.xlabel(\"Position from the centre of the family (au)\")\n",
    "        plt.legend()\n",
    "        plt.axis([min(ApAll-CenterOfFamily)-0.01, max(ApAll-CenterOfFamily)+0.01, -2, max(100*densityK/densityA)])\n",
    "        plt.close(fig1)\n",
    "        \n",
    "        #Export plot as HTML image\n",
    "        buf1 = io.BytesIO()\n",
    "        fig1.savefig(buf1, format='png')\n",
    "        buf1.seek(0)\n",
    "        string1 = base64.b64encode(buf1.read())\n",
    "        \n",
    "        uri1 = 'data:image/png;base64,' + urllib.parse.quote(string1)\n",
    "        FIG1.value='<img src = \"%s\"/>' % uri1\n",
    "        \n",
    "        \n",
    "        #DASHBOARD_1F (without progressbar, include taable and graph)\n",
    "        clear_output(wait=True)\n",
    "        progress.layout.visibility = 'hidden'\n",
    "        #DashBoard1f=ipw.VBox([Title, Copyright, Contact, InsertTab, FindAsteroids])\n",
    "        #display(DashBoard1)\n",
    "        #display(HTML(TargetList.to_html(index=False)))\n",
    "        ObservableAsteroids_text.layout.display = ''\n",
    "        TableList.value=TargetList.to_html(index=False, classes='table table-striped', justify='justify-all')\n",
    "        #display(HTML('<p><span style=\"font-size: 10px; font-family: Helvetica;\">The ephimeris data were computed using the <a href=\"https://ssd.jpl.nasa.gov/?horizons\" rel=\"noopener noreferrer\" target=\"_blank\">JPL&apos;s HORIZON system</a>, physical parametrers have been retrieved from the <a href=\"https://mp3c.oca.eu/catalogue/index.htm\" rel=\"noopener noreferrer\" target=\"_blank\">Minor Planet Physical Properties Catalogue</a> (Delbo, Walsh, Bolin, Avdellidou, Morbidelli. 2017. &quot;Identification of a primordial asteroid family constrains the original planetesimal population&quot;. Science. 357, 1026 - 1029) and the spin pole data from the <a href=\"https://astro.troja.mff.cuni.cz/projects/damit/\" rel=\"noopener noreferrer\" target=\"_blank\">DAMIT</a> (Ďurech, J., Sidorin, V., &amp; Kaasalainen, M. (2010). DAMIT: a database of asteroid models. Astronomy &amp; Astrophysics, 513, A46).</span></p>'))\n",
    "        #display(TableFootNote)\n",
    "        TableFootNote.layout.display = ''\n",
    "        #display(FIG1)\n",
    "        \n",
    "\n",
    "        #DASHBOARD_2\n",
    "        IDtarget_dd.options=OptionTargets #IDtargetList\n",
    "        DashBoard2.layout.display = ''\n",
    "        #display(DashBoard2)\n",
    "\n",
    "\n",
    "Find_btn=ipw.Button(description=\"Find the Observable Asteroids\", \n",
    "                    layout=ipw.Layout(justify_content='center', width='50%', height='50px'))\n",
    "Find_btn.style.button_color = 'mediumpurple'\n",
    "Find_btn.style.font_weight='bold'\n",
    "\n",
    "Find_btn.on_click(find_asteroids)\n",
    "\n",
    "OutputMessage4=ipw.HTML(\"\")\n",
    "\n",
    "FindAsteroids=ipw.HBox([Find_btn, OutputMessage4])\n",
    "FindAsteroids.layout.align_items = 'center'\n",
    "\n",
    "\n",
    "ObservableAsteroids_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 22px; font-family: Helvetica;\">The observable asteroids</span></strong></p>')\n",
    "ObservableAsteroids_text.layout.display = 'None'\n",
    "\n",
    "TableList=ipw.HTMLMath(\"\")\n",
    "TableFootNote=ipw.HTML('<p><span style=\"font-size: 12px; font-family: Helvetica;\">The ephimeris data were computed using the <a href=\"https://ssd.jpl.nasa.gov/?horizons\" rel=\"noopener noreferrer\" target=\"_blank\">JPL&apos;s HORIZON system</a>, physical parametrers have been retrieved from the <a href=\"https://mp3c.oca.eu/catalogue/index.htm\" rel=\"noopener noreferrer\" target=\"_blank\">Minor Planet Physical Properties Catalogue</a> (Delbo, Walsh, Bolin, Avdellidou, Morbidelli. 2017. &quot;Identification of a primordial asteroid family constrains the original planetesimal population&quot;. Science. 357, 1026 - 1029) and the spin pole data from the <a href=\"https://astro.troja.mff.cuni.cz/projects/damit/\" rel=\"noopener noreferrer\" target=\"_blank\">DAMIT</a> (Ďurech, J., Sidorin, V., &amp; Kaasalainen, M. (2010). DAMIT: a database of asteroid models. Astronomy &amp; Astrophysics, 513, A46).</span></p>')\n",
    "TableFootNote.layout.display = 'None'\n",
    "\n",
    "FIG1=ipw.HTML(\"\")\n",
    "\n",
    "#DASHBOARD_1\n",
    "progress.layout.visibility = 'hidden'\n",
    "DashBoard1=ipw.VBox([Title,  Copyright, Contact, InsertTab, FindAsteroids, progress, ObservableAsteroids_text, TableList, TableFootNote, FIG1])\n",
    "display(DashBoard1)\n",
    "#OUT1=ipw.Output()\n",
    "\n",
    "################################################ B PART ####################################################\n",
    "VOID=ipw.HTML(\"\")\n",
    "VOID.layout.height = '20px'\n",
    "\n",
    "#Please type the ID of the target\n",
    "TargetSelection_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 22px; font-family: Helvetica;\">Select your target</span></strong></p>')\n",
    "IDtarget_dd=ipw.Dropdown(description=\"ID Target:\", options=IDtargetList, style={'description_width': '100px'})\n",
    "FIG2=ipw.HTML(\"\")\n",
    "@output.capture(clear_output=False,wait=True)\n",
    "def targetSelection_input(b=None):\n",
    "    time.sleep(0.1)\n",
    "    DashBoard3.layout.display = 'None'\n",
    "    \n",
    "    Message1.value=\"\"\n",
    "    Message2.value=\"\"\n",
    "    TableOP.value=\"\"\n",
    "    FIG3.value=\"\"\n",
    "    \n",
    "    global IDtarget, ct1\n",
    "    ct1=0 #This is for the final button\n",
    "    \n",
    "    IDtarget=IDtarget_dd.value\n",
    "    cin=Family[Family[\"Number\"]==IDtarget].index\n",
    "    #clear_output(wait=True)\n",
    "    #plt.close(fig)\n",
    "    #display(TargetSelection)\n",
    "    #Show where is located\n",
    "    \n",
    "    #plt.close()\n",
    "    \n",
    "    #clear_output(True)\n",
    "    \n",
    "    # Turn interactive plotting off\n",
    "    plt.ioff()\n",
    "    cm = plt.cm.get_cmap('jet')\n",
    "    fig2,ax2 = plt.subplots(1,1, figsize=(10,5))\n",
    "\n",
    "    #ax2.clear()\n",
    "    plt.scatter(Family[\"a_p (au)\"], Family[\"1/D (km^-1)\"], color='k',s=10, label=\"Primordial Family members\")\n",
    "    sc1 = plt.scatter(KnownFamily[\"a_p (au)\"], KnownFamily[\"1/D (km^-1)\"], c=KnownFamily[\"beta 1 (deg)\"], edgecolor='black', vmin=-90, vmax=90, s=50, cmap=cm, label=\"Members with known spin\")\n",
    "    plt.plot([0,2.366], [1.39594, 0], '-.k')\n",
    "    plt.text(2.22,0.05,\"Primordial Family \\n (left border)\", ha='center', va='center')\n",
    "    plt.plot([0,2.42], [2.662, 0], '--k')\n",
    "    plt.plot([2.42, 3], [0, 3*1.1-2.662], '--k')\n",
    "    plt.text(2.42,0.125,\"Polana Family \\n borders\", ha='center', va='center')\n",
    "    plt.scatter(Family[\"a_p (au)\"][cin], Family[\"1/D (km^-1)\"][cin], s=100, marker='x', color='deeppink', label=\"Selected Target\")\n",
    "    \n",
    "    plt.xlabel(\"Proper semimajor axis (au)\")\n",
    "    plt.ylabel(\"$1/D \\ (km^{-1})$\")\n",
    "    plt.xlim([2.12,2.5])\n",
    "    plt.ylim([0, 0.28])\n",
    "    #plt.yticks(np.arange(0, 0.200001, step=0.05))\n",
    "    v = np.linspace(-90, 90, 10, endpoint=True)\n",
    "    \n",
    "    #if ct==0:\n",
    "    cb=plt.colorbar(sc1, ticks=v)\n",
    "    cb.set_label(\"$ \\\\beta \\ (^o)$\")\n",
    "    \n",
    "    plt.legend(loc='upper right', framealpha=0.95)\n",
    "    plt.close(fig2)\n",
    "    \n",
    "    #Export plot as HTML image\n",
    "    buf2 = io.BytesIO()\n",
    "    fig2.savefig(buf2, format='png')\n",
    "    buf2.seek(0)\n",
    "    string2 = base64.b64encode(buf2.read())\n",
    "        \n",
    "    uri2 = 'data:image/png;base64,' + urllib.parse.quote(string2)\n",
    "    #time.sleep(1)\n",
    "    FIG2.value = '<img src = \"%s\"/>' % uri2  \n",
    "    \n",
    "    clear_output(wait=True)\n",
    "    \n",
    "    #Previous Output\n",
    "    #display(DashBoard1)\n",
    "    #display(HTML(TargetList.to_html(index=False)))\n",
    "    #display(TableFootNote)\n",
    "    #display(FIG1)\n",
    "    #display(DashBoard2)\n",
    "    \n",
    "    #display(FIG2)\n",
    "    \n",
    "    #DASHBOARD_3\n",
    "    #if ct==0:\n",
    "    DashBoard3.layout.display = ''\n",
    "    #display(DashBoard3)\n",
    "    \n",
    "    #ct=ct+1\n",
    "\n",
    "\n",
    "#IDtarget_dd=ipw.interact(description=\"ID Target:\", options=IDtargetList, style={'description_width': '100px'})\n",
    "\n",
    "Insert_btn4=ipw.Button(description=\"Insert\", button_style='Success')\n",
    "#out = ipw.Output()\n",
    "#out\n",
    "\n",
    "TargetSelection=ipw.VBox([TargetSelection_text, ipw.HBox([IDtarget_dd, Insert_btn4])])\n",
    "#clear_output(wait=True)\n",
    "#clear_output()\n",
    "#display(TargetSelection)\n",
    "\n",
    "#Click Times\n",
    "ct=0;\n",
    "\n",
    "Insert_btn4.on_click(targetSelection_input)\n",
    "\n",
    "#DASHBOARD_2\n",
    "DashBoard2=ipw.VBox([VOID, TargetSelection, FIG2])\n",
    "DashBoard2.layout.display = 'None'\n",
    "display(DashBoard2)\n",
    "\n",
    "\n",
    "################################################ C PART ####################################################\n",
    "\n",
    "ObservationPlan_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 22px; font-family: Helvetica;\">The Observation Plan</span></strong></p>')\n",
    "\n",
    "\n",
    "TimeStep=ipw.FloatText(description=\"Time step:\", value=10, layout=Layout(width='15%'))\n",
    "TimeUnit=ipw.Dropdown(options=[\"m\", \"h\"], value=\"m\", layout=Layout(width='5%'))\n",
    "\n",
    "TimeStepUnit=ipw.HBox([TimeStep, TimeUnit])\n",
    "TableOP=ipw.HTML(\"\")\n",
    "FIG3=ipw.HTML(\"\")\n",
    "\n",
    "ct1=0\n",
    "@output.capture(clear_output=False,wait=True)\n",
    "def create_obsplan(b=None):\n",
    "    \n",
    "    time.sleep(0.1) #Refresh\n",
    "    OutputMessage5.value = '<p><strong><span style=\"color: blue;\">Please wait...</span></strong></p>'\n",
    "    \n",
    "    global FIG3, NightTime, LocalDate, airmass\n",
    "    \n",
    "    Message1.value=\"\"\n",
    "    Message2.value=\"\"\n",
    "    TableOP.value=\"\"\n",
    "    FIG3.value=\"\"\n",
    "    \n",
    "    step=TimeStep.value\n",
    "    units=TimeUnit.value\n",
    "    \n",
    "    ObsDays={'start':START+' 12:00', 'stop':STOP+' 12:00','step':(str(int(step))+units)}\n",
    "    obj = Horizons(id=IDtarget, location=Observatory, epochs=ObsDays)\n",
    "    eph = obj.ephemerides()\n",
    "    TargetName=eph['targetname'][0]\n",
    "    \n",
    "    if units=='h':\n",
    "        step=step\n",
    "    elif units=='m':\n",
    "        step=step/60\n",
    "    #elif units=='d':\n",
    "    #    step=step*24\n",
    "    else:\n",
    "        print('WARNING!: Please insert again a time-step with units in minutes(\"m\") or hours (\"h\")!')\n",
    "    \n",
    "    eph=eph.to_pandas()\n",
    "    airmass=eph['airmass']\n",
    "    Date=eph['datetime_str']\n",
    "    DateT=pd.to_datetime(Date)\n",
    "    \n",
    "    if (TimeZone=='UTC') | (TimeZone=='Europe/London'):\n",
    "        LocalDate=DateT\n",
    "    else:\n",
    "        #LocalDate = DateT.dt.tz_localize('UTC').dt.tz_convert(TimeZone)\n",
    "        LocalDate = DateT\n",
    "    \n",
    "    Sun=eph['solar_presence']\n",
    "    NightTime=[]\n",
    "    n=0\n",
    "    for x in range(len(Sun)):\n",
    "        if (Sun[x]==\"\"):\n",
    "            #print(LocalDate[x])\n",
    "            NightTime.append(LocalDate[x])\n",
    "            n=+1\n",
    "    \n",
    "    #The target at nightime\n",
    "    Target=eph[['datetime_str', 'solar_presence','V', 'airmass','RA','DEC','alpha_true','PABLon','PABLat']].copy(deep=True) #Suppress the SettingWithCopyWarning\n",
    "    Target['datetime_str']=pd.to_datetime(Target['datetime_str'])\n",
    "    pd.options.mode.chained_assignment = None  # default='warn'\n",
    "    Target.rename(columns={'datetime_str':'datetime'}, inplace=True)\n",
    "    Target=Target[(Target['airmass']<=3) & (Target['solar_presence']==\"\")]\n",
    "    Target.reset_index(drop=True, inplace=True)\n",
    "    \n",
    "    # Turn interactive plotting off\n",
    "    plt.ioff()\n",
    "    \n",
    "    #if ct1==0:\n",
    "    \n",
    "    fig3, ax3 = plt.subplots(2,1, figsize=(18,5), sharex=True)\n",
    "    #ax3[i].clear()\n",
    "    ax3[0].plot(LocalDate,airmass)\n",
    "    #ax3[i].fill_between(NightTime, [1]*len(NightTime), 0, facecolor='green', alpha=0.5)\n",
    "    #plt.fill_between(NightTime,[1]*len(NightTime))  \n",
    "    \n",
    "    #ax3[0].set_xlabel(\"Epoch (UTC)\")\n",
    "    ax3[0].set_ylabel(\"Airmass\")\n",
    "    ax3[0].set_title(\"Observing Window for \" + TargetName + \"\\n\" +\"Location: (\" + str(LAT) + \",\" + str(LON) + \")\")\n",
    "    ax3[0].set_xlim([LocalDate[0],LocalDate[len(LocalDate)-1]])\n",
    "    ax3[0].set_ylim([1, 3])\n",
    "    ax3[0].invert_yaxis()\n",
    "    \n",
    "    ax3[1].plot(LocalDate,eph[\"V\"]) \n",
    "    \n",
    "    ax3[1].set_xlabel(\"Epoch (UTC)\")\n",
    "    ax3[1].set_ylabel(\"Apparent Magnitude (mag)\")\n",
    "    ax3[1].set_title(\"Apparent Magnitude for \" + TargetName)\n",
    "    ax3[1].set_xlim([LocalDate[0],LocalDate[len(LocalDate)-1]])\n",
    "    ax3[1].invert_yaxis()\n",
    "    \n",
    "    for i, ax in enumerate(ax3.flatten()):\n", 
    "        for j in range(len(NightTime)-1):  #highlight from j to j+1\n",
    "            if abs(NightTime[j]-NightTime[j+1])<=dt.timedelta(hours=step):\n",
    "                pn3=ax.axvspan(NightTime[j], NightTime[j+1], facecolor='green', edgecolor='none', alpha=0.5)\n",
    "        pn3.set_label('Nighttime')\n",
    "        ax.legend(loc='upper right') \n",
    "        \n",
    "        if (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=90):\n",
    "            ml = MultipleLocator(1)\n",
    "            ax3[i].xaxis.set_major_locator(ml)\n",
    "            ax3[i].grid(which='major', axis='x', linestyle='-', linewidth=1.5,  color='#000000')\n",
    "            ax3[i].xaxis.set_major_locator(mpl.dates.DayLocator())\n",
    "            ax3[i].xaxis.set_major_formatter(mpl.dates.DateFormatter('%d %b %Y'))\n",
    "            ax3[i].xaxis.set_tick_params(which='major', pad=40)\n",
    "        \n",
    "        if (LocalDate[len(LocalDate)-1]-LocalDate[0])>=dt.timedelta(days=20):\n",
    "            plt.gcf().autofmt_xdate()\n",
    "        \n",
    "        if (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=1):\n",
    "            ml = MultipleLocator(1/24)\n",
    "            ax3[i].xaxis.set_minor_locator(ml)\n",
    "            ax3[i].grid(which='minor', axis='x', linestyle='-', linewidth=0.5)\n",
    "            ax3[i].xaxis.set_minor_formatter(mpl.dates.DateFormatter('%H:%M')) \n",
    "            plt.setp(ax3[i].xaxis.get_minorticklabels(), rotation=90)\n",
    "        elif (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=5):\n",
    "            ml = MultipleLocator(1/12)\n",
    "            ax3[i].xaxis.set_minor_locator(ml)\n",
    "            ax3[i].grid(which='minor', axis='x', linestyle='-', linewidth=0.5)\n",
    "            ax3[i].xaxis.set_minor_formatter(mpl.dates.DateFormatter('%H:%M')) \n",
    "            plt.setp(ax3[i].xaxis.get_minorticklabels(), rotation=90)\n",
    "        elif (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=10):\n",
    "            ml = MultipleLocator(1/8)\n",
    "            ax3[i].xaxis.set_minor_locator(ml)\n",
    "            ax3[i].grid(which='minor', axis='x', linestyle='-', linewidth=0.5)\n",
    "            ax3[i].xaxis.set_minor_formatter(mpl.dates.DateFormatter('%H:%M')) \n",
    "            plt.setp(ax3[i].xaxis.get_minorticklabels(), rotation=90)\n",
    "        elif (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=15):\n",
    "            ml = MultipleLocator(1/6)\n",
    "            ax3[i].xaxis.set_minor_locator(ml)\n",
    "            ax3[i].grid(which='minor', axis='x', linestyle='-', linewidth=0.5)\n",
    "            ax3[i].xaxis.set_minor_formatter(mpl.dates.DateFormatter('%H:%M')) \n",
    "            plt.setp(ax3[i].xaxis.get_minorticklabels(), rotation=90)\n",
    "        elif (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=20):\n",
    "            ml = MultipleLocator(1/3)\n",
    "            ax3[i].xaxis.set_minor_locator(ml)\n",
    "            ax3[i].grid(which='minor', axis='x', linestyle='-', linewidth=0.5)\n",
    "            ax3[i].xaxis.set_minor_formatter(mpl.dates.DateFormatter('%H:%M')) \n",
    "            plt.setp(ax3[i].xaxis.get_minorticklabels(), rotation=90)\n",
    "        elif (LocalDate[len(LocalDate)-1]-LocalDate[0])<=dt.timedelta(days=30):\n",
    "            ml = MultipleLocator(1/2)\n",
    "            ax3[i].xaxis.set_minor_locator(ml)\n",
    "            ax3[i].grid(which='minor', axis='x', linestyle='-', linewidth=0.5)\n",
    "            ax3[i].xaxis.set_minor_formatter(mpl.dates.DateFormatter('%H:%M')) \n",
    "            plt.setp(ax3[i].xaxis.get_minorticklabels(), rotation=90)\n",
    "    \n",
    "    plt.tight_layout()\n",
    "    #plt.draw()\n",
    "    plt.close(fig3)\n",
    "        \n",
    "    #Export plot as HTML image\n",
    "    buf3 = io.BytesIO()\n",
    "    fig3.savefig(buf3, format='png')\n",
    "    buf3.seek(0)\n",
    "    string3 = base64.b64encode(buf3.read())\n",
    "        \n",
    "    uri3 = 'data:image/png;base64,' + urllib.parse.quote(string3)\n",
    "    FIG3.value='<img src = \"%s\"/>' % uri3\n",
    "    \n",
    "    # Turn interactive plotting off\n",
    "    #plt.ioff()\n",
    "    \n",
    "    #if ct1==0:\n",
    "    #fig4, ax4 = plt.subplots(1,1, figsize=(20,4))\n",
    "    #ax4.clear()    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    #CREATE THE OBSERVATION PLAN TABLE\n",
    "    DurationT=[] \n",
    "    if abs(Target[\"datetime\"][0]-Target[\"datetime\"][len(Target[\"datetime\"])-1])<dt.timedelta(hours=24):\n",
    "        #check=1\n",
    "        c1=0\n",
    "        c2=len(Target[\"datetime\"])-1\n",
    "        if (pd.Timestamp.date(Target[\"datetime\"][c1])==pd.Timestamp.date(Target[\"datetime\"][c2])):\n",
    "            D0=(pd.Timestamp.date(Target[\"datetime\"][c2])+dt.timedelta(days=1))\n",
    "        else:\n",
    "            D0=(pd.Timestamp.date(Target[\"datetime\"][c2]))\n",
    "        dd=Target[\"datetime\"][c2]-Target[\"datetime\"][c1]\n",
    "        D1=(dd.total_seconds()/3600)\n",
    "        D2=(Target[\"V\"][int((c2+c1)/2)])\n",
    "        if (TimeZone=='UTC') | (TimeZone=='Europe/London'):\n",
    "            D3=(Target[\"datetime\"][c1])\n",
    "            D4=(Target[\"datetime\"][c2])\n",
    "        else:\n",
    "            D3=(Target[\"datetime\"][c1].tz_localize(TimeZone))\n",
    "            D4=(Target[\"datetime\"][c2].tz_localize(TimeZone))\n",
    "        D5=([Target[\"RA\"][c1],Target[\"DEC\"][c1]])\n",
    "        D6=([Target[\"RA\"][c2],Target[\"DEC\"][c2]])\n",
    "        D7=(Target[\"alpha_true\"][int((c2+c1)/2)])\n",
    "        D8=(Target[\"PABLon\"][int((c2+c1)/2)])\n",
    "        D9=(Target[\"PABLat\"][int((c2+c1)/2)])\n",
    "        DurationT.append([D0,D1,D2,D3,D4,D5,D6,D7,D8,D9])\n",
    "    else:\n",
    "        c1=0\n",
    "        c2=0\n",
    "        check=0\n",
    "        for i in range(len(Target[\"datetime\"])-1):\n",
    "            if  abs(Target[\"datetime\"][i]-Target[\"datetime\"][i+1])>dt.timedelta(hours=step) or i==len(Target[\"datetime\"])-2:\n",
    "                #check=1\n",
    "                c2=i\n",
    "                if (pd.Timestamp.date(Target[\"datetime\"][c1])==pd.Timestamp.date(Target[\"datetime\"][c2])):\n",
    "                    D0=(pd.Timestamp.date(Target[\"datetime\"][c2])+dt.timedelta(days=1))\n",
    "                else:\n",
    "                    D0=(pd.Timestamp.date(Target[\"datetime\"][c2]))\n",
    "                dd=Target[\"datetime\"][c2]-Target[\"datetime\"][c1]\n",
    "                D1=(dd.total_seconds()/3600)\n",
    "                D2=(Target[\"V\"][int((c2+c1)/2)])\n",
    "                if (TimeZone=='UTC') | (TimeZone=='Europe/London'):\n",
    "                    D3=(Target[\"datetime\"][c1])\n",
    "                    D4=(Target[\"datetime\"][c2])\n",
    "                else:\n",
    "                    D3=(Target[\"datetime\"][c1].tz_localize(TimeZone))\n",
    "                    D4=(Target[\"datetime\"][c2].tz_localize(TimeZone))\n",
    "                D5=([Target[\"RA\"][c1],Target[\"DEC\"][c1]])\n",
    "                D6=([Target[\"RA\"][c2],Target[\"DEC\"][c2]])\n",
    "                D7=(Target[\"alpha_true\"][int((c2+c1)/2)])\n",
    "                D8=(Target[\"PABLon\"][int((c2+c1)/2)])\n",
    "                D9=(Target[\"PABLat\"][int((c2+c1)/2)])\n",
    "                DurationT.append([D0,D1,D2,D3,D4,D5,D6,D7,D8,D9])\n",
    "                c1=c2+1\n",
    "    \n",
    "    col=['Mid-Date', 'Duration (h)', 'V (mag)', 'Observation Start (UTC + Timezone)', 'Observation End (UTC + Timezone)', 'Start Point [RA,DEC]', 'End Point [RA,DEC]', 'Solar Phase (\\u00b0)','PABLon (\\u00b0)','PABLat (\\u00b0)']\n",
    "    \n",
    "    DurationT=pd.DataFrame(DurationT, columns=col)\n",
    "    \n",
    "    clear_output(wait=True)\n",
    "    OutputMessage5.value=\"\"\n",
    "    \n",
    "    #Previous Output\n",
    "    #display(DashBoard1)\n",
    "    #display(HTML(TargetList.to_html(index=False)))\n",
    "    #display(TableFootNote)\n",
    "    #display(FIG1)\n",
    "    #display(DashBoard2)\n",
    "    #display(FIG2)\n",
    "    #display(DashBoard3)\n",
    "\n",
    "    \n",
    "    #DASHBOARD_4\n",
    "    #TimeStepUnit.layout.display = 'None' \n",
    "    #ObPlan_btn.layout.display = 'None'\n",
    "    Message1.value = '<p style=\"text-align: right;\"> </p>' + \"Location: (\" + str(LAT) + \",\" + str(LON) + \") \\t Timezone: \" +TimeZone\n",
    "    Message2.value = 'The target '+TargetName+ ' is observable for these dates:'\n",
    "    \n",
    "    #display(ipw.HTML(value=Message2))\n",
    "    #display(HTML(DurationT.to_html(index=False)))\n",
    "    TableOP.value=DurationT.to_html(index=False, classes='table table-striped', justify='justify-all')\n",
    "    #display(FIG3)\n",
    "    #display(FIG4)\n",
    "    \n",
    "    #display(VOID)\n",
    "    #AsteroidFinder_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 18px; font-family: Helvetica;\">Here you can find the path of the asteroid in the sky:</span></strong></p>')\n",
    "    #%%html\n",
    "    #<iframe src=\"https://asteroid.lowell.edu/astfinder/\" frameborder=\"0\" width=\"1800\" height=\"1000\"></iframe>\n",
    "    \n",
    "    #display(VOID)\n",
    "    #AsteroidFinder_text=ipw.HTML('<p style=\"text-align: left;\"><strong><span style=\"font-size: 18px; font-family: Helvetica;\">Here you can check archival lightcurves:</span></strong></p>')\n",
    "    #%%html\n",
    "    #<iframe src=\"http://alcdef.org/\" frameborder=\"0\" width=\"1800\" height=\"1000\"></iframe>\n",
    "    \n",
    "    \n",
    "    #ct1=ct1+1\n",
    "\n",
    "\n",
    "\n",
    "ObPlan_btn=ipw.Button(description=\"Create the Observation Plan\", \n",
    "                    layout=ipw.Layout(justify_content='center', width='50%', height='50px'))\n",
    "ObPlan_btn.style.button_color = 'mediumslateblue'\n",
    "ObPlan_btn.style.font_weight='bold'\n",
    "\n",
    "ObPlan_btn.on_click(create_obsplan)\n",
    "\n",
    "OutputMessage5=ipw.HTML(\"\")\n",
    "\n",
    "Message1= ipw.HTML(\"\")\n",
    "Message2= ipw.HTML(\"\")\n",
    "    \n",
    "FrontTableOP=ipw.VBox([Message1, Message2])\n",
    "\n",
    "ObPlan=ipw.VBox([ObPlan_btn, OutputMessage5, FrontTableOP, TableOP, FIG3])\n",
    "ObPlan.layout.align_items = 'center'\n",
    "\n",
    "\n",
    "#DASHBOARD_3\n",
    "DashBoard3=ipw.VBox([VOID,ObservationPlan_text, TimeStepUnit, ObPlan,])\n",
    "DashBoard3.layout.display = 'None'\n",
    "display(DashBoard3)\n",
    "output"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}