<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr">Hi Carolin,<div><br></div><div>The first argument to ProjectionPlot should be a dataset object and not a data container. Can you try this instead:</div><div><span style="font-size:12.8000001907349px">pl=yt.ProjectionPlot(</span><span style="font-size:12.8000001907349px"><b>mcyt</b></span><span style="font-size:12.8000001907349px">, 'x', 'y')</span><br style="font-size:12.8000001907349px"></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">Additionally, if you want to project just the values in the sphere, you can try this:</span></div><div><div><span style="font-size:12.8000001907349px">pl=yt.ProjectionPlot(</span><span style="font-size:12.8000001907349px"><b>mcyt</b></span><span style="font-size:12.8000001907349px">, 'x', 'y', data_source=circle)</span><br style="font-size:12.8000001907349px"></div></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">You're a pioneer in unexplored territory here, so please let us know how it goes!</span></div><div><span style="font-size:12.8000001907349px"><br></span></div><div><span style="font-size:12.8000001907349px">Britton</span></div><div><span style="font-size:12.8000001907349px"><br></span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 23, 2015 at 10:56 AM, Carolin Villforth <span dir="ltr"><<a href="mailto:cv21@st-andrews.ac.uk" target="_blank">cv21@st-andrews.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello yt-users,<br>
<br>
I saw a yt presentation at the Python Astronomy Meeting and am interested to see if yt could be used for non-spatial data. My specific problem is MCMC chains in 19 (or up to 23) dimensional parameter space in which none of the parameters are nuisance parameters. I thought yt could be used to visualize and analyze such data. I tried to load the chain as particles and managed to do so, but am having problems plotting. I pasted the code below.<br>
<br>
Has anyone used yt for non-spatial data?<br>
<br>
Cheers,<br>
<br>
Carolin<br>
<br>
<br>
-----------<br>
<br>
In [81]: mcchain.keys() ## a dictionary containing dimensions of the chain<br>
Out[81]:<br>
['particle_position_z',<br>
'particle_mass',<br>
'particle_position_x',<br>
'particle_position_y',<br>
'qso']<br>
<br>
In [82]: mcyt = yt.load_particles(mcchain)<br>
yt : [INFO ] 2015-04-23 09:47:39,235 Parameters: current_time = 0.0<br>
yt : [INFO ] 2015-04-23 09:47:39,235 Parameters: domain_dimensions = [2 2 2]<br>
yt : [INFO ] 2015-04-23 09:47:39,237 Parameters: domain_left_edge = [ 0. 0. 0.]<br>
yt : [INFO ] 2015-04-23 09:47:39,237 Parameters: domain_right_edge = [ 1. 1. 1.]<br>
yt : [INFO ] 2015-04-23 09:47:39,238 Parameters: cosmological_simulation = 0.0<br>
<br>
In [83]: circle = mcyt.sphere([0.5, 0.5, 0.5], 0.25)<br>
yt : [INFO ] 2015-04-23 09:47:45,741 Allocating for 1.000e+04 particles<br>
yt : [INFO ] 2015-04-23 09:47:45,763 Identified 7.130e+02 octs<br>
yt : [INFO ] 2015-04-23 09:47:45,897 Loading field plugins.<br>
yt : [INFO ] 2015-04-23 09:47:45,898 Loaded angular_momentum (8 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,898 Loaded astro (15 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,898 Loaded cosmology (22 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,899 Loaded fluid (64 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,900 Loaded fluid_vector (96 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,900 Loaded geometric (112 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,900 Loaded local (112 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,901 Loaded magnetic_field (120 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,901 Loaded my_plugins (120 new fields)<br>
yt : [INFO ] 2015-04-23 09:47:45,901 Loaded species (122 new fields)<br>
<br>
In [84]: pl=yt.ProjectionPlot(circle, 'x', 'y')<br>
---------------------------------------------------------------------------<br>
AttributeError Traceback (most recent call last)<br>
<ipython-input-84-685c784708ec> in <module>()<br>
----> 1 pl=yt.ProjectionPlot(circle, 'x', 'y')<br>
<br>
/Users/cvillforth/Library/Enthought/Canopy_64bit/User/lib/python2.7/site-packages/yt/visualization/plot_window.pyc in __init__(self, ds, axis, fields, center, width, axes_unit, weight_field, max_level, origin, fontsize, field_parameters, data_source, method, proj_style, window_size, aspect)<br>
1212 self.ts = ts<br>
1213 ds = self.ds = ts[0]<br>
-> 1214 axis = fix_axis(axis, ds)<br>
1215 # proj_style is deprecated, but if someone specifies then it trumps<br>
1216 # method.<br>
<br>
/Users/cvillforth/Library/Enthought/Canopy_64bit/User/lib/python2.7/site-packages/yt/funcs.pyc in fix_axis(axis, ds)<br>
649<br>
650 def fix_axis(axis, ds):<br>
--> 651 return ds.coordinates.axis_id.get(axis, axis)<br>
652<br>
653 def get_image_suffix(name):<br>
<br>
AttributeError: 'YTSphere' object has no attribute 'coordinates'<br>
<br>
<br>
--------------------------------------------------------------------------------<br>
Dr Carolin Villforth<br>
SUPA Advanced Fellow<br>
University of St Andrews (SUPA)<br>
School of Physics & Astronomy<br>
North Haugh, St Andrews, Fife<br>
KY16 9SS<br>
UK<br>
<a href="mailto:cv21@st-andrews.ac.uk">cv21@st-andrews.ac.uk</a><br>
Office: 272<br>
phone <a href="tel:%2B44%20%280%29%201334%2046%201612" value="+441334461612">+44 (0) 1334 46 1612</a><br>
<a href="http://www-star.st-and.ac.uk/~cv21/cv21/Home.html" target="_blank">http://www-star.st-and.ac.uk/~cv21/cv21/Home.html</a><br>
--------------------------------------------------------------------------------<br>
<br>
<br>
<br>
<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>
</blockquote></div><br></div>