Next Article in Journal
Application of a Convolutional Neural Network for the Detection of Sea Ice Leads
Next Article in Special Issue
Performance Evaluation of Single-Frequency Precise Point Positioning and Its Use in the Android Smartphone
Previous Article in Journal
InSAR Monitoring of Arctic Landfast Sea Ice Deformation Using L-Band ALOS-2, C-Band Radarsat-2 and Sentinel-1
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preliminary Results on Tropospheric ZTD Estimation by Smartphone

1
Department of Informatics Bioengineering Robotics and Systems Engineering (DIBRIS), University of Genoa, via Dodecaneso 35, 16146 Genoa, Italy
2
Gter srl, Innovation in Geomatics GNSS and GIS, via Jacopo Ruffini 9/1A, 16123 Genoa, Italy
3
Department of Environment, Land, and Infrastructure Engineering (DIATI), Politecnico di Torino, corso Duca degli Abruzzi 24, 10129 Turin, Italy
4
Department of Civil, Chemical and Environmental Engineering (DICCA), University of Genoa, via Montallegro 1, 16145 Genoa, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(22), 4567; https://doi.org/10.3390/rs13224567
Submission received: 4 October 2021 / Revised: 9 November 2021 / Accepted: 11 November 2021 / Published: 13 November 2021

Abstract

:
The Global Navigation Satellite System (GNSS) receiver is one of the many sensors embedded in smartphones. The early versions of the Android operating system could only access limited information from the GNSS, allowing the related Application Program Interface (API) to obtain only the location. With the development of the Android 7.0 (Nougat) operating system in May 2016, raw measurements from the internal GNSS sensor installed in the smartphone could be accessed. This work aims to show an initial analysis regarding the feasibility of Zenith Total Delay (ZTD) estimation by GNSS measurements extracted from smartphones, evaluating the accuracy of estimation to open a new window on troposphere local monitoring. Two different test sites have been considered, and two different types of software for data processing have been used. ZTDs have been estimated from both a dual-frequency and a multi-constellation receiver embedded in the smartphone, and from a GNSS Continuously Operating Reference Station (CORS). The results have shown interesting performances in terms of ZTD estimation from the smartphone in respect of the estimations obtained with a geodetic receiver.

1. Introduction

The applications exploiting Global Navigation Satellite System (GNSS) positioning have increased in recent years [1,2]. The first GNSS-derived Zenith Total Delay (ZTD) estimates were carried out for scientific purposes using a post-processing strategy; the influence of ZTD was then considered in real-time applications to guarantee high-quality positioning, and it is now tackled in a Precise Point Positioning (PPP) strategy [3,4].
The development of networks of Continuously Operating Reference Stations (CORSs) for Network Real-Time Kinematic (NRTK) positioning and the increasing diffusion of PPP strategies allow one to obtain a centimetric accuracy in a much shorter time with respect to relative post-processing. The first methodology (NRTK) has been investigated and considered for many purposes, where precision farming [5], autonomous navigation, maritime survey [6], and meteorological monitoring [7,8,9,10,11] are only a few examples. This positioning technique has allowed the attainment of impressive accuracy, even for single-frequency GNSS receivers [12], and it has permitted the increased performance of the low-cost ones [13,14,15,16] since the implementation of studies into the potentialities of the GNSS system embedded in the smartphones [17]. Thus, the spread of low-cost technologies such as smartphones and tablets, with their rapid evolution in terms of the quality of their installed sensors, has increased the interest in these systems and the management of emergency scenarios. One primary advantages of using portable devices, such as those previously cited, is the possibility of completing a rapid survey and exploiting their embedded sensors, which results in them being useful for multidisciplinary teams cooperating in a coordinated manner on a common task, e.g., reconnaissance, inspection, and the survey of unstable structures.
Starting in 2016, with the Android Nougat 7.0 Operating System (OS) development, Google has permitted direct access to the GNSS chipset raw measurements mounted on some Android-based smartphones. The possibility to manage pseudo-range and carrier-phase measurements from the GNSS chipset installed on smartphones and tablets with an Android OS has changed the concept of precise positioning with portable devices. Several studies have been conducted to verify the feasibility [18] and positioning accuracy [19,20] of smartphones for different purposes, from urban [21,22,23,24,25] to pedestrian positioning applications [26,27], always facing the problems related to the high-level Application Programming Interface (API) and the filtered measurements provided by the GNSS chipset.
In [28], the authors demonstrated that it is possible to reach a decimeter level of accuracy in terms of positioning performances following the post-processing approach, made by double differencing raw smartphone observations. Meanwhile, the authors of [29] first focused their attention on single-base RTK positioning and then demonstrated the possibility of obtaining a centimeter-level accuracy through the use of NRTK corrections [30].
Recently, attention has been moved to PPP for obtaining the absolute position of a single receiver without the use of corrections or base stations [31]. However, to the best of our knowledge, nobody has focused their attention on the derivate products of GNSS positioning, such as the estimation of ionospheric or tropospheric delays.
Recently, ZTD has been estimated using the PPP strategy [32,33], but it was initially estimated within the GNSS station network adjustment to allow it to extend the local estimates to the entire area covered by the network through interpolation.
In this context, this paper aims to investigate the possibility of retrieving ZTD estimates from smartphones and to evaluate their accuracy. To evaluate their quality, ZTD values obtained from smartphones have been compared to those estimated with a geodetic GNSS receiver settled a few meters away from the smartphone. Two different types of software (one online, the CSRS-PPP [34], and one open-source, RTKLIB 2.4.3 b33 [35]) and two different test sites have been considered to make the analyses as general as possible. The positioning solutions have been computed to verify that the improvements in the ZTD estimations do not negatively affect the coordinate estimations, because both have been estimated from the same observations.
The remainder of this paper is organized as follows. Section 2.1 describes the ZTD estimation from GNSS observations, while Section 2.2 presents the case studies. Section 3 highlights the results of ZTD estimations regarding both smartphones and geodetic receivers. Section 4 provides the discussion on the research outcomes, and Section 5 reports the conclusions and gives some ideas for future research activities and investigations.

2. Materials and Methods

2.1. ZTD Estimation from GNSS

