/usr/share/doc/python-mpltoolkits.basemap-doc/html/users/lcc.html is in python-mpltoolkits.basemap-doc 1.0.7+dfsg-3build2.
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 | <!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>Lambert Conformal 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="Lambert Azimuthal Equal Area Projection" href="laea.html" />
<link rel="prev" title="Cylindrial Equal-Area Projection" href="cea.html" />
<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="lambert-conformal-projection">
<span id="lcc"></span><h1>Lambert Conformal Projection<a class="headerlink" href="#lambert-conformal-projection" title="Permalink to this headline">ΒΆ</a></h1>
<p>A conformal projection. The green shapes drawn on the map are equal-area
circles on the surface of the earth. Known as “Tissot’s indicatrix”,
they can be used to show the angular and areal distortion of a map projection.
On a conformal projection, the shape of the circles is preserved, but the
area is not. On a equal-area projection, the area is preserved but the
shape is not.</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 lambert conformal basemap.</span>
<span class="c"># lat_1 is first standard parallel.</span>
<span class="c"># lat_2 is second standard parallel (defaults to lat_1).</span>
<span class="c"># lon_0,lat_0 is central point.</span>
<span class="c"># rsphere=(6378137.00,6356752.3142) specifies WGS4 ellipsoid</span>
<span class="c"># area_thresh=1000 means don't plot coastline features less</span>
<span class="c"># than 1000 km^2 in area.</span>
<span class="n">m</span> <span class="o">=</span> <span class="n">Basemap</span><span class="p">(</span><span class="n">width</span><span class="o">=</span><span class="mi">12000000</span><span class="p">,</span><span class="n">height</span><span class="o">=</span><span class="mi">9000000</span><span class="p">,</span>
<span class="n">rsphere</span><span class="o">=</span><span class="p">(</span><span class="mf">6378137.00</span><span class="p">,</span><span class="mf">6356752.3142</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">area_thresh</span><span class="o">=</span><span class="mf">1000.</span><span class="p">,</span><span class="n">projection</span><span class="o">=</span><span class="s">'lcc'</span><span class="p">,</span>\
<span class="n">lat_1</span><span class="o">=</span><span class="mf">45.</span><span class="p">,</span><span class="n">lat_2</span><span class="o">=</span><span class="mi">55</span><span class="p">,</span><span class="n">lat_0</span><span class="o">=</span><span class="mi">50</span><span class="p">,</span><span class="n">lon_0</span><span class="o">=-</span><span class="mf">107.</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">9</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">12</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">1.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">"Lambert Conformal 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/lcc.py">Source code</a>)</p>
<div class="figure">
<img alt="../_images/lcc.png" src="../_images/lcc.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="cea.html" title="previous chapter">Cylindrial Equal-Area Projection</a></li>
<li>Next: <a href="laea.html" title="next chapter">Lambert Azimuthal Equal Area 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/lcc.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="footer">
©2011, Jeffrey Whitaker.
|
Powered by <a href="http://sphinx-doc.org/">Sphinx 1.3.4</a>
& <a href="https://github.com/bitprophet/alabaster">Alabaster 0.7.6</a>
|
<a href="../_sources/users/lcc.txt"
rel="nofollow">Page source</a>
</div>
</body>
</html>
|