/usr/share/doc/python-pyresample-doc/html/preproc.html is in python-pyresample-doc 1.1.6-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 | <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title>Preprocessing of grids — pyresample 1.1.6 documentation</title>
<link rel="stylesheet" href="_static/classic.css" type="text/css" />
<link rel="stylesheet" href="_static/pygments.css" type="text/css" />
<script type="text/javascript">
var DOCUMENTATION_OPTIONS = {
URL_ROOT: './',
VERSION: '1.1.6',
COLLAPSE_INDEX: false,
FILE_SUFFIX: '.html',
HAS_SOURCE: true
};
</script>
<script type="text/javascript" src="_static/jquery.js"></script>
<script type="text/javascript" src="_static/underscore.js"></script>
<script type="text/javascript" src="_static/doctools.js"></script>
<link rel="top" title="pyresample 1.1.6 documentation" href="index.html" />
<link rel="next" title="Plotting with pyresample and Basemap" href="plot.html" />
<link rel="prev" title="Using multiple processor cores" href="multi.html" />
</head>
<body role="document">
<div class="related" role="navigation" aria-label="related navigation">
<h3>Navigation</h3>
<ul>
<li class="right" style="margin-right: 10px">
<a href="genindex.html" title="General Index"
accesskey="I">index</a></li>
<li class="right" >
<a href="py-modindex.html" title="Python Module Index"
>modules</a> |</li>
<li class="right" >
<a href="plot.html" title="Plotting with pyresample and Basemap"
accesskey="N">next</a> |</li>
<li class="right" >
<a href="multi.html" title="Using multiple processor cores"
accesskey="P">previous</a> |</li>
<li class="nav-item nav-item-0"><a href="index.html">pyresample 1.1.6 documentation</a> »</li>
</ul>
</div>
<div class="document">
<div class="documentwrapper">
<div class="bodywrapper">
<div class="body" role="main">
<div class="section" id="preprocessing-of-grids">
<span id="preproc"></span><h1>Preprocessing of grids<a class="headerlink" href="#preprocessing-of-grids" title="Permalink to this headline">¶</a></h1>
<p>When resampling is performed repeatedly to the same grid significant execution time can be save by
preprocessing grid information.</p>
<div class="section" id="preprocessing-for-grid-resampling">
<h2>Preprocessing for grid resampling<a class="headerlink" href="#preprocessing-for-grid-resampling" title="Permalink to this headline">¶</a></h2>
<p>Using the function <strong>generate_quick_linesample_arrays</strong> or
<strong>generate_nearest_neighbour_linesample_arrays</strong> from <strong>pyresample.utils</strong> arrays containing
the rows and cols indices used to calculate the result in <strong>image.resample_area_quick</strong> or
<strong>resample_area_nearest_neighbour</strong> can be obtained. These can be fed to the method
<strong>get_array_from_linesample</strong> of an <strong>ImageContainer</strong> object to obtain the resample result.</p>
<div class="highlight-python"><div class="highlight"><pre><span class="gp">>>> </span><span class="kn">import</span> <span class="nn">numpy</span>
<span class="gp">>>> </span><span class="kn">from</span> <span class="nn">pyresample</span> <span class="kn">import</span> <span class="n">utils</span><span class="p">,</span> <span class="n">image</span><span class="p">,</span> <span class="n">geometry</span>
<span class="gp">>>> </span><span class="n">area_def</span> <span class="o">=</span> <span class="n">geometry</span><span class="o">.</span><span class="n">AreaDefinition</span><span class="p">(</span><span class="s1">'areaD'</span><span class="p">,</span> <span class="s1">'Europe (3km, HRV, VTC)'</span><span class="p">,</span> <span class="s1">'areaD'</span><span class="p">,</span>
<span class="gp">... </span> <span class="p">{</span><span class="s1">'a'</span><span class="p">:</span> <span class="s1">'6378144.0'</span><span class="p">,</span> <span class="s1">'b'</span><span class="p">:</span> <span class="s1">'6356759.0'</span><span class="p">,</span>
<span class="gp">... </span> <span class="s1">'lat_0'</span><span class="p">:</span> <span class="s1">'50.00'</span><span class="p">,</span> <span class="s1">'lat_ts'</span><span class="p">:</span> <span class="s1">'50.00'</span><span class="p">,</span>
<span class="gp">... </span> <span class="s1">'lon_0'</span><span class="p">:</span> <span class="s1">'8.00'</span><span class="p">,</span> <span class="s1">'proj'</span><span class="p">:</span> <span class="s1">'stere'</span><span class="p">},</span>
<span class="gp">... </span> <span class="mi">800</span><span class="p">,</span> <span class="mi">800</span><span class="p">,</span>
<span class="gp">... </span> <span class="p">[</span><span class="o">-</span><span class="mf">1370912.72</span><span class="p">,</span> <span class="o">-</span><span class="mf">909968.64</span><span class="p">,</span>
<span class="gp">... </span> <span class="mf">1029087.28</span><span class="p">,</span> <span class="mf">1490031.36</span><span class="p">])</span>
<span class="gp">>>> </span><span class="n">msg_area</span> <span class="o">=</span> <span class="n">geometry</span><span class="o">.</span><span class="n">AreaDefinition</span><span class="p">(</span><span class="s1">'msg_full'</span><span class="p">,</span> <span class="s1">'Full globe MSG image 0 degrees'</span><span class="p">,</span>
<span class="gp">... </span> <span class="s1">'msg_full'</span><span class="p">,</span>
<span class="gp">... </span> <span class="p">{</span><span class="s1">'a'</span><span class="p">:</span> <span class="s1">'6378169.0'</span><span class="p">,</span> <span class="s1">'b'</span><span class="p">:</span> <span class="s1">'6356584.0'</span><span class="p">,</span>
<span class="gp">... </span> <span class="s1">'h'</span><span class="p">:</span> <span class="s1">'35785831.0'</span><span class="p">,</span> <span class="s1">'lon_0'</span><span class="p">:</span> <span class="s1">'0'</span><span class="p">,</span>
<span class="gp">... </span> <span class="s1">'proj'</span><span class="p">:</span> <span class="s1">'geos'</span><span class="p">},</span>
<span class="gp">... </span> <span class="mi">3712</span><span class="p">,</span> <span class="mi">3712</span><span class="p">,</span>
<span class="gp">... </span> <span class="p">[</span><span class="o">-</span><span class="mf">5568742.4</span><span class="p">,</span> <span class="o">-</span><span class="mf">5568742.4</span><span class="p">,</span>
<span class="gp">... </span> <span class="mf">5568742.4</span><span class="p">,</span> <span class="mf">5568742.4</span><span class="p">])</span>
<span class="gp">>>> </span><span class="n">data</span> <span class="o">=</span> <span class="n">numpy</span><span class="o">.</span><span class="n">ones</span><span class="p">((</span><span class="mi">3712</span><span class="p">,</span> <span class="mi">3712</span><span class="p">))</span>
<span class="gp">>>> </span><span class="n">msg_con</span> <span class="o">=</span> <span class="n">image</span><span class="o">.</span><span class="n">ImageContainer</span><span class="p">(</span><span class="n">data</span><span class="p">,</span> <span class="n">msg_area</span><span class="p">)</span>
<span class="gp">>>> </span><span class="n">row_indices</span><span class="p">,</span> <span class="n">col_indices</span> <span class="o">=</span> \
<span class="gp">... </span> <span class="n">utils</span><span class="o">.</span><span class="n">generate_nearest_neighbour_linesample_arrays</span><span class="p">(</span><span class="n">msg_area</span><span class="p">,</span> <span class="n">area_def</span><span class="p">,</span> <span class="mi">50000</span><span class="p">)</span>
<span class="gp">>>> </span><span class="n">result</span> <span class="o">=</span> <span class="n">msg_con</span><span class="o">.</span><span class="n">get_array_from_linesample</span><span class="p">(</span><span class="n">row_indices</span><span class="p">,</span> <span class="n">col_indices</span><span class="p">)</span>
</pre></div>
</div>
<p>The numpy arrays returned by <strong>generate_*_linesample_arrays</strong> can be and used with the
<strong>ImageContainer.get_array_from_linesample</strong> method when the same resampling is to be performed
again thus eliminating the need for calculating the reprojection.</p>
<p>Numpy arrays can be saved and loaded using <strong>numpy.save</strong> and <strong>numpy.load</strong>.</p>
</div>
</div>
</div>
</div>
</div>
<div class="sphinxsidebar" role="navigation" aria-label="main navigation">
<div class="sphinxsidebarwrapper">
<h3><a href="index.html">Table Of Contents</a></h3>
<ul>
<li><a class="reference internal" href="#">Preprocessing of grids</a><ul>
<li><a class="reference internal" href="#preprocessing-for-grid-resampling">Preprocessing for grid resampling</a></li>
</ul>
</li>
</ul>
<h4>Previous topic</h4>
<p class="topless"><a href="multi.html"
title="previous chapter">Using multiple processor cores</a></p>
<h4>Next topic</h4>
<p class="topless"><a href="plot.html"
title="next chapter">Plotting with pyresample and Basemap</a></p>
<div role="note" aria-label="source link">
<h3>This Page</h3>
<ul class="this-page-menu">
<li><a href="_sources/preproc.txt"
rel="nofollow">Show Source</a></li>
</ul>
</div>
<div id="searchbox" style="display: none" role="search">
<h3>Quick search</h3>
<form class="search" action="search.html" method="get">
<input type="text" name="q" />
<input type="submit" value="Go" />
<input type="hidden" name="check_keywords" value="yes" />
<input type="hidden" name="area" value="default" />
</form>
<p class="searchtip" style="font-size: 90%">
Enter search terms or a module, class or function name.
</p>
</div>
<script type="text/javascript">$('#searchbox').show(0);</script>
</div>
</div>
<div class="clearer"></div>
</div>
<div class="related" role="navigation" aria-label="related navigation">
<h3>Navigation</h3>
<ul>
<li class="right" style="margin-right: 10px">
<a href="genindex.html" title="General Index"
>index</a></li>
<li class="right" >
<a href="py-modindex.html" title="Python Module Index"
>modules</a> |</li>
<li class="right" >
<a href="plot.html" title="Plotting with pyresample and Basemap"
>next</a> |</li>
<li class="right" >
<a href="multi.html" title="Using multiple processor cores"
>previous</a> |</li>
<li class="nav-item nav-item-0"><a href="index.html">pyresample 1.1.6 documentation</a> »</li>
</ul>
</div>
<div class="footer" role="contentinfo">
© Copyright 2013, Esben S. Nielsen.
Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.3.6.
</div>
</body>
</html>
|