It is widely understood that GNSS was initially designed to determine the position, expressed by cartesian coordinates (X, Y, Z) of a receiver by means of pseudo-range or carrier-phase measurements. Both these approaches are affected by biases, primarily produced by the Earth atmosphere, which increases the real distance between satellite and the receiver, and by instrumental (both receiver and satellites) and site-specific causes: receiver and satellites clock errors, multipath, etc.
The basic GNSS carrier-phase observable, in units of length, can be written as follows [36]:
φ R S ( t ) =     ρ R S ( t ) + c     ( δ R ( t ) δ S ( t ) ) +   λ N R S ( t ) I R S ( t ) +   T R S ( t ) + MP ( t ) + ε
where λ is the wavelength of the carrier-phase, ρ R S is the geometric satellite-receiver distance (obtained by the satellite and the receiver coordinates), c is the speed of light, δS and δR are the satellite and receiver clock errors, respectively, N R S is the unknown initial phase ambiguity, I R S is the ionospheric (slant) delay, T R S is the tropospheric (slant) delay, MP is the multipath effect, and ε contains all the remaining unmodelled errors.
Most of the terms in Equation (1) can be neglected by means of appropriate precautions during the installation of the instrumentation and expedients when the GNSS receiver is operative. For example, ρ R S estimation can be improved through accurately known satellite and receiver positions, δS and δR can be eliminated using the double differences strategy or can be modelled in PPP, N R S is resolved in the initialization phase, and I R S can be eliminated using the iono-free dual frequencies combination for double-frequency receivers, or it can be modelled for single-frequency receivers. The multipath effect MP can be mitigated by choosing an antenna position that is far from reflecting objects, installing special antennas (e.g., choke-ring antennas), masking low elevation signals, or by using techniques and models designed to recognize the reflected signals and discard them [37].
The presence of the atmosphere introduces effects on the signal crossing that influence the antenna positioning. At the same time, these effects may be a starting point to monitor the troposphere itself, representing a useful knowledge base for meteorological purposes. Since the early 1990s, GNSS has also been used for meteorological purposes [38] because of its high accuracy, all-weather operation capability, high temporal resolution, and compatibility with other observation systems, e.g., radiosondes, water vapor radiometers, and lately, environmental satellites. GNSS meteorology exploits the delay of electromagnetic GNSS signal transmission in the atmosphere due to water vapor, dry gases, hydrometeors, and other particulates [39,40]. As already stated, GNSS observations are affected by different sources of bias, one of which is caused by tropospheric refractions during the crossing of the atmosphere by the electromagnetic GNSS signal. The tropospheric effect is frequency-independent, and it cannot be reduced, but it produces a bias in each satellite-receiver observation that can be related to the so-called ZTD. The estimation of ZTD helps to enhance the positioning precision and represents a contribution to meteorological studies.
ZTD can be estimated by GNSS CORSs network adjustments to correct the computed delay, thus exploiting a tropospheric model that can be adopted in atmospheric conditions that are not standard. Several approaches have been carried out: the investigation of the vertical column of the atmosphere over a single station [41], the exploitation of existing national GNSS networks [42,43,44,45], and the implementation of specifically designed GNSS networks [46,47,48,49]. One element of innovation in this field was introduced by [8], who proposed to use existing regional, national, and international GNSS CORSs for the retrieval of ZTD, to be used in meteorological and climatological applications and studies.

2.2. Case Study: Materials, Methods and Processing

The Broadcom 47755 dual-frequency GNSS chip embedded in Xiaomi Mi8 smartphones was considered for these tests. It is capable of tracking GPS L1 C/A, GLONASS L1, BeiDou (BDS) B1, QZSS L1, Galileo (GAL) E1, GPS L5, Galileo E5a, and QZSS L5 signals. Regarding the geodetic receivers, two GNSS CORSs were considered to guarantee the highest quality in terms of collected signals and frequencies, as well as the number of tracked satellites.
Two different types of software have been employed and considered for this work: the CSRS-PPP and a modified version of RTKLIB 2.4.3 b33.
CSRS-PPP is an online application for GNSS data post-processing, allowing users to compute higher accuracy positions from their raw observation data. CSRS-PPP uses precise GNSS satellite orbit ephemerides to produce the corrected coordinates of a user located on a generic point, regardless of proximity to available base stations. The software can process RINEX observation data from single or dual-frequency receivers operating in static or kinematic mode. CSRS-PPP uses the best available ephemerides, and it allows the users to select final (+/−2 cm, combined weekly and available 13–15 days after the end of the week), rapid (+/−5 cm, available the next day), or ultra-rapid (+/−15 cm, available every 90 min) ones, as it is possible to see in the specific IGS website section [50]. With the new version, the software includes PPP with ambiguity resolution (PPP-AR) for data collected on or after 1 January 2018. The output of the processing is sent to the user in a compressed folder. In particular, a tropospheric zenith delay file (with .tro extension) is created, containing hydrostatic and wet zenith path delays and tropospheric gradients for each processed epoch. For all measurement campaigns considered in this research activity and all the employed receivers, the following parameters have been selected for the data processing, considering the CSRS-PPP software: a cut-off angle equal to 7.5°, the Vienna Mapping function as the a priori tropospheric model, and the a posteriori variance factor used to scale the covariance matrix equal to 1.
RTKLIB is an open-source software widely used by both the academic community and other interested parties. The version considered in this paper is based on an Extended Kalman Filter (EKF) employing Zero-Difference (ZD) measurement equations, such as single point positioning, but also considering the receiver and antenna Phase Center Variations (PCVs), to take into account the Phase Center Offset (PCO), defined as the relative position of the receiver antenna phase center with respect to the antenna reference point (ARP). In this research activity, the ANTEX format for the antenna model, including PCO and PCV data, has been considered and used. Even for the RTKLIB software, the cut-off angle was defined as 7.5°; by using the typical EKF formulation, the unknown parameters, including the receiver position and velocity, the receiver clock bias, the troposphere parameters and the ionosphere-free LC, and carrier-phase biases are estimated. For this work, the RTKLIB version, solid earth tides, Ocean Tide Loading (OTL), and pole tides are modelled and properly considered, the same as for the CSRS-PPP software. In both cases, precise ephemerides and clocks have been used considering MGEX products [51,52], as well as the PPP-AR algorithms as the ambiguity resolution method. To weight the observations, SNR values have been selected. For both types of software, the final solutions have been considered in the ITRF realization to be compliant with the GNSS processing products (e.g., ephemerides, ocean tide loading), considering a forward solution based on EKF.
Table 1 summarizes the processing settings for both of the considered software. It is important to highlight the fact that the duty cycle was disabled for all the data collected in the two campaigns.

