[yt-users] Projection Question

[Nathan Goldbaum]

Hi Dave,

It looks like the plot is crashing because of some issue with the color bar - something like all the values in the plot being zero, causing the logarithmic color bar to crash.  I'd need to know more about the dataset you're trying to plot to be able to figure out the exact issue.

This should in principle work, both with the data source and without.  One possible issue with the data source is that you're creating the region and then re-loading the parameter file.  Are you sure you want to do that?

I think a more modern script would look like this:

pf = load(thing)
reg = pf.h.region(center,left,right)
proj =  pf.h.proj(0,'Density',center=Center, source=reg)
prooj = proj.to_pw()
prooj.save()

This leverages the to_pw() function, which is defined for projection, slice, and off-axis slice data objects.  The function returns a plot window plot that uses the user-defined data object to select data.  By default plot window plots use a region that covers the whole dataset to select data.

-Nathan

On Apr 8, 2013, at 4:07 PM, david collins <dcollins4096 at gmail.com> wrote:

> Hi, Everybody!
> 
> I have a probably dumb question.  I keep getting the following error when projecting.  I get this from a variety of methods, but what I'm trying to do is project a region.  What I used to do was something like this, but it's no longer working (with or without the source.)
> 
> reg = pf.h.region(center,left,right)
> pf = load(thing)
> pc = PlotCollection(pf)
> proj =  pf.h.proj(0,'Density',center=Center) #,source=reg)
> prooj = pc.add_projection('Density','x') #,data_source=reg)
> 
> So my question is: with the recent developments to yt, what's the best mechanism to do this?  
> I'm on a self install that was done fresh pretty recently, the hash is 34b95297062b
> 
> Thanks,
> d.
> 
> 
>     prooj = pc.add_projection('Density','x',data_source=reg)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/visualization/plot_collection.py", line 786, in add_projection
>     size=fig_size, periodic=periodic))
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/visualization/plot_types.py", line 326, in __init__
>     self.__init_temp_image(use_colorbar)
>   File "/nics/b/home/collins/local_2013_03_12/src/yt-hg/yt/visualization/plot_types.py", line 358, in __init_temp_image
>     shrink=0.95)
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/figure.py", line 1387, in colorbar
>     cb = cbar.colorbar_factory(cax, mappable, **kw)
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 1173, in colorbar_factory
>     cb = Colorbar(cax, mappable, **kwargs)
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 854, in __init__
>     ColorbarBase.__init__(self, ax, **kw)
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 298, in __init__
>     self.draw_all()
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 323, in draw_all
>     self._config_axes(X, Y)
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 409, in _config_axes
>     self.update_ticks()
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 345, in update_ticks
>     ticks, ticklabels, offset_string = self._ticker()
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 560, in _ticker
>     ticks = self._locate(b)
>   File "/nics/b/home/collins/local_2013_03_12/lib/python2.7/site-packages/matplotlib/colorbar.py", line 808, in _locate
>     z = np.take(y, i0) + (xn-np.take(b,i0))*dy/db
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org