[yt-users] Projection of 2D cylindrical FLASH data
Ken Shen
kenshen at astro.berkeley.edu
Fri May 15 10:29:35 PDT 2015
Thanks for the suggestions, I'll give them a try!
On Fri, May 15, 2015 at 8:18 AM, Roman Yurchak <roman.yurchak at crans.org>
wrote:
> If your observation direction is orthogonal to the revolution axis, you
> could also map your 2D cylindrical data to a 2D numpy array (i.e.
> re-sample on a uniform grid) and then do a forward Abel transform (the
> integral transform to do this particular projection), implemented for
> instance in https://github.com/luli/hedp/blob/master/hedp/math/abel.py#L19
> Although, volume rendering approaches would probably be more general and
> better integrated with yt.
> --
> Roman
>
>
> On 15/05/15 03:27, Suoqing JI wrote:
> > One way (imperfect but definitely doable) you might want to try is to
> > map 2D cylindrical data into a 3D numpy array, which could be loaded
> > by yt.load_uniform_grid and passed into the volume rendering module.
> >
> > See here: http://yt-project.org/doc/examining/generic_array_data.html
> >
> > Best wishes,
> > --
> > Suoqing JI
> > Ph.D Student
> > Department of Physics
> > University of California, Santa Barbara
> > http://web.physics.ucsb.edu/~suoqing
> >
> >> On May 14, 2015, at 6:01 PM, Nathan Goldbaum <nathan12343 at gmail.com
> >> <mailto:nathan12343 at gmail.com>> wrote:
> >>
> >>
> >>
> >> On Thu, May 14, 2015 at 5:35 PM, Ken Shen <kenshen at astro.berkeley.edu
> >> <mailto:kenshen at astro.berkeley.edu>> wrote:
> >>
> >> I'm really sorry, I can't seem to explain this very clearly in
> words!
> >>
> >> I have 2D axisymmetric data. Let's suppose that I map this to a
> >> 3D Cartesian dataset. I'd like to make a 2D plot from this 3D
> >> data that consists of integrating the density along rays from the
> >> observer through the dataset onto a plane on the other side of the
> >> data. Ideally, I'd like to do this with the observer at an
> >> arbitrary position with respect to the axes.
> >>
> >> I haven't done it, but it seems like this is built in to yt's
> >> functionality for 3D datasets. But I haven't found how to do it
> >> for a 2D axisymmetric dataset.
> >>
> >>
> >> In yt we would call this a volume rendering using rays that are
> >> aligned with the cartesian axes and a projection transfer function.
> >>
> >> I don't think this has been implemented yet, but volume renderings of
> >> data with non-cartesian geometries or with unstructured meshes is
> >> something we would like to do.
> >>
> >>
> >>
> >> Sorry for all the messages!
> >> Ken
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >> On Thu, May 14, 2015 at 5:28 PM, Michael Zingale
> >> <michael.zingale at stonybrook.edu
> >> <mailto:michael.zingale at stonybrook.edu>> wrote:
> >>
> >> to be clear here, what Flash calls cylindrical is 2-d
> >> axisymmetric, r-z. If I understand what you are asking, you
> >> basically want to still have a 2-d plot where each cell is the
> >> integral over an annular region as if that cell were rotated
> >> about the symmetry axis -- is that right?
> >>
> >> On Thu, May 14, 2015 at 8:21 PM, Ken Shen
> >> <kenshen at astro.berkeley.edu
> >> <mailto:kenshen at astro.berkeley.edu>> wrote:
> >>
> >> Hi, sorry I wasn't very clear before. I'd like to make 2D
> >> column density plots of my 2D cylindrical data as if my
> >> data were 3D cylindrical data with no phi-dependence.
> >>
> >> (The ProjectionPlot fails for axis = 1 or 3. For axis =
> >> 2, it essentially recreates the regular density slice plot.)
> >>
> >>
> >> Thanks,
> >> Ken
> >>
> >>
> >>
> >>
> >> On Thu, May 14, 2015 at 5:00 PM, Matthew Turk
> >> <matthewturk at gmail.com <mailto:matthewturk at gmail.com>>
> wrote:
> >>
> >> Hi Ken,
> >>
> >> (I think Suoqing's answer may work as well.) If you
> >> want to project
> >> along r, z, or theta, you can do a ProjectionPlot (or
> >> "ds.proj") and
> >> specify "r", "z", or "theta", and the result will be
> >> along that axis
> >> with the correct path lengths. I think that the "r"
> >> direction may be
> >> what you're looking for here, i.e., a soup-can label.
> >>
> >> -Matt
> >>
> >> On Thu, May 14, 2015 at 6:31 PM, Ken Shen
> >> <kenshen at astro.berkeley.edu
> >> <mailto:kenshen at astro.berkeley.edu>> wrote:
> >> > Hi all, I'd like to make column density plots of my
> >> 2D cylindrical FLASH
> >> > data. Ideally, these would be at arbitrary angles
> >> to the z-axis. I'm new
> >> > to yt, so I could very well have missed it, but I
> >> couldn't seem to find this
> >> > in the documentation. Any hints?
> >> >
> >> >
> >> > Thanks,
> >> > Ken
> >> >
> >> > _______________________________________________
> >> > yt-users mailing list
> >> > yt-users at lists.spacepope.org
> >> <mailto:yt-users at lists.spacepope.org>
> >> >
> >>
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >> >
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org
> >> <mailto:yt-users at lists.spacepope.org>
> >>
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >>
> >>
> >>
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org
> >> <mailto:yt-users at lists.spacepope.org>
> >>
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >>
> >>
> >>
> >>
> >> --
> >> Michael Zingale
> >> Associate Professor
> >>
> >> Dept. of Physics & Astronomy • Stony Brook University • Stony
> >> Brook, NY 11794-3800
> >> /phone/: 631-632-8225 <tel:631-632-8225>
> >> /e-mail/: Michael.Zingale at stonybrook.edu
> >> <mailto:Michael.Zingale at stonybrook.edu>
> >> /web/: http://www.astro.sunysb.edu/mzingale
> >>
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org <mailto:
> yt-users at lists.spacepope.org>
> >> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >>
> >>
> >>
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> >> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >>
> >>
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org <mailto:yt-users at lists.spacepope.org>
> >> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> >
> >
> > _______________________________________________
> > yt-users mailing list
> > yt-users at lists.spacepope.org
> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20150515/3c8b69a4/attachment.htm>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
More information about the yt-users
mailing list