2.2.1. Case Study 1

The acquisition test was performed on 11 December 2018; the smartphone acquisition lasted just over one hour and a half, starting from 11:49 a.m. until 1:30 p.m. UTC. This case study’s chosen location was the rooftop of the Politecnico di Torino, where the TORI GNSS CORS, which belongs to the SPIN3 GNSS network, is installed. A few meters away from that location, a Xiaomi Mi 8 smartphone (referred to as TOSM for this case study) was placed. TORI, used as a reference, is a geodetic-level multi-frequency (L1, L2, and L5) and multi-constellations receiver. The GNSS receiver embedded in the Xiaomi Mi 8 (Broadcom BCM47755) can record multi-frequency (L1 and L5) and multi-constellations GNSS signals. The pseudo-ranges and carrier-phase measurements were collected via the GEO++ RINEX logger app (available on Google Play) and then processed using the previously described software.

2.2.2. Case Study 2

The acquisition test was performed on 25 June 2020, and lasted seven hours, starting from 7:40 a.m. until 2:40 p.m. UTC. The chosen location for this case study was the Genoa University rooftop, and the GNSS receivers involved were GENU, which belongs to the Regione Liguria GNSS CORSs network, and a Xiaomi Mi 8 (referred to as GESM for this case study) smartphone. GENU, used as a reference, is a geodetic-level multi-frequency (L1 and L2) and multi-constellations receiver. The smartphone used in this case study was the same model as the one used for Case study 1. The GEO++ RINEX logger app was also employed for the processing of this dataset.
Figure 1 shows the GNSS receiver layout used in this case study: the GENU CORS is on the left upper corner of the rooftop, while the smartphone is placed a few meters away.

3. Results

In this section, the main results obtained after the data processing phase are shown. Firstly, a comparison between the TORI GNSS station and TOSM (Xiaomi Mi8 smartphone close to TORI) is presented, to show the different ZTD estimations obtained considering both CSRS-PPP and RTKLIB software and these two receivers. Particular attention is paid to the positioning solutions, to not only verify the tropospheric estimation; knowing that the implemented algorithms are based on Kalman filter, this choice has been made to exclude the possibility that a good ZTD estimation produces the detriment of positioning solutions. For this reason, the Up component of the positioning solution has also been considered for investigating this aspect. Considering these stations, the performances obtained using the RTKLIB software are presented, again making the same consideration on the positioning solutions as well as on the ZTD. Moreover, considering only the TOSM station, the obtained results with the two types of software have been compared to highlight the differences in ZTD estimations and make a statistical analysis of these differences. Finally, a second case study has been investigated over a longer time span, considering again both CSRS-PPP and RTKLIB.

3.1. Case Study 1

The elaboration results carried out for Case study 1 are reported below, starting with the one obtained with CSRS-PPP. In the CSRS-PPP solution, after the initialization phase, lasting about 20–30 min, both the positioning solutions (Figure 2) and the ZTD estimations (Figure 3) become more stable, with a standard deviation of a few millimeters in this last case. As it is possible to see from Figure 2 and Figure 3, there is a lack of measurements of about 13 min at about 12:45 UTC, due to unexpected events. Despite that, after a second initialization phase, the solution converges quicker than the first one.
It is worth highlighting that, considering TOSM, CSRS-PPP is not able to estimate ZTD, but it provides the values coming from the Vienna Mapping Function (VMF) model due to the poor quality of the observation and noisy data collected and saved in the RINEX file, which contains long gaps. To check the behavior of different software, RTKLIB was considered on the same data.
As expected, RTKLIB highlights a poor quality of the solution, obtained from the smartphone observations (TOSM), for both Up positioning component and ZTD (Figure 4 and Figure 5, respectively). The residuals of the three positioning components of TORI and TOSM with respect to the corresponding reference positions, i.e., the convergence solution, are reported in Figure 6 and Figure 7, respectively. Table 2 summarizes the statistics deriving from Figure 6 and Figure 7.
The solution instability in the initial time span (approximately the starting 20–30 min) is particularly evident in Figure 4. This corresponds to the solution needing time to converge. After this phase, the positioning solution becomes more stable, and it is coherent with those available in the literature. The average difference between TORI and TOSM is about 0.8 m in the vertical component, while the ZTD estimation is at the order of 0.2 m. In this case, RTKLIB is able to estimate ZTD, for both TORI and TOSM, without introducing the VMF model as CRSR-PPP did.
To analyze the differences of the ZTD estimations with the two types of considered software and to verify the independence of the software used, Figure 8 shows the comparison of tropospheric delay values obtained with the CSRS-PPP (in blue) and RTKLIB (in orange) software. Figure 9 shows the differences in ZTD values for the TOSM station, considering the CSRS-PPP and RTKLIB software.
Concerning the smartphone receiver, the two results are quite different; as already mentioned, if the online software is considered (CSRS-PPP), the solution is smoothed and stable over time because the software is not able to produce an estimation and it provides model-derived ZTD values, whereas considering the open-source one (RTKLIB), the solution is less stable, even if the values are quite reasonable. As shown in Figure 10, the differences between these two solutions vary from −0.04 m to 0.3 m, with a mean value of 0.2 m and a standard deviation of 0.1 m, as summarized in Table 3.
Considering the TORI GNSS station, the results are definitely better than in the previous case: again, CSRS-PPP software provides more stable and smoothed ZTD solutions if compared to the RTKLIB ones (Figure 11 and Figure 12), but the ZTD residuals are smaller than for TOSM because they are derived from an estimation process rather than from a model. The range of difference is from −0.01 m to 0.03 m, as seen from Figure 13, with a mean value of 0.007 m and a standard deviation of 0.01 m (Table 3).
For the comparison shown in Figure 13, it should be noted that solutions falling in the convergence time have been excluded. The considered time interval, therefore, is from 11:00 a.m. to 1:00 p.m. UTC. This improvement is undoubtedly due to the quality of raw measurements, which are better in the case of TORI than TOSM. Nevertheless, even using the smartphone, it is possible to reach feasible results and interesting values for ZTD using RTKLIB.

