-
Notifications
You must be signed in to change notification settings - Fork 48
/
QReadAnt.h
144 lines (131 loc) · 7.69 KB
/
QReadAnt.h
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
139
140
141
142
143
144
/*************************************************************************
**
** MG-APPS----Multi-GNSS-Automatic Precise Positioning Software
** Copyright (C) 2016-2019 XiaoGongWei
** This file is part of MG-APPS.
**
** GNU Lesser General Public License Usage
** Alternatively, this file may be used under the terms of the GNU Lesser
** General Public License version 3 as published by the Free Software
** Foundation and appearing in the file LICENSE.LGPL3 included in the
** packaging of this file. Please review the following information to
** ensure the GNU Lesser General Public License version 3 requirements
** will be met: https://www.gnu.org/licenses/lgpl-3.0.html.
**
** MPL License Usage
** This Source Code Form is subject to the terms of the Mozilla Public
** License, v. 2.0. If a copy of the MPL was not distributed with this
** file, You can obtain one at http://mozilla.org/MPL/2.0/.
**
** GPLv3.0 License Usage
** This program is free software: you can redistribute it and/or modify
** it under the terms of the GNU General Public License as published by
** the Free Software Foundation, either version 3 of the License, or
** (at your option) any later version.
**
** This program 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 General Public License for more details.
**
** You should have received a copy of the GNU General Public License
** along with this program. If not, see http://www.gnu.org/licenses/.
**
**************************************************************************
** Author: XiaoGongWei
** Website/Contact: http://github.com/xiaogongwei
** Date: 26.04.2019
****************************************************************************/
#ifndef QREADANT_H
#define QREADANT_H
#include "QCmpGPST.h"
/* program description
If the "receiver" antenna type and Julian obvJD at the "start" moment of observation cannot be passed in during initialization, setObsJD function must be set
Antenna type and Julian day otherwise cannot match the data of satellite start effective time antenna
1. Class initialization
QString AntFileName = "D:\\antmod.atx", AntType = "ASH700936B_M SNOW";//AntFileName: antenna file path (empty will go to the current directory to open "antmod.atx") AntType:O file (observation file) to obtain the antenna type
double ObsJD = 2455296.50;// start time (approximately) Julian day
QReadAnt readAnt(AntFileName,AntType,ObsJD);
or
QReadAnt readAnt(AntFileName);
readAnt.setObsJD(AntType,ObsJD);
2. Read data files
Read data only once at initialization (time taken)
readAnt.getAllData();// read all data from satellites and receivers
3. Function call method
double E = 0.7,A = 0.3;// altitude Angle and azimuth Angle
double L1Offset = 0, L2Offset = 0;// save the correction of L1 and L2 line of sight (unit m)
readAnt.getRecvL12(E,A,&L1Offset,&L2Offset);// calculate the correction of L1 and L2 line of sight (unit m)
int Year = 2010,Month = 4,Day = 10,Hours = 0,Minuts = 0,PRN = 32;
double Seconds = 0.0;
double StaPos[3] = {9999,999,9999},RecPos[3] = {9999,9999,9999};
double L12OffSet = 0;
// time is UTC, month, year, day, hour and second, StaPos: satellite XYZ coordinate, RecPos: receiver XYZ coordinate, SatAntH: satellite antenna correction (since the PCO and PCV of each satellite antenna frequency are the same, it only needs to return a correction distance) to return m
gL12OffSet = readAnt.getSatOffSet(Year,Month,Day,Hours,Minuts,Seconds,PRN,StaPos,RecPos);
*/
// define the storage antenna data structure. Only commonly used data is saved here
typedef struct _FrqunceData
{// define a PCO and PCV data format for the frequency band
QString FrqFlag;//START OF FREQUENCY
double PCO[3];
//PCV data
short int Hang;// PCVAZI line number
short int Lie;// number of PCVAZI and PCVNoAZI columns
QVector< double > PCVNoAZI;
QVector< double > PCVAZI;// the matrix is converted to a 1-dimensional vector
}FrqData;
typedef struct _AntDataType
{
QString StrAntType;// antenna type
QString SatliCNN;//CNN satellite Numbers G,R,E...
double DAZI;// store azimuth interval 360/DAZI
double ZEN_12N[3];// store high Angle range and interval (ZEN_12N [1] - ZEN_12N [0])/ZEN_12N [2]
short NumFrqu;// number of frequency bands of satellites
double ValidJD;// valid start time
double EndJD;// valid end time
QVector <FrqData> PCOPCV;// PCO and PCV data for multiple frequency bands are stored
bool IsSat;//true represents satellite false represents receiver data
bool isFrqData;// if all antenna files are not found, mark false; // initialization hypothesis can be found that the antenna is true
}AntDataType;
class QReadAnt:public QBaseObject
{
// function part
public:
QReadAnt(QString AntFileName = "",QString AnTypeName = "",double ObvJD = 0);
void setAntFileName(QString AntFileName = "",QString AnTypeName = "",double ObvJD = 0);
~QReadAnt(void);
bool getRecvL12(double E, double A, char SatType, double &L1Offset, double &L2Offset, QVector<QString> FrqFlag);//E: satellite altitude Angle; A: satellite azimuth (unit radian), L1Offset and L2Offset wavelength correction (unit m)
bool getSatOffSet(int Year,int Month,int Day,int Hours,int Minuts,double Seconds,int PRN,char SatType,double *StaPos,double *RecPos, double *L12Offset,
QVector<QString> FrqFlag);// StaPos: satellite XYZ coordinate, RecPos: receiver XYZ coordinate,L12Offset:FrqFlag first two frequency value correction return m
void setObsJD(QString AnTypeName,double ObsJD);// a Julian date for observation must be set
bool getAllData();// read valid satellite antenna data and receiver antenna data
private:
void initVar();// initializes variables
bool openFiles(QString AntTypeName);// open the file
bool readFileAntHead();// read the header file
bool readAntFile(QString AntTypeName);// read the antenna file
bool readSatliData();// read the data needed by the satellite
bool readRecvData();// read the documents needed by the satellite
double computeJD(int Year,int Month,int Day,int HoursInt=0,int Minutes=0,int Seconds=0);// calculate Julian days
bool getPCV(const AntDataType &tempAntData,char SatType,double *PCV12,double Ztop,double AZI, QVector<QString> FrqFlag);// calculate the PCV; Z: fixed Angle of satellite = PI /2 - altitude Angle; A: satellite azimuth (unit radians)
void findPCO(double *pco, QVector<FrqData> &PCOPCV, char SatType = 'G', QString frqStr = "C1C");// add by xiaogongwei 2019.04.12
// variable section
public:
double m_sunpos[3],m_moonpos[3],m_gmst;// each epoch is updated once to provide coordinates for other programs, reducing repeated calculations
QCmpGPST m_CmpClass;// compute library public for compute sun and moon position
private:
QString m_AntFileName;// save the antenna data file name
QString m_AntType;// save the receiver antenna type name
QFile m_ReadFileClass;// read the file class
bool isReadAllData;// whether to read all data flags
bool isReadHead;// whether to read header file
bool isReadSatData;// whether to read satellite antenna data
bool isReadRecvData;// whether to read receiver antenna data
QString m_tempLine;// cache read row data
double m_ObeservTime;// Julian days to store observation O files to match the effective time of the satellite antenna
AntDataType m_RecvData;// data of receiver antenna PCO and PCV are saved
QVector< AntDataType > m_SatData;// data of multiple satellite antenna PCO and PCV are saved
double m_pi;// PI radian
double m_sunSecFlag;// if the reference epoch changes, the solar coordinates need to be recalculated to store the reference epoch in seconds
};
#endif