Locating mobile terminals

[0001] The present invention relates to a method and an apparatus for locating a mobile terminal in a communications network, particularly but not exclusively using multiple sources of information.

[0002] The ability to pinpoint the location of mobile terminals is a desirable feature of a mobile telephone network. This is because of the need to provide customer services which rely on knowing the whereabouts of users of these services. For example, up-to-date local traffic information can be provided to enable a user to avoid nearby traffic jams. A user may also wish to know, for example, how to get to the nearest pub or restaurant from their present location. Clearly the location of the user must be ascertained to within even a few metres for this type of service to work.

[0003] Another reason for wishing to know the location of a mobile terminal is so that emergency services can locate a caller who is unable to provide an accurate personal location themselves.

[0004] It is known in a GSM mobile network to provide the location of a mobile telephone in terms of the cell of the network in which the telephone is located. Each cell contains one base station and a telephone is only ever in traffic communication with one base station at a given time. Hence the location of the telephone to an accuracy of the cell area can be determined simply by ascertaining with which base station the telephone is communicating. Such methods are known as cell-based location methods. Other methods can be combined with the cell identity (CI), such as a triangulation system, in which the location of a particular mobile phone is calculated using control signals from at least the three base stations closest to it (two of which are located in adjacent cells to the cell in which the mobile telephone is located). This system uses the assumption that the distance of the phone from a base station is proportional to the strength of the signal which the base station receives from it, or the time taken for the signal to travel between the phone and the respective base station. Thus the position of the phone can be determined by comparing the relative strengths or travel times of received signals between the three base stations and thus assessing the distance of the user from each base station. The actual location of the user is then obtainable geometrically since the location of the base stations is known and fixed.

[0005] In a 3GPP (3^rd Generation Partnership Project) network using a Wideband Code Division Multiple Access (W-CDMA) signalling system, it is possible for a mobile terminal to be in active communication with more than one base station at any one time. This situation is known as “soft handover” and differs from (hard) handover in a GSM system, in which a mobile terminal is “handed over” from one base station to another as it moves between cells of the network. Because of the nature of the soft handover, the above-described cell-based mobile location procedures suitable for GSM can not always be used in a W-CDMA type signalling system. It has therefore been necessary to develop more reliable ways of locating a mobile terminal in this type of signalling system.

[0006] In W-CDMA a “softer handover” is defined as well. In the case of “softer handover” the antennas of the base stations with which the mobile station is communicating are co-located (e.g. they are installed at the same physical location or site). In the remainder of this document, the term “soft handover” will be used also to cover the case of “softer handover”, and it will be understood by those skilled in the art that the invention and the described embodiments thereof are applicable to a softer handover situation as well as a soft handover situation.

[0007] Location techniques are known which are directed towards calculation of a location estimate and a “confidence region” associated with the location estimate of a User Equipment (UE).

[0008] One such technique is a Cell Identity and Round Trip Time (CI+RTT) location method.

[0009] The CI+RTT location method in UMTS relies on the availability of Round Trip Time (RTT) and Rx−Tx Time Difference (RxTxTD) measurements. RTT and RxTxTD measurements are introduced in UMTS FDD (Frequency Division Duplex) to allow the implementation of the CI+RTT location method.

[0010] The RTT is defined as RTT=T_Rx^UL−T_TX^DL, where TD_TX^DL is the time of transmission of the beginning of a downlink dedicated physical channel (DPCH) frame to a User Equipment (UE) and T_RX^UL is the time of reception of the beginning (the first detected path, in time) of the corresponding uplink DPCCH (Dedicated Physical Control Channel)/DPDCH (Dedicated Physical Data Channel) frame from the UE.

[0011] The RTxTxTD=T_Tx^UL−T_Rx^DL is the difference in time between the UE uplink DPCCH/DPDCH frame transmission (T_Tx^UL) and the first detected path (in time) of the downlink DPCH frame from the measured radio link (T_Rx^DL).

[0012] RTTs are measured by the base stations, RxTxTDs are measured by the UE.

[0013] By combining a pair of RTT and RxTxTD measurements referred to the same base station the distance between the UE and that base station can be estimated. Such distance estimate is analogous to the distance estimate that can be obtained from one Timing-Advance (TA) in GSM. In this sense, the CI+RTT location method corresponds to the Cell Identity+Timing Advance (CI+TA) location method in GSM. However, two particular features of UMTS FDD make the CI+RTT method potentially more accurate than the CI+TA method in GSM:

• [0014] 1. The much shorter UMTS chip period as compared to the GSM bit period affects the resolution with which a distance estimate can be determined from a TA in GSM or from an (RTT, RxTxTD) pair in UMTS. One GSM bit period is equivalent to approximately 1100 meters while one UMTS chip period is equivalent to approximately 80 meters, thus the distance measurements resolution in UMTS is better than in GSM.
• [0015] 2. In UMTS a UE can be in soft handover. UMTS standards require that RTTs and RxTxTDs are measured for each active radio link, thus multiple distance estimates can be potentially available for locating one UE in UMTS. In GSM this is not possible since the TA is available only for the unique serving cell.

[0016] In the CI+RTT location method the unknown geographical coordinates of the UE whose position it is desired to determine are estimated by combining absolute distance measurements between the UE and the base stations in the active set. Each absolute distance measurement is calculated from each (RTT, RxTxTD) pair.

[0017] Existing CI+RTT location methods can be considered to fall broadly into two general categories, single-site location methods and multi-site location methods. There are a number of good, reasonably accurate algorithms for implementing multi-site location methods when the UE have active radio links with cells located at three or more sites. However, such methods often fail when two or less sites are available for location calculation purposes.

[0018] One of the difficulties with existing location systems is that a location system method is selected to generate a location estimate which may then fail for one or more reason. It then becomes necessary to try a different method, which may or may not be implementable at the site where the first location method was carried out.

[0019] It is an aim of the present invention to obviate or at least mitigate this difficulty.

[0020] One aspect of the invention provides a method of providing a location estimate for a mobile station in a wireless communications network, the method comprising: implementing in parallel a plurality of different location estimate generation methods, each method generating a location estimate; and selecting one of the location estimates based on predetermined criteria.

[0021] Another aspect of the invention provides in a cellular communications network, a network entity adapted to provide a location estimate for a mobile station in the network, the network entity comprising: means for implementing in parallel a plurality of different location estimate generation methods, each method generating a location estimate; and means for selecting one of the location estimates based on predetermined criteria.

[0022] A first location estimate generation method is a weighted mass centre (WMC) method which uses location data for respective sites in radio communication with the mobile station.

[0023] Second and third methods are WSC and WSS methods which involve: using first measurement data relating to at least one radio link between the mobile station and a first site in a location calculation method with radio network data of the first site to provide a first preliminary location estimate; using second measurement data relating to at least one radio link between the mobile station and a second site in the location calculation method with radio network data of the second site to provide a second preliminary location estimate; and generating the location estimate as the weighted mass centre of the first and second preliminary location estimates.

[0024] In the preferred embodiment, the steps of providing the first and second preliminary location estimates comprise generating first and second respective distance estimates of the distance between the mobile station from the first and second site respectively.

[0025] According to the second method, where the first site supports a plurality of cells, the first measurement data relates to radio links of at least two of said cells and a single first distance estimate is generated using the first measurement data of the two radio links.

[0026] In the third method, where the first site supports at least first and second located cells, the first measurement data relates to a first radio link from the first cell, and further measurement data is provided which relates to a second link from the second cell. In that case, a further preliminary location estimate is generated which is used with the first and second preliminary location estimates to generate the weighted mass centre of the location estimate of the mobile station. In this case, respective first, second and further distance estimates are preferably generated in the step of providing the first, second and further preliminary location estimates respectively.

[0027] The predetermined criteria can be based on distances which are determined between each site and the location estimates. Where distance estimates are obtained as part of the estimate generation location methods. The distance estimates can be used in combination with these distances to minimise an error function, on the basis of which the location estimate to be used can be selected.

[0028] For a better understanding of the present invention and to show how the same may be carried into effect, reference will now be made to the accompanying drawings.

[0039] Before describing methods in accordance with embodiments of the present invention, a detailed description of one possibility for calculating a distance estimate will firstly be described with reference to FIG. 1, because this is an important component of the location methods discussed in the following. It will be appreciated however that other distance estimation techniques could be used.