3.2. Case Study 2

As previously stated, for this case study, the analyses have been made considering the use of both CSRS-PPP and RTKLIB software for ZTD estimation. In this case, the two considered stations are GENU and GESM, the geodetic CORS and the smartphone, respectively. As highlighted in Figure 14, the behavior of the two types of software is analogous to what already emerged from the previous case study: CSRS-PPP produces modelled ZTD for GESM whereas it correctly estimates ZTD for GENU.
Similar to Case study 1, GENU and GESM data were processed also using RTKLIB. The ZTD estimations obtained with this elaboration are shown in Figure 15. The positioning outputs for this elaboration are shown in Figure 16, while Figure 17 highlights the differences between the estimated and reference coordinates for GENU and GESM estimated with RTKLIB.
Table 4 shows the statistics of the solution in Figure 17. For the two solutions, the statistics are computed in respect of the reference position, i.e., the convergence solution, and excluding the convergence time (i.e., from 8:15 UTC onwards). The residuals of the three positioning components of GENU and GESM with respect to the corresponding reference positions are reported in Figure 18 and Figure 19, respectively. Finally, Figure 20 shows the ZTD differences between GENU and GESM as computed with RTKLIB software
In order to check the coherence of the ZTD values with official products, the obtained values, both estimated and modelled, were compared with those evaluated from EUREF [53] for the GENO CORS [54], which is a few kilometers away from the considered test site. These results are shown in Figure 21, where it is possible to notice that the estimations made using CSRS-PPP are more consistent with the EUREF solutions than with the ones obtained with the modified version of RTKLIB, although this last one provides promising results.
The same comparison with official estimates has been performed considering GESM. Provided that the previous tests demonstrated the lack of performances of CSRS-PPP in estimating ZTD from smartphone data, only the comparison between the estimates computed with RTKLIB and EUREF is presented in Figure 22.
This demonstrates the software capability to process GNSS measurements extracted from portable devices in a good way, with differences in the order of 5 cm from three hours later than the convergence time (i.e., from about 12:00 UTC onward), and maximum differences in the order of 10 cm in the whole considered time span. This highlights the possibility to exploit raw GNSS measurements obtained from smartphones for ZTD estimations, at the same time guaranteeing the positioning accuracy. This aspect opens new frontiers, not only in positioning solutions with portable devices, but also for atmospheric monitoring.

4. Discussion

ZTD estimations were analyzed for two different case studies, considering two different applications: an online application (CSRS-PPP) and a free and open-source software (RTKLIB 2.4.3 b33), modified by the authors. Two different types of GNSS receiver were employed for each test site: a smartphone and a geodetic receiver. This latter was considered as a comparison to test and verify if the results obtained from the software could be regarded as reliable. Firstly, attention was focused on the Turin test site (TORI and TOSM GNSS stations); the differences of ZTD estimations were not negligible if the open-source software was considered, even for the TORI GNSS station. This behavior was even reflected in the coordinate estimations, reaching differences of about 1 m for each component after the initialization phase. It is important to highlight the fact that all analyses and comparisons were made after the convergence of the solution, which was reached in approximately 20 min. Considering the CSRS-PPP software, the performances were better than the previous case considering the positioning solutions, whereas nothing can be confirmed concerning the ZTD values because CSRS-PPP was not able to estimate ZTD and it provided ZTD values from the VMF model. Starting from these analyses, it seems that the processing software plays a crucial role not only in the positioning but also in the ZTD estimations. For this reason, another case study was selected, which also extended the duration of the measurement campaign, considering a time interval of about 7 h. Considering the CSRS-PPP software, also in the present case, the ZTD estimations for GESM (smartphone) were not estimated by the software due to the poor quality of observations collected in the RINEX file derived from the smartphone, whereas the software provided ZTD estimates for GENU (CORS). Considering the RTKLIB software, the performances were quite different; in this case, the maximum span between GENU and GESM ZTD estimates was about 15 cm (Figure 15), even excluding the convergence time. In order to better understand how the processing software impacts the solutions, a deep analysis was made, comparing the ZTD estimations obtained for GENU and GESM with RTKLIB (Figure 18), respectively. In Figure 21, a comparison between the ZTD values estimated by these two types of software and those computed by the EUREF service were made for GENU; from this, it was possible to notice that the estimations made using CSRS-PPP were more in accordance with the EUREF ones with respect to the estimates obtained with the modified version of the RTKLIB software. This does not mean the one software is better than the other, because it could be interesting to deeply understand which are the processing strategies for both types of software. If, for RTKLIB, it is possible to make changes in the processing parameters thanks to its open-source state, considering the online one (CSRS-PPP) only a few settings can be modified. Thus, it is not possible to deeply investigate the constraints applied. However, this research aims to verify if a ZTD estimation made using smartphone devices and free software is possible and feasible, and not to evaluate what is the best software in terms of performances. Our work was intended to show what can be obtained with these two free types of software and to test their performances dealing with smartphone data. In both cases, the differences were about 20 cm, even if, in the latter case, the trend was less stable than in the first one. Thus, even if the quality of the data collected and the characteristics of the receivers have great importance for the ZTD estimations, it has been demonstrated that the processing software plays a crucial role and can provide quite different results. Despite that, the ZTD estimations obtained from smartphones can be considered feasible and useful.

5. Conclusions

