Systematic Evaluation of Applying Space-Filling Curves to Automotive Maneuver Detection
Abstract
Identifying driving maneuvers plays an essential role on-board vehicles to monitor driving and driver states, as well as off-board to train and evaluate machine learning algorithms for automated driving for example. Maneuvers can be characterized by vehicle kinematics or data from its surroundings including other traffic participants. Extracting relevant maneuvers therefore requires analyzing time-series of (i) structured, multi-dimensional kinematic data, and (ii) unstructured, large data samples for video, radar, or LiDAR sensors. However, such data analysis requires scalable and computationally efficient approaches, especially for non-annotated data. In this paper, we are presenting a maneuver detection approach based on two variants of space-filling curves (Z-order and Hilbert) to detect maneuvers when passing roundabouts that do not use GPS data. We systematically evaluate their respective performance by including permutations of selections of kinematic signals at varying frequencies and compare them with two alternative baselines: All manually identified roundabouts, and roundabouts that are marked by geofences. We find that encoding just longitudinal and lateral accelerations sampled at 10 Hz using a Hilbert space-filling curve is already successfully identifying roundabout maneuvers, which allows to avoid the use of potentially sensitive signals such as GPS locations to comply with data protection and privacy regulations like GDPR.
Index Terms:
space-filling curve, manoeuvre detection, Z-order curve, Hilbert curve, Morton codes, roundaboutI Introduction
The development and evaluation of Advanced Driver Assistance Systems (ADAS) and Automated Driving Systems (ADS) is a data-driven process that covers systematic and reproducible validation & verification (V&V) using closed-loop simulations, controlled experimentation at confined test sites using prototypical vehicles, and field data collections with larger vehicle fleets for open-loop replay or scenario reconstruction. While virtual environments provide systematic and cost-effective annotations by design, the degree of fidelity in tests using real sensors in real scenarios is still higher, as subtle effects that may not be part of the models used for the virtual sensors or the environment model can be captured. In addition, a much broader variety of situations can be potentially recorded that would require too much effort to be modeled manually in virtual test environments.
Large-scale field data collections apparently have a greater value during the system development and for V&V because of the representativeness of the captured scenarios and because of the data originates from the real sensor configuration. Furthermore, such field data can grow quickly by equipping hundreds of vehicles with data loggers that stream their data to a centralized storage, where it can be analyzed further. While the data collection itself can be conducted rather cost-effectively, the actual data processing afterwards to make use of the value of such data is becoming a growing challenge: Firstly, field data is usually non-annotated right after collection and relevant labels need to be added and verified (semi-)automatically before queries for maneuvers are possible; secondly, data privacy and protection regulations such as General Data Protection Regulation (GDPR) may not only require the anonymization of video data (for instance, blurring faces or license plates) but also potentially the removal of signal types such as GPS locations of the vehicle that collected the data to prevent the identification of individuals later.
The removal of certain signal types, though, may render the implementation and use of some maneuver detection and querying approaches infeasible. For instance, finding maneuvers when a vehicle is passing through a roundabout by defining geofences and checking when a vehicle is entering such geofences of interest would require potentially GDPR-sensitive information. Alternatively, adding labels for entering/leaving roundabouts at the time of collecting data would require the existence of a taxonomy of scenarios to label upfront, and added elements in such a taxonomy would then only be applicable to newly collected data but not to existing data from the past.
Motivation: There is a need for an approach to detect maneuvers and events that (a) does not require the definition and maintenance of an annotation taxonomy before data is collected, (b) does not depend on
(semi-)automatically added labels, and (c) is scalable to cope with growing amounts of data.
Berger and Birkemeyer (cf. [1]) introduced using a Z-order Space-Filling Curve (SFC) as a novel approach to explore automotive data in a computationally efficient way for event detection. In this paper, we are conducting a qualitative experiment with two different types of SFCs, namely Z-order and Hilbert curves, to answer the following research questions:
RQ-1:
To what extent can Z-order and Hilbert SFC be exploited to identify automotive maneuvers in an efficient way without relying on GDPR-sensitive data?
RQ-2:
How do different experiment configurations (e.g., used signals, sampling frequency, etc.) compare to manual annotations of automotive maneuvers?
We limit our systematic experiment on the detection of driving through roundabouts with a passenger vehicle (Volvo XC90). While the parametrization is adjusted for the experimental vehicle, the fundamental concept is transferable to other vehicles and scenarios. However, the best-performing combination of kinematic signals, though, may be affected by the chosen automotive maneuver in our experiments and may differ for other maneuvers and/or geographic regions with different lane layouts.
II Related Work
One fundamental issue with event detection based on kinematic signals originates from their multi-dimensionality, which is posing a computational challenge at scale. Brute-force-based approaches like the one suggested by Perez et al. [2] are on the one hand traceable and allow to understand contributing factors for a positive hit, but they scale proportionally with a data set’s size on the other hand and hence, are not practical for large data sets.
Alternative approaches that explicitly address multi-dimensionality are based on SFCs. An SFC enables a single-dimensional representation of multi-dimensional data (cf. Bader, [3]) and hence, effectively allows for dimensionality reduction. Intuitively, an SFC passes through each data point from a multi-dimensional space in a recursive or repetitive manner. Certain types of SFCs even preserve locality properties, ie., data points that are close to each other in the multi-dimensional data space will be located near to each after being mapped onto an SFC. This property in turn can be exploited for efficient data retrieval as operations on an SFC are usually faster and computationally cheaper on single dimensions than on multiple dimensions simultaneously. SFCs have been explored by Hulbert et al. [4] for example who evaluated the performance of SFCs on large data. Similarly, Dai et al. [5] compared both, Z-order and Hilbert curves applied to two-dimensional data, and Moon et al. [6] investigated the Hilbert SFC for indexing multi-dimensional data that can be used for data queries.
There are also complementary studies that employ SFCs for event detection across diverse application fields. For instance, Safia et al. [7] propose an algorithm to detect environmental events like air pollution that can be used for monitoring by sensor networks. Liu et al. [8] suggest a method called HGST that uses Hilbert curves to improve the efficiency of spatiotemporal range queries. The method was tested on a real taxi trajectory data set. Nair et al. [9] present their study on dynamic planning of evasive maneuvers in an environment, which represents obstacles and holes by using Hilbert curves.
Locations of roundabouts:
| latitude | longitude | |
|---|---|---|
| RA1 | ||
| RA2 | ||
| RA3 | ||
| RA4 | ||
| RA5 | ||
| RA6 | ||
| RA7 |
We introduced the novel pattern detection approach named Z-order Curve-based Event Retrieval Approach (ZEBRA) based on Z-order curves that we exploited for computationally efficient event detection in the automotive context (cf. Berger and Birkemeyer, [1]). We leveraged a Z-order SFC to reduce multi-dimensional data into their corresponding single-dimensional representations (vertical stripes) resulting in Characteristic Stripe Patterns (CSPs) as shown in Fig. 3(b). We showed that the temporal occurrence, distribution, and spread of such stripes correlate with certain maneuvers in the multi-dimensional space. The locality preservation property of SFCs in combination with operating just on a single dimension allowed us to design computationally efficient operations to detect maneuvers. We showed in [1] that maneuver detection can be conducted orders of magnitude faster in comparison to directly processing time-series data in the original data space. The key idea behind this performance gain is motivated by the fact that the application context restricts the value range of the underlying SFC, on which operations like inserting or looking up entries can be realized efficiently as a consequence.
In this study, we are using two types of SFCs, namely Z-order and Hilbert curves: A Z-order curve is exhibiting a Z-like pattern when repetitively traversing a multi-dimensional space as shown in Fig. 1(a); we chose the Z-order curve as the mapping from the multi-dimensional space to the single-dimensional representation can be computed efficiently. A Hilbert curve is a recursive traversal through a multi-dimensional data space and its pattern is depicted in Fig. 1(b); we selected the Hilbert curve as it exhibits even better locality preservation (cf. Bader, [3]) for data points that are in close proximity to each other.
III Methodology
| signal | |||||
| GPS position | - | - | - | - | |
| heading, | - | - | |||
| steering wheel angle, | - | - | |||
| lateral acceleration, | |||||
| longitudinal acceleration, | |||||
| vertical acceleration, | |||||
| yaw rate, | |||||
| pitch rate, | |||||
| roll rate, |
The goal of the study is to compare the performance (in terms of precision, accuracy, and recall) of maneuver detectors to find potential passing maneuvers when driving through a roundabout by exploiting multi-dimensional kinematic signals as input for the SFCs. We use a Volvo XC90 from the vehicle laboratory Chalmers Revere that has been instrumented with an Applanix GPS/IMU system to obtain a baseline reference. In addition, on-board vehicle data as listed in Tab. I is accessible via CAN for data logging. The vehicle was used for a manual data collection of about for a trip of in an urban area within Gothenburg, Sweden as shown in Fig. 2. The trip covered in total 17 roundabout scenarios by visiting 7 different roundabouts on the route. The baseline for the experiments as described in the following is defined two-fold: (a) manual annotations were added to extract the time-points when entering a roundabout, information about roundabout layout (ie., number of exits and chosen exit), and the time-points when leaving the roundabout, and (b) by defining geofences around all possible roundabouts along the trip as marked in Fig. 2.
accelerations over time.
The geofences serve as input for an automatic GPS position-based maneuver detector that encodes the GDPR-sensitive vehicle location onto an SFC curve as follows: The two dimensions (latitude, longitude) were converted to single-dimensional Morton codes (cf. Morton, [10] and Bader, [3]) resulting in a positive integer number. This number was then used as sorted index to the corresponding time of a GPS location resulting in a tuple .
The maneuver look-up is then a range query where the bottom/left and top/right GPS location of a geofence are the lower and upper limits (, ) on the single-dimensional space. All Morton indices between such lower and upper limits correspond to potential candidate locations where the vehicle potentially was passing through a roundabout. Finally, all such potential pairs , with between the lower/upper limits as returned from the range query, are converted back to the original tuple (latitude, longitude) and filtered subsequently using a point-in-polygon test for the corresponding geofence to compensate for potential false-positive Morton candidates resulting from sudden “jumps” in the traversal of the Z-order curve through the multi-dimensional space (see Fig. 1(a), right).
In contrast to the aforementioned maneuver detector, which is using a signal that can be considered potentially being GDPR-sensitive, we setup a systematic experiment for different maneuver detectors based on the two variants (a) Z-order curve-based SFC using Morton codes (abbreviated as M in Fig. 4), and (b) Hilbert curve (abbreviated as H in Fig. 4) to answer both research questions. Both variants of SFCs convert multi-dimensional data samples into their respective single-dimensional representations similar to the aforementioned description for the geofence-based detector. Specific CSP as shown in Fig. 3(b), in which the stripes occur in dedicated ranges over time, correlate with specific automotive maneuvers in the original multi-dimensional space. We use a CSP from a known reference maneuver when passing through a roundabout () to obtain the distribution of SFC values as shown in Fig. 3. To detect subsequent passing maneuvers through roundabouts, a dataset’s SFC values are traversed to check whether respective CSP ranges and value distributions in the stripes are within the thresholds extracted from the reference passing maneuver that we used for calibrating the maneuver detector.
We systematically created permutations as listed in Fig. 4 for our experimentation using the following configuration parameters to evaluate the respective performance in finding drive-through maneuvers at roundabouts:
- 1.
Select Morton or Hilbert for SFC conversion.
- 2.
Select sampling frequency between 5 and .
- 3.
Select combinations of minimum two signals from the set: .
The motivation for the first option is to systematically evaluate what SFC variant is showing a superior performance for detecting a passing maneuver. The second option is to evaluate how much data is essential for a successful maneuver detection; the motivation here is that data minimization is not only reducing processing time and storage usage, but also obeying the data minimization principle according to GDPR. Finally, systematic combinations of various signals is evaluating to what degree certain kinematic parameters are influencing a successful maneuver detection.
IV Results
The results from our systematic experiments are aggregated and presented in Fig. 4. The first column lists the specific configuration for a single experiment, characterized by the chosen parameters: An SFC method for encoding the multi-dimensional values, a combination of kinematic signals from the multi-dimensional space, and a data sampling frequency. Other columns depict each roundabout pass-through in the order of occurrence. Therefore, the x-axis shows the temporal dimension: the length of each box refers to the event’s duration. The top two rows show both baselines: The manually annotated roundabout events and the results obtained from the geofence-based maneuver detector. All subsequent rows show the results from a specific configuration of an SFC-based maneuver detector.
The green boxes represent correctly identified passes through a roundabout (true-positive, TP) and the red boxes indicate when a maneuver detector returned time ranges from a maneuver that it believed was a passing maneuver through a roundabout but it was not (false-positive, FP). The color was assigned to each box by conducting a binary test whether the time range returned by a detector is partially overlapping with the manual annotations. Thus, green indicates an overlap, and red no overlap. In our experiments, we consider a minimal overlap being sufficient for a TP hit. We determine such overlapping time ranges by calculating the intersection-over-union (IoU) of the manually annotated ground truth with the returned time ranges of an experiment. Each finding is set to green (a TP) if it has an IoU value greater than . In this setting, a 1-second-long overlap with a passing maneuver, a 100% overlap with a passing maneuver, or even ten different short detections that overlap with a passing maneuver are all considered TPs, which we refer to as Boolean-identified passes.
Alternatively, successful passes can also be defined as the amount of intersecting seconds, FPs are then defined as the detected seconds without the intersecting seconds, false-negative (FN) is defined as the ground truth time windows without the intersections, and true-negative (TN) is the remaining time (neither passing in the ground truth nor according to a detector). Fig. 4 is ordered using these time-based definitions by the F1-score of the respective experiment because we prefer configurations for maneuver detectors that find as many TPs as possible without being penalized by too many FPs; intuitively speaking, we want as many green ranges for the costs of the least red ranges as possible. This is also reflected by having the manual annotations at the top (only TPs with no FPs), followed by the geofence-based detector that is using GDPR-sensitive signals. It can be noted, though, that while this detector correctly identifies all passing maneuvers as expected, the dimensioning of the geofences may unintentionally reach into nearby streets resulting in FPs (two in our experiment).
Based on the analysis of the results, we choose three exemplary experiment configurations with different success rates and discuss their strengths: (i) Hilbert with , to which we will refer as experiment from now on; (ii) Morton with , here forth ; and (iii) Morton with , here forth . The first two configurations use normalized data.
Out of the three selected experiments, has the greater F1 score, which we consider to indicate a good balance between finding the maximum number of roundabout maneuvers (TPs) and minimizing false detections (FPs). Nevertheless, only detects 6 out of the 17 manually annotated passes and two of the roundabouts (RA1 and RA2) are never detected with it. On the other hand, is able to detect 14 of the passes at the expense of a higher number of false detections (54). In fact, is one of the 31% of experiments to detect pass 2. is also able to detect most of the passes and only missing 5. However, unlike , most of the predictions (63%) using overlap with the manually annotated events.
Out of the 17 passes as depicted in Tab. II, 5 correspond to roundabout maneuvers where the first exit was taken, and 5 to maneuvers where the second exit was taken. These passes are only detected by 15% and 32% of the experiment configurations, respectively. In contrast, maneuvers where the third exit was taken are detected by 60% of the experiment configurations. This shows that the longer an actual passing maneuver through a roundabout takes, the higher is the maximum accuracy from all experiments as well as more detector configurations will spot such a maneuver. However, as depicted in Fig. 4, the exit taken does not necessarily affect the length of the maneuver, yet maneuvers that consist of a full turn (or more) are detected by more than 94% of the experiment configurations, including , , and .
Another aspect to consider are multiple detections of a single roundabout. This can be seen in Fig. 4 as multiple green boxes overlapping with the same pass in the manually annotated data. For instance, pass 15 was detected twice, and pass 17 was detected three times using ; and when using and , each of the 4 passes that went through RA5 was detected twice. This raises the question as discussed in Sec. V about whether a part of the roundabout in SFC space is enough to claim that a roundabout maneuver was successfully detected. Moreover, experiment configurations that break down single maneuvers into multiple detections also have a higher number of false predictions as can be seen for in Fig. 4.
| Pass | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Roundabout | RA1 | RA2 | RA3 | RA3 | RA4 | RA5 | RA4 | RA5 | RA6 | RA6 | RA7 | RA7 | RA7 | RA5 | RA4 | RA5 | RA6 |
| Exit taken | 3/4 | 2/4 | 6/5 | 2/5 | 2/4 | 4/4 | 1/4 | 2/4 | 1/4 | 1/4 | 1/4 | 4/4 | 1/4 | 2/4 | 8/4 | 5/4 | 14/4 |
| Hits (over all) | 60% | 31% | 96% | - | 40% | 81% | 15% | 46% | 38% | - | - | 85% | 25% | 42% | 94% | 96% | 98% |
| Max. overlap (%) | 83% | 83% | 86% | - | 82% | 86% | 69% | 83% | 82% | - | - | 86% | 83% | 83% | 86% | 86% | 86% |
The selected experiments described before therefore trigger three important discussions for Sec. V: (i) the balance of metrics such as precision and recall to determine success, (ii) the trade-off between the data inputs and the amount of correctly identified passes, and (iii) whether fractions of a roundabout (small time-window), can be considered a TP as indicated above.
V Analysis & Discussion
As expected, the automatic geofence-based detector, which exploits the GPS location to detect roundabouts, is scoring perfect results in terms of both precision and recall. However, as discussed before, GPS information may be considered as a GDPR-sensitive information. An alternative detection strategy that does not rely on such information is therefore preferred.
Answering RQ-1, all experiment configurations were able to detect at least one roundabout passing maneuver, relying on non-GDPR-sensitive data only. However, three passes were not detected by any of the experiment configurations.
According to the ordering as shown in Fig. 4, the most successful performing configuration is Morton using at (marked as in Fig. 4) returning 14 out of 17 maneuvers. Nevertheless, this high detection rate often comes at the expense of a low precision due to an increasing number of FPs. Moreover, the length in time of the detections is shorter than the manually annotated passes as the green ranges that are overlapping between the manual annotations and the experiment’s results become smaller. Hence, in order to decide what experiment can be considered more successful over others is requiring a definition of a correctly identified passing maneuver (TP). The first option is Boolean detection, ie., caring about how many intersections with the ground truth are present. Therefore, the actual width of the green ranges is of secondary interest but rather the amount of distinct yet correct TPs is preferred. The second option is caring about which percentage of the manually annotated time windows corresponding to roundabout passes is detected by each of the experiments.
Answering RQ-2, we discussed a number of experiment configurations that are able to detect roundabout passing maneuvers successfully. However, metrics for “success” might depend on the application of the proposed technique.
It is important to note that three passing maneuvers are not detected by any of the configurations listed Fig. 4: Passes 4, 10, and 11 as described in Tab. II. These passes may be particularly challenging: In the former, the second exit was taken, which may be appearing an evasive maneuver; while in the latter the first exit was taken, making it similar to a simple right turn maneuver.
Experiment configurations that do not use steering wheel angle information like and detect on average a higher percentage of the manually annotated maneuvers. At the same time, configurations that do not use steering wheel angle information were only able to detect a maximum of 9 out of the 17 passes. Similarly, configurations that do not use the yaw rate to detect roundabout events also ranked higher yet none of them was able to detect passes 2, 4, 10, or 11. Interestingly, the passes which benefitted the most from not using steering wheel angle or yaw rate, were maneuvers when the roundabout was passed for at least one complete round.
In order to decide what experiment was more successful, the relative importance of TPs, FPs and FNs needs to be determined independently on whether we use Boolean or continuous definitions for correctly identified passes. In the context of Advanced Driver Assistance System / Automated Driving (ADAS/AD), FPs are typically less problematic than FNs, so we ordered the experiments in Fig. 4 based on the F1-score computed based on overlapping time windows. It is therefore important before using our technique to automatically detect maneuvers without relying on GDPR-sensitive data to decide what is preferable: A small intersection with the ground truth might be enough to trigger a roundabout-specific ADAS/AD behavior, ie., making the Boolean approach useful depending on such detections shall be post-processed. In a similar way, a different balance between TP, FP, and FN might be needed depending on the application context otherwise.
We report threats to validity following Feldt and Magazinius [11]. The data used in the experiments detailed in Section 4 correspond to a limited data collection and showcase a selected number of events corresponding to roundabout pass-through maneuvers, which may pose an internal threat to validity. Even though the discussed input signals as reported in Tab. I are typically used for maneuver identification, we cannot claim generalizability of these results over alternative selections of signals or newly collected data at other locations; for example, many roundabouts along the trip included also a road bump to reduce speed, which on the one hand would even allow to improve the identification due to this particularity visible in ; however, this aspect is on the other not necessarily standard in other locations or countries and hence, we excluded this parameter.
VI Conclusion and Future Work
We are presenting a systematic evaluation of SFC to be used for automotive maneuver detection on the example of finding maneuvers for passing a roundabout. SFCs allow for mapping multi-dimensional data into their single-dimensional representations and hence, such data can be processed computationally efficiently. This becomes important when working with ever-growing data lakes collected from large vehicle fleets, or when using automatic maneuver detectors that shall be used in real-time on the vehicle.
Uniquely for our study is that we are interested in the performance of the maneuver detectors when dropping signals that may be GDPR-sensitive such as GPS locations. Changing the sampling rate or the SFC encoding can improve the number of true-positives but at an expense of finding more FPs as well, rendering such detectors impractical. In order to make decisions on which experiment configuration is more successful, application domain details need to be considered as discussed in Sec. IV using three examples with different time overlap and number of detections.
Future work may refine the design of the maneuver detectors further to lower the false-positive rates. In addition, the particularly challenging maneuver when leaving a roundabout already at the first exit in comparison to a plain turning right maneuver needs to be addressed further to better separate different automotive maneuvers when using non-GDPR-sensitive signals. As we conducted our experiments with data from Gothenburg, Sweden, also the transferability of our findings needs to be explored further when looking at data from other regions.
References
- [1] C. Berger and L. Birkemeyer, “ZEBRA: Z-order Curve-based Event Retrieval Approach to Efficiently Explore Automotive Data,” in 2023 IEEE Intelligent Vehicles Symposium (IV), 2023, pp. 1–7. [Online]. Available: https://ieeexplore.ieee.org/document/10186770
- [2] M. A. Perez, J. D. Sudweeks, E. Sears, J. Antin, S. Lee, J. M. Hankey, and T. A. Dingus, “Performance of basic kinematic thresholds in the identification of crash and near-crash events within naturalistic driving data,” Accident Analysis & Prevention, vol. 103, pp. 10–19, 2017. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0001457517301033
- [3] M. Bader, Space-Filling Curves. Springer Berlin, Heidelberg, 2013. [Online]. Available: https://link.springer.com/book/10.1007/978-3-642-31046-1
- [4] A. Hulbert, T. Kunicki, J. N. Hughes, A. D. Fox, and C. N. Eichelberger, “An experimental study of big spatial data systems,” in 2016 IEEE International Conference on Big Data (Big Data). Washington, DC, USA: IEEE, 2016, pp. 2664–2671. [Online]. Available: https://ieeexplore.ieee.org/abstract/document/7840909
- [5] H. Dai and H.-C. Su, “Clustering Analyses of Two-Dimensional Space-Filling Curves: Hilbert and z-Order Curves,” SpringerNature Computer Science, vol. 4, no. 1, p. 8, 2022. [Online]. Available: https://link.springer.com/article/10.1007/s42979-022-01320-9
- [6] B. Moon, H. V. Jagadish, C. Faloutsos, and J. H. Saltz, “Analysis of the clustering properties of the Hilbert space-filling curve,” IEEE Transactions on Knowledge and Data Engineering, vol. 13, no. 1, pp. 124–141, 2001. [Online]. Available: https://ieeexplore.ieee.org/abstract/document/908985
- [7] A. A. Safia, Z. A. Aghbari, and I. Kamel, “Distributed Environmental Event Monitoring using Mobile Wireless Sensor Network,” Procedia Computer Science, vol. 155, pp. 335–342, 2019. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S1877050919309627
- [8] H. Liu, J. Yan, J. Wang, B. Chen, M. Chen, and X. Huang, “HGST: A Hilbert-GeoSOT Spatio-Temporal Meshing and Coding Method for Efficient Spatio-Temporal Range Query on Massive Trajectory Data,” ISPRS International Journal of Geo-Information, vol. 12, no. 3, 2023. [Online]. Available: https://www.mdpi.com/2220-9964/12/3/113
- [9] S. H. Nair, A. Sinha, and L. Vachhani, “Hilbert’s space-filling curve for regions with holes,” in 2017 IEEE 56th Annual Conference on Decision and Control (CDC), 2017, pp. 313–319. [Online]. Available: https://ieeexplore.ieee.org/abstract/document/8263684
- [10] G. M. Morton, “A Computer Oriented Geodetic Data Base and a New Technique in File Sequencing,” 1966. [Online]. Available: https://dominoweb.draco.res.ibm.com/reports/Morton1966.pdf
- [11] R. Feldt and A. Magazinius, “Validity Threats in Empirical Software Engineering Research - An Initial Survey,” in SEKE, 2010. [Online]. Available: https://www.cse.chalmers.se/~feldt/publications/feldt_2010_validity_threats_in_ese_initial_survey.pdf