[yt-users] Projection Question

[Nathan Goldbaum]

In analogy to the to_pw() function, there is also to_frb(), which will return a FixedResolutionBuffer object that you can use to take a direct look at the data: http://paste.yt-project.org/show/3349/

Note this uses the enzo IsolatedGalaxy test dataset available on yt-project.org/data.  On my machine the output from this script looks like this: http://paste.yt-project.org/show/3347/

Hope that's useful and that you're able to figure out what's going wrong.

-Nathan


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

> 
> 
> 
> 
> 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.
> 
> 
> Yeah, it looks like what happens when I have a bad dataset, like Density = 0 everywhere.  However, I've plotted (and published) this dataset many times, so unless it's been corrupted on disk the set itself is ok.  To check that I did 
> for grid in pf.h.grids:
>  print grid['Density'].min(), grid['Density'].max()
> 
> and it didn't barf-- I'll double check those values, though.
> 
> I think a more modern script would look like this:
> 
> 
> Thanks a ton for the update.  I'll try this, and ensure the data isn't borked.
> 
> Thanks a ton!
> d.
>  
> 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
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20130408/a868b0c8/attachment.html>