The paper introduces an essentially new algorithm for calculating the GNSS position as an alternative to the least-square method. The proposed approach can be widely applied to any positioning method that uses multiple position lines for position calculation and is an example of how using a numerical solution can improve position accuracy without access to historical data. In essence, the method is based on the adaptation of the median filtering method widely used in the field of image processing, while at the same time applying a combinatorial approach and order statistics. The proposed solution makes it possible to improve on and assess the credibility of a single measurement. The article highlights the differences between the proposed and currently used approaches, as well as their advantages and disadvantages. The algorithm has been extensively tested under various environmental and weather conditions. The tests were carried out in typical and also in very demanding conditions, thus taking into account the real application context, i.e., pedestrian and car navigation in densely urbanized areas.