Skip to main content

geographics - Generating different GeoPaths for multiple equivalent segments


Is there an efficient and elegant method so that multiple calls of GeoPath for the same segment (path) yield visibly separate segments on a map?


For instance, this code calls two "equivalent" GeoPaths (from China to Brazil):


GeoGraphics[
{Red, Thick,
GeoPath[#] & /@
{{Entity["Country", "China"], Entity["Country", "Brazil"]},
{Entity["Country", "China"], Entity["Country", "Brazil"]}}},

GeoRange -> "World"]

Two GeoPaths overlapping


Because the GeoPaths are the same, the two overlap and appear as a single segment or path.


The same problem arises if one uses specific locations, such as:


{CityData["BuenosAires", "Coordinates"], CityData["Beijing", "Coordinates"]}

I'd like to displace each path slightly, so that the two paths (or in general $n$ paths) appear separate (even if $n-1$ paths are not geodesics). Note that such displacement is done automatically in renderings of edges in Multigraphs. (Background: I want to plot a complicated multi-city itinerary in which several legs of the trip were between the same pair of cities.)


I could do this "by hand," for each such set of overlapping paths (either by displacing, or using dashes and coloring, etc.), but I'm seeking a far more elegant and automatic method.


Ideally, I think Wolfram should implement an option for GeoGraphics such as MultiSegment -> True which draws the $n$ separate paths, even if they are specified in a non-sequential order such as:



{{Entity["Country", "India"], Entity["Country", "China"]}, 
{Entity["Country", "China"], Entity["Country", "Brazil"]},
{Entity["Country", "Brazil"], Entity["Country", "Chile"]},
{Entity["Country", "Chile"], Entity["Country", "Brazil"]},
{Entity["Country", "Brazil"], Entity["Country", "China"]}}

Answer



A completely different answer, this time ignoring Geo methods and focusing on Graph connectivity: just place the Graph object on the map.


Take your entity pairs:


ents = {
{Entity["Country", "India"], Entity["Country", "China"]},

{Entity["Country", "China"], Entity["Country", "Brazil"]},
{Entity["Country", "Brazil"], Entity["Country", "Chile"]},
{Entity["Country", "Chile"], Entity["Country", "Brazil"]},
{Entity["Country", "Brazil"], Entity["Country", "China"]}
};

Draw a map with those entities. This time we use the Mollweide projection and a vector background:


map = GeoGraphics[ents, GeoProjection -> "Mollweide", GeoBackground -> "CountryBorders"];

Extract the actual projection used by GeoGraphics. Recall that GeoGraphics adds parameters to the projection to adapt it to the situation, unless those parameters are explicit:



proj = GeoProjection /. Options[map, GeoProjection]
(* {"Mollweide", "Centering" -> GeoPosition[{-1.07727, 29.5607}]} *)

Now, eliminate duplicate entities and find their respective projected coordinates:


uents = Union[Flatten[ents]]
(* {Entity["Country", "Brazil"], Entity["Country", "Chile"],
Entity["Country", "China"], Entity["Country", "India"]} *)

coords = GeoGridPosition[GeoPosition[uents], proj][[1]]
(* {{-1.31625, -0.19348}, {-1.44548, -0.571304},

{1.04746, 0.66214}, {0.717329, 0.384687}} *)

Construct a Graph object using those coordinates:


graph = Graph[uents, UndirectedEdge @@@ ents, VertexCoordinates -> coords]

And finally show the map and the graph together (Show should be able to merge directly the GeoGraphics and Graph objects, but currently we need to extract the Graphics from GeoGraphics):


Show[map[[1]], graph]

enter image description here


Comments

Popular posts from this blog

plotting - Plot 4D data with color as 4th dimension

I have a list of 4D data (x position, y position, amplitude, wavelength). I want to plot x, y, and amplitude on a 3D plot and have the color of the points correspond to the wavelength. I have seen many examples using functions to define color but my wavelength cannot be expressed by an analytic function. Is there a simple way to do this? Answer Here a another possible way to visualize 4D data: data = Flatten[Table[{x, y, x^2 + y^2, Sin[x - y]}, {x, -Pi, Pi,Pi/10}, {y,-Pi,Pi, Pi/10}], 1]; You can use the function Point along with VertexColors . Now the points are places using the first three elements and the color is determined by the fourth. In this case I used Hue, but you can use whatever you prefer. Graphics3D[ Point[data[[All, 1 ;; 3]], VertexColors -> Hue /@ data[[All, 4]]], Axes -> True, BoxRatios -> {1, 1, 1/GoldenRatio}]

plotting - Filling between two spheres in SphericalPlot3D

Manipulate[ SphericalPlot3D[{1, 2 - n}, {θ, 0, Pi}, {ϕ, 0, 1.5 Pi}, Mesh -> None, PlotPoints -> 15, PlotRange -> {-2.2, 2.2}], {n, 0, 1}] I cant' seem to be able to make a filling between two spheres. I've already tried the obvious Filling -> {1 -> {2}} but Mathematica doesn't seem to like that option. Is there any easy way around this or ... Answer There is no built-in filling in SphericalPlot3D . One option is to use ParametricPlot3D to draw the surfaces between the two shells: Manipulate[ Show[SphericalPlot3D[{1, 2 - n}, {θ, 0, Pi}, {ϕ, 0, 1.5 Pi}, PlotPoints -> 15, PlotRange -> {-2.2, 2.2}], ParametricPlot3D[{ r {Sin[t] Cos[1.5 Pi], Sin[t] Sin[1.5 Pi], Cos[t]}, r {Sin[t] Cos[0 Pi], Sin[t] Sin[0 Pi], Cos[t]}}, {r, 1, 2 - n}, {t, 0, Pi}, PlotStyle -> Yellow, Mesh -> {2, 15}]], {n, 0, 1}]

plotting - Mathematica: 3D plot based on combined 2D graphs

I have several sigmoidal fits to 3 different datasets, with mean fit predictions plus the 95% confidence limits (not symmetrical around the mean) and the actual data. I would now like to show these different 2D plots projected in 3D as in but then using proper perspective. In the link here they give some solutions to combine the plots using isometric perspective, but I would like to use proper 3 point perspective. Any thoughts? Also any way to show the mean points per time point for each series plus or minus the standard error on the mean would be cool too, either using points+vertical bars, or using spheres plus tubes. Below are some test data and the fit function I am using. Note that I am working on a logit(proportion) scale and that the final vertical scale is Log10(percentage). (* some test data *) data = Table[Null, {i, 4}]; data[[1]] = {{1, -5.8}, {2, -5.4}, {3, -0.8}, {4, -0.2}, {5, 4.6}, {1, -6.4}, {2, -5.6}, {3, -0.7}, {4, 0.04}, {5, 1.0}, {1, -6.8}, {2, -4.7}, {3, -1....