1. Introduction
In the modern era of ubiquitous computing where the demand for Location Based Services (LBS) is increasing exponentially, positioning requires two fundamental characteristics: pervasive coverage and infrastructure independence. Precise location determination plays a very crucial role, not only for LBS but also for health care, assets monitoring, emergency search, rescue response, etc. Outdoor localization is mostly based on the GPS technology as it can provide an estimated location as accurate as 5 m. However, in many situations where the user location is surrounded by high-rise buildings or cloudy environments, the positioning accuracy can drop drastically to as low as 50 to 76 meters [
1,
2]. GPS, as well as Russian GLONASS constellation, provide similar positioning accuracy [
3]. However, for the indoor environment, it does not fulfill the first requirement of being pervasive due to the signal attenuation caused by reflection and refraction, and lower signal-to-noise ratio (SNR).
The researchers spotlighted other emerging technologies that can be used for positioning, especially in the indoor, and GPS denied environments. The primary focus of such research is to utilize the modern smartphones which are pervasive and infrastructure independent. Unlike RFID (Radio Frequency ID), IR (Infrared), iBeacons, Ultrasound, vision, and Wi-Fi based indoor positioning systems, smartphone based systems use built-in sensors of smartphones [
4]. RFID, Wi-Fi, and Bluetooth based systems are called infrastructure dependent systems as they require additional hardware to perform localization. Infrastructure dependence in this paper refers to the hardware infrastructure without which the positioning and localization systems cannot work.
The smartphone based systems are cheap as the installation and maintenance of additional infrastructure is not needed. One limitation of these systems is their dependence on GPS to know the place/building information before indoor positioning can be performed. Many such techniques consider that the starting position is known as a priori or consider a naive approach that the GPS calculated position of a building is correct where the user currently is, before localizing it indoors. The building identification is an important element of localization systems which perform both indoor and outdoor localization. GPS accuracy to determine the user’s current building is limited as stated above. Many current indoor localization systems operate on the assumption that the correct building information is already known. Thus, the localization is performed with respect to the known building, e.g., in the case of the Wi-Fi based system, the map of the current building and the database is loaded to perform the localization. The correct building information of the user plays a pivotal role in such localization systems. Thus, many research works have proposed the integration of GPS with other sensors [
5,
6]. However, the problem is the uncertainty of GPS satellites, which is caused mainly by the signal fading. The signal fading on account of the canyons and closely constructed multi-structural buildings can hinder GPS positioning accuracy. Similarly, semi-outdoor environment affects GPS performance and can wrongly localize a pedestrian to a wrong building. Additionally, GPS localization is unreliable and can be inaccurate when the signal becomes weak e.g., in the proximity of tall buildings [
7]. Moreover, GPS sensors become erroneous and uncertain in the case of sporadic unavailability of the satellite signals [
8]. Many indoor localization schemes simply assume that a clear pre-knowledge on the building information has been known, but such assumptions hardly remain valid in real scenarios.
We propose the use of the magnetic field data to identify the building where the user is currently walking. The magnetic sensor in the smartphone provides magnetic
x,
y, and
z values that can be used for this purpose. The earth’s uniform magnetic field within a small constrained area is disturbed by the presence of ferromagnetic materials like elevators, door, and concrete structures and causes direction anomalies. However, such disturbances are researched to serve as a unique signature also called the fingerprint [
9] and can be leveraged to identify different buildings. The magnetic field has some inherent limitations, e.g., varying magnetic attitude where different devices show different values for the same magnetic field [
10] as well as different values by the same device during the different time of the day. Similarly, low discernibility (strength of magnetic field is very weak, i.e., in few
T) and sensitivity (magnetic field reading is affected by indoor infrastructure) also need to be considered. However, in spite of such limitations, many localization systems [
11,
12,
13,
14,
15,
16] have been proposed, which utilize the magnetic field for localization.
This research presents a novel technique which makes the use of the smartphone magnetometer for building identification and, to the best of authors’ knowledge, it is the very first attempt to utilize the magnetic field to identify buildings. The contributions of this research include an infrastructure free building identification model which solely works on smartphone sensors and a deep learning based ANN model which can predict a user’s activity of walking and stationary. Additionally, we use the Binary Grid (BG) to prepare the fingerprint database which stores the patterns formed by the magnetic field instead of magnetic magnitude itself. This overcomes the limitation of Dynamic Time Warping (DTW) whose performance is affected by the varying magnitude of magnetic data collected over different time.
The rest of the paper is organized in five sections.
Section 2 describes the related work.
Section 3 illustrates the proposed method and its working mechanism. The experiment setup and results are discussed in
Section 4 while discussions and the conclusion are given in
Section 5 and
Section 6, respectively.
2. Related Work
The building identification is not a new problem and has been studied for quite a long time. GPS is considered the chief technology for building localization and it works quite fine. Besides GPS, other assisting technologies including cameras and most recently barometers are used to achieve a higher accuracy.
The authors in [
5] presented a method to identify a building using a smartphone camera, local 2D Google maps, and GPS. A partial top view of the building, orientation, and distance relative to the camera are fed as inputs to a probability model which is based on kernel density estimation. In addition, it also utilizes smartphone accelerometer and magnetic sensors to estimate the camera ray directions as input features. The results indicate that the proposed system achieves an accuracy of 92.7% to correctly identify the buildings with a processing time of approximately six seconds for 87% of the cases.
The authors in [
17] propose an approach for large area classification using a WiFi fingerprint based method. The approach determines whether the user to be localized is indoors in the area for which a radio map is already built or outdoors in the vicinity. The study includes the investigation of the accuracy of several techniques for classification including support vector data description (SVDD), self-organizing map (SOM), minimax probability machine (MPM) and principal component analysis (PCA). The results show that PCA can achieve a classification accuracy of 95.69%. The authors in [
18] present a method which utilizes the SNR of GPS to locate a user near the building entrances. The SNR of GPS drops closer to the building walls while GPS provides the accurate location information for a short period of time before entering into the buildings. This transition time is used to find the most probable building location of the user. The technique makes the use of GPS satellites behind the user and at high altitudes to this end and asserts that such satellites can provide more reliable information. The research work also confirms that the accuracy of the GPS localization is affected near the buildings, especially, when LOS (line of sight) with GPS satellites is lost.
The authors in [
19] propose the fusion of GPS signals with smartphone camera images to complement the deficiencies of each technique. The approach is based on a bags of words (BOW) approach which utilizes the important features selected from camera images to rectify the positioning error caused by GPS noisy signals. Both scenarios, where GPS signals correct the visual localization and vision-based localization to improve the GPS localization, are studied. The results show that the proposed method can perform the localization within a 7.5 m distance. The authors in [
20] present a method to mitigate the error of GPS position estimates in the urban areas. The proposed approach reduces the effect of signal reflection from NLOS (non-line-of-sight) satellites. The approach utilizes the publicly available 3D models but re-calculates the heights using the smartphone camera. Using 3D models and height information, the model estimates the path inflation caused by the signal reflection and re-estimates the calculated location. The adopted method is able to reduce the localization error by 6 to 8 m.
Several similar research works focus on IO (indoor-outdoor) environment detection as it is one of the important tasks in the active field of indoor positioning and localization. The majority of existing research works on indoor localization operates on the assumption that the user is in an indoor environment already when they start to locate him indoors. However, this poses challenges in scenarios where the localization environment contains both the indoor and outdoor environments, like the open covered space between two buildings and over bridge open corridors, etc. Thus, building identification is also very important in such situations where the positioning system needs contextual status to switch between indoor and outdoor positioning modules. The authors in [
21] propose a system for indoor–outdoor detection by leveraging the onboard Bluetooth module. Sparse low power Bluetooth devices are installed at the specified landmarks at intermediate regions e.g., covered corridors and entry/exit points. The conducted experiments reveal that the system achieves higher accuracy especially for semi-outdoor environments and consumes less battery than the existing schemes. The system achieves a localization accuracy of 10 m at 90 percent.
The authors in [
22] use Bluetooth signals to detect the user indoor and outdoor location. The underlying fact about the GPS accuracy that it is good outdoors and worse indoors, is used in an inverse way for the Bluetooth. Thus, when the distance estimated from Bluetooth is larger, the user is outdoors and vice versa. The experimental results show good indoor–outdoor switching as compared to the GPS. The authors in [
23] present an IO detector which is based on four various components including a signal serving cell tower, Wi-Fi, activity recognition, and light intensity captured from the smartphone. It makes the use of these modules and can discriminate four different areas of indoor, outdoor, urban, and complex environments. The experimental results demonstrate that the proposed system outperforms GPS in all four of the test environments. The discussed research works perform well in certain specific environments; however, they are limited in functionality on account of certain factors. First, most of the listed methods rely heavily on GPS, so the inaccessibility of GPS severely limits their functionality. Secondly, other similar systems pertain to additional installation, e.g., Bluetooth or Wi-Fi based systems that make them vulnerable and incur additional cost as well. Additionally, GPS based systems are power-hungry and consume more batteries. This research aims to present a solution which neither utilizes GPS nor necessitates any additional infrastructure. It rather makes use of the built-in sensors which are already made available in the smartphones.
3. The Proposed Method
The proposed building identification method utilizes the built-in sensors of the smartphone.
Figure 1 shows the block diagram of the proposed system. Two built-in sensors including accelerometer and magnetometer are utilized in the proposed method. The accelerometer data is used to identify whether the user is walking or stationary. On the other hand, the magnetometer data is utilized to identify the building where the user currently is.
The fingerprint database is built using the magnetometer data for this purpose. The concept of the fingerprint is to connect a fingerprint/signature to relevant locations in a given indoor environment. The fingerprint represents a feature like signal strength, a pattern of a signal or any other similar characteristic which makes it possible to distinguish between various locations. Analogous to the Wi-Fi fingerprinting where RSSI (Received Signal Strength Indicator) [
10] serves as a fingerprint, in magnetic fingerprinting, the magnetic flux density, also referred to as the magnetic field strength, is used as a fingerprint. The area is divided into grid points separated by a distance of 1 meter. A total of 100 samples are collected at each point and their normalized value serves as the fingerprint for the associated location. There are two possibilities for a user to be in a specific building:
- (i)
The user is entering a specific building now.
- (ii)
The user is already in a building.
For the first option, the building identification can be performed by making a fingerprint database of only the portion close to the entrances. Since we use only 4 sec data, we need to collect the fingerprints of only a 4 to 6 meter area, starting from the entrances. However, the second scenario requires building the database for the complete building, as we do not know where the user is exactly right now. Thus, keeping these two possibilities in mind, we make the fingerprint database of the whole building. Since during the building identification the magnetic data is collected while the user is walking, we need the fingerprints between the selected location points as well. The spline interpolation is applied to generate the intermittent values between the grid points. The fingerprint collection process and the resulting interpolated values are shown in
Figure 2.
Traditional magnetic fingerprint-based systems use the magnetic intensity as the fingerprint and thus have two major limitations. First, the magnetic intensity value may be different during different collection times even with the same device, which limits the performance. Second, various devices show different magnetic intensity even for the same location, which leads to device dependence. These issues are handled by using the magnetic field patterns in [
24]. The magnetic field intensity is transformed into magnetic patterns for this purpose. The patterns are stored as a matrix called Binary Grid (BG) and Algorithm 1 is used for that purpose. The BG contains 1’s and 0’s only; 1 where a magnetic field value is present and 0 otherwise. The pattern forming procedure is added here for completeness.
Algorithm 1 Generate the magnetic pattern |
- 1:
- 2:
- 3:
- 4:
- 5:
- 6:
fordo - 7:
- 8:
- 9:
end for
|
Let’s assume that the magnetic values shown in
Figure 3 are to be transformed into a magnetic pattern. First, a BG is declared with columns equal to the length of the magnetic data. The number of rows in the grid is the difference between the maximum and minimum magnitude of the given magnetic data. First, a matrix containing zeros of
size is declared where
N shows the number of rows and
M is the number of columns. Then, each magnetic data point is taken and converted to a row number by subtracting it from the maximum magnitude of data and rounding it. Afterwards, a 1 is added to the calculated row and this process is repeated for all the data values. The BG shown in
Figure 3a is for the first ten values only of the graph shown in
Figure 3b.
The matching between the user collected samples and magnetic pattern database is performed using the Dynamic Time Warping (DTW). The DTW was first proposed by Berndt and Clifford in 1994 [
25]. It was originally designed for speech recognition systems; however, it was later adopted by other domains as well. The DTW is a similarity measure that is used to calculate the distance between two time series. A time series is a sequence of measurement data over time. The DTW is used to measure the similarity for time series where the varying speed of the data collection source may result in a different length of the measured data. Similarly, the shape (magnitude) of the measured data may also change due to the speed, which makes it difficult for simple distance measuring methods to find the similarity. Simple distance measuring techniques are one-to-one linear alignment techniques while the DTW is one-to-many non-linear alignment, which makes it possible to find the similarity between the samples of different length. Consider that the
and
are two time series that we want to match. A warping path
W shows a contiguous set of matrix elements which defines a mapping between the given elements of
X and
Y. The warping path is built by arranging the elements of
X and
Y on the sides of a grid as shown in
Figure 4. Both sequences start from the bottom left of the grid.
Each cell of the grid shows the distance between the elements of two sequences. The warping path is through the cells that show the minimum distance. The overall distance represents the minimum of the sum of distances between the individual elements of given sequences. The DTW technique is good at comparing the sequences of varying lengths, yet it compares the magnitude of the given sequences. This limitation severely confines the systems which are based on the magnetic field data. The reason is the varying magnitude of the magnetic data. For example,
Figure 5 shows the magnetic data samples collected during the different time of the day with the same smartphone and the calculated distance using the DTW. It shows the Euclidean distance of 11.94 and 26.92 for the samples of the same place. Thus, although the samples are for the same location taken during different collection time, the distance is very different.
We already mentioned that the magnetic data for a given location may show different magnitudes every time we collect the data. Thus, if we use traditional DTW, it will produce erroneous results. Experiments reveal that the magnetic data magnitude may change; however, the patterns of magnetic data remain almost similar. Thus, we utilize the patterns formed by the magnetic data.
Figure 6 shows the distance calculated by using the DTW for the magnetic patterns and the fingerprint database. It is also important to mention that the fingerprint database also contains the magnetic patterns now. The magnetic pattern fingerprint database is built using Algorithm 1 discussed above.
The Euclidean distance calculated using the magnetic pattern of the given samples are very similar. Contrary to the previously calculated distance of 11.94 and 26.92, now the distance is 11.31 and 12.72 for the given samples. It shows that, when the DTW is used with the BG of magnetic patterns, it becomes more accurate and less prone to the changes in the magnitude of magnetic field data. The data collected from the magnetic sensor is used to make the fingerprint database while the building identification process uses the data from the accelerometer and magnetometer. The data from these sensors are collected over 4 sec and split into four frames each containing 1 sec data. Each frame is processed separately and the final decision is made after the processing of four frames is completed.
3.1. The Accelerometer
The accelerometer is used to identify the user’s state of walking or stationary in a given time frame
f. The size of a frame is 1 sec and it contains 10 measured values of acceleration data along
x,
y, and
z-axis as the data is collected at a sampling rate of 10 Hz (new measurement after 100 ms). Using the given three axes of accelerometer, the total acceleration is calculated as:
where
i shows the measurement for the given frame
f. We want to use the acceleration magnitude to detect if the user is walking. The acceleration data as shown in
Figure 7 is very different for the user’s activities of walking and being stationary.
A simple threshold-based movement detection can also be used to identify if the user is walking or standing stationary. However, the limitation of using the threshold-based heuristic approach is its low accuracy in many cases. Occasionally, phone shaking during the stationary mode generates the acceleration and the data crosses the threshold set to detect the user walking state. The ANN (Artificial Neural Network) identifies the sequence patterns formed while walking and its accuracy is higher than that of a simple threshold-based movement detection. The ANN has shown good performance even in such cases. We use an ANN for the classification of user activities as many types of research [
26,
27,
28] have utilized ANN for similar tasks and showed good classification results. The correct classification of user’s activity of walking is very important for the proposed method. The magnetic patterns need to be collected while the user is walking as when standing stationary there is no change in the magnetic field data. Thus, the magnetic pattern is a straight line as a result which will affect the overall accuracy of building identification.
The ANN is the mathematical modeling of human brain which imitates the function of brain neurons. Similar to a neuron, it is comprised of three layers: an input layer, an output layer and hidden layer(s). Edges are used to connect these layers and each edge bears a weight. Edge weights are determined during the process called back-propagation. An important process in ANN modeling is the selection of the appropriate feature set. Thus, a feature importance analysis is desirable before features are used for training. The Repetitive Feature Elimination (RFE) is performed for the current study to select the most suitable features. The RFE recursively eliminates less important features, and builds the prediction model on those features only, which remain present after the specified iterations. We use a set of eight features including x, y, z, and total acceleration, and variance of x, y, z, and total acceleration. The variance of acceleration in x, y, z, and total acceleration have higher importance so they are used in the ANN model. The Soft-max activation function is used for ANN. A total of 3000 samples of each feature are used for supervised training. The data is collected at a sampling frequency of 10 Hz and a low pass filter is used in the pre-processing phase to reduce the effect of noise. The ANN achieves an accuracy of 92.67% with 10,000 samples tested.
3.2. The Magnetometer Module
The pivotal part in identifying a building is performed by the magnetometer module as the identification process is totally performed using the data collected by the magnetometer. The correct building is found by matching the collected data against the fingerprint database which is built during the offline phase of the system. The proposed algorithm BLocate is used for localization.
Table 1 shows the notations used in algorithm 2 which identifies the buildings.
Algorithm 2BLocate Building identification with sensors’ data. |
- 1:
if User is walking then - 2:
for do - 3:
- 4:
for do - 5:
- 6:
- 7:
end for - 8:
- 9:
- 10:
end for - 11:
end if - 12:
;
|
The algorithm follows these steps to predict the building where the user is currently walking.
Step 1: The data is collected using magnetometer and accelerometer with a smartphone. First of all, the accelerometer data is utilized to classify the user’s mode of walking or stationary. If the user is in a walking state, then the magnetic data is used for further processing.
Step 2: The magnetic data is transformed into magnetic patterns. As it is previously mentioned that the magnetic fingerprint database is made up of magnetic patterns during the offline phase, we need to transform the magnetic data to magnetic patterns during the online phase to match them against the fingerprint database.
Step 3: During the third step, the online magnetic field pattern is matched against the database of all the floors of all buildings. Since the user’s walking speed varies, it results in different magnetic patterns depending upon the walking speed. We use the DTW for matching as it serves as a tool to compare two time-dependent sequences of varying length [
25]. The DTW can become computationally expensive when long sequences are matched. Thus, we take a database sample equal to user sample (1 sec data) and match them to calculate the Euclidean distance:
Step 4: Once the Euclidean distance between the user sample and the fingerprint database is calculated, we need to find the point where the error/distance is the minimum as that point is the most suitable candidate. Thus, each floor’s candidate is the lowest error estimate.
Step 5: Once the minimum error for each floor of each building is calculated, these errors are to be normalized. Since each building has a different number of floors and each floor may have a different error depending upon which floor the user is actually on. The magnetic field shows higher similarity to floors of the building where they are taken from. It infers that the magnetic data collected from a floor of a building A should show similar matching results to the other floors of building A as well. Thus, comparing the error value of a single floor is not appropriate. Instead, if we normalize the error with the total number of floors that would give more accurate results. The normalization is done using:
where
is the minimum Euclidean distance for a floor of the given building and
represents the number of floors in the given building.
Step 6: After normalization, the minimum error value represents the building where the user is right now. As mentioned earlier, the total user collected data is split into four consecutive frames, so this process of building identification is repeated for the four frames. Each frame possesses a vote and the final decision is based upon the votes of the four frames. After this process is complete, we have four predicted buildings which are called building candidates. In certain situations, it is also possible that the predicted buildings are very different so we need to correct the predictions. Since we use four frames data so the output candidates are four buildings and can have one of the following forms:
The prediction set has only one building as output.
The prediction set has three predictions for one building and one prediction for the other.
The prediction set has two predictions each for two different buildings.
The prediction set has one prediction for four different buildings.
The first set is very perfect. As the prediction is only one building, there are no outliers. In the second case, the building that has three predictions in the prediction set is selected as the building where the user is positioned. For the third option, where we get two various buildings, the building with lower error is selected as the target buildings. However, if the error is the same as well, then the building with the higher number of floors is the target building. The same procedure is adopted for the fourth case and first priority is given to the candidate with the lower error, and then to the building with a higher number of floors, when each building candidate possesses the same value of error:
3.3. Accuracy Measures
We need to understand four basic notations used in the classification before we move forward to calculate the accuracy measures of a classifier [
29].
True Positives (TP): These represent the positive tuples of a class that are correctly labeled by the classifier.
True Negatives (TN): These represent the negative tuples that are correctly labeled by the classifier.
False Positives (FP): These are the negative tuples that are incorrectly labeled as positive by the classifier.
False Negatives (FN): These are the positive tuples that are incorrectly labeled as negative by the classifier.
Accuracy is a measure widely used to evaluate the performance of a classifier and it is defined in terms of the percentage of the buildings predicted correctly [
29]. Using the above described terms of
TP,
TN,
FP, and
FN, accuracy is calculated using:
Precision and recall are two measures which are commonly used to evaluate the performance of a classifier. Precision is a measure which tells what percentage of all tuples are labeled positive that are actually positive. Precision is calculated as:
The recall is often referred to as the measure of completeness and it represents the percentage of true positive tuples which are labeled as such. It is also called the sensitivity and calculated as:
The
F score is a statistical analysis measure of classification, which considers both precision and recall of the classifier and computes a score between 0 and 1. The values closer to 1 represent the higher performance of the classifier and
F reaches 1 when the classifier performs 100% correct classification. The
F score reflects the average effect of both precision and recall. It is calculated as:
The Cohen’s kappa coefficient [
30] also represented as
k is a statistic that is used to measure the true inter-rater agreement [
31]. It defines the agreement between two evaluators, whereas we define it as the automatic prediction of the building and already defined label of the buildings. It takes the following mathematical form:
where
represents the observed agreements while
is the portion of the agreement that is expected by chance. There are many interpretations of Kappa value, although no values are universally accepted upon. The authors in [
29] state that if the Kappa value is higher than 0.80 it represents the higher strength of agreement and shows that the classification model has high classification accuracy. The values between 0.61 and 0.80 represent good agreement, between 0.41 to 0.60 moderate, between 0.21 to 0.40 fair, while the value lower than 0.20 shows very poor agreement.
5. Discussion
This study presents an approach to identify the building where a pedestrian is present now. The proposed approach is totally infrastructure free and works without any additional assisting technology. The technique is founded on the built-in sensors of the smartphone and leverages the accelerometer and magnetometer data only. The user collected data is matched against the fingerprint databases of each building of the university departments. The fingerprint databases are pre-built during the surveying phase. However, contrary to a traditional fingerprinting approach which stores magnetic magnitude as the fingerprint, we adopt the Binary Grid (BG) approach which transforms the magnetic magnitude into magnetic patterns and then stores them. The BG approach overcomes the limitation of Dynamic Time Warping (DTW) as well. The DTW approach shows different results for the data of the same place collected during different times. The BG approach followed in this paper shows very similar results for the data collected over different times. The data is collected at points which are 1 m apart and then interpolation is applied to generate the intermittent values between these points. The two important factors associated with a successful model i.e., response time and accuracy are investigated using two scenarios.
The first scenario involves the localization with four buildings which are geographically well separated. The proposed method performs very well both in terms of accuracy and response time and obtains an accuracy of 90.54% and a response time of 0.21 sec, which involves processing 4 sec data gathered from the accelerometer and magnetometer. The second scenario is more realistic where the experiment involves eight buildings with 29 floors in total. In spite of the fact that with an increased number of databases the accuracy can be severely affected, the overall accuracy of the proposed method is 86.37% and the average response time is 0.32 sec only. The accuracy of only two buildings is depraved and lowers under 70%, while other six buildings sustain the accuracy over 80%. In addition, Cohen’s Kappa value for classification with eight buildings is 0.84, which shows the strong inter-rater agreement. Supplementary experiments are carried out to investigate the factors which reduce the accuracy. The fundamental cause of the reduced accuracy is the higher similarity between the fingerprint databases which reduces the classification accuracy in essence. Despite the fact that the results illustrate that the proposed method has huge potential and can be utilized to work in GPS denied environments.
One criticism of the traditional process of fingerprinting approach is the laborious survey that requires both cost and time. Traditionally, during surveys, the data are collected at specified points separated by 1 m; we followed the same procedure. In addition, we, however, performed experiments with a database built using the data collected at points separated by 2 m. The results demonstrate that the overall accuracy is reduced to 81.47% from 86.37%. The elementary reason is the incorrectness of the interpolated values that are generated using the interpolation. The higher the distance between the points used for the interpolation, the more smoothed the interpolated values are that affect the uniqueness of the fingerprints. Thus, the labor for surveys can be reduced, but there is a trade-off between the surveying cost and accuracy.