Skip to main content

plotting - How to embed graphic-element-specific metadata within a Graphics object?


One of the most crucial requirements for rich, interactive scientific graphics is being able to "annotate" individual graphic elements (e.g. the data points of a scatterplot, the individual curves in a plot of multiple curve fits to data, individual subregions of a density plot) with additional information. I'm using the term "annotate" very broadly here to stand for "associate additional information with". Minimally, one would like to be able to assign unique identifiers to each feature of interest, that can be used as keys in auxiliary data structures.


(For example, one way to make a scatterplot interactive is to equip each data point in the scatterplot with a tooltip that will display additional information about that point when the user hovers the cursor over the point. Another way would be to give the user the ability to "light up" a subset of the data points based on some shared metadata value.)


Can the Mathematica Graphics object accommodate such annotations within it?


(NB: I don't doubt that it would be possible to implement this sort of capability "on top of" Mathematica. At the moment I am interested only in whether such functionality is already built into Mathematica.)


EDIT: added one more example of possible modes of interactivity.



Answer



I suggest the dictionary data structure, to create semantic and maintainable code.



The most common approach that I've seen on here for handling related data is to stick it all in a list, so people end up with code like this:


{city[[54, 1]], Tooltip[Disk[Reverse@city[[54, 2]], 0.1], city[[54, 3]] <> " (" <> ToString@city[[54, 4]] <> ")"]}

Which is utterly nonsense to the rest of us, who haven't got a clue about what those indices represents. With the dictionary data structure on the other hand we can associate with each object several key/value pairs. For example, to get the population of a certain city we may write:


city["Goteborg"]["Population"]

This is a lot easier to remember than if we had some index for the city and some index for population:


city[[45,3]]

To implement the dictionary data structure, all you have to do is to load the package that I've posted previously here.



<< "~/Documents/Mathematica/datadictionary.m"

Here's an example I've created. It uses tooltips to display the name of a city and its population on mouse over, and it colors each city according to which region it belongs to. It also uses MouseAnnotation to color the cities in the same region as the city you're hovering over white. So it's meant to both exemplify how one can do the things mentioned in the post, and at the same time the dictionary data structure.


getKeys[symbol_] := DownValues[symbol][[All, 1, 1, 1]];

colors = MapIndexed[# -> ColorData[3][First@#2] &, CountryData["Sweden", "Regions"]];

addCity[{name_, region_, country_}] := (
  city[name] = makeDictionary[];
  dictStore[city[name], "Name", name];

  dictStore[city[name], "Region", region];
  dictStore[city[name], "Coordinates", CityData[{name, country}, "Coordinates"]];
  dictStore[city[name], "Population", CityData[{name, country}, "Population"]];
  dictStore[city[name], "Color", region /. colors]
  )

addCity /@ CityData[{All, "Sweden"}];

Deploy@Dynamic@Graphics[{
    RGBColor[0.896`, 0.8878`, 0.8548`], EdgeForm[GrayLevel[0]],

    CountryData["Sweden", "FullPolygon"], {
       If[
        MouseAnnotation[] === city[#]["Region"],
        White,
        city[#]["Color"]],
       Annotation[
        Tooltip[Disk[Reverse@city[#]["Coordinates"], 0.1], 
         city[#]["Name"] <> " (" <> ToString@city[#]["Population"] <> 
          ")"]
        , city[#]["Region"], "Mouse"]

       } & /@ getKeys[city]
    }, AspectRatio -> Full]

map


As you may have noticed, the first example of how horribly unreadable code with related data may become if you keep the data in a list is just a converted version of the code I'm actually using to draw the disks. Only with the dictionary data structure it looks like this:


{city["Goteborg"]["Color"], Tooltip[Disk[Reverse@city["Goteborg"]["Coordinates"], 0.1], city["Goteborg"]["Name"] <> " (" <> ToString@city["Goteborg"]["Population"] <> ")"]}

Which is much more intuitive.


Comments

Post a Comment

Popular posts from this blog

mathematical optimization - Minimizing using indices, error: Part::pkspec1: The expression cannot be used as a part specification

I want to use Minimize where the variables to minimize are indices pointing into an array. Here a MWE that hopefully shows what my problem is. vars = u@# & /@ Range[3]; cons = Flatten@ { Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; Minimize[{Total@((vec1[[#]] - vec2[[u[#]]])^2 & /@ Range[1, 3]), cons}, vars, Integers] The error I get: Part::pkspec1: The expression u[1] cannot be used as a part specification. >> Answer Ok, it seems that one can get around Mathematica trying to evaluate vec2[[u[1]]] too early by using the function Indexed[vec2,u[1]] . The working MWE would then look like the following: vars = u@# & /@ Range[3]; cons = Flatten@{ Table[(u[j] != #) & /@ vars[[j + 1 ;; -1]], {j, 1, 3 - 1}], 1 vec1 = {1, 2, 3}; vec2 = {1, 2, 3}; NMinimize[ {Total@((vec1[[#]] - Indexed[vec2, u[#]])^2 & /@ R...
Post a Comment
Read more

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

What is and isn't a valid variable specification for Manipulate?

I have an expression whose terms have arguments (representing subscripts), like this: myExpr = A[0] + V[1,T] I would like to put it inside a Manipulate to see its value as I move around the parameters. (The goal is eventually to plot it wrt one of the variables inside.) However, Mathematica complains when I set V[1,T] as a manipulated variable: Manipulate[Evaluate[myExpr], {A[0], 0, 1}, {V[1, T], 0, 1}] (*Manipulate::vsform: Manipulate argument {V[1,T],0,1} does not have the correct form for a variable specification. >> *) As a workaround, if I get rid of the symbol T inside the argument, it works fine: Manipulate[ Evaluate[myExpr /. T -> 15], {A[0], 0, 1}, {V[1, 15], 0, 1}] Why this behavior? Can anyone point me to the documentation that says what counts as a valid variable? And is there a way to get Manpiulate to accept an expression with a symbolic argument as a variable? Investigations I've done so far: I tried using variableQ from this answer , but it says V[1...
Post a Comment
Read more