ECE 150 Fundamentals of programming

In UTM, the world is divided into 100 km × 100 km squares and locally, within these squares, the Earth is approximated as if it were flat. Due to the Earth being a sphere (although the validity of this fact appears to be in dispute among those who would prefer to see the human species as the center of the universe), these squares do not perfectly line up, but for most applications, navigation occurs within one of these squares.

Given two points, a 10-figure grid reference are two pairs of five decimal-digit numbers, the first specifying the number of meters from the bottom left corner of the square in the easterly direction, and the next specifying the numbers of meters from the bottom left corner of the square in the northerly direction.

Aside: to help remember that the first number is the east, and the second number is north, the Canadian Forces use the aide-memoire "In the door, up the stairs."

Given two points within the same 100 km square, it may be useful to know the distance from one point to the other. For a 10-figure grid reference, the distance simply uses the Pythagorean theorem.

For example, the distance between 0523941178 and 0768243596 is

$\sqrt{(05239 - 07682)^2 + (41178 - 43596)^2} = 3437.291520950762$

or approximately 3.44 km. You can see this in Figure 1.

Figure 1. The distance between the points 0523941178 and 0768243596.

Normally, grid references are written as a 10-digit number, but we will use an array of capacity two where the first number is the distance east, and the second is the distance north.

Write a function that takes two pairs of UTM grid references and returns the distance from the first to the second:

double distance_utm( unsigned int source[2], unsigned int destination[2] );