groops-devs / groops

A software toolkit for gravity field recovery and GNSS processing
GNU General Public License v3.0
185 stars 107 forks source link

Maybe there is something wrong in the gnssResiduals2Skyplot. #156

Closed lambert87041 closed 1 year ago

lambert87041 commented 1 year ago

Description

I have found something wrong about the function gnssResiduals2Skyplot. When I set the parameters ‘type’ to ‘G’, ‘typeTransmitter’ to ‘G05’, which means that I just want to draw only one satellite in skyplot. But the output image is like this:

skyplot_2022-01-01 WUH2

After I checked the code and added the content about PRN matching in the gnssResiduals2Skyplot, I found that the output image became normal.

skyplot_2022-01-01 WUH2 - 副本

It seems that there is still a problem with the function of using ‘typeTransmitter’ to match a satellite. Of course, I am not sure whether it is the problem of my configuration file. I am using the GROOPS compiled by Visual Studio 2017.

GROOPS version

main (latest commit)

Operating systems

Log output

The log seems normal
2023-02-25 15:34:10 Status  === Starting GROOPS ===
2023-02-25 15:34:10 Info    Config file: <D:\data\2021\scenarioGnssPPP\skyPlot.xml>
2023-02-25 15:34:10 Status  --- FileCreateDirectories (D:/data/2021/scenarioGnssPPP/skyplot) ---
2023-02-25 15:34:10 Status  Create directory <D:/data/2021/scenarioGnssPPP/skyplot>
2023-02-25 15:34:10 Status  --- LoopPrograms (daily) ---
2023-02-25 15:34:10 Status  === 1. loop ===
2023-02-25 15:34:10 Status  --- FileCreateDirectories (D:/data/2021/scenarioGnssPPP/tmp/2022-01-01) ---
2023-02-25 15:34:10 Status  Create directory <D:/data/2021/scenarioGnssPPP/tmp/2022-01-01>
2023-02-25 15:34:10 Status  --- LoopPrograms (positions) ---
2023-02-25 15:34:10 Status  === 1. loop ===
2023-02-25 15:34:10 Status  --- GnssResiduals2Skyplot (D:/data/2021/scenarioGnssPPP/tmp/2022-01-01/grid_2022-01-01.WUH2.gps.dat) ---
2023-02-25 15:34:10 Status  read GNSS residuals <D:/data/2021/scenarioGnssPPP/output/2022-01-01/residuals_2022-01-01.WUH2.dat>
2023-02-25 15:34:14 Status  save values to file <D:/data/2021/scenarioGnssPPP/tmp/2022-01-01/grid_2022-01-01.WUH2.gps.dat>
2023-02-25 15:34:15 Info    grid statistics
2023-02-25 15:34:15 Info      count: 23516
2023-02-25 15:34:15 Info      area:  187134% of Earth's surface (9.56646e+11 km^2)
2023-02-25 15:34:15 Info    data statistics
2023-02-25 15:34:15 Info      rms:  40.2402      
2023-02-25 15:34:15 Info      avg:  34.0066      
2023-02-25 15:34:15 Info      min:  3.65687      
2023-02-25 15:34:15 Info      max:  178.719      
2023-02-25 15:34:15 Info      mean: 34.0066      
2023-02-25 15:34:15 Status  --- PlotMap (D:/data/2021/scenarioGnssPPP/skyplot/skyplot_2022-01-01.WUH2.png) ---
2023-02-25 15:34:15 Status  create temporary data files
2023-02-25 15:34:15 Status  create scriptfile
2023-02-25 15:34:15 Status  run scriptfile
2023-02-25 15:34:16 Status  1 loops in time 00:00:06
2023-02-25 15:34:16 Status  1 loops in time 00:00:06
2023-02-25 15:34:16 Status  === Finished GROOPS ===
tmayerguerr commented 1 year ago

Hi,

From the documentation:

If with typeTransmitter (e.g. '***G05') a single transmitter is selected, the azimuth and elevation are computed from the transmitter point of view.

The plot is correct: it shows the motion of the WHU station as seen from the G05 satellite (azimuth and elevation of the station in the satellite body frame). If you want to have the skyplot as seen from WHU station for a single satellite, set the parameters ‘type’ instead of typeTransmitter to ‘***G05'.

regards Torsten