[0040] FIG. 1 is a schematic block diagram illustrating the blocks involved in a distance estimate calculation obtained from a set of measurement pairs of the round trip time (RTT) and Rx−Tx time difference (TD). FIG. 2 illustrates what these distances represent in a mobile communications network. FIG. 2 illustrates a user equipment UE which can be in the form of a mobile telephone or other mobile terminal and a Node B, which represents a base station in the universal telecommunications radio access network (UTRAN) system. A downlink dedicated physical channel signal labelled DL DPCH in FIG. 2 dispatched from node B takes a certain amount of time to reach the user equipment UE. There is a time delay TD at the user equipment before an uplink response to the downlink signal is returned from the user equipment to Node B. This response is labelled ULDP(C/D)CH in FIG. 2. The total circuit time measured from dispatch of the downlink signal to receipt of the corresponding uplink signal is the round trip time RTT. That is, it represents the difference between the time of transmission of the beginning of a downlink dedicated physical channel (DPCH) frame to a user equipment UE and the time of reception of the beginning of the corresponding uplink UL dedicated physical control channel (DPCCH)/dedicated physical data channel (DPDCH) frame from the user equipment.

[0041] The time difference TD is, as explained, the difference in time between the uplink ULDPCCH/DPDCH frame transmission and the first detected part of the downlink DPCH frame. This delay is in a typical case nominally equal to 1024 chip periods. A measurement pair of RTT_i, TD_i values is supplied to the distance estimation equipment from each i-th active radio link in the mobile communications network. The manner in which measurements for co-located cells (that is cells served from the same site) are dealt with differs depending on the location method which is to be implemented, as will become clear in the following. For the purpose of describing FIG. 1, it is assumed that one distance estimate is obtained from a set of RTT, TD measurement pairs measured from active radio links supplied by co-located cells.

[0042] The measurement pairs are denoted by arrow 2 in FIG. 1, each pair being supplied to a respective probability function determination block 4. Each probability determination block 4 comprises a double range distance estimation calculator 6 which calculates a double range distance estimate d_DR from each RTT and TD measurement pair. d_DR is an estimate expressed in metres of the distance between the user equipment and the location of node B at the site of the serving cell. It is calculated by converting into metres an integer number, the so-called double range (DR) which is obtained by removing the TD value contribution from the RTT value. A statistical data block 8 holds a priori statistical data that has been gathered which gives an estimate of the likely measurement error in the double range distance estimate. That is, given the exact distance between the user equipment UE and the serving cell, d, it is possible to calculate the double range measurement error, x=d−d_DR. Such an error can arise as a result of several factors: multipath and non-line of sight propagation, measurement algorithms inaccuracy, granularity introduced by the measurements' finite resolution, etc. By premeasuring data of this kind, a probability density function of the double range measurement area can be determined and this is supplied to a PDF determination (probability density function determination) block 10 which combines the PDF of DR measurement error f_x|DR (x|DR) with the DR distance estimate d_DR, to obtain the PDF of the UE distance from each co-located serving cell f_d|DR(d|DR)=f_x|DR(x=d−d_DR|DR).

[0043] The probability density functions of the distance from all the co-located serving cells is then combined into one probability density function, f_d(d) representing the probability density function of the distance of the user equipment UE from the site where all the serving cells' antennas are installed. This is carried out in the combine PDF block 12. Once the combined PDF is available, the distance estimation is performed in the distance estimation block 14. The result of the distance estimation is represented by three values:

• • [0044] a distance estimate, d_EST≧0
  • [0045] an inner radius, 0≦R₁≦d_EST
  • [0046] an uncertainty radius, R₂≧0

[0047] The distance estimate d_EST is calculated as the median of the combined distance distribution from the PDFs. The inner and uncertainty radii define a confidence interval [R₁,R₁+R₂] within which the exact distance d of the user equipment UB is estimated to fall with a predefined confidence coefficient, ξ. The value of the confidence coefficient depends on the target confidence level requested by the user of the location algorithm, and is input into the distance estimation block 14 as shown in FIG. 1.

[0048] A description of a technique for estimating the location of a user equipment UE in accordance with an embodiment of the invention will now be described. FIG. 3 illustrates the physical situation in which this method may be used.

