Skip to main content

plotting - How do I make framed plots the same size?


Suppose that I have two plots, plot1 and plot2. The plots have different ranges and different axes. Here is a fictitious example that makes no scientific sense, but demonstrates my issue with the size of framed plots:


imgSize = 475;

plot1 = Plot[x^2, {x, 0, 1}, Frame -> True,
FrameLabel -> {"x (nm)", "y (nm)"},
BaseStyle -> {FontFamily -> "Arial", 20}, ImageSize -> imgSize];

plot2 = Plot[x^2, {x, -10, 10}, PlotRange -> {-10, 1000},

Frame -> True,
FrameLabel -> {"\!\(\*SuperscriptBox[SubscriptBox[\"x\", \"a\"], \
\"2\"]\) (\!\(\*SuperscriptBox[\"nm\", \"2\"]\))",
"\!\(\*SuperscriptBox[SubscriptBox[\"y\", \"b\"], \"2\"]\) (\!\(\
\*SuperscriptBox[\"nm\", \"2\"]\))"},
BaseStyle -> {FontFamily -> "Arial", 20}, ImageSize -> imgSize];

Grid[{{plot1, plot2}}]

I get this output:



plotsize


Notice that plot1 and plot2 are not really the same size -- in terms of the size of the outer frame of each. In particular, the outer frame of plot1 has greater height than that of plot2. In addition, I think that the outer frame of plot1 has greater width than that of plot2. I think this is because of the different ranges and different axis labels of the two plots.


Observed separately, one probably could not discern a difference in size between plot1 and plot2. But when they are next to each other, as in a Grid, plot1 looks noticeably "larger" than plot2. This would look rather poor in a publication, like an article in a scientific journal.


Is there any way that I can make the outer frames the same size, such that the plots look like they are the same size?



Answer



As Jagra said, the usual solution is to manually specify the ImagePadding values. The problem is that if the plot or frame labels are too large, a fixed ImagePadding may cut them off.


But can we automate this?


Ideally what we would do is:





  • Create the two plots of fixed vertical size and retrieve their vertical ImagePadding




  • Change the ImagePadding of both to use the larger value. This will ensure that they have the same size while no labels are cut off.




So how do we measure the ImagePadding of an existing plot? This is unfortunately tricky as the value depends on the size of the plot. Since the plots will be stacked horizontally, we need to fix the vertical image size before trying to retrieve the padding. But here's a pretty useful solution (that I use regularly):



First note that I'm fixing the vertical size instead of the horizontal one. This allows the graphics to have differing horizontal sizes if necessary while still aligning perfectly when stacked in a row.


verticalSize = 250;


plot1 = Plot[x^2, {x, 0, 1}, Frame -> True,
FrameLabel -> {"x (nm)", "y (nm)"},
BaseStyle -> {FontFamily -> "Arial", 20},
ImageSize -> {Automatic, verticalSize}];

plot2 = Plot[x^2, {x, -10, 10}, PlotRange -> {-10, 1000},
Frame -> True,
FrameLabel -> {"\!\(\*SuperscriptBox[SubscriptBox[\"x\", \"a\"], \
\"2\"]\) (\!\(\*SuperscriptBox[\"nm\", \"2\"]\))",

"\!\(\*SuperscriptBox[SubscriptBox[\"y\", \"b\"], \"2\"]\) \
(\!\(\*SuperscriptBox[\"nm\", \"2\"]\))"},
BaseStyle -> {FontFamily -> "Arial", 20},
ImageSize -> {Automatic, verticalSize}];

This function measures the padding (based on @Heike's code):


getPadding[g_] := Module[{im},
im = Image[Show[g, LabelStyle -> White, Background -> White]];
BorderDimensions[im]
]


Now let's choose the larger one of both the top and bottom paddings of the two figures. This will give us the minimum image padding that still does not cut off labels.


{p1h, p1v} = getPadding[plot1];

{p2h, p2v} = getPadding[plot2];

verticalPadding = Max /@ Transpose[{p1v, p2v}]

Row[{
Show[plot1, ImagePadding -> {p1h, verticalPadding}],

Show[plot2, ImagePadding -> {p2h, verticalPadding}]
}]

Mathematica graphics


The problem with this approach is that often one would wish to fix the horizontal size of the whole graphic (to fit the text width of the document). I admit that when I had the same problem I did this by iteratively refining the sizes, which is not very elegant, but produces good results automatically.


Comments

Popular posts from this blog

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 - 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 - Adding a thick curve to a regionplot

Suppose we have the following simple RegionPlot: f[x_] := 1 - x^2 g[x_] := 1 - 0.5 x^2 RegionPlot[{y < f[x], f[x] < y < g[x], y > g[x]}, {x, 0, 2}, {y, 0, 2}] Now I'm trying to change the curve defined by $y=g[x]$ into a thick black curve, while leaving all other boundaries in the plot unchanged. I've tried adding the region $y=g[x]$ and playing with the plotstyle, which didn't work, and I've tried BoundaryStyle, which changed all the boundaries in the plot. Now I'm kinda out of ideas... Any help would be appreciated! Answer With f[x_] := 1 - x^2 g[x_] := 1 - 0.5 x^2 You can use Epilog to add the thick line: RegionPlot[{y < f[x], f[x] < y < g[x], y > g[x]}, {x, 0, 2}, {y, 0, 2}, PlotPoints -> 50, Epilog -> (Plot[g[x], {x, 0, 2}, PlotStyle -> {Black, Thick}][[1]]), PlotStyle -> {Directive[Yellow, Opacity[0.4]], Directive[Pink, Opacity[0.4]],