Hi Elizabeth,

There is probably a better way to do this, but since Orion doesn't have particle support at the moment, I have a routine for drawing the particle at position "r" onto a volume rendered image with arbitrary orientation.

import numpy as na

def add_star_projected(r,ax,cam,ms):
    """ Adds the particle at position r to Axes ax. Needs information on the orientation of the image plane contained in cam """
    north_vector = cam.unit_vectors[0]
    east_vector = cam.unit_vectors[1]
    N1 = cam.resolution[0]
    N2 = cam.resolution[1]
    c = cam.center
    r = r - c
    W = cam.width[0]
    y_pixels = floor(na.dot(r,north_vector) * (N2/W) + N2/2.0)
    x_pixels = floor(na.dot(r,east_vector) * (N1/W) + N1/2.0)
    ax.hold(True)
    ax.plot([x_pixels],[y_pixels],'wo', markersize = ms)
    ax.hold(False)
    ax.axis((0,N1,0,N2))

Assuming you're using the Camera interface, that contains all the information you need to do the conversion. If you can get a list of particles w/ their positions (we have our own way of doing this that probably doesn't apply), you can then loop through that and add them at the correct position on the off-axis projection.

Andrew M

On Thu, Jan 6, 2011 at 6:15 AM, Elizabeth Harper-Clark <h-clark@astro.utoronto.ca> wrote:
Hi all,

Thanks for all your help over the last couple of days. One more question:
    -  Can I plot particles on a volume rendered image?
I have stars and I want to show where they are!

Thanks,

Libby

--
Elizabeth Harper-Clark MA MSci
PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT
Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT

www.astro.utoronto.ca/~h-clark
h-clark@cita.utoronto.ca
Astronomy office phone: +1-416-978-5759

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