[0049] FIG. 3 represents part of a cellular communications network illustrating a first base station at Site A and a second base station at Site B. Each base station serves three cells with accordingly oriented antennas. The cells are labelled C₁, C₂ etc and the antennas are labelled a₁, a₂ etc. The orientation of the antennas at each site is known, and the assumption is made that the cell identity CI provides network data including for example antenna orientation data. A user equipment UE is shown in softer handover at Site A, that is with a first radio link RL₁ established with antenna a₁, and a second radio link RL₃ established with antenna a₃. The user equipment UE also has a radio link RL₅ established with antenna a₅ of the base station at Site B.

[0050] FIG. 3A is a schematic diagram showing how a selecting method for location methods (SELM) can be implemented in the network as illustrated in FIG. 3. FIG. 3A shows the two base stations labelled as Site A and Site B with their respective antennas which are indicated diagrammatically by arrows in approximately the antenna directions. According to the SELM method, a plurality of different location estimates are generated by implementing in parallel a number of different location estimate generation methods. The estimates which have been generated as a result of three such methods are labelled in FIG. 3A Est WMC, Est WSS and Est WCS. Est WSS is a location estimate generated by a weighted single site (WSS) method. Est WCS is a location estimate generated by a weighted single cell (WSC) method. Est WMC is an estimate generated by a weighted mass centre (WMC) method. The WSS and WCS methods will be described in more detail in the following. The WMC method is described in our Patent Application Publication No. (Nokia Ref. NC19511). The weighted mass centre method uses a weighted average of the site's coordinates. The weight applied to each site location is proportional to the reciprocal of the estimated distance between the mobile station and the site. The distance can be estimated using the technique described above with reference to FIG. 1, or in any other way.

[0051] The distances from the sites to the generated location estimates is then calculated. These are referred to herein as r_mn, where m denotes the method which has been used (by an integer 1, 2 . . . m . . . M) and n denotes the site. As an example, in FIG. 3A r1A and r1B refer to the distances from the WSC estimate to the Sites A and B respectively.

[0052] A measure of “goodness” of the estimate can be considered as e_mn=d_EST,n−r_mn, where d_EST,n is the distance estimate obtained from the location method itself, for example by the technique described above with reference to FIG. 1. In the schematic diagram of FIG. 3A, Site A has two distance estimates associated with it (being the inner and outer radii of the circles centred at Site A), and Site B has a single distance estimate associated with it, being the radius of the circle centred at Site B.

[0053] The estimate which is finally selected is based on the estimate which minimises the cost function f_m=sum_n{e_mn²}, or f_m=sum_n{|e_mn|}.

[0054] FIG. 3B is a schematic block diagram of the circuitry in terms of functional blocks for implementing the SELM method described above.

[0055] The WSS and WSC methods will now be described in more detail with reference to FIGS. 4 and 5.

[0056] FIG. 4 is a schematic diagram illustrating how a location estimate can be obtained using the WSS method in this situation. A first estimated location which is labelled Est 4 in FIG. 4 is obtained in the following way. RTT, TD pairs from each of the radio links RL₁ and RL₃ are combined in the distance estimation equipment illustrated in FIG. 1 and an estimated distance d_EST is produced, which represents the radius of the circle which is shown in full and labelled 20 in FIG. 4. Note that the dotted circles just inside and outside circle 20 represents different measurements for the estimated distance, which will be discussed more fully in relation to the WSS method discussed below.

[0057] A preliminary single site CI+RTT location estimate is obtained using the distance estimate obtained using the technique described above with reference to FIG. 1 and an angular estimation arrived at using a single site CI+RTT location based on the network data from the CI. An “averaged” orientation of the two antennas generating the radio links, in this case a₁ and a₃ is produced. Each preliminary location estimate is calculated in a polar reference system centred at each cell's site, using the single site CI+RTT method discussed later. The first location estimate Est 4 can thus be seen to lie at a distance d_EST4 along an averaged orientation of the antennas a₁, a₃. It will readily be appreciated that the orientation value which is used, being based on an average of the two antennas' orientation is a simple but primitive measurement. If more accurate orientation information relating to the user equipment UE is available, that could also be used to determine the first location estimate Est 4.

