<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 22, 2016 at 11:23 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-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Wed, Jun 22, 2016 at 2:03 PM, Hansheng Chen <span dir="ltr"><<a href="mailto:jasonhs1221@gmail.com" target="_blank">jasonhs1221@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear yt-users,<div><br></div><div>I got a image of projected gas density of a galaxy, and from this projected image, I could get a image array of projected gas density with a resolution of (800, 800). Now, I would like to set the center at (400, 400) and to get a profile that shows the relationship between the projected gas density and radius in 2-D. Could some one give me any help in getting a radial profile from the image array?</div></div></blockquote><div><br></div></span><div>There isn't anything in yt 3.2.3 that does this directly, although one can straightforwardly accomplish this using numpy's histogram function if you bin pixels according to radius.</div></div></div></div></blockquote><div><br></div><div>Although, if you're willing to install the development version of yt, you can use the ytdata frontend to save a projection data object and reload it as a dataset:</div><div><br></div><div><a href="http://yt-project.org/docs/dev/analyzing/saving_data.html#spatial-plots">http://yt-project.org/docs/dev/analyzing/saving_data.html#spatial-plots</a><br></div><div><br></div><div>Here's a small example that illustrates this workflow using a test dataset from <a href="http://yt-project.org/data">yt-project.org/data</a>:</div><div><br></div><div><div> import yt</div><div><br></div><div> ds = yt.load('IsolatedGalaxy/galaxy0030/galaxy0030')</div><div> prj = ds.proj(['density', 'radius'], 0)</div><div> prj.save_as_dataset()</div><div><br></div><div> ds2 = yt.load('galaxy0030_proj.h5')</div><div> plot = yt.ProfilePlot(ds2.all_data(), 'radius', 'density', weight_field=None)</div></div><div> plot.save()</div><div><br></div><div>Which produces the following image:</div><div><br></div><div><a href="http://i.imgur.com/L0nxEyy.png">http://i.imgur.com/L0nxEyy.png</a><br></div><div><br></div><div>-Nathan</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>Thank you!</div><div><br></div><div>Sincerely,</div><div>Hansheng</div></div>
<br>_______________________________________________<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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div></div>
</blockquote></div><br></div></div>