The spread of low-cost technologies such as smartphones and tablets, with their rapid evolution in terms of the quality of their installed sensors, has increased interest in these systems and their employment for positioning purposes. Of course, because they are not primarily developed for these applications, their performances are not comparable with those obtainable with geodetic receivers. However, starting from 2016 with the advent of the Android Nougat 7.0 OS, their use as positioning tools has been investigated, and attempts have been made to reach a sub-meter accuracy by exploiting different techniques. In this work, one of these devices was employed for the ZTD estimation, comparing the obtained values with those obtainable with geodetic receivers, commonly used for these applications. Two different types of software (one online and one open-source) and two different case studies have been considered to obtain solutions independently from the used processing software and the test site.
The obtained results are encouraging: the ZTD estimations obtained with the smartphone receiver are comparable with those obtained with the geodetic one, especially if RTKLIB is used. Indeed, the online tool (CSRS-PPP) incapable of producing ZTD estimates for smartphone-derived RINEX files, whereas it produces consistent ZTD estimates for geodetic GNSS CORSs, without losing the quality of the positioning solutions. The measurement campaigns considered in this work do not cover an entire day due to the smartphone battery capacity. In the future, we plan to find solutions to extend the session lengths to more than 24 h, mainly by exploiting an external power supply. It is quite challenging to perform measurement campaigns considering different consecutive days because we have not had the possibility to provide a power supply to the smartphones up to now. Furthermore, another improvement will be considering different smartphones to verify the results obtained with the model considered in this research activity. Because there are few GNSS receivers installed inside smartphones, we believe that interesting results could also be obtained with other devices, and our goal is to prove that in the future.
This study represents initial analysis regarding the capability of online software to positively process GNSS measurements extracted from portable devices, exploiting the PPP strategy. Moreover, it has also shown the possibility to exploit raw GNSS measurements obtained from smartphones for ZTD estimations. This aspect opens new frontiers, not only in the positioning solutions with portable devices, but also for the contribution that smartphones could provide for atmospheric monitoring; if used correctly, they could become an interesting and widespread tool that is useful for creating dense sensor networks, even ones that are deployable for atmospheric monitoring. This aspect will be investigated in the future, and it is one of the next steps that the authors want to pursue.

Author Contributions

Conceptualization, data interpretation, and analysis, writing—original draft preparation, all authors; GENU/GESM GNSS tests, L.B. and D.S.; TORI/TOSM GNSS tests, P.D.; GNSS processing, L.B. and P.D.; writing—review and editing, I.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available upon reasonable request.

Acknowledgments

