Mailman 3 Projection of 2D cylindrical FLASH data - yt-users

newer
Re: [yt-users] Projection of 2D...

Projection of 2D cylindrical FLASH data

older
Cython Build Errors in Development...

Ken Shen

14 May 2015 14 May '15

11:31 p.m.

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 http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Attachments:

attachment.html (text/html — 425 bytes)

Show replies by date

Suoqing JI

14 May 14 May

11:47 p.m.

Hi Ken, Just to clarify: when you referred to the column density plot, do you mean the 1D plot of column density which is accumulated along z direction? If so, you could something like the following: profile = BinnedProfile1D(mydisk, Nbin, 'cylindrical_r', rmin, rmax, log_space=True, lazy_reader=True, end_collect=False) profile.add_fields('cell_mass', weight=None) R = profile['cylindrical_r’] Sigma = profile[‘cell_mass'] R_edge = np.logspace(np.log10(rmin), np.log10(rmax), num=Nbin+1) for i in range(0, Nbin): Sigma[i] = Sigma[i] / (np.pi*(R_edge[i+1]**2 - R_edge[i]**2.)) 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 4:31 PM, Ken Shen 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 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

Matthew Turk

15 May 15 May

midnight

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 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 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

Ken Shen

12:21 a.m.

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 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

...
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

On Thu, May 14, 2015 at 6:31 PM, Ken Shen wrote: this

...
in the documentation. Any hints?

Thanks, Ken

_______________________________________________ yt-users mailing list 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

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Michael Zingale

12:28 a.m.

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 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 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

...
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

On Thu, May 14, 2015 at 6:31 PM, Ken Shen wrote: this

...
in the documentation. Any hints?

Thanks, Ken

_______________________________________________ yt-users mailing list 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

_______________________________________________ yt-users mailing list 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 *e-mail*: Michael.Zingale@stonybrook.edu *web*: http://www.astro.sunysb.edu/mzingale _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Ken Shen

12:35 a.m.

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. Sorry for all the messages! Ken On Thu, May 14, 2015 at 5:28 PM, Michael Zingale < 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 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 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 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 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

_______________________________________________ yt-users mailing list 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 *e-mail*: Michael.Zingale@stonybrook.edu *web*: http://www.astro.sunysb.edu/mzingale

_______________________________________________ yt-users mailing list 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

Nathan Goldbaum

1:01 a.m.

On Thu, May 14, 2015 at 5:35 PM, Ken Shen 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> 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 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 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 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 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

_______________________________________________ yt-users mailing list 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 *e-mail*: Michael.Zingale@stonybrook.edu *web*: http://www.astro.sunysb.edu/mzingale

_______________________________________________ yt-users mailing list 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

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Suoqing JI

1:27 a.m.

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 wrote:

On Thu, May 14, 2015 at 5:35 PM, Ken Shen 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 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 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 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 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 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 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 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 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 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 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

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Roman Yurchak

3:24 p.m.

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 mailto:nathan12343@gmail.com> wrote:

On Thu, May 14, 2015 at 5:35 PM, Ken Shen 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 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 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 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 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

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Michael Zingale

5:19 p.m.

that's really a nice solution. On Fri, May 15, 2015 at 11:24 AM, Roman Yurchak 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 mailto:nathan12343@gmail.com> wrote:

On Thu, May 14, 2015 at 5:35 PM, Ken Shen 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 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 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 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 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
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

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

Matthew Turk

12:34 a.m.

Hi Ken, I've just now tested out on some data that I have here and I'm getting odd things for projections in 2D cylindrical coordinates, which I'll check and add tests for to the test suite. On Thu, May 14, 2015 at 7:21 PM, Ken Shen 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 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 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 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

_______________________________________________ yt-users mailing list 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

3263

Age (days ago)

3264

Last active (days ago)

List overview

Download

10 comments

6 participants

participants (6)

Ken Shen
Matthew Turk
Michael Zingale
Nathan Goldbaum
Roman Yurchak
Suoqing JI

Projection of 2D cylindrical FLASH data

Ken Shen

Suoqing JI

Matthew Turk

Ken Shen

Michael Zingale

Ken Shen

Nathan Goldbaum

Suoqing JI

Roman Yurchak

Michael Zingale

Matthew Turk

tags

participants (6)