[0058] A second location estimate is obtained for Site B, in this case applying a single RTT, TD measurement pair to the distance estimation equipment of FIG. 1 and generating an estimated distance d_EST1 which in this case represents the radius of the circle denoted 24 in FIG. 4. Note that when a single RTT, TD pair is input into the distance estimation equipment of FIG. 1, there is no need to perform the combined PDF step of block 12, because only a single PDF is determined which is supplied directly to the distance estimation block 14. In this case, there is a single radio link RL5 and so the location estimate which is labelled Est 1 in FIG. 4 is generated by taking the estimated distance d_EST1 with the orientation of that antenna a₅, is the simplest case. Again, the angular estimation for the antenna orientation can be generated using the single site CI+RTT location method.

[0059] The final location estimate which is labelled Est WSS is obtained as the weighted mass centre of the preliminary location estimates, Est 1, Est 4. That is, the coordinates of the preliminary location estimate are taken (on a polar or Cartesian reference frame), weighted and then averaged.

[0060] In a Cartesian reference system, an algorithm which can be used in embodiments of the present invention to obtain the estimated x (and respectively, y) UE coordinates as a weighted average of the x (and respectively, y) site coordinates, the signal of which is received by the MS is given below:

[0061] $\hat{x}=\frac{\sum _{i=1}^N{w}^i{x}^i}{\sum _{i=1}^N{w}^i};\hat{y}=\frac{\sum _{i=1}^N{w}^i{y}^i}{\sum _{i=1}^N{w}^i};\left(\hat{x},\hat{y}\right)\in 𝒟$
where w¹, . . . , w^N are suitable weights assigned to each one of the N estimates involved, {circumflex over (x)},ŷ are the coordinates of the final estimate and xiⁱ, yⁱ are the coordinates of each preliminary location estimate.

[0062] One suitable definition for the i-th weight used in embodiments of the present invention is to use the inverse of the estimated distance, d_ESTi associated with each preliminary location estimate:

[0063] $w^i=\frac{1}{d_{\mathrm{ESTi}}}$
An alternative is to apply constant weights for all cells. A further alternative is to use weights equal to the inverse of the area of an arc-shaped confidence region which is determined as a by-product to the preliminary single site CI+RTT location calculation (discussed later). Another possible criterion for determining the weights is to assign heigher weights to preliminary location estimates which are believed to originate from more accurate distance measurements.

[0064] The Weighted Single Cell (WSC) method will now be described with reference to FIG. 5. This method is similar to the method WSS, the main difference between the two consisting in the way that measurements from co-located cells are treated. Note that the physical situation in which this embodiment can be used is the same as that for the WSS method, namely the situation of FIG. 3. In the WSC method each cell is considered separately from other cells. Thus, the number of preliminary location estimates is equal to the number of RTT, TD measurement pairs. The method is implemented as follows.

[0065] Firstly, a distance estimate for each of the serving cells is obtained using the distance estimation equipment of FIG. 1 using RTT, TD measurement pairs formed on each of the three active radio links RL1, RL2, RL3. These distance estimates are the respective radii of the circles labelled 21, 23 and 24 (as before). Note that the circles 21 and 23 are those which were illustrated dotted in the embodiment of FIG. 4.

[0066] A separate preliminary single site CI+RTT location estimate is then obtained for each distance estimate, calculated in a polar reference system centred at each cell site. Each preliminary location estimate has a distance from the cell site equal to the distance estimate, and an orientation based on the orientation of the particular antenna for that radio link (unless more accurate direction information is available). The preliminary location estimates are labelled Est 1, Est 2 and Est 3 in FIG. 5. The antenna orientations can be established using the single site CI+RTT location method.

[0067] The final location estimate which is labelled EST WSC is the weighted mass centre of the preliminary location estimates Est 1, Est 2, Est 3. Obtaining the weighted mass centre, and the weighting values is the same as for the WSS method described above.

[0068] In the above methods, reference is made to a CI+RTT location method. One example of such location method will now be described with reference to FIG. 6. The aim of the method is to estimate the UE geographical coordinates and a confidence region where the true UE location is supposedly located with a given degree of confidence ξ. The estimated UE coordinates are distance ρand orientation θof the UE from the serving site, represented in the polar reference system. Note that ρ—in FIG. 6 corresponds to d_EST in FIG. 1. In the following, there is some duplicate description relating to the distance estimation already discussed above, but in this case there is no discussion of how to deal with multiple sites. Note too that like numerals in FIG. 6 denote like parts in FIG. 1.