The authors wish to thank the two anonymous reviewers whose precise and valuable comments helped to significantly improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, X.; Montillet, J.P.; Fernandes, R.; Bos, M.; Yu, K.; Hua, X.; Jiang, W. Review of current GPS methodologies for producing accurate time series and their error sources. J. Geodyn. 2017, 106, 12–29. [Google Scholar] [CrossRef]
  2. Bock, Y.; Melgar, D. Physical applications of GPS geodesy: A review. Rep. Prog. Phys. 2016, 79, 106801. [Google Scholar] [CrossRef]
  3. Angrisano, A.; Dardanelli, G.; Innac, A.; Pisciotta, A.; Pipitone, C.; Gaglione, S. Performance Assessment of PPP Surveys with Open Source Software Using the GNSS GPS–GLONASS–Galileo Constellations. Appl. Sci. 2020, 10, 5420. [Google Scholar] [CrossRef]
  4. Mendez Astudillo, J.; Lau, L.; Tang, Y.-T.; Moore, T. Analysing the Zenith Tropospheric Delay Estimates in On-line Precise Point Positioning (PPP) Services and PPP Software Packages. Sensors 2018, 18, 580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhang, N.; Wang, M.; Wang, N. Precision agriculture—A worldwide overview. Comput. Electron. Agric. 2002, 36, 113–132. [Google Scholar] [CrossRef]
  6. Moore, T.; Hill, C.; Norris, A.; Hide, C.; Park, D.; Ward, N. The Potential Impact of GNSS/INS Integration on Maritime Navigation. J. Navig. 2008, 61, 221–237. [Google Scholar] [CrossRef]
  7. Tsuda, T.; Sato, K.; Realini, E.; Oigawa, M.; Iwaki, Y.; Shoji, Y.; Seko, H. A real-time monitoring system of Precipitable Water Vapor (PWV) using a dense GNSS receiver network. J. Disaster Res. 2013, 8, 155–156. [Google Scholar]
  8. Sguerso, D.; Labbouz, L.; Walpersdorf, A. 14 years of GPS tropospheric delays in the French-Italian border region: A data base for meteorological and climatological analyses. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2013, XL-5/W3, 7–14. [Google Scholar] [CrossRef] [Green Version]
  9. Sguerso, D.; Labbouz, L.; Walpersdorf, A. 14 years of GPS tropospheric delays in the French–Italian border region: Comparisons and first application in a case study. Appl. Geomat. 2016, 8, 13–25. [Google Scholar] [CrossRef]
  10. Ferrando, I.; Federici, B.; Sguerso, D. 2D PWV monitoring of a wide and orographically complex area with a low dense GNSS network. Earth Planets Space 2018, 70, 54. [Google Scholar] [CrossRef]
  11. Benvenuto, L.; Ferrando, I.; Federici, B.; Sguerso, D. The GNSS for meteorology (G4M) procedure and its application to four significant weather events. In IX Hotine-Marussi Symposium on Mathematical Geodesy. International Association of Geodesy Symposia International; Novák, P., Crespi, M., Sneeuw, N., Sansò, F., Eds.; Springer: Cham, Switzerland, 2020; Volume 151. [Google Scholar] [CrossRef]
  12. Dabove, P.; Manzino, A.M. GPS & GLONASS mass-market receivers: Positioning performances and peculiarities. Sensors 2014, 14, 22159–22179. [Google Scholar] [CrossRef] [Green Version]
  13. Cina, A.; Dabove, P.; Manzino, A.M.; Piras, M. Augmented Positioning with CORSs Network Services Using GNSS Mass-market Receivers. In Proceedings of the IEEE/ION Position, Location and Navigation Symposium (PLANS), Monterey, CA, USA, 5-8 May 2014; pp. 359–366. [Google Scholar] [CrossRef]
  14. Cina, A.; Piras, M. Performance of low-cost GNSS receiver for landslides monitoring: Test and results. Geomat. Nat. Hazards Risk 2015, 6, 497–514. [Google Scholar] [CrossRef]
  15. Lachapelle, G.; Gratton, P.; Horrelt, J.; Lemieux, E.; Broumandan, A. Evaluation of a low cost hand held unit with GNSS raw data capability and comparison with an Android smartphone. Sensors 2018, 18, 4185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Odolinski, R.; Teunissen, P.J.G. An assessment of smartphone and low-cost multi-GNSS single-frequency RTK positioning for low, medium and high ionospheric disturbance periods. J. Geod. 2018, 93, 701–722. [Google Scholar] [CrossRef] [Green Version]
  17. Paziewski, J. Recent advances and perspectives for positioning and applications with smartphone GNSS observations. Meas. Sci. Technol. 2020, 31, 091001. [Google Scholar] [CrossRef]
  18. Humphreys, T.E.; Murrian, M.; van Diggelen, F.; Podshivalov, S.; Pesyna, K.M. On the feasibility of cm-accurate positioning via a smartphone’s antenna and GNSS chip. In Proceedings of the 2016 IEEE/ION Position, Location and Navigation Symposium (PLANS), Savannah, GA, USA, 11–14 April 2016; pp. 232–242. [Google Scholar] [CrossRef] [Green Version]
  19. Pesyna Jr, K.M.; Heath Jr, R.W.; Humphreys, T.E. Centimeter positioning with a smartphone-quality GNSS antenna. In Proceedings of the 27th International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS+ 2014), Tampa, FL, USA, 8–12 September 2014; pp. 1568–1577. [Google Scholar]
  20. Zhang, X.; Tao, X.; Zhu, F.; Shi, X.; Wang, F. Quality assessment of GNSS observations from an Android N smartphone and positioning performance analysis using time-differenced filtering approach. GPS Solut. 2018, 22, 70. [Google Scholar] [CrossRef]
  21. Masiero, A.; Guarnieri, A.; Pirotti, F.; Vettore, A. A particle filter for smartphone-based indoor pedestrian navigation. Micromachines 2014, 5, 1012–1033. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, L.; Groves, P.D.; Ziebart, M.K. Smartphone shadow matching for better cross-street GNSS positioning in urban environments. J. Navig. 2015, 68, 411–433. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, L.; Li, Z.; Zhao, J.; Zhou, K.; Wang, Z.; Yuan, H. Smart device-supported BDS/GNSS real-time kinematic positioning for sub-meter-level accuracy in urban location-based services. Sensors 2016, 16, 2201. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, L.; Li, Z.; Yuan, H.; Zhao, J.; Zhou, K.; Yuan, C. Influence of the time-delay of correction for BDS and GPS combined real-time differential positioning. Electron. Lett. 2016, 52, 1063–1065. [Google Scholar] [CrossRef]
  25. Adjrad, M.; Groves, P.D. Intelligent urban positioning: Integration of shadow matching with 3D-mapping-aided GNSS ranging. J. Navig. 2017, 71, 1–20. [Google Scholar] [CrossRef] [Green Version]
  26. Al-Azizi, J.I.; Shafri, H.Z.M. Performance evaluation of pedestrian locations based on contemporary smartphones. Int. J. Navig. Obs. 2017, 2017, 1–10. [Google Scholar] [CrossRef] [Green Version]
  27. Fissore, F.; Masiero, A.; Piragnolo, M.; Pirotti, F.; Guarnieri, A.; Vettore, A. Towards surveying with a smartphone. In New Advanced GNSS and 3D Spatial Techniques. Lecture Notes in Geoinformation and Cartography; Cefalo, R., Zieliński, J., Barbarella, M., Eds.; Springer: Cham, Switzerland, 2018. [Google Scholar] [CrossRef]
  28. Realini, E.; Caldera, S.; Pertusini, L.; Sampietro, D. Precise GNSS positioning using smart devices. Sensors 2017, 17, 2434. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Dabove, P.; Di Pietra, V. Single-baseline RTK positioning using dual-frequency GNSS receivers inside smartphones. Sensors 2019, 19, 4302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Dabove, P.; Di Pietra, V. Towards high accuracy GNSS real-time positioning with smartphones. Adv. Space Res. 2019, 63, 94–102. [Google Scholar] [CrossRef]
  31. Aggrey, J.; Bisnath, S.; Naciri, N.; Shinghal, G.; Yang, S. Multi-GNSS precise point positioning with next-generation smartphone measurements. J. Spat. Sci. 2020, 65, 79–98. [Google Scholar] [CrossRef]
  32. Hadas, T.; Hobiger, T.; Hordyniec, P. Considering different recent advancements in GNSS on real-time zenith troposphere estimates. GPS Solut. 2020, 24, 1–14. [Google Scholar] [CrossRef]
  33. Zhou, F.; Cao, X.; Ge, Y.; Li, W. Assessment of the positioning performance and tropospheric delay retrieval with precise point positioning using products from different analysis centers. GPS Solut. 2020, 24, 12. [Google Scholar] [CrossRef]
  34. CSRS-PPP. Available online: https://webapp.geod.nrcan.gc.ca/geod/tools-outils/ppp.php (accessed on 1 October 2021).
  35. RTKLIB 2.4.3. Available online: https://github.com/tomojitakasu/RTKLIB_bin/tree/rtklib_2.4.3 (accessed on 1 October 2021).
  36. Teunissen, P.; Montenbruck, O. (Eds.) Springer Handbook of Global Navigation Satellite Systems; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  37. Marzocchi, R.; Leotta, M.; Federici, B.; Delzanno, G. The NARVALO project: Real time collision avoidance system in a GIS environment based on precise GNSS positioning. Geoing. Ambient. E Min. 2017, 151, 33–38. [Google Scholar]
  38. Bevis, M.; Businger, S.; Herring, T.; Rocken, C.; Anthes, R.A.; Ware, R.H. GPS meteorology: Remote sensing of atmospheric water vapor using the Global Positioning System. J. Geophys. Res.: Atmos. 1992, 97, 15787–15801. [Google Scholar] [CrossRef]
  39. Niell, A.E. Global mapping functions for the atmosphere delay at radio wavelengths. J. Geophys. Res. 1996, 101, 3227–3246. [Google Scholar] [CrossRef]
  40. Solheim, F.S.; Vivekanandan, J.; Ware, R.H.; Rocken, C. Propagation delays induced in GPS signals by dry air, water vapor, hydrometeors, and other particulates. J. Geophys. Res. 1999, 104, 9663–9670. [Google Scholar] [CrossRef]
  41. Rocken, C.; Van Hove, T.; Ware, R. Near real-time GPS sensing of atmospheric water vapor. Geophys. Res. Lett. 1997, 24, 3221–3224. [Google Scholar] [CrossRef] [Green Version]
  42. Seko, H.; Shoji, Y.; Fujibe, F. Evolution and air flow structure of a Kanto thunderstorm on July 21 1999 (the Nerima Heavy Rainfall Event). J. Meteorol. Soc. Jpn. 2007, 85, 455–477. [Google Scholar] [CrossRef] [Green Version]
  43. Inoue, H.Y.; Inoue, T. Characteristics of the water-vapor field over the Kanto district associated with summer thunderstorm activities. Sci. Online Lett. Atmos. 2007, 3, 101–104. [Google Scholar] [CrossRef] [Green Version]
  44. Chen, B.; Dai, W.; Liu, Z.; Wu, L.; Kuang, C.; Ao, M. Constructing a precipitable water vapor map from regional GNSS network observations without collocated meteorological data for weather forecasting. Atmos. Meas. Tech. 2018, 11, 5153–5166. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, H.; Yuan, Y.; Li, W.; Zhang, B. A Real-Time Precipitable Water Vapor Monitoring System Using the National GNSS Network of China: Method and Preliminary Results. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 1587–1598. [Google Scholar] [CrossRef]
  46. Zhang, C.; Kuo, Y.; Dai, L.; Chu, Y.; Braun, J.; Zhang, J.; Li, Q.; Chen, M. The design and application of network of ground-based GPS water vapor monitoring stations to improve precipitation prediction in the Greater Beijing metropolitan area. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, 37, 517–522. [Google Scholar]
  47. Oigawa, M.; Realini, E.; Tsuda, T. Study of Water Vapor Variations Associated with Meso-γ Scale Convection: Comparison between GNSS and Non-Hydrostatic Model Data. SOLA 2015, 11, 27–30. [Google Scholar] [CrossRef] [Green Version]
  48. Barindelli, S.; Realini, E.; Venuti, G.; Fermi, A.; Gatti, A. Detection of water vapor time variations associated with heavy rain in northern Italy by geodetic and low-cost GNSS receivers. Earth Planets Space 2018, 70, 28. [Google Scholar] [CrossRef] [Green Version]
  49. Zhao, C.; Zhang, B.; Li, W.; Yuan, Y.; Li, M. Simultaneous retrieval of precipitable water vapor (PWV) and vertical total electron content (VTEC) by low-cost multi-GNSS single-frequency receivers. Earth Space Sci. 2019, 6, 1694–1709. [Google Scholar] [CrossRef] [Green Version]
  50. IGS Products. Available online: https://www.igs.org/products/ (accessed on 1 October 2021).
  51. Montenbruck, O.; Steigenberger, P.; Prange, L.; Deng, Z.; Zhao, Q.; Perosanz, F.; Romero, I.; Noll, C.; Stürze, A.; Weber, G.; et al. The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS)—Achievements, prospects and challenges. Adv. Space Res. 2017, 59, 1671–1697. [Google Scholar] [CrossRef]
  52. MGEX Products. Available online: https://igs.org/mgex/data-products/#products (accessed on 20 October 2021).
  53. Bruyninx, C.; Legrand, J.; Fabian, A.; Pottiaux, E. GNSS metadata and data validation in the EUREF Permanent Network. GPS Solut. 2019, 23, 1–14. [Google Scholar] [CrossRef]
  54. EUREF Permanent GNSS Network Tropospheric Delays. Available online: http://www.epncb.oma.be/_productsservices/troposphere/ (accessed on 4 November 2021).
