<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 4, 2016 at 3:26 PM, Scott Feister <span dir="ltr"><<a href="mailto:sfeister@gmail.com" target="_blank">sfeister@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Hi Nathan,<br><br></div>Thanks a lot for the clarifications. I got confused on just one of your sentences:<span class="gmail-"><br><br><i>It will average all the cells along the line of sight at that x,y value. Those cells might have different spatial resolutions.</i><br><br></span>So, a projection doesn't factor in every cell in the simulation, but it takes some subset? So, if there is a single cell in the grid that happens to have an absolutely enormous field value, it's quite possible that this will be under-factored or over-factored in a projection (depending on if the projection rays "misses" or "hits" it)?<br></div></div></blockquote><div><br></div><div>For each pixel, it takes into account all cells that contribute along the line of sight for that pixel.</div><div><br></div><div>If a single cell in the grid has an enormous value, it will be factored in at the same level of weighting as other cells at the same resolution level.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br>When you say the cells might have different resolution, do you mean they would have <i>sub-cells </i>or that they would have different lengths/widths/heights?<br></div></div></blockquote><div><br></div><div>I'm not quite sure what you mean when you refer to "sub-cells" here. Let me see if I can explain in a little bit more detail.</div><div><br></div><div>For a single pixel in a projection, yt figures out a "stack" of cells in your simulation that are at different Z values but also enclose the X and Y values of your pixel.</div><div><br></div><div>To figure out the projection at that pixel, yt computes the sum of cell values times dZ in the "stack" of cells at that X,Y location. Optionally, there might be a weight field here as well, so for a weighted projection it would be the sum of the cell values multiplied by the weight field values multiplied by dZ.</div><div><br></div><div>There's a little bit of text in the docs describing how different projection types are formulated mathematically:</div><div><br></div><div><a href="http://yt-project.org/docs/dev/visualizing/plots.html#types-of-projections">http://yt-project.org/docs/dev/visualizing/plots.html#types-of-projections</a><br></div><div><br></div><div>This is also covered in the yt method paper (section 3.2):</div><div><br></div><div><a href="http://adsabs.harvard.edu/abs/2011ApJS..192....9T">http://adsabs.harvard.edu/abs/2011ApJS..192....9T</a></div><div><br></div><div>Hope that explains it,</div><div><br></div><div>Nathan</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>Thanks a lot!<br><br></div><div>Best,<br><br></div><div>Scott<br></div><div><br><br></div><div class="gmail_extra"><span class="gmail-"><br clear="all"><div><div class="gmail-m_7286382021884814447m_258656256167903431gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><br>Scott Feister, Ph.D.<br>Postdoctoral Researcher, Flash Center for Computational Science<br></div><div>University of Chicago, Department of Astronomy and Astrophysics<br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
<br></span><div><div class="gmail-h5"><div class="gmail_quote">On Tue, Oct 4, 2016 at 2:17 PM, Nathan Goldbaum <span dir="ltr"><<a href="mailto:nathan12343@gmail.com" target="_blank">nathan12343@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span>On Tue, Oct 4, 2016 at 11:26 AM, Scott Feister <span dir="ltr"><<a href="mailto:sfeister@gmail.com" target="_blank">sfeister@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Hi all,<br><br>1) Is there a way to do a "mean" ProjectionPlot? As far as I can tell, the options are "integrate", "mip", and "sum", so I am guessing the answer is no. I know you can do this with regions by reg.mean(), and make your own plot in matplotlib.<br></div></div></div></blockquote><div><br></div></span><div>If you use a weight field you'll get back a mean value (weighted by that field). If you want an unweighted mean, you can weight by the "ones" field (although this will likely not be what you want for AMR data). If you want a mass-weighted mean, you can weight by the cell_mass field. If you want a volume weighted mean, use cell_volume.</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div></div><div>2) If I do a Z projection on a domain with high X,Y spatial resolution, will it average cells along X and Y (within the pixel limits) at each Z depth as it projects? Or, just pick a single cell at X,Y for each pixel?<br></div></div></div></blockquote><div><br></div></span><div>It will average all the cells along the line of sight at that x,y value. Those cells might have different spatial resolutions.</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div></div><div dir="ltr"><br></div><div>Also, does anyone know if the reg.mean() function weights the cells equally in the mean if there is varied spatial resolution? For example, I know the reg.integrate() function multiplies by path length, but reg.sum() simply adds up cell values.<br></div></div></div></blockquote><div><br></div></span><div>The mean() function takes a field to weight by:</div><div><br></div><div><a href="http://yt-project.org/docs/dev/reference/api/generated/yt.data_objects.selection_data_containers.YTRegion.mean.html#yt.data_objects.selection_data_containers.YTRegion.mean" target="_blank">http://yt-project.org/docs/dev<wbr>/reference/api/generated/yt.da<wbr>ta_objects.selection_data_cont<wbr>ainers.YTRegion.mean.html#yt.<wbr>data_objects.selection_data_<wbr>containers.YTRegion.mean</a><br></div><div><br></div><div>By default it uses "ones" or "particle_ones", depending on whether you're averaging a particle or mesh field.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span><div dir="ltr"><div><div><br></div><div>Thanks!<br><br></div><div dir="ltr">Scott<br><br>Scott Feister, Ph.D.<br>Postdoctoral Researcher, Flash Center for Computational Science<br></div></div><div><div><div><div><div class="gmail-m_7286382021884814447m_258656256167903431m_5604617467800397975gmail-m_669731357677844511gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div>University of Chicago, Department of Astronomy and Astrophysics<br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div></div>
<br></span>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br></blockquote></div><br></div></div>
<br>______________________________<wbr>_________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/lis<wbr>tinfo.cgi/yt-users-spacepope.o<wbr>rg</a><br>
<br></blockquote></div><br></div></div></div></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br></blockquote></div><br></div></div>