# Geometric Median

In normal median, we find a point that has minimum sum of distances. Similar concept applies in 2-D space.

Given **N** points in **2-D space**, the task is to find out a single point **(x, y)** from which the **sum of distances** to the input points are **minimized** (also known as the center of minimum distance).

**Examples:**

Input:(1, 1), (3, 3)

Output:Geometric Median = (2, 2) with minimum distance = 2.82843

Input:(0, 0), (0, 0), (0, 12)

Output:Geometric Median = (0, 0) with minimum distance = 12

**Approach:**

At first thought, it seems that the problem asks us to find the Mid point or the Geometric Center point (in other words, centroid) of the given input points. Since it is the “center” point of the input, sum of distances from the center to all the given input points should automatically be minimized. This process is analogous to finding the Center of Gravity of discrete Mass particles. The first example test case even gives the correct answer. But what happens when we apply the same logic to the second example?

We can clearly see that the Geometric Center, or the Centroid of is at . So according to the Euclidean Distance formula, the total distance to travel from Centroid to all 3 of the input points is But the optimal point should be , giving us a total distance of So, where are we wrong?

Intuitively, you can think that Centroid of input points gives us the Arithmetic Mean of the input points. But what we require is the Central Tendency of the input points such that the cost to reach that central tendency (or in other words, the Euclidean Distance) is minimized. This is called the Geometric Median of a set of points.It is kind of like how conceptually, a Median is drastically different from Mean of given inputs.

There isn’t any defined correct algorithm for finding the Geometric Median. What we do to approach this kind of problems is approximating a solution and determining whether our solution is indeed the Geometric Median or not.

**Algorithm **

There are two important variables :

- current_point – stores the x and y coordinates of the point which could be the Geometric Median.
- minimum_distance – stores the sum of Euclidean distances from current_point to all input points.

After every approximation, if we find a new point from which the sum of distances is lower, then we update both the values of current_point and minimum_distance to the new point and new distance.

First, we find the Centroid of the given points, take it as the current_point (or the median) and store the sum of distances in minimum_distance. Then, we iterate over the given input points, by turn assuming each input point to be the median, and then calculating the distance to other points. If this distance is lower than the minimum_distance, then we update the old values of current_point and minimum_distance to the new values. Else, the old values remains the same.

Then we enter a while loop. Inside that loop, we move a distance of test_distance (we assume a test_distance of 1000 for this example) from the current_point in all directions (left, up, right, down). Hence we get new points. Then we calculate the distance from these new points to the given input points. If this sum of distances is lower than the previous minimum_distance then we update the old values of current_point and minimum_distance to the new values and repeat the while loop. Else, we divide the test_distance by and then repeat the while loop.

The terminating condition for the while loop is a certain value called the “lower_limit”. Lower the value, higher the accuracy of our approximation. Loop terminates when lower_limit exceeds the test_distance.

Below is the implementation of the above approach:

`// C++ implementation of the approach ` `#include <bits/stdc++.h> ` `using` `namespace` `std; ` ` ` `// To store a point in 2-D space ` `struct` `Point { ` ` ` `double` `x, y; ` `}; ` ` ` `// Test points. These points are the left, ` `// up, right and down relative neighbours ` `// (arranged circularly) to the ` `// current_point at a distance of ` `// test_distance from current_point ` `Point test_point[] = { { -1.0, 0.0 }, ` ` ` `{ 0.0, 1.0 }, ` ` ` `{ 1.0, 0.0 }, ` ` ` `{ 0.0, -1.0 } }; ` ` ` `// Lowest Limit till which we are going ` `// to run the main while loop ` `// Lower the Limit higher the accuracy ` `double` `lower_limit = 0.01; ` ` ` `// Function to return the sum of Euclidean ` `// Distances ` `double` `distSum(Point p, ` ` ` `Point arr[], ` `int` `n) ` `{ ` ` ` `double` `sum = 0; ` ` ` `for` `(` `int` `i = 0; i < n; i++) { ` ` ` `double` `distx = ` `abs` `(arr[i].x - p.x); ` ` ` `double` `disty = ` `abs` `(arr[i].y - p.y); ` ` ` `sum += ` `sqrt` `((distx * distx) + (disty * disty)); ` ` ` `} ` ` ` ` ` `// Return the sum of Euclidean Distances ` ` ` `return` `sum; ` `} ` ` ` `// Function to calculate the required ` `// geometric median ` `void` `geometricMedian(Point arr[], ` `int` `n) ` `{ ` ` ` ` ` `// Current x coordinate and y coordinate ` ` ` `Point current_point; ` ` ` ` ` `for` `(` `int` `i = 0; i < n; i++) { ` ` ` `current_point.x += arr[i].x; ` ` ` `current_point.y += arr[i].y; ` ` ` `} ` ` ` ` ` `// Here current_point becomes the ` ` ` `// Geographic MidPoint ` ` ` `// Or Center of Gravity of equal ` ` ` `// discrete mass distributions ` ` ` `current_point.x /= n; ` ` ` `current_point.y /= n; ` ` ` ` ` `// minimum_distance becomes sum of ` ` ` `// all distances from MidPoint to ` ` ` `// all given points ` ` ` `double` `minimum_distance = ` ` ` `distSum(current_point, arr, n); ` ` ` ` ` `int` `k = 0; ` ` ` `while` `(k < n) { ` ` ` `for` `(` `int` `i = 0; i < n, i != k; i++) { ` ` ` `Point newpoint; ` ` ` `newpoint.x = arr[i].x; ` ` ` `newpoint.y = arr[i].y; ` ` ` `double` `newd = ` ` ` `distSum(newpoint, arr, n); ` ` ` `if` `(newd < minimum_distance) { ` ` ` `minimum_distance = newd; ` ` ` `current_point.x = newpoint.x; ` ` ` `current_point.y = newpoint.y; ` ` ` `} ` ` ` `} ` ` ` `k++; ` ` ` `} ` ` ` ` ` `// Assume test_distance to be 1000 ` ` ` `double` `test_distance = 1000; ` ` ` `int` `flag = 0; ` ` ` ` ` `// Test loop for approximation starts here ` ` ` `while` `(test_distance > lower_limit) { ` ` ` ` ` `flag = 0; ` ` ` ` ` `// Loop for iterating over all 4 neighbours ` ` ` `for` `(` `int` `i = 0; i < 4; i++) { ` ` ` ` ` `// Finding Neighbours done ` ` ` `Point newpoint; ` ` ` `newpoint.x = current_point.x ` ` ` `+ (` `double` `)test_distance * test_point[i].x; ` ` ` `newpoint.y = current_point.y ` ` ` `+ (` `double` `)test_distance * test_point[i].y; ` ` ` ` ` `// New sum of Euclidean distances ` ` ` `// from the neighbor to the given ` ` ` `// data points ` ` ` `double` `newd = distSum(newpoint, arr, n); ` ` ` ` ` `if` `(newd < minimum_distance) { ` ` ` ` ` `// Approximating and changing ` ` ` `// current_point ` ` ` `minimum_distance = newd; ` ` ` `current_point.x = newpoint.x; ` ` ` `current_point.y = newpoint.y; ` ` ` `flag = 1; ` ` ` `break` `; ` ` ` `} ` ` ` `} ` ` ` ` ` `// This means none of the 4 neighbours ` ` ` `// has the new minimum distance, hence ` ` ` `// we divide by 2 and reiterate while ` ` ` `// loop for better approximation ` ` ` `if` `(flag == 0) ` ` ` `test_distance /= 2; ` ` ` `} ` ` ` ` ` `cout << ` `"Geometric Median = ("` ` ` `<< current_point.x << ` `", "` ` ` `<< current_point.y << ` `")"` `; ` ` ` `cout << ` `" with minimum distance = "` ` ` `<< minimum_distance; ` `} ` ` ` `// Driver code ` `int` `main() ` `{ ` ` ` ` ` `int` `n = 2; ` ` ` `Point arr[n]; ` ` ` `arr[0].x = 1; ` ` ` `arr[0].y = 1; ` ` ` `arr[1].x = 3; ` ` ` `arr[1].y = 3; ` ` ` `geometricMedian(arr, n); ` ` ` ` ` `return` `0; ` `} ` |

*chevron_right*

*filter_none*

**Output:**

Geometric Median = (2, 2) with minimum distance = 2.82843

References: Geometric Median, Center of minimum distance

## Recommended Posts:

- Number of terms in Geometric Series with given conditions
- Draw geometric shapes on images using OpenCV
- Median
- Mean and Median of a matrix
- Maximize the median of an array
- Median after K additional integers
- Median Of Running Stream of Numbers - (using Set)
- Program for Mean and median of an unsorted array
- Median and Mode using Counting Sort
- Find median of BST in O(n) time and O(1) space
- Median of two sorted arrays with different sizes in O(log(min(n, m)))
- Median of two sorted arrays of same size
- Randomized Algorithms | Set 3 (1/2 Approximate Median)
- Median of two sorted arrays of different sizes
- Finding Median in a Sorted Linked List

If you like GeeksforGeeks and would like to contribute, you can also write an article using contribute.geeksforgeeks.org or mail your article to contribute@geeksforgeeks.org. See your article appearing on the GeeksforGeeks main page and help other Geeks.

Please Improve this article if you find anything incorrect by clicking on the "Improve Article" button below.