/usr/share/doc/python-mpltoolkits.basemap-doc/html/users/paeqd.html is in python-mpltoolkits.basemap-doc 1.0.7+dfsg-4.
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 165 166 167 | <!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>Polar Azimuthal Equidistant Projection — Basemap Matplotlib Toolkit 1.0.7 documentation</title>
<link rel="stylesheet" href="../_static/alabaster.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.0.7',
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="search" type="application/opensearchdescription+xml"
title="Search within Basemap Matplotlib Toolkit 1.0.7 documentation"
href="../_static/opensearch.xml"/>
<link rel="top" title="Basemap Matplotlib Toolkit 1.0.7 documentation" href="../index.html" />
<link rel="up" title="Setting up the map" href="mapsetup.html" />
<link rel="next" title="van der Grinten Projection" href="vandg.html" />
<link rel="prev" title="Polar Lambert Azimuthal Projection" href="plaea.html" />
<link rel="stylesheet" href="../_static/custom.css" type="text/css" />
<meta name="viewport" content="width=device-width, initial-scale=0.9, maximum-scale=0.9" />
</head>
<body role="document">
<div class="document">
<div class="documentwrapper">
<div class="bodywrapper">
<div class="body" role="main">
<div class="section" id="polar-azimuthal-equidistant-projection">
<span id="paeqd"></span><h1>Polar Azimuthal Equidistant Projection<a class="headerlink" href="#polar-azimuthal-equidistant-projection" title="Permalink to this headline">ΒΆ</a></h1>
<p>For convenience, the projections <code class="docutils literal"><span class="pre">npaeqd</span></code> and <code class="docutils literal"><span class="pre">spaeqd</span></code> are provided
for easy access to the polar aspect of the azimuthal equidistant projection.</p>
<div class="highlight-python"><div class="highlight"><pre><span></span><span class="kn">from</span> <span class="nn">mpl_toolkits.basemap</span> <span class="kn">import</span> <span class="n">Basemap</span>
<span class="kn">import</span> <span class="nn">numpy</span> <span class="kn">as</span> <span class="nn">np</span>
<span class="kn">import</span> <span class="nn">matplotlib.pyplot</span> <span class="kn">as</span> <span class="nn">plt</span>
<span class="kn">from</span> <span class="nn">matplotlib.patches</span> <span class="kn">import</span> <span class="n">Polygon</span>
<span class="c1"># setup north polar aimuthal equidistant basemap.</span>
<span class="c1"># The longitude lon_0 is at 6-o'clock, and the</span>
<span class="c1"># latitude circle boundinglat is tangent to the edge</span>
<span class="c1"># of the map at lon_0.</span>
<span class="n">m</span> <span class="o">=</span> <span class="n">Basemap</span><span class="p">(</span><span class="n">projection</span><span class="o">=</span><span class="s1">'npaeqd'</span><span class="p">,</span><span class="n">boundinglat</span><span class="o">=</span><span class="mi">10</span><span class="p">,</span><span class="n">lon_0</span><span class="o">=</span><span class="mi">270</span><span class="p">,</span><span class="n">resolution</span><span class="o">=</span><span class="s1">'l'</span><span class="p">)</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawcoastlines</span><span class="p">()</span>
<span class="n">m</span><span class="o">.</span><span class="n">fillcontinents</span><span class="p">(</span><span class="n">color</span><span class="o">=</span><span class="s1">'coral'</span><span class="p">,</span><span class="n">lake_color</span><span class="o">=</span><span class="s1">'aqua'</span><span class="p">)</span>
<span class="c1"># draw parallels and meridians.</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawparallels</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="o">-</span><span class="mf">80.</span><span class="p">,</span><span class="mf">81.</span><span class="p">,</span><span class="mf">20.</span><span class="p">))</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawmeridians</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="o">-</span><span class="mf">180.</span><span class="p">,</span><span class="mf">181.</span><span class="p">,</span><span class="mf">20.</span><span class="p">))</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawmapboundary</span><span class="p">(</span><span class="n">fill_color</span><span class="o">=</span><span class="s1">'aqua'</span><span class="p">)</span>
<span class="c1"># draw tissot's indicatrix to show distortion.</span>
<span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">gca</span><span class="p">()</span>
<span class="k">for</span> <span class="n">y</span> <span class="ow">in</span> <span class="n">np</span><span class="o">.</span><span class="n">linspace</span><span class="p">(</span><span class="n">m</span><span class="o">.</span><span class="n">ymax</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="mi">19</span><span class="o">*</span><span class="n">m</span><span class="o">.</span><span class="n">ymax</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="mi">10</span><span class="p">):</span>
<span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">np</span><span class="o">.</span><span class="n">linspace</span><span class="p">(</span><span class="n">m</span><span class="o">.</span><span class="n">xmax</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="mi">19</span><span class="o">*</span><span class="n">m</span><span class="o">.</span><span class="n">xmax</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="mi">10</span><span class="p">):</span>
<span class="n">lon</span><span class="p">,</span> <span class="n">lat</span> <span class="o">=</span> <span class="n">m</span><span class="p">(</span><span class="n">x</span><span class="p">,</span><span class="n">y</span><span class="p">,</span><span class="n">inverse</span><span class="o">=</span><span class="bp">True</span><span class="p">)</span>
<span class="n">poly</span> <span class="o">=</span> <span class="n">m</span><span class="o">.</span><span class="n">tissot</span><span class="p">(</span><span class="n">lon</span><span class="p">,</span><span class="n">lat</span><span class="p">,</span><span class="mf">2.5</span><span class="p">,</span><span class="mi">100</span><span class="p">,</span>\
<span class="n">facecolor</span><span class="o">=</span><span class="s1">'green'</span><span class="p">,</span><span class="n">zorder</span><span class="o">=</span><span class="mi">10</span><span class="p">,</span><span class="n">alpha</span><span class="o">=</span><span class="mf">0.5</span><span class="p">)</span>
<span class="n">plt</span><span class="o">.</span><span class="n">title</span><span class="p">(</span><span class="s2">"North Polar Azimuthal Equidistant Projection"</span><span class="p">)</span>
<span class="n">plt</span><span class="o">.</span><span class="n">show</span><span class="p">()</span>
</pre></div>
</div>
<p>(<a class="reference external" href="../users/figures/npaeqd.py">Source code</a>)</p>
<div class="figure">
<img alt="../_images/npaeqd.png" src="../_images/npaeqd.png" />
</div>
<div class="highlight-python"><div class="highlight"><pre><span></span><span class="kn">from</span> <span class="nn">mpl_toolkits.basemap</span> <span class="kn">import</span> <span class="n">Basemap</span>
<span class="kn">import</span> <span class="nn">numpy</span> <span class="kn">as</span> <span class="nn">np</span>
<span class="kn">import</span> <span class="nn">matplotlib.pyplot</span> <span class="kn">as</span> <span class="nn">plt</span>
<span class="c1"># setup south polar aimuthal equidistant basemap.</span>
<span class="c1"># The longitude lon_0 is at 6-o'clock, and the</span>
<span class="c1"># latitude circle boundinglat is tangent to the edge</span>
<span class="c1"># of the map at lon_0.</span>
<span class="n">m</span> <span class="o">=</span> <span class="n">Basemap</span><span class="p">(</span><span class="n">projection</span><span class="o">=</span><span class="s1">'spaeqd'</span><span class="p">,</span><span class="n">boundinglat</span><span class="o">=-</span><span class="mi">10</span><span class="p">,</span><span class="n">lon_0</span><span class="o">=</span><span class="mi">90</span><span class="p">,</span><span class="n">resolution</span><span class="o">=</span><span class="s1">'l'</span><span class="p">)</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawcoastlines</span><span class="p">()</span>
<span class="n">m</span><span class="o">.</span><span class="n">fillcontinents</span><span class="p">(</span><span class="n">color</span><span class="o">=</span><span class="s1">'coral'</span><span class="p">,</span><span class="n">lake_color</span><span class="o">=</span><span class="s1">'aqua'</span><span class="p">)</span>
<span class="c1"># draw parallels and meridians.</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawparallels</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="o">-</span><span class="mf">80.</span><span class="p">,</span><span class="mf">81.</span><span class="p">,</span><span class="mf">20.</span><span class="p">))</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawmeridians</span><span class="p">(</span><span class="n">np</span><span class="o">.</span><span class="n">arange</span><span class="p">(</span><span class="o">-</span><span class="mf">180.</span><span class="p">,</span><span class="mf">181.</span><span class="p">,</span><span class="mf">20.</span><span class="p">))</span>
<span class="n">m</span><span class="o">.</span><span class="n">drawmapboundary</span><span class="p">(</span><span class="n">fill_color</span><span class="o">=</span><span class="s1">'aqua'</span><span class="p">)</span>
<span class="c1"># draw tissot's indicatrix to show distortion.</span>
<span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">gca</span><span class="p">()</span>
<span class="k">for</span> <span class="n">y</span> <span class="ow">in</span> <span class="n">np</span><span class="o">.</span><span class="n">linspace</span><span class="p">(</span><span class="mi">19</span><span class="o">*</span><span class="n">m</span><span class="o">.</span><span class="n">ymin</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="n">m</span><span class="o">.</span><span class="n">ymin</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="mi">10</span><span class="p">):</span>
<span class="k">for</span> <span class="n">x</span> <span class="ow">in</span> <span class="n">np</span><span class="o">.</span><span class="n">linspace</span><span class="p">(</span><span class="mi">19</span><span class="o">*</span><span class="n">m</span><span class="o">.</span><span class="n">xmin</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="n">m</span><span class="o">.</span><span class="n">xmin</span><span class="o">/</span><span class="mi">20</span><span class="p">,</span><span class="mi">10</span><span class="p">):</span>
<span class="n">lon</span><span class="p">,</span> <span class="n">lat</span> <span class="o">=</span> <span class="n">m</span><span class="p">(</span><span class="n">x</span><span class="p">,</span><span class="n">y</span><span class="p">,</span><span class="n">inverse</span><span class="o">=</span><span class="bp">True</span><span class="p">)</span>
<span class="n">poly</span> <span class="o">=</span> <span class="n">m</span><span class="o">.</span><span class="n">tissot</span><span class="p">(</span><span class="n">lon</span><span class="p">,</span><span class="n">lat</span><span class="p">,</span><span class="mf">2.5</span><span class="p">,</span><span class="mi">100</span><span class="p">,</span>\
<span class="n">facecolor</span><span class="o">=</span><span class="s1">'green'</span><span class="p">,</span><span class="n">zorder</span><span class="o">=</span><span class="mi">10</span><span class="p">,</span><span class="n">alpha</span><span class="o">=</span><span class="mf">0.5</span><span class="p">)</span>
<span class="n">plt</span><span class="o">.</span><span class="n">title</span><span class="p">(</span><span class="s2">"South Polar Azimuthal Equidistant Projection"</span><span class="p">)</span>
<span class="n">plt</span><span class="o">.</span><span class="n">show</span><span class="p">()</span>
</pre></div>
</div>
<p>(<a class="reference external" href="../users/figures/spaeqd.py">Source code</a>)</p>
<div class="figure">
<img alt="../_images/spaeqd.png" src="../_images/spaeqd.png" />
</div>
</div>
</div>
</div>
</div>
<div class="sphinxsidebar" role="navigation" aria-label="main navigation">
<div class="sphinxsidebarwrapper"><div class="relations">
<h3>Related Topics</h3>
<ul>
<li><a href="../index.html">Documentation overview</a><ul>
<li><a href="index.html">The Matplotlib Basemap Toolkit User’s Guide</a><ul>
<li><a href="mapsetup.html">Setting up the map</a><ul>
<li>Previous: <a href="plaea.html" title="previous chapter">Polar Lambert Azimuthal Projection</a></li>
<li>Next: <a href="vandg.html" title="next chapter">van der Grinten Projection</a></li>
</ul></li>
</ul></li>
</ul></li>
</ul>
</div>
<div role="note" aria-label="source link">
<h3>This Page</h3>
<ul class="this-page-menu">
<li><a href="../_sources/users/paeqd.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">
<div><input type="text" name="q" /></div>
<div><input type="submit" value="Go" /></div>
<input type="hidden" name="check_keywords" value="yes" />
<input type="hidden" name="area" value="default" />
</form>
</div>
<script type="text/javascript">$('#searchbox').show(0);</script>
</div>
</div>
<div class="clearer"></div>
</div>
<div class="footer">
©2016, Jeffrey Whitaker.
|
Powered by <a href="http://sphinx-doc.org/">Sphinx 1.4.5</a>
& <a href="https://github.com/bitprophet/alabaster">Alabaster 0.7.8</a>
|
<a href="../_sources/users/paeqd.txt"
rel="nofollow">Page source</a>
</div>
</body>
</html>
|