[0069] The estimates are obtained from a probability distribution of the UE geographical location, which is determined by combining observations and a-priori statistical information. A-priori information includes statistical information about the accuracy of DR as a distance measurement and distribution of the UE orientation within cells border.

[0070] Referring to FIG. 6, the location procedure can be seen as a combination of different processing steps, briefly described in the following. The quantization rules at the basis of RTT and Rx−Tx Time Difference definition in UMTS technical specifications, are used to calculate from each {RTT,TD} pair a Double Range Value, DR, and Double Range distance estimate, d_DR in block 6. The Double Range is an integer number indicatively equal to the difference RTT−TD; thus it represents an estimate of twice the absolute distance between target UE and (co-located) cells providing the active radio link.

[0071] As a result of the nominal distance estimation, the observation a=[RI,DR] (block 7) becomes available. RI is a radio link identifier, for example a cell identifier CI, providing network data.

[0072] For example, radio network data for all cells having an active radio link with the user equipment of interest, can include the following:

• • [0073] antenna coordinates
  • [0074] bearing of the node B's antenna measured clockwise from x-axis in radian φ
  • [0075] half power beam width (HPBW) of the node B's antenna, Δφ
  • [0076] R_F the maximum front radius to the service cell
  • [0077] R_B the maximum back radius to the service cell

[0078] The Double Range distance estimate is a raw conversion of RTT and TD, reported from UMTS network elements, in a distance estimate. This conversion does not take into account effects of measurement process impairments (multipath and non-line-of-sight propagation over the mobile radio channel, measurement errors, quantization error, etc.) that contribute to determining the so-called Double range measurement error, χ=ρ−d_DR, as defined as the difference between the exact distance between UE and node B, ρ, and Double Range distance estimate, d_DR. One way to take into account the DR measurement error, that has proved to be very successful in GSM applications, is to supply the location algorithms with a-priori information describing the statistical behaviour of Double Range measurement error, χ as in block 8.

[0079] The behaviour of −χ depends on many factors, among them environment and distance between UE and node B. For this reason, the PDF of −χ conditioned by the observation, ƒ(χ|a), is used as a-priori PDF from which the statistical behaviour of the UE distance from the serving site is derived. Once the PDF ƒ(χ|a) is available, the DR distance estimate, d_DR, is used to obtain the probability density function (PDF) (block 10) of the UE distance in the cell providing the active radio link, ƒ(ρ|a).

[0080] An analogous but simpler procedure is used to determine the PDF of the UE orientation in the cell providing the active radio link, ƒ(θ|a). The additional simplicity is due to the fact that an orientation measurement is not available, thus statistical distribution of the UE angular distribution can be defined only by using a-priori cell-based orientation information (block 9) such as cell orientation, width, etc.

[0081] The PDFs of distance conditioned by the observation and of orientation conditioned by the observations are then combined with each other (block 11) to determine the joint PDF of UE distance and orientation within the cell providing each active radio link ƒ(ρ,θa). The PDFs available at this stage describe the geographical distribution of users within the borders of each single cell served by node B providing each active radio link.

[0082] Once the joint PDFs of distance and orientation within all cells providing active radio links are available (blocks 4′) they are combined (block 12′) to determine the final joint PDF of distance and orientation, ƒ(ρθa). Such distribution includes all information gathered from observations (radio link identifiers and measurements) and a-priori statistics; thus it can be finally used to estimate UE coordinates ({circumflex over (ρ)},{circumflex over (θ)}) and confidence region R (block 14′). The priming of blocks 4′, 12′, 14′ denotes that they are equivalent to blocks 4, 12 and 14 of FIG. 1, but providing a more complex function.

[0083] Thus, the CI+RTT location algorithm provides as outputs:

• • [0084] location calculation results being a set of parameters identifying the coordinates of a location estimate and the parameters of a certain confidence region (a geographical region where the exact UE location is estimated to be with a certain probability)
  • [0085] QoS data which is an indication of whether the location calculation results met the accuracy requirements set in the input quality of service parameters.
    The CI+RTT location method involves the two main steps of:
  • [0086] 1. estimating the location of the UE in terms of x-y coordinates and
  • [0087] 2. calculating a confidence region for this location estimate.

