Friday 29 January 2016

distance - Finding most remote spot in Eastern United States?


I have this hunch that I could travel no more than 50 miles from the most remote spot in the Eastern United States (east of the Mississippi River), in the direction of the nearest road, and find a road.


Definitions:
Most Remote: Spot furthest from a road.
Road: Google Maps definition of a road.


How could I prove or disprove this claim (i.e where is the most remote spot in the Eastern US)?



Answer




A fast and informative way is to create a distance grid based on the roads. This is usually done in a projected coordinate system, which necessarily introduces some error, but by choosing a good coordinate system the error will not be too great (and can be corrected).


The following example defines a "road" as a US Interstate or US or state highway of comparable magnitude. These roads are shown as red polylines. It uses a Lambert Conformal Conic projection. Although its metric distortion can readily be corrected in terms of latitude, that's not really necessary in this example because the distortion is less than 0.6% except in Florida, where it grows to 2.3%: good enough for this illustration.


Road distance map


The distances are color coded from dark cyan (short) through yellow (long) and hillshaded to emphasize the local maxima. A glance shows the greatest distances are attained in central Wisconsin and the North Carolina coast. The GIS tells me the maximum distances attained are 194 km and 180 km, respectively. (The maximum attained in Michigan is 120 km, less even than the maximum in central Mississippi, 137 km.)


Using any raster GIS (such as ArcGIS, GRASS, Manifold, etc.) one can perform a similar computation using any roads layer desired (such as Census TIGER streets features). Straightforward post-processing will find all local maxima of the distance grid (seen as peaks on this map), thereby identifying all points that locally are as far from a road as you can get. Very simple post-processing will identify all points exceeding a distance threshold such as 50 miles (about 80 km).


A variant uses a "costdistance" calculation, instead of Euclidean distance (as a proxy for spherical distance), to determine points that are (say) a maximum travel time from the nearest road. This is not an onerous task: typical computation times are a few seconds (at most) at the 1 km resolution used here.


No comments:

Post a Comment

arcpy - Changing output name when exporting data driven pages to JPG?

Is there a way to save the output JPG, changing the output file name to the page name, instead of page number? I mean changing the script fo...