Utilice un vuelo Levy para encontrar Stonehenge en Inglaterra 

Utilice los llamados vuelos Levy, los cuales son óptimos para búsquedas, para encontrar Stonehenge en Inglaterra.

Obtenga el polígono de Inglaterra y las coordenadas de Stonehenge.

In[1]:=

X polygon = Entity["AdministrativeDivision", {"England", "UnitedKingdom"}][ EntityProperty["AdministrativeDivision", "Polygon"]]; stonehenge = Entity["Building", "Stonehenge"][ EntityProperty["Building", "Position"]];

Asumiendo que Stonehenge es visible desde una distancia de 2 km, defina un predicado que compruebe si Stonehenge se puede ver desde un lugar determinado.

In[2]:=

X stonehengeVisibleQ[start_, {dist_, dir_}] := With[{m = Floor[dist/1000]}, Min[(GeoDistance[GeoPosition[#1], stonehenge] &) /@ GeoDestination[start, GeoDisplacement[{(Range[0, m] dist)/m, dir}]][[1]]] < Quantity[2, "Kilometers"]]

Empezando desde el Big Ben, realice un trayecto Levy hasta que Stonehenge se pueda ver.

In[3]:=

X searchpath = {Entity["Building", "ClockTower"][ EntityProperty["Building", "Position"]]}; mr = DiscretizeGraphics[Polygon[Reverse /@ polygon[[1, 1, -1]]]]; While[ {d, \[Alpha]} = { 1000 RandomVariate[StableDistribution[1, 1, 1, 10, 2]], RandomReal[{0, 360}]}; gd = GeoDestination[Last[searchpath], GeoDisplacement[{d, \[Alpha]}]]; If[sd = SignedRegionDistance[mr, Reverse[First[ gd = GeoDestination[Last[searchpath], GeoDisplacement[{d, \[Alpha]}]] ]]] < 0, ! stonehengeVisibleQ[Last[searchpath], {d, \[Alpha]}], True], If[sd, AppendTo[searchpath, gd]]]

Muestre la ruta de búsqueda.

In[4]:=

X GeoGraphics[{GeoStyling["StreetMap"], polygon , PointSize[Medium], RGBColor[1, 0, 0], Point[stonehenge], RGBColor[ 0.5, 0, 0.5], Thickness[Small], GeoPath[searchpath]}]

Out[4]=