User Portlet User Portlet

Francisco Rodríguez
Discussions
You should take advantage of `GeoDistance` listability. In[2]:= level1$dists = GeoDistance[level1, level1, DistanceFunction -> "Center"]; // AbsoluteTiming Out[2]= {0.91751, Null} Also with the default option...
I couldn't find it either, but I thought it was documented, anyway I will see where should it be added.
In `GeoGraphics` the recommended way is using `GeoStyling`, which depeding on what you want as output, you can use `"GeoImage"` or `"Image"`: GeoGraphics[{GeoStyling[{"GeoImage", Rasterize[CountryData["UnitedStates", "Flag"],...
I think the problem comes when the _lcp_ image is resized (because the transparency). Doing the composition in other way lets you have the nice image we all expect: ImageResize[ ImageCompose[ GeoGraphics[{}, ImageSize -> 2000, ...
Just as an idea, something that also considers repetitions in the list of searched items: deleteCasesOnce[list_List, cases_List] := Module[{countq}, countq[_] := 0; Scan[(countq[#] = countq[#] + 1;) &, cases]; ...
Very nice pictures!!! By the way, you can also compute how far is Odesa from the partial eclipse border: GeoDistance[Entity["City", {"Odesa", "Odesa", "Ukraine"}], SolarEclipse[DateObject[{2015, 3, 20}], "PartialPhasePolygon", ...
Even the partial view is very far from Peru! GeoDistance[ SolarEclipse[DateObject[{2015, 3, 20}], "PartialPhasePolygon", EclipseType -> "Total"], Entity["Country", "Peru"]] > Quantity[3627.84, "Kilometers"] And you can...
If you see documentation: > `Abs[z]` is left unevaluated if z is not a numeric quantity. So it has always to be forced with `ComplexExpand`, which also works without any option by the way.
Given that `s1` and `s2` are two matrix with the same dimensions, you can do it in many ways, just the two first I can think of: Join[s1, s2] ArrayFlatten[{{s1}, {s2}}] If you want to join `s2` at the right of `s2` instead of below: ...
Give it a try now, Wolfram|Alpha has updated its data and now your examples work as expected: regions = CountryData["Guinea", "AdministrativeDivisions"] prefectures = AdministrativeDivisionData[#, "Subdivisions"] & /@ regions ...