[yt-users] projection depth

[John Wise]

Hi,

I would think that you'd have to change the region boundaries, if you 
wanted to change the projection depth.  To have L=0.2, you have to 
change your region boundaries to reflect that change, i.e. [0.4,0.4,0.4] 
-> [0.6,0.6,0.6].


If I wanted to make projections with a given width and depth, I'd do the 
following.

L = 0.1  # projection width
d = 0.2  # projection depth
center = na.array([0.5]*3)
dim = 0

left_edge = center - 0.5*L
right_edge = center + 0.5*L

# For the projection dimension, use the depth
left_edge[dim] = center - 0.5*d
right_edge[dim] = center + 0.5*d

pf = get_pf()
reg = pf.h.region(center, left_edge, right_edge)

pc = PlotCollection(reg, center=center)
pc.add_projection("Density", 0, "Density")
pc.set_width(L, '1')
pc.save("test", override=True)

Cheers,
John

On 04/13/2010 10:21 AM, Luigi Iapichino wrote:
> Dear all,
>
> by comparing volume renderings and projections in my simulation I realised
> that I don't have control on the projection depth, in the latter case.
> This is the relevant part of my script:
>
> ...
>
> L = 0.1 # projection width
>
> pf = load(data_dir+fn) # load data
>
> reg = pf.h.region([0.5, 0.5, 0.5],[0.45,0.45,0.45],[0.55,0.55,0.55]) # a small
> region
> #reg = pf.h.region([0.5, 0.5, 0.5],[0,0,0],[1,1,1]) # the whole box
>
> pc = PlotCollection(reg, center=[0.5,0.5,0.5])
>
> p = pc.add_projection("Density", 0, weight_field="Density",fig_size=(7.5,6))
>
> p.modify["contour"]("Density", ncont=5,take_log=True)
>
> pc.set_width(L, '1')
>
> pc.save("%s_test-proj" % pf,override=True)
>
>
> If I set L = 0.2, I don't see in the central region of the projection more
> structures than with L=0.1 (only the plot width is, of course, larger).
> Moreover, the result do not depend on the region that I select (either third
> or fourth line). It seems to me that, in all cases, the projection is
> actually performed along the whole simulation box, and not on the region that
> I select. In other words, I can control the plot width, but not the
> projection depth. In the volume rendering, on the other hand, the plot width
> sets the depth, too. Is anything wrong or missing in my script?
>
> If needed, I can send you a few explanatory plots and the whole sample script
> for reproducing my issue.
> Looking forward to your suggestions,
>
>    Luigi
>