Skip to main content

pdf format - Export MatrixPlot graphics to PDF


I am trying to export a plot generated by MatrixPlot to PDF but am having problems with the way the exported PDF is rendered.



  mPlot = MatrixPlot[length, 
  Epilog -> {Cyan, 
    MapIndexed[If[#1 != 0, Text[#1, Reverse[#2 - 1/2]]] &, 
     Reverse[number], {2}]}, Mesh -> True, 
  FrameTicks -> {{Table[i, {i, 9}], 
     MapIndexed[{First[#2], regionCodes[[#1]]} &, 
      Table[i, {i, 9}]]}, {Table[i, {i, 9}], 
     MapIndexed[{First[#2], Rotate[regionCodes[[#1]], 90 Degree]} &, 
      Table[i, {i, 9}]]}}, ColorFunctionScaling -> False, 
  PlotRange -> {Min[length], Max[length]}]


The matrices length (plotted as shades of gray in the figure below), number (shown as the cyan numbers in the figure below), and regionCodes have the following values


length = {{0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 0.079774, 0.0545131, 
   0.0484267, 0.0614324, 0.0627509, 0.0572262, 0.0424639, 0}, {0, 0, 
   0.0600802, 0.100418, 0.0615561, 0.0737833, 0.0732888, 0.0624553, 
   0}, {0, 0, 0, 0.0740108, 0.0559982, 0.0666054, 0.0749018, 
   0.0493431, 0}, {0, 0, 0, 0, 0.0774839, 0.0554172, 0.0564699, 
   0.0491287, 0}, {0, 0, 0, 0, 0, 0.0798436, 0.0643044, 0.0606039, 
   0}, {0, 0, 0, 0, 0, 0, 0.107009, 0.0713126, 0}, {0, 0, 0, 0, 0, 0, 
   0, 0.0951722, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}}

number = {{0, 0, 0, 0, 0, 0, 0, 0, 0}, {0, 37, 53, 21, 778, 1064, 228,
    15, 0}, {0, 0, 27, 40, 761, 1164, 272, 21, 0}, {0, 0, 0, 7, 331, 
   507, 125, 9, 0}, {0, 0, 0, 0, 5935, 12430, 2917, 245, 0}, {0, 0, 0,
    0, 0, 10665, 4323, 393, 0}, {0, 0, 0, 0, 0, 0, 664, 90, 0}, {0, 0,
    0, 0, 0, 0, 0, 7, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0}}
regionCodes = {"New England", "Middle Atlantic", "East North Central",
   "West North Central", "South Atlantic", "East South Central", 
  "West South Central", "Mountain", "Pacific"}

This is length//MatrixForm



enter image description here


And this is number//MatrixForm


enter image description here


The plot looks like this when shown in Mathematica enter image description here


But the exported PDF is rendered as follows (shown in OS X Preview) enter image description here


There is a lot of shading and (and anti-aliasing?) going on (it shouldn't!). The only solution I have come up with is to rasterize the plot before exporting it. It does the job, but I don't consider it very elegant.


Export["regionIBD_MatrixPlot.pdf", Rasterize[mPlot, ImageResolution -> 600], ImageSize -> Full]

The funny thing is that if I export it as TIFF, it renders correctly. Is there a way to export the plot as a vector graphics and not as a bitmap?



Answer




To long for a comment, but looks great on "10.0 for Mac OS X x86 (64-bit) (June 29, 2014)" with AR, Preview and Skim. Shut-down and restart your system and try with an alternative viewer ...


enter image description here


Edit


@Jens response is fantastic, right? The procedure improves the resolution dramatically and even reduced the size.


enter image description here


I was able to test both algorithms on iMac & MacBookRetina (Home) and iMac (work). No difference, still the same results. The PDF's I opened on W7 & W8 with AR without significant differences between Mac and W be seen.


I am also "R" user and therefore use X11 and Cairo. On my Mac the Cairo Engine is the PDF Creator (I guess?). I wanted to share this information with the community.


Maybe this Information leads to a new question with even better answers.


Comments

Post a Comment

Popular posts from this blog

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

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

How to remap graph properties?

Graph objects support both custom properties, which do not have special meanings, and standard properties, which may be used by some functions. When importing from formats such as GraphML, we usually get a result with custom properties. What is the simplest way to remap one property to another, e.g. to remap a custom property to a standard one so it can be used with various functions? Example: Let's get Zachary's karate club network with edge weights and vertex names from here: http://nexus.igraph.org/api/dataset_info?id=1&format=html g = Import[ "http://nexus.igraph.org/api/dataset?id=1&format=GraphML", {"ZIP", "karate.GraphML"}] I can remap "name" to VertexLabels and "weights" to EdgeWeight like this: sp[prop_][g_] := SetProperty[g, prop] g2 = g // sp[EdgeWeight -> (PropertyValue[{g, #}, "weight"] & /@ EdgeList[g])] // sp[VertexLabels -> (# -> PropertyValue[{g, #}, "name"]...
Post a Comment
Read more