<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div><br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">On Fri, Jul 13, 2018 at 1:27 PM Macon Abernathy <<a href="mailto:maber001@ucr.edu">maber001@ucr.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hey all, <div><br></div><div>Has anyone had success loading data from larger pixel arrays into Larch? I have data collected on a 100 element Ge detector that results in Athena crashing upon import (running on Mac OS 10.13.4, Demeter v9.25). I figure this is a good task for Larch, but am having a difficult time figuring out how to import all 100 channels for each data set simultaneously. </div><div><br></div><div>Thanks in advance for any tips, </div></div></blockquote><div><br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">What format is the data in?  ​Larch can read very large files, including HDF5 and netCDF files. I don't know anyone using those formats for XAFS data, but I use them all the time for fluorescence maps, where the files can easily get to be 10s of Gb.<br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default"><br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">But most XAFS data can be held in ASCII files, even for 100 element detectors.  If the data is in an ASCII file, something like</div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default"><br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">    fdata = read_ascii('MyDataFile.txt')</div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default"><br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">might just work, and will give you a 2D array in 'fdata.data' that has a shape of  (ncolumns, n_energy_points).  Read_ascii() will also try to parse the column labels, which may work for some data files, but might not work great for all data files. With the full 2D data table you might be able to do something like (totally making up the array indices):<br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">    fdata.energy  = fdata.data[0, :]<br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">    fdata.i0  = fdata[2, :]</div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">    fdata.fluor = fdata.data[4:104, :].sum(axis=0)  # sum columns 5 to 104</div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">    <a href="http://fdata.mu">fdata.mu</a> = fdata.fluor / fdata.i0</div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default"><br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">If something like that doesn't work, send a data file and we'll figure it out.<br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default"><br></div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default">--Matt</div><div style="font-family:arial,helvetica,sans-serif" class="gmail_default"><br></div></div></div>