/usr/share/doc/python-mpltoolkits.basemap-doc/html/users/lcc.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 | <!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/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="Lambert Azimuthal Equal Area Projection" href="laea.html" />
<link rel="prev" title="Cylindrial Equal-Area Projection" href="cea.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="laea.html" title="Lambert Azimuthal Equal Area Projection"
accesskey="N">next</a> |</li>
<li class="right" >
<a href="cea.html" title="Cylindrial Equal-Area 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="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>
</div>
</div>
</div>
<div class="sphinxsidebar">
<div class="sphinxsidebarwrapper">
<h4>Previous topic</h4>
<p class="topless"><a href="cea.html"
title="previous chapter">Cylindrial Equal-Area Projection</a></p>
<h4>Next topic</h4>
<p class="topless"><a href="laea.html"
title="next chapter">Lambert Azimuthal Equal Area Projection</a></p>
<h3>This Page</h3>
<ul class="this-page-menu">
<li><a href="../_sources/users/lcc.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="laea.html" title="Lambert Azimuthal Equal Area Projection"
>next</a> |</li>
<li class="right" >
<a href="cea.html" title="Cylindrial Equal-Area 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>
|