Compare mapas históricos y contemporáneos  

Compare el plan original L'Enfant de Washington D.C. de 1792 con un mapa de calles actual.

muestre la entrada completa de Wolfram Languageoculte la entrada

In[1]:=

X DC1792 = Import[ "http://upload.wikimedia.org/wikipedia/en/9/98/Ellicott_plan.jpg"];

In[2]:=

X DC1792 = ColorConvert[Sharpen[ ImageTake[ImageTake[DC1792, {28, 600}], {20, 600}], 20], "Grayscale"];

In[3]:=

X DCToday = GeoGraphics[ Entity["City", {"Washington", "DistrictOfColumbia", "UnitedStates"}], GeoProjection -> "Equirectangular", GeoRange -> {{38.85, 38.93}, {-77.08, -76.96}}, GeoRangePadding -> None, GeoZoomLevel -> 14];

In[4]:=

X landmarks = {{GeoPosition[{38.8843, -76.9773}], {0.865707, 0.426307}}, {GeoPosition[{38.9172, -77.0367}], {0.373301, 0.92497}}, {GeoPosition[{38.8989, -77.0538}], {0.227818, 0.650244}}};

In[5]:=

X map[{x_, y_}] = {ax x + ay y + a , bx x + by y + b} /. Solve[Flatten[ {{ax, ay}.#2 + a == #1[[1, 2]], {bx, by}.#2 + b == #1[[1, 1]]} & @@@ landmarks ], {a, ax, ay, b, bx, by}][[1]];

In[6]:=

X overlayGraphics[{{xmin_, ymin_}, {xmax_, ymax_}}, color_, opacity_] := Graphics[Raster[Reverse@ImageData@ With[{colim = Image[ConstantArray[1. List @@ ColorConvert[color, "RGB"], Reverse@ImageDimensions[DC1792]]]}, ColorCombine[ {colim, ImageApply[opacity (1 - #)^(1/3) &, DC1792]}, "RGB"] ], {{xmin, ymin}, {xmax, ymax}}], Frame -> True]

In[7]:=

X Show[{First[DCToday], overlayGraphics[{map[{0, 0}], map[{1, 1}]}, Purple, 0.8]}]

Out[7]=