Hi Matt<div><br></div><div>I have until now:</div><meta http-equiv="content-type" content="text/html; charset=utf-8"><div>* Projection plugin: it shows the image in paraview.</div><div>* Slice plugin: it's also an image, but it' placed inside the volume, so you can see the image as a common slice in 3D.</div>
<div>* Halo finder: halos information are recovered from YT and represented in Paraview using spheres.</div><div><br></div><div>Now, I'm trying to run a simple volume rendering, and display the image in paraview. It runs. but the image is not</div>
<div>the same as loading data directly with yt. I think the problem is the conversion factor. because using the</div><div>StreamHandler, the data is saved without conversion. </div><div>My script is:</div><div><div><br></div>
<div><font class="Apple-style-span" face="'courier new', monospace">pf = stream.StreamStaticOutput(handler)</font></div></div><div><div><font class="Apple-style-span" face="'courier new', monospace">L = [0.5, 0.2, 0.7]</font></div>
<div><font class="Apple-style-span" face="'courier new', monospace">W = 1.0</font></div><div><font class="Apple-style-span" face="'courier new', monospace">c = [0.5, 0.5, 0.5]</font></div><div><font class="Apple-style-span" face="'courier new', monospace">Nvec = 512</font></div>
<div><font class="Apple-style-span" face="'courier new', monospace">dd = pf.h.all_data()</font></div><div><font class="Apple-style-span" face="'courier new', monospace">mi, ma = dd.quantities["Extrema"]("Density")[0]</font></div>
<div><font class="Apple-style-span" face="'courier new', monospace">tf = ColorTransferFunction((na.log10(mi), na.log10(ma)))</font></div><div><font class="Apple-style-span" face="'courier new', monospace">tf.add_layers(6, w=0.02)</font></div>
<div><font class="Apple-style-span" face="'courier new', monospace">vp = pf.h.volume_rendering(L, W, c, Nvec, tf, whole_box=True)</font></div><div><font class="Apple-style-span" face="'courier new', monospace">vp.ray_cast()</font></div>
</div><div><br></div><div>What would be the best way to solve it?</div><div><br></div><div><br></div><div>Regards</div><div><br></div><div><br></div><div><br></div><div>Halo Finder is working with paraview, halos are represented in paraview using spheres.<br>
<br><div class="gmail_quote">On Fri, Jun 17, 2011 at 1:39 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>
<div class="im"><br>
On Fri, Jun 17, 2011 at 9:51 AM, Jorge Poco <<a href="mailto:jorge.poco@kitware.com">jorge.poco@kitware.com</a>> wrote:<br>
> Hi Matt<br>
> I will try to use HaloFinder from YT, the script would be:<br>
<br>
</div>Wow, cool!<br>
<div class="im"><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 spheres<br>
> which<br>
><br>
> could be visualized using paraview.<br>
><br>
> Questions:<br>
><br>
> * HaloFinder needs the particle information right?<br>
<br>
</div>Yup, that's right.<br>
<div class="im"><br>
><br>
> * Current stream.StreamHandler support the particle data?<br>
<br>
</div>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>
</blockquote></div><br></div>