gap of current data, and furthermore increase the depth and quality of observational data sets. In this paper, we will use the model-independent method Gaussian processes (GP) to reconstruct one-dimensional function of the angular diameter distance ratios, with fixed lens (or source) redshift. An obvious benefit of this approach is that GP allow one to reconstruct a function from data directly without any parametric assumption, which has been widely used in various studies [@Seikel12a; @Seikel12b; @Cai; @Yennapureddy; @Melia18b]. The first (to our knowledge) formulations of this approach can be traced back to @Yennapureddy, which revisited the most recent and significantly improved observations of early-type gravitational lenses (158 combined systems) to distinguish $\Lambda$CDM another Friedmann-Robertson-Walker (FRW) cosmology known as the $R_h=ct$ universe. Their results showed that, the probability of $R_h=ct$ (which is characterized by a total equation of state $w=-1/3$) being the correct cosmology is higher than that of $\Lambda$CDM, with a degree of significance that grows with the number of sources considered. Therefore, although the differentiation of competing cosmologies is already quite competitive compared with those from other methods, it still suffer from the small number of lenses in the statistical sample.

In the near future, the next generation of 