<div dir="ltr">I tried this, but I still get a blank black figure:<div><br></div><div><div>import numpy as np</div><div>from yt.mods import *</div><div>from yt.analysis_modules.halo_finding.api import *</div><div>Npixels=512</div>
<div>L=[1,1,1]</div><div>c_x =np.array( [ 0.5696202374 , 0.5330663554 , 0.4554251862 ])</div><div>c_X =np.array( [ 0.5657944233 , 0.524785043 , 0.4587818409 ])</div><div>for i in range(42,41,-1):</div><div> t1="RD00%s"%i</div>
<div> t2="RedshiftOutput00%s"%i</div><div> pf=load(t1+'/'+t2)</div><div><br></div><div>W=1/pf["Mpc"]</div><div>c=(c_x+c_X)/2.</div><div>le=c-np.array(3*[2])/pf['Mpc']</div>
<div>re=c+np.array(3*[2])/pf['Mpc']</div><div><br></div><div>region=pf.h.region(c, le,re)</div><div><br></div><div>region=pf.h.region(c, le,re,fields="particle_density" )</div><div>mi,ma=region.quantities["Extrema"]("particle_density")[0]</div>
<div>mi=-24</div><div>ma=np.log10(ma)</div><div>tf = ColorTransferFunction((mi,ma))</div><div>tf.add_layers(5, w=0.02, colormap="spectral")</div><div>cam=region.pf.h.camera(c, L, W, Npixels, tf,fields=["particle_density"],log_fields=True)</div>
<div>cam.snapshot("%s_volume_rendered.png" % pf, clip_ratio=8.0)</div></div><div><br></div><div style>-Mohammad</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Sep 4, 2013 at 1:13 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com" target="_blank">matthewturk@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Mohammad,<br>
<br>
You have to tell the transfer function that you want some levels to<br>
correspond to various colors. An easy way to do that is with<br>
tf.add_layers(). Otherwise, the transfer function is blank, and it<br>
won't update the image as the rays traverse the domain.<br>
<br>
-Matt<br>
<div class="HOEnZb"><div class="h5"><br>
On Wed, Sep 4, 2013 at 1:04 PM, mohammad safarzadeh<br>
<<a href="mailto:mtsafarzadeh@gmail.com">mtsafarzadeh@gmail.com</a>> wrote:<br>
> Hi,<br>
> I am trying to do VR on a halo in my zoom-in region which is a dark matter<br>
> only sim, but I get a black image when I try this:<br>
><br>
> L=[1,1,1]<br>
><br>
> W=1/pf["Mpc"]<br>
><br>
> n [16]: mi,ma=region.quantities["Extrema"]("particle_density")[0]<br>
><br>
><br>
> In [17]: mi<br>
><br>
> Out[17]: 0.0<br>
><br>
><br>
> In [18]: ma<br>
><br>
> Out[18]: 4.5282895435060503e-24<br>
><br>
><br>
> In [21]: np.log10(ma)<br>
><br>
> Out[21]: -23.344065811694637<br>
><br>
><br>
> mi=-24<br>
><br>
> ma=np.log10(ma)<br>
><br>
> tf = ColorTransferFunction((mi,ma))<br>
><br>
> cam=region.pf.h.camera(c, L, W, Npixels,<br>
> tf,fields=["particle_density"],log_fields=True)<br>
><br>
><br>
> ImageArray([[[ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> ...,<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.]],<br>
><br>
><br>
> [[ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> ...,<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.],<br>
><br>
> [ 0., 0., 0., 1.]],….<br>
><br>
><br>
> what is the mistake that I am doing?<br>
><br>
><br>
> -Regards,<br>
><br>
> Mohammad<br>
><br>
><br>
</div></div><div class="HOEnZb"><div class="h5">> _______________________________________________<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>
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>
</div></div></blockquote></div><br></div>