Re: Polar(List)DensityPlot
- To: mathgroup at smc.vnet.net
- Subject: [mg21289] Re: Polar(List)DensityPlot
- From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
- Date: Wed, 22 Dec 1999 00:34:43 -0500 (EST)
- Organization: Universitaet Leipzig
- References: <83nf1r$5al@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Hi,
I assume you have a 2d array {r,phi,f[r,phi]}
data = Table[ N[{r, phi, r^2*Sin[phi]}], {r, 0, 4, 0.5}, {phi, 0, 2Pi,
Pi/9}];
than:
makePolarMesh[vl_List] :=
Module[{l = vl, l1 = Map[RotateLeft, vl], mesh, f},
f = Function[{r, phi, z}, {r*Cos[phi], r*Sin[phi], z}];
mesh = {l, l1, RotateLeft[l1], RotateLeft[l]};
mesh = Map[Drop[#, -1] &, mesh, {1}];
mesh = Map[Drop[#, -1] &, mesh, {2}];
mesh = Map[
Polygon[Reverse[#]] &,
Transpose[
Flatten[#, 1] & /@
mesh
]
] // Flatten;
mesh = mesh /. Polygon[pnts_] :> Polygon[f @@@ pnts]
]
and
Show[Graphics3D[makePolarMesh[data]]]
will work.
Hope that helps
Jens
"W.F.J.C. Meyboom" wrote:
>
> Hi,
>
> I'd like to make something best described as a PolarDensityPlot or a
> PolarListDensityPlot.
>
> Does anyone know how to do that ?
>
> Greets,
>
> Wilfred Meyboom