[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