/usr/share/doc/python-mpltoolkits.basemap-doc/html/users/paeqd.html is in python-mpltoolkits.basemap-doc 1.0.7+dfsg-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 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 | <!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/default.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" />
</head>
<body>
<div class="related">
<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="vandg.html" title="van der Grinten Projection"
accesskey="N">next</a> |</li>
<li class="right" >
<a href="plaea.html" title="Polar Lambert Azimuthal Projection"
accesskey="P">previous</a> |</li>
<li><a href="../index.html">Basemap Matplotlib Toolkit 1.0.7 documentation</a> »</li>
<li><a href="index.html" >The Matplotlib Basemap Toolkit User’s Guide</a> »</li>
<li><a href="mapsetup.html" accesskey="U">Setting up the map</a> »</li>
</ul>
</div>
<div class="document">
<div class="documentwrapper">
<div class="bodywrapper">
<div class="body">
<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 <tt class="docutils literal"><span class="pre">npaeqd</span></tt> and <tt class="docutils literal"><span class="pre">spaeqd</span></tt> are provided
for easy access to the polar aspect of the azimuthal equidistant projection.</p>
<div class="highlight-python"><div class="highlight"><pre><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="c"># setup north polar aimuthal equidistant basemap.</span>
<span class="c"># The longitude lon_0 is at 6-o'clock, and the</span>
<span class="c"># latitude circle boundinglat is tangent to the edge</span>
<span class="c"># 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="s">'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="s">'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="s">'coral'</span><span class="p">,</span><span class="n">lake_color</span><span class="o">=</span><span class="s">'aqua'</span><span class="p">)</span>
<span class="c"># 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="s">'aqua'</span><span class="p">)</span>
<span class="c"># 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="s">'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="s">"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="highlight-python"><div class="highlight"><pre><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="c"># setup south polar aimuthal equidistant basemap.</span>
<span class="c"># The longitude lon_0 is at 6-o'clock, and the</span>
<span class="c"># latitude circle boundinglat is tangent to the edge</span>
<span class="c"># 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="s">'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="s">'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="s">'coral'</span><span class="p">,</span><span class="n">lake_color</span><span class="o">=</span><span class="s">'aqua'</span><span class="p">)</span>
<span class="c"># 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="s">'aqua'</span><span class="p">)</span>
<span class="c"># 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="s">'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="s">"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>
</div>
</div>
</div>
<div class="sphinxsidebar">
<div class="sphinxsidebarwrapper">
<h4>Previous topic</h4>
<p class="topless"><a href="plaea.html"
title="previous chapter">Polar Lambert Azimuthal Projection</a></p>
<h4>Next topic</h4>
<p class="topless"><a href="vandg.html"
title="next chapter">van der Grinten Projection</a></p>
<h3>This Page</h3>
<ul class="this-page-menu">
<li><a href="../_sources/users/paeqd.txt"
rel="nofollow">Show Source</a></li>
</ul>
<div id="searchbox" style="display: none">
<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">
<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="vandg.html" title="van der Grinten Projection"
>next</a> |</li>
<li class="right" >
<a href="plaea.html" title="Polar Lambert Azimuthal Projection"
>previous</a> |</li>
<li><a href="../index.html">Basemap Matplotlib Toolkit 1.0.7 documentation</a> »</li>
<li><a href="index.html" >The Matplotlib Basemap Toolkit User’s Guide</a> »</li>
<li><a href="mapsetup.html" >Setting up the map</a> »</li>
</ul>
</div>
<div class="footer">
© Copyright 2011, Jeffrey Whitaker.
Last updated on Jan 11, 2014.
Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.
</div>
</body>
</html>
|