Hi Matt<div><br></div><div>I found the problem.</div><div><div><font class="Apple-style-span" face="'courier new', monospace"><br></font></div><div><font class="Apple-style-span" face="'courier new', monospace">pf = stream.StreamStaticOutput(handler)</font></div>
<div><font class="Apple-style-span" face="'courier new', monospace">pf.h</font></div><div><font class="Apple-style-span" face="'courier new', monospace"><span class="Apple-style-span" style="background-color: rgb(255, 255, 102);">pf.field_info['Density'].take_log = True</span></font></div>
<div><br></div><div>pf.field_info['Density'].take_log was False and should be True, because Transfer function was </div><div>working in log scale.</div><div> </div><div>Regards</div><div><br></div><div class="gmail_quote">
On Wed, Jun 22, 2011 at 12:22 PM, Matthew Turk <span dir="ltr"><<a href="mailto:matthewturk@gmail.com">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 Jorge,<br>
<br>
Sure -- I'm going to be away from my desk for a little while, but I<br>
will get back to you this afternoon (pacific time, unfortunately).<br>
<br>
-Matt<br>
<div><div></div><div class="h5"><br>
On Wed, Jun 22, 2011 at 9:21 AM, Jorge Poco <<a href="mailto:jorge.poco@kitware.com">jorge.poco@kitware.com</a>> wrote:<br>
> Hi Matt<br>
><br>
> I have tested the volume rendering using your first script to convert from<br>
> Enzo loader to StreamHandler<br>
> and it's giving me bad results. Could you check the results with that<br>
> script, please?<br>
><br>
> Regards<br>
><br>
> On Wed, Jun 22, 2011 at 11:35 AM, Jorge Poco <<a href="mailto:jorge.poco@kitware.com">jorge.poco@kitware.com</a>> wrote:<br>
>><br>
>> Hi Matt<br>
>> images are attached.<br>
>><br>
>> On Wed, Jun 22, 2011 at 11:24 AM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
>> wrote:<br>
>>><br>
>>> Hi Jorge,<br>
>>><br>
>>> It's not entirely clear to me. I think what you have done should be<br>
>>> immune to a difference in multiplicative factor. Would it be possible<br>
>>> to see the image that comes out of PV and compare to the image<br>
>>> directly out of yt?<br>
>>><br>
>>> -Matt<br>
>>><br>
>>> On Wed, Jun 22, 2011 at 8:07 AM, Jorge Poco <<a href="mailto:jorge.poco@kitware.com">jorge.poco@kitware.com</a>><br>
>>> wrote:<br>
>>> > Hi Matt<br>
>>> > I have until now:<br>
>>> > * Projection plugin: it shows the image in paraview.<br>
>>> > * Slice plugin: it's also an image, but it' placed inside the volume,<br>
>>> > so you<br>
>>> > can see the image as a common slice in 3D.<br>
>>> > * Halo finder: halos information are recovered from YT and represented<br>
>>> > in<br>
>>> > Paraview using spheres.<br>
>>> > Now, I'm trying to run a simple volume rendering, and display the image<br>
>>> > in<br>
>>> > paraview. It runs. but the image is not<br>
>>> > the same as loading data directly with yt. I think the problem is the<br>
>>> > conversion factor. because using the<br>
>>> > StreamHandler, the data is saved without conversion.<br>
>>> > My script is:<br>
>>> > pf = stream.StreamStaticOutput(handler)<br>
>>> > L = [0.5, 0.2, 0.7]<br>
>>> > W = 1.0<br>
>>> > c = [0.5, 0.5, 0.5]<br>
>>> > Nvec = 512<br>
>>> > dd = pf.h.all_data()<br>
>>> > mi, ma = dd.quantities["Extrema"]("Density")[0]<br>
>>> > tf = ColorTransferFunction((na.log10(mi), na.log10(ma)))<br>
>>> > tf.add_layers(6, w=0.02)<br>
>>> > vp = pf.h.volume_rendering(L, W, c, Nvec, tf, whole_box=True)<br>
>>> > vp.ray_cast()<br>
>>> > What would be the best way to solve it?<br>
>>> ><br>
>>> > Regards<br>
>>> ><br>
>>> ><br>
>>> > Halo Finder is working with paraview, halos are represented in paraview<br>
>>> > using spheres.<br>
>>> ><br>
>>> > On Fri, Jun 17, 2011 at 1:39 PM, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>><br>
>>> > wrote:<br>
>>> >><br>
>>> >> Hi Jorge,<br>
>>> >><br>
>>> >> On Fri, Jun 17, 2011 at 9:51 AM, Jorge Poco <<a href="mailto:jorge.poco@kitware.com">jorge.poco@kitware.com</a>><br>
>>> >> wrote:<br>
>>> >> > Hi Matt<br>
>>> >> > I will try to use HaloFinder from YT, the script would be:<br>
>>> >><br>
>>> >> Wow, cool!<br>
>>> >><br>
>>> >> > # create Paraview Handler<br>
>>> >> > pf = stream.StreamStaticOutput(handler)<br>
>>> >> > # halos will be added to this polydata<br>
>>> >> > spheres = vtkAppendPolyData()<br>
>>> >> > halos = HaloFinder(pf)<br>
>>> >> > for halo in halos:<br>
>>> >> > # create a vtkSphere using halo information<br>
>>> >> > s = vtkSphereSource()<br>
>>> >> > s.SetCenter(halo.get_sphere().center)<br>
>>> >> > s.SetRadius(halo.get_sphere().radius)<br>
>>> >> > s.Update()<br>
>>> >> > spheres.AddInput(s.GetOutput())<br>
>>> >> > spheres.Update()<br>
>>> >> ><br>
>>> >> > output.SetPoints(s1.GetOutput().GetPoints())<br>
>>> >> ><br>
>>> >> > output.SetPolys(s1.GetOutput().GetPolys())<br>
>>> >> ><br>
>>> >> > Using this script we will be able to recover the halos as a set of<br>
>>> >> > spheres<br>
>>> >> > which<br>
>>> >> ><br>
>>> >> > could be visualized using paraview.<br>
>>> >> ><br>
>>> >> > Questions:<br>
>>> >> ><br>
>>> >> > * HaloFinder needs the particle information right?<br>
>>> >><br>
>>> >> Yup, that's right.<br>
>>> >><br>
>>> >> ><br>
>>> >> > * Current stream.StreamHandler support the particle data?<br>
>>> >><br>
>>> >> Yup, it should, although in the future we can make this work more<br>
>>> >> seamlessly. The easiest way to supply this would be to push them<br>
>>> >> through the stream IO handler, and supply them as (1D) fields called<br>
>>> >> "particle_mass", "particle_position_[xyz]", "particle_velocity_[xyz]".<br>
>>> >> I think this should work just fine, but there may be other hangups.<br>
>>> >> If you run into any problems, let me know.<br>
>>> >><br>
>>> >> Best,<br>
>>> >><br>
>>> >> -Matt<br>
>>> >><br>
>>> >> ><br>
>>> >> > Regards,<br>
>>> >> ><br>
>>> >> > Att. Jorge Poco<br>
>>> >> ><br>
>>> >> > _______________________________________________<br>
>>> >> > Amr mailing list<br>
>>> >> > <a href="mailto:Amr@public.kitware.com">Amr@public.kitware.com</a><br>
>>> >> > <a href="http://public.kitware.com/cgi-bin/mailman/listinfo/amr" target="_blank">http://public.kitware.com/cgi-bin/mailman/listinfo/amr</a><br>
>>> >> ><br>
>>> >> ><br>
>>> ><br>
>>> ><br>
>><br>
><br>
><br>
</div></div></blockquote></div><br></div>