[0088] A confidence region is a geographical region where the exact UE location is estimated to be with a given probability, referred to as the confidence coefficient 0<ξ≦1.

The location calculation method is implemented by location calculation algorithms. Two classes of location calculation algorithms are used:

• • [0089] Location Estimate Calculation Algorithms to implement step 1
  • [0090] Confidence Region Calculation Algorithms to implement step 2

[0091] The combination of location estimate and confidence region parameters is referred to as “shape”. The shape definitions supported by the location calculation algorithms described above are:

• • [0092] (i) Point Shape (i.e. including only the location estimate)
  • [0093] (ii) Point with Uncertainty Ellipse Shape (where the confidence region is an ellipse)
  • [0094] (iii) Point with Uncertainty Polygon Shape (where the confidence region is a polygon)
  • [0095] (iv) Point with Uncertainty Arc Shape (where the confidence region is a circular arc)

[0096] The only confidence region which will be described herein is the Arc-Shape shown in FIG. 7.

This has the following features:

• • [0097] Coordinates of a location estimate ({circumflex over (x)},ŷ)
  • [0098] Coordinates of an origin x_o and y_o
  • [0099] Inner radius R1 and uncertainty radius R2
  • [0100] Offset (orientation) angle α and included angle β
  • [0101] Confidence coefficient value ξ

[0102] In the context of the Single-Site CI+RTT method, the polar coordinates can be used to introduce the following formal definition of R:
PR((ρ,θ)εR|a₁, . . . , a_L)=ξ
The Arc is the “natural” shape for the Single-Site CI+RTT confidence region. The origin of the arc is at the antenna coordinates, common to all cells providing an active radio link to the UE being located:

[0103] $\{\begin{array}{c}{x}_0=x_S\\ y_0=y_S\end{array}\hspace{1em}$
Inner radius R₁, uncertainty radius R₂, offset angle and included angle (see FIG. 7) are calculated by means of the above that, when the confidence region is arc shaped and the PDFs of interest are not conditioned by the observations, can be rewritten as:

[0104] $\begin{array}{c}\mathrm{Pr}\left(\left(\rho ,\theta \right)\varepsilon R\right)=\int {\int }_{R}f\left(\rho ,\theta \right)d\rho d\theta \\ ={\int }_{\alpha }^{\alpha =\beta }d\theta {\int }_{R_1}^{R_1+R}f\left(\rho ,\theta \right)dp\\ =\xi \end{array}$
The total probability theorem can be invoked to write ƒ(ρ,θ) as:

[0105] $f\left(\rho ,\theta \right)=\sum _{m=1}^M{p}_{m}f\left(\rho ,\theta ❘a_m\right)\mathrm{hence}$ $\mathrm{Pr}\left(\left(\rho ,\theta \right)\varepsilon R\right)={\int }_{\alpha }^{\alpha +\beta }d\theta {\int }_{R_1}^{R_1+R_2}\left[\sum _{m+1}^M{p}_{m}f\left(\rho ,\theta ❘a_m\right)\right]d\rho$

[0106] The confidence region parameters R₁,R₂, α and β can be uniquely determined if appropriate additional constraints are introduced, according to requirements of a particular implementation.

[0107] FIG. 8 is a schematic diagram of the architecture of a network to illustrate where the location methods described above are implemented. A single base station BTS is shown, which can for example be the base station at site A or site B. In fact, base stations at both site A and site B will be in communication with a radio network controller RNC as illustrated for the single base station in FIG. 8. The radio network controller is in communication with a mobile switching MSC which is in turn in communication with a gateway mobile location center GMLC. The gateway mobile location center can receive an LIF request from commercial applications being implemented in the network. The gateway mobile location center dispatches a location request to the mobile switching center which in turn contacts the radio network controller to obtain the necessary network data and return it to the mobile switching center. A response is thus returned to the gateway mobile location center which supplies CI+RTT, TD pairs to a serving mobile location center SMLC. The location methods described hereinabove are implemented at the serving mobile location center and the x, y data of the selected estimate is returned to the gateway mobile location center. The gateway mobile location center returns this x, y data together with the confidence region to the commercial application which requested the location data.

[0108] Alternatively, the location methods can be implemented at the RNC or in a stand-alone module elsewhere in the network.