File indexing completed on 2024-04-28 15:11:24

 /*
     SPDX-FileCopyrightText: 2012 Samikshan Bairagya <samikshan@gmail.com>
 
     SPDX-License-Identifier: GPL-2.0-or-later
 */
 
 #include "obsconditions.h"
 
 #include <QDebug>
 
 #include <cmath>
 
 ObsConditions::ObsConditions(int bortle, double aperture, Equipment equip, TelescopeType telType)
     : m_BortleClass(bortle), m_Equip(equip), m_TelType(telType), m_Aperture(aperture)
 {
     // 't' parameter
     switch (m_TelType)
     {
         case Reflector:
             m_tParam = 0.7;
             break;
         case Refractor:
             m_tParam = 0.9;
             break;
         case Invalid:
             m_tParam = -1; //invalid value
             break;
     }
     setLimMagnitude();
 
     qDebug() << Q_FUNC_INFO << "Aperture value being used:" << m_Aperture;
 }
 
 QMap<int, double> ObsConditions::setLMMap()
 {
     QMap<int, double> LMMap;
     LMMap.insert(1, 7.8);
     LMMap.insert(2, 7.3);
     LMMap.insert(3, 6.8);
     LMMap.insert(4, 6.3);
     LMMap.insert(5, 5.8);
     LMMap.insert(6, 5.3);
     LMMap.insert(7, 4.8);
     LMMap.insert(8, 4.3);
     LMMap.insert(9, 3.8);
 
     return LMMap;
 }
 
 const QMap<int, double> ObsConditions::m_LMMap = setLMMap();
 
 void ObsConditions::setLimMagnitude()
 {
     m_LM = m_LMMap[m_BortleClass];
 }
 
 double ObsConditions::getOptimumMAG()
 {
     double power = (2.81 + 2.814 * m_LM - 0.3694 * pow(m_LM, 2)) / 5;
     return 0.1333 * m_Aperture * sqrt(m_tParam) * pow(power, 10);
 }
 
 double ObsConditions::getTrueMagLim()
 {
     //     qDebug()<< (4.12 + 2.5 * log10( pow(aperture,2)*t ));
     //     return 4.12 + 2.5 * log10( pow(aperture,2)*t ); //Taking optimum magnification into consideration
 
     ///If there is no equipment available then return limiting magnitude for naked-eye
     if (m_Equip == None || m_Aperture == -1)
         return m_LM;
 
     /**
      * This is a more traditional formula which does not take the
      * 't' parameter into account. It also does not take into account
      * the magnification being used. The formula used is:
      *
      * TLM_trad = LM + 5*log10(aperture/7.5)
      *
      * The calculation is just based on the calculation of the
      * telescope's aperture to eye's pupil surface ratio.
      */
 
     return m_LM + 5 * log10(m_Aperture / 7.5);
 }
 
 bool ObsConditions::isVisible(GeoLocation *geo, dms *lst, SkyObject *so)
 {
     if (so->type() == SkyObject::SATELLITE)
     {
         return so->alt().Degrees() > 6.0;
     }
     KStarsDateTime ut = geo->LTtoUT(KStarsDateTime(QDateTime::currentDateTime().toLocalTime()));
     SkyPoint sp       = so->recomputeCoords(ut, geo);
 
     //check altitude of object at this time.
     sp.EquatorialToHorizontal(lst, geo->lat());
 
     return (sp.alt().Degrees() > 6.0 && so->mag() < getTrueMagLim());
 }
 
 void ObsConditions::setObsConditions(int bortle, double aperture, ObsConditions::Equipment equip,
                                      ObsConditions::TelescopeType telType)
 {
     m_BortleClass = bortle;
     setLimMagnitude();
     m_Aperture = aperture;
     m_Equip    = equip;
     m_TelType  = telType;
 
     qDebug() << Q_FUNC_INFO << "Aperture value being used:" << m_Aperture;
 }