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@gmail.com
>> <mailto:nathan12343@gmail.com>> wrote:
>>
>>
>>
>> On Thu, May 14, 2015 at 5:35 PM, Ken Shen <kenshen@astro.berkeley.edu
>> <mailto:kenshen@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@stonybrook.edu
>> <mailto:michael.zingale@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@astro.berkeley.edu
>> <mailto:kenshen@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@gmail.com <mailto:matthewturk@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@astro.berkeley.edu
>> <mailto:kenshen@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@lists.spacepope.org
>> <mailto:yt-users@lists.spacepope.org>
>> >
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>> >
>> _______________________________________________
>> yt-users mailing list
>> yt-users@lists.spacepope.org
>> <mailto:yt-users@lists.spacepope.org>
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users@lists.spacepope.org
>> <mailto:yt-users@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@stonybrook.edu
>> <mailto:Michael.Zingale@stonybrook.edu>
>> /web/: http://www.astro.sunysb.edu/mzingale
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users@lists.spacepope.org <mailto:yt-users@lists.spacepope.org>
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users@lists.spacepope.org <mailto:yt-users@lists.spacepope.org>
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users@lists.spacepope.org <mailto:yt-users@lists.spacepope.org>
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users@lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>