1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
|
/*******************************************************************************
Copyright(c) 2021 Jasem Mutlaq. All rights reserved.
AstroTrac Mount Driver.
This library is free software; you can redistribute it and/or
modify it under the terms of the GNU Library General Public
License version 2 as published by the Free Software Foundation.
.
This library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Library General Public License for more details.
.
You should have received a copy of the GNU Library General Public License
along with this library; see the file COPYING.LIB. If not, write to
the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
Boston, MA 02110-1301, USA.
*******************************************************************************/
#pragma once
#include "AlignmentSubsystemForMathPlugins.h"
#include "ConvexHull.h"
namespace INDI
{
namespace AlignmentSubsystem
{
/*!
* \struct AlignmentDatabaseEntry
* \brief Entry in the in memory alignment database
*
*/
struct ExtendedAlignmentDatabaseEntry : public AlignmentDatabaseEntry
{
/// \brief Default constructor
ExtendedAlignmentDatabaseEntry()
: AlignmentDatabaseEntry(),
CelestialAzimuth(0),
CelestialAltitude(0),
TelescopeAzimuth(0),
TelescopeAltitude(0) {}
/// \brief Copy constructor
ExtendedAlignmentDatabaseEntry(const ExtendedAlignmentDatabaseEntry &Source)
: AlignmentDatabaseEntry(Source),
CelestialAzimuth(Source.CelestialAzimuth),
CelestialAltitude(Source.CelestialAltitude),
TelescopeAzimuth(Source.TelescopeAzimuth),
TelescopeAltitude(Source.TelescopeAltitude)
{
}
/// Override the assignment operator to provide a const version
inline const ExtendedAlignmentDatabaseEntry &operator=(const ExtendedAlignmentDatabaseEntry &RHS)
{
ObservationJulianDate = RHS.ObservationJulianDate;
RightAscension = RHS.RightAscension;
Declination = RHS.Declination;
TelescopeDirection = RHS.TelescopeDirection;
PrivateDataSize = RHS.PrivateDataSize;
CelestialAzimuth = RHS.CelestialAzimuth;
CelestialAltitude = RHS.CelestialAltitude;
TelescopeAzimuth = RHS.TelescopeAzimuth;
TelescopeAltitude = RHS.TelescopeAltitude;
if (0 != PrivateDataSize)
{
PrivateData.reset(new unsigned char[PrivateDataSize]);
memcpy(PrivateData.get(), RHS.PrivateData.get(), PrivateDataSize);
}
return *this;
}
/**
* @brief Celestial Azimuth of Sync Point at the time it was added to the list.
*/
double CelestialAzimuth;
/**
* @brief Celestial Altitude of the Sync point at the time it was added to the list.
*/
double CelestialAltitude;
/**
* @brief Telescope Azimuth of Sync Point at the time it was added to the list.
*/
double TelescopeAzimuth;
/**
* @brief Telescope Altitude of the Sync point at the time it was added to the list.
*/
double TelescopeAltitude;
};
class NearestMathPlugin : public AlignmentSubsystemForMathPlugins
{
public:
NearestMathPlugin();
virtual ~NearestMathPlugin();
virtual bool Initialise(InMemoryDatabase *pInMemoryDatabase);
virtual bool TransformCelestialToTelescope(const double RightAscension, const double Declination,
double JulianOffset,
TelescopeDirectionVector &ApparentTelescopeDirectionVector);
virtual bool TransformTelescopeToCelestial(const TelescopeDirectionVector &ApparentTelescopeDirectionVector,
double &RightAscension, double &Declination);
private:
std::vector<ExtendedAlignmentDatabaseEntry> ExtendedAlignmentPoints;
/**
* @brief SphereUnitDistance Get distance between two points on a sphere.
* @param theta1 latitudal angle of object 1
* @param theta2 latitudal angle of object 2
* @param phi1 longitudal angle of object 1
* @param phi2 longitudal angle of object 2
* @return disatnce in degrees.
*/
double SphereUnitDistance(double theta1, double theta2, double phi1, double phi2);
/**
* @brief GetNearestPoint Traverses the ExtendedAlignmentPoints to find the closest point in horizontal coordinates on
* a sphere.
* @param Azimuth Object azimuth in degrees.
* @param Altitude Object altitude in degrees.
* @param isCelestial If true, compute difference between Celestial coords, otherwise compute using Telescope coords.
* @return Closest point in data set.
*/
ExtendedAlignmentDatabaseEntry GetNearestPoint(const double Azimuth, const double Altitude, bool isCelestial);
};
} // namespace AlignmentSubsystem
} // namespace INDI
|