Figure 1. Layout for Case study 2.
Figure 1. Layout for Case study 2.
Remotesensing 13 04567 g001
Figure 2. TORI (blue) and TOSM (Xiaomi Mi 8—orange) Up component estimated with CSRS-PPP.
Figure 2. TORI (blue) and TOSM (Xiaomi Mi 8—orange) Up component estimated with CSRS-PPP.
Remotesensing 13 04567 g002
Figure 3. ZTD values for TORI (estimated, blue) and TOSM (Xiaomi Mi 8—modelled, orange) obtained with CSRS-PPP.
Figure 3. ZTD values for TORI (estimated, blue) and TOSM (Xiaomi Mi 8—modelled, orange) obtained with CSRS-PPP.
Remotesensing 13 04567 g003
Figure 4. Differences between estimated and reference coordinates for TORI (blue) and TOSM (Xiaomi Mi 8—orange) receivers, obtained with RTKLIB.
Figure 4. Differences between estimated and reference coordinates for TORI (blue) and TOSM (Xiaomi Mi 8—orange) receivers, obtained with RTKLIB.
Remotesensing 13 04567 g004
Figure 5. ZTD estimates for TORI (blue) and TOSM (Xiaomi Mi8—orange) obtained with RTKLIB.
Figure 5. ZTD estimates for TORI (blue) and TOSM (Xiaomi Mi8—orange) obtained with RTKLIB.
Remotesensing 13 04567 g005
Figure 6. Residuals of positioning components for TORI with respect to the reference position, as obtained with RTKLIB.
Figure 6. Residuals of positioning components for TORI with respect to the reference position, as obtained with RTKLIB.
Remotesensing 13 04567 g006
Figure 7. Residuals of positioning components for TOSM with respect to the reference position, as obtained with RTKLIB.
Figure 7. Residuals of positioning components for TOSM with respect to the reference position, as obtained with RTKLIB.
Remotesensing 13 04567 g007
Figure 8. Comparison of ZTD values for TOSM considering CSRS-PPP (modelled, blue) and RTKLIB (estimated, orange).
Figure 8. Comparison of ZTD values for TOSM considering CSRS-PPP (modelled, blue) and RTKLIB (estimated, orange).
Remotesensing 13 04567 g008
Figure 9. Differences in ZTD values for TOSM with CSRS-PPP and RTKLIB.
Figure 9. Differences in ZTD values for TOSM with CSRS-PPP and RTKLIB.
Remotesensing 13 04567 g009
Figure 10. ZTD residuals estimated on TOSM by CSRS-PPP and RTKLIB software.
Figure 10. ZTD residuals estimated on TOSM by CSRS-PPP and RTKLIB software.
Remotesensing 13 04567 g010
Figure 11. ZTD estimations comparison for TORI with CSRS-PPP (blue) and RTKLIB (orange).
Figure 11. ZTD estimations comparison for TORI with CSRS-PPP (blue) and RTKLIB (orange).
Remotesensing 13 04567 g011
Figure 12. ZTD difference between TORI considering CSRS-PPP and RTKLIB.
Figure 12. ZTD difference between TORI considering CSRS-PPP and RTKLIB.
Remotesensing 13 04567 g012
Figure 13. ZTD residuals estimated on TORI by CSRS-PPP and RTKLIB.
Figure 13. ZTD residuals estimated on TORI by CSRS-PPP and RTKLIB.
Remotesensing 13 04567 g013
Figure 14. ZTD values for GENU (estimated, blue) and GESM (Xiaomi Mi8—modelled, orange) obtained with CSRS-PPP software.
Figure 14. ZTD values for GENU (estimated, blue) and GESM (Xiaomi Mi8—modelled, orange) obtained with CSRS-PPP software.
Remotesensing 13 04567 g014
Figure 15. ZTD estimates for GENU (blue) and GESM (Xiaomi Mi8—orange) obtained with RTKLIB.
Figure 15. ZTD estimates for GENU (blue) and GESM (Xiaomi Mi8—orange) obtained with RTKLIB.
Remotesensing 13 04567 g015
Figure 16. Differences between estimated and reference coordinates for GENU (blue) and GESM (Xiaomi Mi 8—orange) receivers, obtained with RTKLIB.
Figure 16. Differences between estimated and reference coordinates for GENU (blue) and GESM (Xiaomi Mi 8—orange) receivers, obtained with RTKLIB.
Remotesensing 13 04567 g016
Figure 17. Differences between estimated and reference coordinates for GENU and GESM estimated with RTKLIB.
Figure 17. Differences between estimated and reference coordinates for GENU and GESM estimated with RTKLIB.
Remotesensing 13 04567 g017
Figure 18. Residuals of positioning components for GENU with respect to the reference position, as obtained with RTKLIB.
Figure 18. Residuals of positioning components for GENU with respect to the reference position, as obtained with RTKLIB.
Remotesensing 13 04567 g018
Figure 19. Residuals of positioning components for GESM with respect to the reference position, as obtained with RTKLIB.
Figure 19. Residuals of positioning components for GESM with respect to the reference position, as obtained with RTKLIB.
Remotesensing 13 04567 g019
Figure 20. ZTD difference between GENU and GESM considering RTKLIB software.
Figure 20. ZTD difference between GENU and GESM considering RTKLIB software.
Remotesensing 13 04567 g020
Figure 21. Comparison of ZTD values for GENU considering CSRS-PPP (orange, modelled) and RTKLIB (blue, estimated) software. The green dots represent GENO CORS ZTD estimates provided by EUREF.
Figure 21. Comparison of ZTD values for GENU considering CSRS-PPP (orange, modelled) and RTKLIB (blue, estimated) software. The green dots represent GENO CORS ZTD estimates provided by EUREF.
Remotesensing 13 04567 g021
Figure 22. Comparison of ZTD estimates for GESM considering RTKLIB (blue). The green dots represent GENO CORS ZTD estimates provided by EUREF.
Figure 22. Comparison of ZTD estimates for GESM considering RTKLIB (blue). The green dots represent GENO CORS ZTD estimates provided by EUREF.
Remotesensing 13 04567 g022
Table 1. Processing parameters for RTKLIB and CSRS-PPP software for both case studies.
Table 1. Processing parameters for RTKLIB and CSRS-PPP software for both case studies.
RINEX formatRINEX 3
Cut-off angle7.5°
EphemerisFinal (MGEX)
ClocksFinal (MGEX)
IonosphereIono-free
OTLBased on Chalmers grid solution
Ambiguity resolutionPPP-AR
Table 2. Statistical parameters for the position outputs for TORI and TOSM with respect to the corresponding reference positions.
Table 2. Statistical parameters for the position outputs for TORI and TOSM with respect to the corresponding reference positions.
Average (m)Standard Deviation (m)Root Mean Square (m)
TORIE-W0.160.290.33
N-S−0.100.200.22
U-D−0.040.140.14
TOSME-W−0.450.560.72
N-S0.300.340.45
U-D0.851.361.60
Table 3. Statistical parameters considering the ZTD residuals on TOSM (CSRS-PPP vs. RTKLIB).
Table 3. Statistical parameters considering the ZTD residuals on TOSM (CSRS-PPP vs. RTKLIB).
ReceiverMean (m)Standard Deviation (m)
TOSM0.1700.114
TORI0.0070.011
Table 4. Statistical parameters for the positions output for GENU and GESM with respect to the corresponding reference positions.
Table 4. Statistical parameters for the positions output for GENU and GESM with respect to the corresponding reference positions.
Average (m)Standard Deviation (m)Root Mean Square (m)
GENUE-W−0.030.020.03
N-S0.020.010.02
U-D−0.080.040.09
GESME-W0.230.350.42
N-S0.070.230.24
U-D0.020.180.18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Benvenuto, L.; Dabove, P.; Ferrando, I.; Sguerso, D. Preliminary Results on Tropospheric ZTD Estimation by Smartphone. Remote Sens. 2021, 13, 4567. https://doi.org/10.3390/rs13224567

AMA Style

Benvenuto L, Dabove P, Ferrando I, Sguerso D. Preliminary Results on Tropospheric ZTD Estimation by Smartphone. Remote Sensing. 2021; 13(22):4567. https://doi.org/10.3390/rs13224567

Chicago/Turabian Style

Benvenuto, Lorenzo, Paolo Dabove, Ilaria Ferrando, and Domenico Sguerso. 2021. "Preliminary Results on Tropospheric ZTD Estimation by Smartphone" Remote Sensing 13, no. 22: 4567. https://doi.org/10.3390/rs13224567

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop