Skip to main content

plotting - More on combining 2d and 3d plots


This is a spin off/follow up to Can 2D and 3D plots be combined so that the 2D plot is the bottom surface of the 3D plot boundary? ...


Following Sjoerd's method, I did this:


image = Plot[x^2, {x, -2, 2}, PlotStyle -> {Thick, Black}]; 
Show[Graphics3D[{EdgeForm[], Texture[image], 
Polygon[{{-1, -1, -1}, {1, -1, -1}, {1, 1, -1}, {-1, 1, -1}}, 
VertexTextureCoordinates -> {{0, 0}, {1, 0}, {1, 1}, {0, 1}}]}, 

Lighting -> "Neutral", BoxRatios -> {1, 1, .1}], Boxed -> True]

Mathematica graphics


Two questions:


(1) When I click on the cell bracket and save as... PDF, the resulting graphic is either empty or only contains a snippet of the intended graphic.


(2) I would like to have Mathematica draw the solid which has the region bounded by the parabola and the x axis as its base, but with cross-sections (perpendicular to the base) to be squares. I am not sure how to go about this using Graphics3D.


Any ideas?


Edit: RegionPlot3D indeed allows me to answer (2) very quickly:


p1 = Plot[x^2, {x, -2, 2}, PlotStyle -> {Thick, Black}];
p2 = RegionPlot3D[-Sqrt[y] <= x <= Sqrt[y] && 

 0 <= z <= 2 Sqrt[y], {x, -2, 2}, {y, 0, 4}, {z, 0, 4}, 
 Mesh -> False, PlotPoints -> 100, PlotStyle -> Opacity[.3], 
 Boxed -> False, Axes -> False]
p3 = Show[
 Graphics3D[{EdgeForm[], Texture[p1], 
 Polygon[{{-1, -1, -1}, {1, -1, -1}, {1, 1, -1}, {-1, 1, -1}}, 
 VertexTextureCoordinates -> {{0, 0}, {1, 0}, {1, 1}, {0, 1}}]}, 
 Lighting -> "Neutral", BoxRatios -> {1, 1, 1}], Boxed -> True]

Now my only issue is how to lay the solid in p2 over the image in p3? I tried



Show[p2, p3, PlotRange -> All]

and other variations, to no avail.


Edit #2: Should have had RegionPlot z range as {z,0,4}. Fixed now. Also, the slices can be visualized via


slices = Table[RegionPlot3D[-Sqrt[y] <= x <= Sqrt[y] && 0 <= z <= 2 Sqrt[y] && 
 j - .5 <= y <= j, {x, -2, 2}, {y, 0, 4}, {z, 0, 4}, 
 Mesh -> False, PlotPoints -> 100, PlotStyle -> Opacity[.3], 
 Boxed -> False, Axes -> False], {j, .5, 4, .5}];
Show[Table[slices[[j]], {j, 1, 8}]]




Answer



I'm guessing you want something like this:


p1 = Plot[x^2, {x, -2, 2}, PlotStyle -> {Thick, Black}];
p2 = RegionPlot3D[-Sqrt[y] <= x <= Sqrt[y] && 
    0 <= z <= 2 Sqrt[y], {x, -2, 2}, {y, 0, 4}, {z, 0, 2}, 
   Mesh -> False, PlotPoints -> 100, PlotStyle -> Opacity[.3], 
   Boxed -> False, Axes -> False];
p3 = Graphics3D[{EdgeForm[], Texture[p1], 
    Polygon[{{-2, 0, -1}, {2, 0, -1}, {2, 4, -1}, {-2, 4, -1}}, 

     VertexTextureCoordinates -> {{0, 0}, {1, 0}, {1, 1}, {0, 1}}]}];

Graphics3D[{p2[[1]], p3[[1]]}, Lighting -> "Neutral", 
 BoxRatios -> {1, 1, 1}, Boxed -> True]

solid


Here I modified your p3 to make the polygon bigger so its boundaries are described by {{-2, 0, -1}, {2, 0, -1}, {2, 4, -1}, {-2, 4, -1}}. The texture on the polygon is of course not perfectly aligned with the 3D RegionPlot because the 2D plot contains additional margins for the axes and labels. But I'm using what you provided in the question.


Finally, the main point is that in order to combine the two Graphics3D objects p2 and p3, I used the fact that the actual 3D graphics objects are always stored as the first element of the Graphics3D, so one can extract them by doing p2[[1]] and p3[[1]]. The results can then be re-packaged into a new Graphics3D as I do on the last line, adding the lighting and frame options that you chose for p3 in your example.


Note


If you do this kind of combination more frequently, you may also want to look into the closely related answer I posted under “Covering up” text in Graphics where I define a function label3D which can be used to place arbitrary 2D objects into a 3D scene, including transparency effects.



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

How to thread a list

I have data in format data = {{a1, a2}, {b1, b2}, {c1, c2}, {d1, d2}} Tableform: I want to thread it to : tdata = {{{a1, b1}, {a2, b2}}, {{a1, c1}, {a2, c2}}, {{a1, d1}, {a2, d2}}} Tableform: And I would like to do better then pseudofunction[n_] := Transpose[{data2[[1]], data2[[n]]}]; SetAttributes[pseudofunction, Listable]; Range[2, 4] // pseudofunction Here is my benchmark data, where data3 is normal sample of real data. data3 = Drop[ExcelWorkBook[[Column1 ;; Column4]], None, 1]; data2 = {a #, b #, c #, d #} & /@ Range[1, 10^5]; data = RandomReal[{0, 1}, {10^6, 4}]; Here is my benchmark code kptnw[list_] := Transpose[{Table[First@#, {Length@# - 1}], Rest@#}, {3, 1, 2}] &@list kptnw2[list_] := Transpose[{ConstantArray[First@#, Length@# - 1], Rest@#}, {3, 1, 2}] &@list OleksandrR[list_] := Flatten[Outer[List, List@First[list], Rest[list], 1], {{2}, {1, 4}}] paradox2[list_] := Partition[Riffle[list[[1]], #], 2] & /@ Drop[list, 1] RM[list_] := FoldList[Transpose[{First@li...
Post a Comment
Read more

front end - keyboard shortcut to invoke Insert new matrix

I frequently need to type in some matrices, and the menu command Insert > Table/Matrix > New... allows matrices with lines drawn between columns and rows, which is very helpful. I would like to make a keyboard shortcut for it, but cannot find the relevant frontend token command (4209405) for it. Since the FullForm[] and InputForm[] of matrices with lines drawn between rows and columns is the same as those without lines, it's hard to do this via 3rd party system-wide text expanders (e.g. autohotkey or atext on mac). How does one assign a keyboard shortcut for the menu item Insert > Table/Matrix > New... , preferably using only mathematica? Thanks! Answer In the MenuSetup.tr (for linux located in the $InstallationDirectory/SystemFiles/FrontEnd/TextResources/X/ directory), I changed the line MenuItem["&New...", "CreateGridBoxDialog"] to read MenuItem["&New...", "CreateGridBoxDialog", MenuKey["m", Modifiers-...
Post a Comment
Read more