source: sasview/_modules/sas/dataloader/readers/danse_reader.html @ a462c6a

gh-pages
Last change on this file since a462c6a was a462c6a, checked in by ajj, 9 years ago

Rebuild to fix index and modules docs
  • Property mode set to 100644
File size: 44.4 KB
RevLine 
[a462c6a]1<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
2  "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
3
4
5<html xmlns="http://www.w3.org/1999/xhtml">
6  <head>
7    <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
8   
9    <title>sas.dataloader.readers.danse_reader &mdash; SasView 3.0.0 documentation</title>
10   
11    <link rel="stylesheet" href="../../../../_static/default.css" type="text/css" />
12    <link rel="stylesheet" href="../../../../_static/pygments.css" type="text/css" />
13   
14    <script type="text/javascript">
15      var DOCUMENTATION_OPTIONS = {
16        URL_ROOT:    '../../../../',
17        VERSION:     '3.0.0',
18        COLLAPSE_INDEX: false,
19        FILE_SUFFIX: '.html',
20        HAS_SOURCE:  true
21      };
22    </script>
23    <script type="text/javascript" src="../../../../_static/jquery.js"></script>
24    <script type="text/javascript" src="../../../../_static/underscore.js"></script>
25    <script type="text/javascript" src="../../../../_static/doctools.js"></script>
26    <script type="text/javascript" src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
27    <link rel="top" title="SasView 3.0.0 documentation" href="../../../../index.html" />
28    <link rel="up" title="sas.dataloader.readers" href="../readers.html" /> 
29  </head>
30  <body>
31    <div class="related">
32      <h3>Navigation</h3>
33      <ul>
34        <li class="right" style="margin-right: 10px">
35          <a href="../../../../genindex.html" title="General Index"
36             accesskey="I">index</a></li>
37        <li class="right" >
38          <a href="../../../../py-modindex.html" title="Python Module Index"
39             >modules</a> |</li>
40        <li><a href="../../../../index.html">SasView 3.0.0 documentation</a> &raquo;</li>
41          <li><a href="../../../index.html" >Module code</a> &raquo;</li>
42          <li><a href="../readers.html" accesskey="U">sas.dataloader.readers</a> &raquo;</li> 
43      </ul>
44    </div> 
45
46    <div class="document">
47      <div class="documentwrapper">
48        <div class="bodywrapper">
49          <div class="body">
50           
51  <h1>Source code for sas.dataloader.readers.danse_reader</h1><div class="highlight"><pre>
52<span class="sd">&quot;&quot;&quot;</span>
53<span class="sd">    DANSE/SANS file reader</span>
54<span class="sd">&quot;&quot;&quot;</span>
55<span class="c">############################################################################</span>
56<span class="c">#This software was developed by the University of Tennessee as part of the</span>
57<span class="c">#Distributed Data Analysis of Neutron Scattering Experiments (DANSE)</span>
58<span class="c">#project funded by the US National Science Foundation. </span>
59<span class="c">#If you use DANSE applications to do scientific research that leads to</span>
60<span class="c">#publication, we ask that you acknowledge the use of the software with the</span>
61<span class="c">#following sentence:</span>
62<span class="c">#This work benefited from DANSE software developed under NSF award DMR-0520547.</span>
63<span class="c">#copyright 2008, University of Tennessee</span>
64<span class="c">#############################################################################</span>
65<span class="kn">import</span> <span class="nn">math</span>
66<span class="kn">import</span> <span class="nn">os</span>
67<span class="kn">import</span> <span class="nn">sys</span>
68<span class="kn">import</span> <span class="nn">numpy</span>
69<span class="kn">import</span> <span class="nn">logging</span>
70<span class="kn">from</span> <span class="nn">sas.dataloader.data_info</span> <span class="kn">import</span> <span class="n">Data2D</span><span class="p">,</span> <span class="n">Detector</span>
71<span class="kn">from</span> <span class="nn">sas.dataloader.manipulations</span> <span class="kn">import</span> <span class="n">reader2D_converter</span>
72
73<span class="c"># Look for unit converter</span>
74<span class="n">has_converter</span> <span class="o">=</span> <span class="bp">True</span>
75<span class="k">try</span><span class="p">:</span>
76    <span class="kn">from</span> <span class="nn">sas.data_util.nxsunit</span> <span class="kn">import</span> <span class="n">Converter</span>
77<span class="k">except</span><span class="p">:</span>
78    <span class="n">has_converter</span> <span class="o">=</span> <span class="bp">False</span>
79
80
81<div class="viewcode-block" id="Reader"><a class="viewcode-back" href="../../../../dev/api/sas.dataloader.readers.html#sas.dataloader.readers.danse_reader.Reader">[docs]</a><span class="k">class</span> <span class="nc">Reader</span><span class="p">:</span>
82    <span class="sd">&quot;&quot;&quot;</span>
83<span class="sd">    Example data manipulation</span>
84<span class="sd">    &quot;&quot;&quot;</span>
85    <span class="c">## File type</span>
86    <span class="n">type_name</span> <span class="o">=</span> <span class="s">&quot;DANSE&quot;</span>
87    <span class="c">## Wildcards</span>
88    <span class="nb">type</span> <span class="o">=</span> <span class="p">[</span><span class="s">&quot;DANSE files (*.sans)|*.sans&quot;</span><span class="p">]</span>
89    <span class="c">## Extension</span>
90    <span class="n">ext</span>  <span class="o">=</span> <span class="p">[</span><span class="s">&#39;.sans&#39;</span><span class="p">,</span> <span class="s">&#39;.SANS&#39;</span><span class="p">]</span>
91       
92<div class="viewcode-block" id="Reader.read"><a class="viewcode-back" href="../../../../dev/api/sas.dataloader.readers.html#sas.dataloader.readers.danse_reader.Reader.read">[docs]</a>    <span class="k">def</span> <span class="nf">read</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">filename</span><span class="o">=</span><span class="bp">None</span><span class="p">):</span>
93        <span class="sd">&quot;&quot;&quot;</span>
94<span class="sd">        Open and read the data in a file</span>
95<span class="sd">        @param file: path of the file</span>
96<span class="sd">        &quot;&quot;&quot;</span>
97       
98        <span class="n">read_it</span> <span class="o">=</span> <span class="bp">False</span>
99        <span class="k">for</span> <span class="n">item</span> <span class="ow">in</span> <span class="bp">self</span><span class="o">.</span><span class="n">ext</span><span class="p">:</span>
100            <span class="k">if</span> <span class="n">filename</span><span class="o">.</span><span class="n">lower</span><span class="p">()</span><span class="o">.</span><span class="n">find</span><span class="p">(</span><span class="n">item</span><span class="p">)</span> <span class="o">&gt;=</span> <span class="mi">0</span><span class="p">:</span>
101                <span class="n">read_it</span> <span class="o">=</span> <span class="bp">True</span>
102               
103        <span class="k">if</span> <span class="n">read_it</span><span class="p">:</span>
104            <span class="k">try</span><span class="p">:</span>
105                <span class="n">datafile</span> <span class="o">=</span> <span class="nb">open</span><span class="p">(</span><span class="n">filename</span><span class="p">,</span> <span class="s">&#39;r&#39;</span><span class="p">)</span>
106            <span class="k">except</span><span class="p">:</span>
107                <span class="k">raise</span>  <span class="ne">RuntimeError</span><span class="p">,</span><span class="s">&quot;danse_reader cannot open </span><span class="si">%s</span><span class="s">&quot;</span> <span class="o">%</span> <span class="p">(</span><span class="n">filename</span><span class="p">)</span>
108       
109            <span class="c"># defaults</span>
110            <span class="c"># wavelength in Angstrom</span>
111            <span class="n">wavelength</span> <span class="o">=</span> <span class="mf">10.0</span>
112            <span class="c"># Distance in meter</span>
113            <span class="n">distance</span>   <span class="o">=</span> <span class="mf">11.0</span>
114            <span class="c"># Pixel number of center in x</span>
115            <span class="n">center_x</span>   <span class="o">=</span> <span class="mi">65</span>
116            <span class="c"># Pixel number of center in y</span>
117            <span class="n">center_y</span>   <span class="o">=</span> <span class="mi">65</span>
118            <span class="c"># Pixel size [mm]</span>
119            <span class="n">pixel</span>      <span class="o">=</span> <span class="mf">5.0</span>
120            <span class="c"># Size in x, in pixels</span>
121            <span class="n">size_x</span>     <span class="o">=</span> <span class="mi">128</span>
122            <span class="c"># Size in y, in pixels</span>
123            <span class="n">size_y</span>     <span class="o">=</span> <span class="mi">128</span>
124            <span class="c"># Format version</span>
125            <span class="n">fversion</span>   <span class="o">=</span> <span class="mf">1.0</span>
126           
127            <span class="n">output</span> <span class="o">=</span> <span class="n">Data2D</span><span class="p">()</span>
128            <span class="n">output</span><span class="o">.</span><span class="n">filename</span> <span class="o">=</span> <span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">basename</span><span class="p">(</span><span class="n">filename</span><span class="p">)</span>
129            <span class="n">detector</span> <span class="o">=</span> <span class="n">Detector</span><span class="p">()</span>
130            <span class="n">output</span><span class="o">.</span><span class="n">detector</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">detector</span><span class="p">)</span>
131           
132            <span class="n">output</span><span class="o">.</span><span class="n">data</span> <span class="o">=</span> <span class="n">numpy</span><span class="o">.</span><span class="n">zeros</span><span class="p">([</span><span class="n">size_x</span><span class="p">,</span><span class="n">size_y</span><span class="p">])</span>
133            <span class="n">output</span><span class="o">.</span><span class="n">err_data</span> <span class="o">=</span> <span class="n">numpy</span><span class="o">.</span><span class="n">zeros</span><span class="p">([</span><span class="n">size_x</span><span class="p">,</span> <span class="n">size_y</span><span class="p">])</span>
134           
135            <span class="n">data_conv_q</span> <span class="o">=</span> <span class="bp">None</span>
136            <span class="n">data_conv_i</span> <span class="o">=</span> <span class="bp">None</span>
137           
138            <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span> <span class="o">!=</span> <span class="s">&#39;1/A&#39;</span><span class="p">:</span>
139                <span class="n">data_conv_q</span> <span class="o">=</span> <span class="n">Converter</span><span class="p">(</span><span class="s">&#39;1/A&#39;</span><span class="p">)</span>
140                <span class="c"># Test it</span>
141                <span class="n">data_conv_q</span><span class="p">(</span><span class="mf">1.0</span><span class="p">,</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
142               
143            <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">output</span><span class="o">.</span><span class="n">I_unit</span> <span class="o">!=</span> <span class="s">&#39;1/cm&#39;</span><span class="p">:</span>
144                <span class="n">data_conv_i</span> <span class="o">=</span> <span class="n">Converter</span><span class="p">(</span><span class="s">&#39;1/cm&#39;</span><span class="p">)</span>
145                <span class="c"># Test it</span>
146                <span class="n">data_conv_i</span><span class="p">(</span><span class="mf">1.0</span><span class="p">,</span> <span class="n">output</span><span class="o">.</span><span class="n">I_unit</span><span class="p">)</span>
147       
148            <span class="n">read_on</span> <span class="o">=</span> <span class="bp">True</span>
149            <span class="k">while</span> <span class="n">read_on</span><span class="p">:</span>
150                <span class="n">line</span> <span class="o">=</span> <span class="n">datafile</span><span class="o">.</span><span class="n">readline</span><span class="p">()</span>
151                <span class="k">if</span> <span class="n">line</span><span class="o">.</span><span class="n">find</span><span class="p">(</span><span class="s">&quot;DATA:&quot;</span><span class="p">)</span> <span class="o">&gt;=</span> <span class="mi">0</span><span class="p">:</span>
152                    <span class="n">read_on</span> <span class="o">=</span> <span class="bp">False</span>
153                    <span class="k">break</span>
154                <span class="n">toks</span> <span class="o">=</span> <span class="n">line</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s">&#39;:&#39;</span><span class="p">)</span>
155                <span class="k">if</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;FORMATVERSION&quot;</span><span class="p">:</span>
156                    <span class="n">fversion</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
157                <span class="k">if</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;WAVELENGTH&quot;</span><span class="p">:</span>
158                    <span class="n">wavelength</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
159                <span class="k">elif</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;DISTANCE&quot;</span><span class="p">:</span>
160                    <span class="n">distance</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
161                <span class="k">elif</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;CENTER_X&quot;</span><span class="p">:</span>
162                    <span class="n">center_x</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
163                <span class="k">elif</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;CENTER_Y&quot;</span><span class="p">:</span>
164                    <span class="n">center_y</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
165                <span class="k">elif</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;PIXELSIZE&quot;</span><span class="p">:</span>
166                    <span class="n">pixel</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
167                <span class="k">elif</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;SIZE_X&quot;</span><span class="p">:</span>
168                    <span class="n">size_x</span> <span class="o">=</span> <span class="nb">int</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
169                <span class="k">elif</span> <span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span> <span class="o">==</span> <span class="s">&quot;SIZE_Y&quot;</span><span class="p">:</span>
170                    <span class="n">size_y</span> <span class="o">=</span> <span class="nb">int</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
171           
172            <span class="c"># Read the data</span>
173            <span class="n">data</span> <span class="o">=</span> <span class="p">[]</span>
174            <span class="n">error</span> <span class="o">=</span> <span class="p">[]</span>
175            <span class="k">if</span> <span class="n">fversion</span> <span class="o">==</span> <span class="mf">1.0</span><span class="p">:</span>
176                <span class="n">data_str</span> <span class="o">=</span> <span class="n">datafile</span><span class="o">.</span><span class="n">readline</span><span class="p">()</span>
177                <span class="n">data</span> <span class="o">=</span> <span class="n">data_str</span><span class="o">.</span><span class="n">split</span><span class="p">(</span><span class="s">&#39; &#39;</span><span class="p">)</span>
178            <span class="k">else</span><span class="p">:</span>
179                <span class="n">read_on</span> <span class="o">=</span> <span class="bp">True</span>
180                <span class="k">while</span> <span class="n">read_on</span><span class="p">:</span>
181                    <span class="n">data_str</span> <span class="o">=</span> <span class="n">datafile</span><span class="o">.</span><span class="n">readline</span><span class="p">()</span>
182                    <span class="k">if</span> <span class="nb">len</span><span class="p">(</span><span class="n">data_str</span><span class="p">)</span> <span class="o">==</span> <span class="mi">0</span><span class="p">:</span>
183                        <span class="n">read_on</span> <span class="o">=</span> <span class="bp">False</span>
184                    <span class="k">else</span><span class="p">:</span>
185                        <span class="n">toks</span> <span class="o">=</span> <span class="n">data_str</span><span class="o">.</span><span class="n">split</span><span class="p">()</span>
186                        <span class="k">try</span><span class="p">:</span>
187                            <span class="n">val</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">0</span><span class="p">])</span>
188                            <span class="n">err</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">toks</span><span class="p">[</span><span class="mi">1</span><span class="p">])</span>
189                            <span class="k">if</span> <span class="n">data_conv_i</span> <span class="ow">is</span> <span class="ow">not</span> <span class="bp">None</span><span class="p">:</span>
190                                <span class="n">val</span> <span class="o">=</span> <span class="n">data_conv_i</span><span class="p">(</span><span class="n">val</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">_yunit</span><span class="p">)</span>
191                                <span class="n">err</span> <span class="o">=</span> <span class="n">data_conv_i</span><span class="p">(</span><span class="n">err</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">_yunit</span><span class="p">)</span>
192                            <span class="n">data</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">val</span><span class="p">)</span>
193                            <span class="n">error</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">err</span><span class="p">)</span>
194                        <span class="k">except</span><span class="p">:</span>
195                            <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s">&quot;Skipping line:</span><span class="si">%s</span><span class="s">,</span><span class="si">%s</span><span class="s">&quot;</span> <span class="o">%</span><span class="p">(</span><span class="n">data_str</span><span class="p">,</span>
196                                                                <span class="n">sys</span><span class="o">.</span><span class="n">exc_value</span><span class="p">))</span>
197           
198            <span class="c"># Initialize</span>
199            <span class="n">x_vals</span> <span class="o">=</span> <span class="p">[]</span>
200            <span class="n">y_vals</span> <span class="o">=</span> <span class="p">[]</span>
201            <span class="n">ymin</span> <span class="o">=</span> <span class="bp">None</span>
202            <span class="n">ymax</span> <span class="o">=</span> <span class="bp">None</span>
203            <span class="n">xmin</span> <span class="o">=</span> <span class="bp">None</span>
204            <span class="n">xmax</span> <span class="o">=</span> <span class="bp">None</span>
205           
206            <span class="c"># Qx and Qy vectors</span>
207            <span class="n">theta</span> <span class="o">=</span> <span class="n">pixel</span> <span class="o">/</span> <span class="n">distance</span> <span class="o">/</span> <span class="mf">100.0</span>
208            <span class="n">stepq</span> <span class="o">=</span> <span class="mf">4.0</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span> <span class="o">/</span> <span class="n">wavelength</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">sin</span><span class="p">(</span><span class="n">theta</span> <span class="o">/</span> <span class="mf">2.0</span><span class="p">)</span>
209            <span class="k">for</span> <span class="n">i_x</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">size_x</span><span class="p">):</span>
210                <span class="n">theta</span> <span class="o">=</span> <span class="p">(</span><span class="n">i_x</span> <span class="o">-</span> <span class="n">center_x</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="n">pixel</span> <span class="o">/</span> <span class="n">distance</span> <span class="o">/</span> <span class="mf">100.0</span>
211                <span class="n">qx</span> <span class="o">=</span> <span class="mf">4.0</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span> <span class="o">/</span> <span class="n">wavelength</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">sin</span><span class="p">(</span><span class="n">theta</span> <span class="o">/</span> <span class="mf">2.0</span><span class="p">)</span>
212               
213                <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span> <span class="o">!=</span> <span class="s">&#39;1/A&#39;</span><span class="p">:</span>
214                    <span class="n">qx</span> <span class="o">=</span> <span class="n">data_conv_q</span><span class="p">(</span><span class="n">qx</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
215               
216                <span class="n">x_vals</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">qx</span><span class="p">)</span>
217                <span class="k">if</span> <span class="n">xmin</span> <span class="o">==</span> <span class="bp">None</span> <span class="ow">or</span> <span class="n">qx</span> <span class="o">&lt;</span> <span class="n">xmin</span><span class="p">:</span>
218                    <span class="n">xmin</span> <span class="o">=</span> <span class="n">qx</span>
219                <span class="k">if</span> <span class="n">xmax</span> <span class="o">==</span> <span class="bp">None</span> <span class="ow">or</span> <span class="n">qx</span> <span class="o">&gt;</span> <span class="n">xmax</span><span class="p">:</span>
220                    <span class="n">xmax</span> <span class="o">=</span> <span class="n">qx</span>
221           
222            <span class="n">ymin</span> <span class="o">=</span> <span class="bp">None</span>
223            <span class="n">ymax</span> <span class="o">=</span> <span class="bp">None</span>
224            <span class="k">for</span> <span class="n">i_y</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">size_y</span><span class="p">):</span>
225                <span class="n">theta</span> <span class="o">=</span> <span class="p">(</span><span class="n">i_y</span> <span class="o">-</span> <span class="n">center_y</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="n">pixel</span> <span class="o">/</span> <span class="n">distance</span> <span class="o">/</span> <span class="mf">100.0</span>
226                <span class="n">qy</span> <span class="o">=</span> <span class="mf">4.0</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span> <span class="o">/</span> <span class="n">wavelength</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">sin</span><span class="p">(</span><span class="n">theta</span><span class="o">/</span><span class="mf">2.0</span><span class="p">)</span>
227               
228                <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span> <span class="o">!=</span> <span class="s">&#39;1/A&#39;</span><span class="p">:</span>
229                    <span class="n">qy</span> <span class="o">=</span> <span class="n">data_conv_q</span><span class="p">(</span><span class="n">qy</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
230               
231                <span class="n">y_vals</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">qy</span><span class="p">)</span>
232                <span class="k">if</span> <span class="n">ymin</span> <span class="o">==</span> <span class="bp">None</span> <span class="ow">or</span> <span class="n">qy</span> <span class="o">&lt;</span> <span class="n">ymin</span><span class="p">:</span>
233                    <span class="n">ymin</span> <span class="o">=</span> <span class="n">qy</span>
234                <span class="k">if</span> <span class="n">ymax</span> <span class="o">==</span> <span class="bp">None</span> <span class="ow">or</span> <span class="n">qy</span> <span class="o">&gt;</span> <span class="n">ymax</span><span class="p">:</span>
235                    <span class="n">ymax</span> <span class="o">=</span> <span class="n">qy</span>
236           
237            <span class="c"># Store the data in the 2D array</span>
238            <span class="n">i_x</span> <span class="o">=</span> <span class="mi">0</span>
239            <span class="n">i_y</span> <span class="o">=</span> <span class="o">-</span><span class="mi">1</span>
240           
241            <span class="k">for</span> <span class="n">i_pt</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="nb">len</span><span class="p">(</span><span class="n">data</span><span class="p">)):</span>
242                <span class="k">try</span><span class="p">:</span>
243                    <span class="n">value</span> <span class="o">=</span> <span class="nb">float</span><span class="p">(</span><span class="n">data</span><span class="p">[</span><span class="n">i_pt</span><span class="p">])</span>
244                <span class="k">except</span><span class="p">:</span>
245                    <span class="c"># For version 1.0, the data were still</span>
246                    <span class="c"># stored as strings at this point.</span>
247                    <span class="n">msg</span> <span class="o">=</span> <span class="s">&quot;Skipping entry (v1.0):</span><span class="si">%s</span><span class="s">,</span><span class="si">%s</span><span class="s">&quot;</span> <span class="o">%</span> <span class="p">(</span><span class="nb">str</span><span class="p">(</span><span class="n">data</span><span class="p">[</span><span class="n">i_pt</span><span class="p">]),</span>
248                                                           <span class="n">sys</span><span class="o">.</span><span class="n">exc_value</span><span class="p">)</span>
249                    <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="n">msg</span><span class="p">)</span>
250               
251                <span class="c"># Get bin number</span>
252                <span class="k">if</span> <span class="n">math</span><span class="o">.</span><span class="n">fmod</span><span class="p">(</span><span class="n">i_pt</span><span class="p">,</span> <span class="n">size_x</span><span class="p">)</span> <span class="o">==</span> <span class="mi">0</span><span class="p">:</span>
253                    <span class="n">i_x</span> <span class="o">=</span> <span class="mi">0</span>
254                    <span class="n">i_y</span> <span class="o">+=</span> <span class="mi">1</span>
255                <span class="k">else</span><span class="p">:</span>
256                    <span class="n">i_x</span> <span class="o">+=</span> <span class="mi">1</span>
257                   
258                <span class="n">output</span><span class="o">.</span><span class="n">data</span><span class="p">[</span><span class="n">i_y</span><span class="p">][</span><span class="n">i_x</span><span class="p">]</span> <span class="o">=</span> <span class="n">value</span>
259                <span class="k">if</span> <span class="n">fversion</span><span class="o">&gt;</span><span class="mf">1.0</span><span class="p">:</span>
260                    <span class="n">output</span><span class="o">.</span><span class="n">err_data</span><span class="p">[</span><span class="n">i_y</span><span class="p">][</span><span class="n">i_x</span><span class="p">]</span> <span class="o">=</span> <span class="n">error</span><span class="p">[</span><span class="n">i_pt</span><span class="p">]</span>
261               
262            <span class="c"># Store all data</span>
263            <span class="c"># Store wavelength</span>
264            <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">output</span><span class="o">.</span><span class="n">source</span><span class="o">.</span><span class="n">wavelength_unit</span> <span class="o">!=</span> <span class="s">&#39;A&#39;</span><span class="p">:</span>
265                <span class="n">conv</span> <span class="o">=</span> <span class="n">Converter</span><span class="p">(</span><span class="s">&#39;A&#39;</span><span class="p">)</span>
266                <span class="n">wavelength</span> <span class="o">=</span> <span class="n">conv</span><span class="p">(</span><span class="n">wavelength</span><span class="p">,</span>
267                                  <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">source</span><span class="o">.</span><span class="n">wavelength_unit</span><span class="p">)</span>
268            <span class="n">output</span><span class="o">.</span><span class="n">source</span><span class="o">.</span><span class="n">wavelength</span> <span class="o">=</span> <span class="n">wavelength</span>
269               
270            <span class="c"># Store distance</span>
271            <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">detector</span><span class="o">.</span><span class="n">distance_unit</span> <span class="o">!=</span> <span class="s">&#39;m&#39;</span><span class="p">:</span>
272                <span class="n">conv</span> <span class="o">=</span> <span class="n">Converter</span><span class="p">(</span><span class="s">&#39;m&#39;</span><span class="p">)</span>
273                <span class="n">distance</span> <span class="o">=</span> <span class="n">conv</span><span class="p">(</span><span class="n">distance</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">detector</span><span class="o">.</span><span class="n">distance_unit</span><span class="p">)</span>
274            <span class="n">detector</span><span class="o">.</span><span class="n">distance</span> <span class="o">=</span> <span class="n">distance</span>
275           
276            <span class="c"># Store pixel size</span>
277            <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">detector</span><span class="o">.</span><span class="n">pixel_size_unit</span> <span class="o">!=</span> <span class="s">&#39;mm&#39;</span><span class="p">:</span>
278                <span class="n">conv</span> <span class="o">=</span> <span class="n">Converter</span><span class="p">(</span><span class="s">&#39;mm&#39;</span><span class="p">)</span>
279                <span class="n">pixel</span> <span class="o">=</span> <span class="n">conv</span><span class="p">(</span><span class="n">pixel</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">detector</span><span class="o">.</span><span class="n">pixel_size_unit</span><span class="p">)</span>
280            <span class="n">detector</span><span class="o">.</span><span class="n">pixel_size</span><span class="o">.</span><span class="n">x</span> <span class="o">=</span> <span class="n">pixel</span>
281            <span class="n">detector</span><span class="o">.</span><span class="n">pixel_size</span><span class="o">.</span><span class="n">y</span> <span class="o">=</span> <span class="n">pixel</span>
282
283            <span class="c"># Store beam center in distance units</span>
284            <span class="n">detector</span><span class="o">.</span><span class="n">beam_center</span><span class="o">.</span><span class="n">x</span> <span class="o">=</span> <span class="n">center_x</span> <span class="o">*</span> <span class="n">pixel</span>
285            <span class="n">detector</span><span class="o">.</span><span class="n">beam_center</span><span class="o">.</span><span class="n">y</span> <span class="o">=</span> <span class="n">center_y</span> <span class="o">*</span> <span class="n">pixel</span>
286           
287            <span class="c"># Store limits of the image (2D array)</span>
288            <span class="n">xmin</span> <span class="o">=</span> <span class="n">xmin</span> <span class="o">-</span> <span class="n">stepq</span> <span class="o">/</span> <span class="mf">2.0</span>
289            <span class="n">xmax</span> <span class="o">=</span> <span class="n">xmax</span> <span class="o">+</span> <span class="n">stepq</span> <span class="o">/</span> <span class="mf">2.0</span>
290            <span class="n">ymin</span> <span class="o">=</span> <span class="n">ymin</span> <span class="o">-</span> <span class="n">stepq</span> <span class="o">/</span><span class="mf">2.0</span>
291            <span class="n">ymax</span> <span class="o">=</span> <span class="n">ymax</span> <span class="o">+</span> <span class="n">stepq</span> <span class="o">/</span> <span class="mf">2.0</span>
292           
293            <span class="k">if</span> <span class="n">has_converter</span> <span class="o">==</span> <span class="bp">True</span> <span class="ow">and</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span> <span class="o">!=</span> <span class="s">&#39;1/A&#39;</span><span class="p">:</span>
294                <span class="n">xmin</span> <span class="o">=</span> <span class="n">data_conv_q</span><span class="p">(</span><span class="n">xmin</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
295                <span class="n">xmax</span> <span class="o">=</span> <span class="n">data_conv_q</span><span class="p">(</span><span class="n">xmax</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
296                <span class="n">ymin</span> <span class="o">=</span> <span class="n">data_conv_q</span><span class="p">(</span><span class="n">ymin</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
297                <span class="n">ymax</span> <span class="o">=</span> <span class="n">data_conv_q</span><span class="p">(</span><span class="n">ymax</span><span class="p">,</span> <span class="n">units</span><span class="o">=</span><span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
298            <span class="n">output</span><span class="o">.</span><span class="n">xmin</span> <span class="o">=</span> <span class="n">xmin</span>
299            <span class="n">output</span><span class="o">.</span><span class="n">xmax</span> <span class="o">=</span> <span class="n">xmax</span>
300            <span class="n">output</span><span class="o">.</span><span class="n">ymin</span> <span class="o">=</span> <span class="n">ymin</span>
301            <span class="n">output</span><span class="o">.</span><span class="n">ymax</span> <span class="o">=</span> <span class="n">ymax</span>
302           
303            <span class="c"># Store x and y axis bin centers</span>
304            <span class="n">output</span><span class="o">.</span><span class="n">x_bins</span> <span class="o">=</span> <span class="n">x_vals</span>
305            <span class="n">output</span><span class="o">.</span><span class="n">y_bins</span> <span class="o">=</span> <span class="n">y_vals</span>
306           
307            <span class="c"># Units</span>
308            <span class="k">if</span> <span class="n">data_conv_q</span> <span class="ow">is</span> <span class="ow">not</span> <span class="bp">None</span><span class="p">:</span>
309                <span class="n">output</span><span class="o">.</span><span class="n">xaxis</span><span class="p">(</span><span class="s">&quot;</span><span class="se">\\</span><span class="s">rm{Q_{x}}&quot;</span><span class="p">,</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
310                <span class="n">output</span><span class="o">.</span><span class="n">yaxis</span><span class="p">(</span><span class="s">&quot;</span><span class="se">\\</span><span class="s">rm{Q_{y}}&quot;</span><span class="p">,</span> <span class="n">output</span><span class="o">.</span><span class="n">Q_unit</span><span class="p">)</span>
311            <span class="k">else</span><span class="p">:</span>
312                <span class="n">output</span><span class="o">.</span><span class="n">xaxis</span><span class="p">(</span><span class="s">&quot;</span><span class="se">\\</span><span class="s">rm{Q_{x}}&quot;</span><span class="p">,</span> <span class="s">&#39;A^{-1}&#39;</span><span class="p">)</span>
313                <span class="n">output</span><span class="o">.</span><span class="n">yaxis</span><span class="p">(</span><span class="s">&quot;</span><span class="se">\\</span><span class="s">rm{Q_{y}}&quot;</span><span class="p">,</span> <span class="s">&#39;A^{-1}&#39;</span><span class="p">)</span>
314               
315            <span class="k">if</span> <span class="n">data_conv_i</span> <span class="ow">is</span> <span class="ow">not</span> <span class="bp">None</span><span class="p">:</span>
316                <span class="n">output</span><span class="o">.</span><span class="n">zaxis</span><span class="p">(</span><span class="s">&quot;</span><span class="se">\\</span><span class="s">rm{Intensity}&quot;</span><span class="p">,</span> <span class="n">output</span><span class="o">.</span><span class="n">I_unit</span><span class="p">)</span>
317            <span class="k">else</span><span class="p">:</span>
318                <span class="n">output</span><span class="o">.</span><span class="n">zaxis</span><span class="p">(</span><span class="s">&quot;</span><span class="se">\\</span><span class="s">rm{Intensity}&quot;</span><span class="p">,</span> <span class="s">&quot;cm^{-1}&quot;</span><span class="p">)</span>
319           
320            <span class="k">if</span> <span class="ow">not</span> <span class="n">fversion</span> <span class="o">&gt;=</span> <span class="mf">1.0</span><span class="p">:</span>
321                <span class="n">msg</span> <span class="o">=</span> <span class="s">&quot;Danse_reader can&#39;t read this file </span><span class="si">%s</span><span class="s">&quot;</span> <span class="o">%</span> <span class="n">filename</span>
322                <span class="k">raise</span> <span class="ne">ValueError</span><span class="p">,</span> <span class="n">msg</span>
323            <span class="k">else</span><span class="p">:</span>
324                <span class="n">logging</span><span class="o">.</span><span class="n">info</span><span class="p">(</span><span class="s">&quot;Danse_reader Reading </span><span class="si">%s</span><span class="s"> </span><span class="se">\n</span><span class="s">&quot;</span> <span class="o">%</span> <span class="n">filename</span><span class="p">)</span>
325           
326            <span class="c"># Store loading process information</span>
327            <span class="n">output</span><span class="o">.</span><span class="n">meta_data</span><span class="p">[</span><span class="s">&#39;loader&#39;</span><span class="p">]</span> <span class="o">=</span> <span class="bp">self</span><span class="o">.</span><span class="n">type_name</span>
328            <span class="n">output</span> <span class="o">=</span> <span class="n">reader2D_converter</span><span class="p">(</span><span class="n">output</span><span class="p">)</span>
329            <span class="k">return</span> <span class="n">output</span>
330       
331        <span class="k">return</span> <span class="bp">None</span></div></div>
332</pre></div>
333
334          </div>
335        </div>
336      </div>
337      <div class="sphinxsidebar">
338        <div class="sphinxsidebarwrapper">
339<div id="searchbox" style="display: none">
340  <h3>Quick search</h3>
341    <form class="search" action="../../../../search.html" method="get">
342      <input type="text" name="q" />
343      <input type="submit" value="Go" />
344      <input type="hidden" name="check_keywords" value="yes" />
345      <input type="hidden" name="area" value="default" />
346    </form>
347    <p class="searchtip" style="font-size: 90%">
348    Enter search terms or a module, class or function name.
349    </p>
350</div>
351<script type="text/javascript">$('#searchbox').show(0);</script>
352        </div>
353      </div>
354      <div class="clearer"></div>
355    </div>
356    <div class="related">
357      <h3>Navigation</h3>
358      <ul>
359        <li class="right" style="margin-right: 10px">
360          <a href="../../../../genindex.html" title="General Index"
361             >index</a></li>
362        <li class="right" >
363          <a href="../../../../py-modindex.html" title="Python Module Index"
364             >modules</a> |</li>
365        <li><a href="../../../../index.html">SasView 3.0.0 documentation</a> &raquo;</li>
366          <li><a href="../../../index.html" >Module code</a> &raquo;</li>
367          <li><a href="../readers.html" >sas.dataloader.readers</a> &raquo;</li> 
368      </ul>
369    </div>
370    <div class="footer">
371        &copy; Copyright 2013, The SasView Project.
372      Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.3.
373    </div>
374  </body>
375</html>
Note: See TracBrowser for help on using the repository browser.