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