Re: Light and surface colors
- To: mathgroup at smc.vnet.net
- Subject: [mg58704] Re: [mg58688] Light and surface colors
- From: "David Park" <djmp at earthlink.net>
- Date: Fri, 15 Jul 2005 03:02:19 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Daniele,
I tried to send you a Mathematica notebook with the code for this but your
server rejected it because of a full mailbox.
To make a shaded plot with regular Mathematica graphics we could use
ParametricPlot3D with the 4th argument in the parametrization. The 4th
argument gives graphics directives for plotting the surface. We don't have
to turn the Lighting off if we give SurfaceColor instead of a raw color.
We redefine the pressure for the parameters of the ellipsoid.
ellipsoid[a_, b_, c_][u_, v_] := {a Cos[v]Cos[u], b Cos[v]Sin[u], c Sin[v]}
pressure[x_, y_, z_] := x^2 + y^2 + Sin[4z];
pressure[u_, v_] =
pressure[x, y, z] /. Thread[{x, y, z} -> ellipsoid[7, 5, 4][u, v]]
49*Cos[u]^2*Cos[v]^2 + 25*Cos[v]^2*Sin[u]^2 + Sin[16*Sin[v]]
Then define the parametrization for the ellipsoid with the function on it.
ellipsoidfunction[a_, b_, c_][u_, v_] := {a Cos[v]Cos[u], b Cos[v]Sin[u],
c Sin[v], {EdgeForm[], SurfaceColor[Hue[pressure[u, v]/50]]}}
Then we can plot.
ParametricPlot3D[ellipsoidfunction[7, 5, 4][u, v], {u, 0, 2*Pi},
{v, -Pi/2, Pi/2},
PlotPoints -> 100,
ImageSize -> 450];
There are three problems with this plot. 1) The regular Mathematica lighting
tends to interfer with the colors on the surface. 2) The Hue color function
is not easy to interpret. What is high pressure and what is low pressure? 3)
It doesn't give a very accurate indication of what the shape of the contours
are. They are not quite as regular as one might expect.
My second solution uses DrawGraphics (from my web site below). It trims a
grid of polygons in the uv-plane to the contour regions and then maps the
resulting grids to the ellipsoidal surface. It uses the TwoColorScale
function to give a color scale that gives shades of blue for values below 25
and shades of red for colors above 25. It is much easier to interpret
two-color plots - they are similar to cartographic maps. The NeutralLighting
option is used to turn down the color saturation of the regular lighting so
it doesn't interfer with the surface colors.
Needs["DrawGraphics`DrawingMaster`"]
colorfun =
TwoColorScale[{-0.29, 25, 49.7}, {ColorMix[RoyalBlue, Black][0.5],
RoyalBlue, ColorMix[RoyalBlue, White][0.8], ColorMix[Red,
White][0.8],
Red, ColorMix[Red, Black][0.5]}];
polygrid = N[MakePolyGrid[{40, 20}, {{0, -Pi/2}, {2*Pi, Pi/2}}]];
contours = Range[-0.29, 49.7, (49.7 + 0.29)/10];
Do[polys[i] = TrimPolygons[pressure[#1, #2] & ,
{contours[[i]], contours[[i + 1]]}][polygrid],
{i, 1, Length[contours] - 1}]
Draw3DItems[
{Table[{SurfaceColor[colorfun[(contours[[i]] + contours[[i + 1]])/2]],
EdgeForm[], polys[i] /. {(u_)?NumberQ, (v_)?NumberQ} ->
ellipsoid[7, 5, 4][u, v]}, {i, 1, Length[contours] - 1}]},
BoxRatios -> Automatic,
NeutralLighting[0.1, 0.5, 0],
Boxed -> False,
Axes -> False,
PlotRegion -> {{-0.1, 1.1}, {-0.3, 1.2}},
ViewPoint -> {1.534, -2.831, 1.04},
Background -> Linen,
ImageSize -> 450];
The contours actually are somewhat irregular and wavy as one can check by
making a ContourPlot of the pressure in the uv-plane using a high number of
plot points.
David Park
djmp at earthlink.net
http://home.earthlink.net/~djmp/
From: Daniele Lupo [mailto:danwolf_no_spam_please_ at libero.it]
To: mathgroup at smc.vnet.net
Hi to everyone.
I've just obtained a Student License of Mathematica, and I'm working in
graphics...
I've created, after a half-day of work, an image representing an ellipsoid,
with a colored surface
When I try to represent it, I must put Ligthing->False to see surface
colors.
I'd like to know if there's a way to obtain a rendered image, in which I
can have both light and surface color, to give a shadow to the figure and
create a more 3D effect.
Thanks for replies.
Below there's the code...
Daniele
--------------------------
<< Graphics`ContourPlot3D`;
a = 7; b = 5; c = 4; r = 1;
ellipsoid = ContourPlot3D[
(x/a)^2 + (y/b)^2 + (z/c)^2 - r, {
x, -10, 10}, {y, -10, 10}, {z, -10, 10},
PlotPoints -> 12,
AspectRatio -> Automatic
];
(*It finds the median point of vertices of a list of a 3D coordinates*)
midpoint[l_] := Apply[Plus, Transpose[l], 2]/Length[l];
(*Function to map in the ellipsoid surface*)
pressure[x_, y_, z_] := x^2 + y^2 + Sin[4z];
(*Color Function*)
color[x_] := Hue[pressure @@ midpoint[x]/50];
(*Rule to change original graphic object*)
colorpolygon = Polygon[coord_] :>
{color[coord],EdgeForm[{}],Polygon[coord]};
(*creation of new, colored ellipsoid*)
coloredellipse = ellipsoid /. colorpolygon;
(*Visualization of the shape*)
(*I'd like to have also the light...*)
Show[coloredellipse, Lighting -> False];