forked from astroML/astroml.github.com
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathplot_crossmatch.html
More file actions
285 lines (228 loc) · 15 KB
/
Copy pathplot_crossmatch.html
File metadata and controls
285 lines (228 loc) · 15 KB
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
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
<!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>Catalog cross-matching — astroML 0.2 documentation</title>
<link rel="stylesheet" href="../../_static/astroMLstyle.css" type="text/css" />
<link rel="stylesheet" href="../../_static/pygments.css" type="text/css" />
<script type="text/javascript">
var DOCUMENTATION_OPTIONS = {
URL_ROOT: '../../',
VERSION: '0.2',
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="shortcut icon" href="../../_static/favicon.ico"/>
<link rel="top" title="astroML 0.2 documentation" href="../../index.html" />
<link rel="up" title="Data Processing Algorithms" href="index.html" />
<link rel="next" title="Linear Sum of Gaussians" href="plot_spectrum_sum_of_norms.html" />
<link rel="prev" title="Bayesian Blocks for Histograms" href="plot_bayesian_blocks.html" />
<!-- Following code is for Google custom search bar -->
<script>
(function() {
var cx = '011400076584591653333:hjd_fbqk1u0';
var gcse = document.createElement('script'); gcse.type = 'text/javascript'; gcse.async = true;
gcse.src = (document.location.protocol == 'https:' ? 'https:' : 'http:') +
'//www.google.com/cse/cse.js?cx=' + cx;
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(gcse, s);
})();
</script>
<!-- Following code is for toggle/glide features -->
<script type="text/javascript" src="http://code.jquery.com/jquery-latest.js"></script>
<script type="text/javascript">
$(document).ready(function(){
$(".toggle_container").hide();
$(".toggle_trigger").click(function(){
$(this).toggleClass("active").next(".toggle_container").slideToggle("fast");
return false; <!-- Prevent the link from being followed -->
});
$(".toggle_trigger#start_open").toggleClass("active").next().show();
});
</script>
<!-- Following code is for Google Analytics -->
<script type="text/javascript">
var _gaq = _gaq || [];
_gaq.push(['_setAccount', 'UA-35748160-1']);
_gaq.push(['_trackPageview']);
(function() {
var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true;
ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s);
})();
</script>
</head>
<body>
<div class="header-wrapper">
<div class="header">
<p class="logo"><a href="../../index.html">
<img src="../../_static/Logo.gif" alt="Logo"/>
</a>
</p><div class="navbar">
<ul>
<li><a href="../../index.html">Home</a></li>
<li><a href="../../user_guide/index.html">User Guide</a></li>
<li><a href="../../book_figures/index.html">Book Figures</a></li>
<li><a href="../index.html">Examples Plots</a></li>
</ul>
<!-- Google custom search. Javascript that enables this is in the header above -->
<div class="search_form">
<gcse:search></gcse:search>
</div>
</div> <!-- end navbar --></div>
</div>
<div class="content-wrapper">
<!-- <div id="blue_tile"></div> -->
<div class="sphinxsidebar">
<div class="rel">
<!-- rellinks[1:] is an ugly hack to avoid link to module
index -->
<div class="rellink">
<a href="plot_bayesian_blocks.html" title="Bayesian Blocks for Histograms"
accesskey="P">Previous
<br/>
<span class="smallrellink">
Bayesian Blocks ...
</span>
<span class="hiddenrellink">
Bayesian Blocks for Histograms
</span>
</a>
</div>
<div class="spacer">
</div>
<div class="rellink">
<a href="plot_spectrum_sum_of_norms.html" title="Linear Sum of Gaussians"
accesskey="N">Next
<br/>
<span class="smallrellink">
Linear Sum of Ga...
</span>
<span class="hiddenrellink">
Linear Sum of Gaussians
</span>
</a>
</div>
<!-- Ad a link to the 'up' page -->
<div class="spacer">
</div>
<div class="rellink">
<a href="index.html" title="Data Processing Algorithms" >
Up
<br/>
<span class="smallrellink">
Data Processing ...
</span>
<span class="hiddenrellink">
Data Processing Algorithms
</span>
</a>
</div>
</div>
<p style="text-align: center">This documentation is
for astroML <strong>version 0.2</strong>
<!-- — <a href="http://scikit-learn.org/stable/support.html#documentation-resources">Other versions</a> -->
</p>
<h3>This page</h3>
<ul>
<li><a class="reference internal" href="#">Catalog cross-matching</a></li>
</ul>
<h3>Links</h3>
<p><a href="https://groups.google.com/forum/#!forum/astroml-general">astroML Mailing List</a></p>
<p><a href="http://github.com/astroML/astroML/issues">GitHub Issue Tracker</a></p>
<h3>Videos</h3>
<p><a href="http://pyvideo.org/video/1218/astroml-data-mining-and-machine-learning-for-ast">Scipy 2012</a> (15 minute talk)</p>
<p><a href="http://pyvideo.org/video/2035/opening-up-astronomy-with-python-and-astroml-sci">Scipy 2013</a> (20 minute talk)</p>
<h3>Citing</h3>
<p>If you use the software, please consider
<!-- hard-coded link here because pathto() adds an extra # -->
<a href="/index.html#citing-astroml">citing astroML</a>.</p>
</div>
<div class="content">
<div class="documentwrapper">
<div class="bodywrapper">
<div class="body">
<div class="section" id="catalog-cross-matching">
<span id="example-algorithms-plot-crossmatch"></span><h1>Catalog cross-matching<a class="headerlink" href="#catalog-cross-matching" title="Permalink to this headline">¶</a></h1>
<p>This plots the cross-matched samples between the SDSS imaging data and
the SDSS Stripe 82 standard stars.</p>
<a class="reference internal image-reference" href="../../_images/plot_crossmatch_1.png"><img alt="../../_images/plot_crossmatch_1.png" class="align-center" src="../../_images/plot_crossmatch_1.png" style="width: 800.0px; height: 600.0px;" /></a>
<div class="toggle_trigger"><a href="#"><p><strong>Code output:</strong></p>
</a></div>
<div class="toggle_container"><div class="highlight-python"><pre>Optimization terminated successfully.
Current function value: -82329.624551
Iterations: 19
Function evaluations: 52
</pre>
</div>
</div>
<div class="toggle_trigger" id="start_open"><a href="#"><p><strong>Python source code:</strong></p>
</a></div>
<div class="toggle_container"><div class="highlight-python"><div class="highlight"><pre><span class="c"># Author: Jake VanderPlas <vanderplas@astro.washington.edu></span>
<span class="c"># License: BSD</span>
<span class="c"># The figure is an example from astroML: see http://astroML.github.com</span>
<span class="kn">import</span> <span class="nn">os</span>
<span class="kn">import</span> <span class="nn">sys</span>
<span class="kn">from</span> <span class="nn">time</span> <span class="kn">import</span> <span class="n">time</span>
<span class="kn">import</span> <span class="nn">numpy</span> <span class="kn">as</span> <span class="nn">np</span>
<span class="kn">from</span> <span class="nn">matplotlib</span> <span class="kn">import</span> <span class="n">pyplot</span> <span class="k">as</span> <span class="n">plt</span>
<span class="kn">from</span> <span class="nn">astroML.datasets</span> <span class="kn">import</span> <span class="n">fetch_imaging_sample</span><span class="p">,</span> <span class="n">fetch_sdss_S82standards</span>
<span class="kn">from</span> <span class="nn">astroML.crossmatch</span> <span class="kn">import</span> <span class="n">crossmatch_angular</span>
<span class="kn">from</span> <span class="nn">astroML.plotting</span> <span class="kn">import</span> <span class="n">hist</span>
<span class="c"># get imaging data</span>
<span class="n">image_data</span> <span class="o">=</span> <span class="n">fetch_imaging_sample</span><span class="p">()</span>
<span class="n">imX</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">((</span><span class="nb">len</span><span class="p">(</span><span class="n">image_data</span><span class="p">),</span> <span class="mi">2</span><span class="p">),</span> <span class="n">dtype</span><span class="o">=</span><span class="n">np</span><span class="o">.</span><span class="n">float64</span><span class="p">)</span>
<span class="n">imX</span><span class="p">[:,</span> <span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="n">image_data</span><span class="p">[</span><span class="s">'ra'</span><span class="p">]</span>
<span class="n">imX</span><span class="p">[:,</span> <span class="mi">1</span><span class="p">]</span> <span class="o">=</span> <span class="n">image_data</span><span class="p">[</span><span class="s">'dec'</span><span class="p">]</span>
<span class="c"># get standard stars</span>
<span class="n">standards_data</span> <span class="o">=</span> <span class="n">fetch_sdss_S82standards</span><span class="p">()</span>
<span class="n">stX</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">empty</span><span class="p">((</span><span class="nb">len</span><span class="p">(</span><span class="n">standards_data</span><span class="p">),</span> <span class="mi">2</span><span class="p">),</span> <span class="n">dtype</span><span class="o">=</span><span class="n">np</span><span class="o">.</span><span class="n">float64</span><span class="p">)</span>
<span class="n">stX</span><span class="p">[:,</span> <span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="n">standards_data</span><span class="p">[</span><span class="s">'RA'</span><span class="p">]</span>
<span class="n">stX</span><span class="p">[:,</span> <span class="mi">1</span><span class="p">]</span> <span class="o">=</span> <span class="n">standards_data</span><span class="p">[</span><span class="s">'DEC'</span><span class="p">]</span>
<span class="c"># crossmatch catalogs</span>
<span class="n">max_radius</span> <span class="o">=</span> <span class="mf">1.</span> <span class="o">/</span> <span class="mi">3600</span> <span class="c"># 1 arcsec</span>
<span class="n">dist</span><span class="p">,</span> <span class="n">ind</span> <span class="o">=</span> <span class="n">crossmatch_angular</span><span class="p">(</span><span class="n">imX</span><span class="p">,</span> <span class="n">stX</span><span class="p">,</span> <span class="n">max_radius</span><span class="p">)</span>
<span class="n">match</span> <span class="o">=</span> <span class="o">~</span><span class="n">np</span><span class="o">.</span><span class="n">isinf</span><span class="p">(</span><span class="n">dist</span><span class="p">)</span>
<span class="n">dist_match</span> <span class="o">=</span> <span class="n">dist</span><span class="p">[</span><span class="n">match</span><span class="p">]</span>
<span class="n">dist_match</span> <span class="o">*=</span> <span class="mi">3600</span>
<span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">axes</span><span class="p">()</span>
<span class="n">hist</span><span class="p">(</span><span class="n">dist_match</span><span class="p">,</span> <span class="n">bins</span><span class="o">=</span><span class="s">'knuth'</span><span class="p">,</span> <span class="n">ax</span><span class="o">=</span><span class="n">ax</span><span class="p">,</span>
<span class="n">histtype</span><span class="o">=</span><span class="s">'stepfilled'</span><span class="p">,</span> <span class="n">ec</span><span class="o">=</span><span class="s">'k'</span><span class="p">,</span> <span class="n">fc</span><span class="o">=</span><span class="s">'#AAAAAA'</span><span class="p">)</span>
<span class="n">ax</span><span class="o">.</span><span class="n">set_xlabel</span><span class="p">(</span><span class="s">'radius of match (arcsec)'</span><span class="p">)</span>
<span class="n">ax</span><span class="o">.</span><span class="n">set_ylabel</span><span class="p">(</span><span class="s">'N(r, r+dr)'</span><span class="p">)</span>
<span class="n">ax</span><span class="o">.</span><span class="n">text</span><span class="p">(</span><span class="mf">0.95</span><span class="p">,</span> <span class="mf">0.95</span><span class="p">,</span>
<span class="s">"Total objects: </span><span class="si">%i</span><span class="se">\n</span><span class="s">Number with match: </span><span class="si">%i</span><span class="s">"</span> <span class="o">%</span> <span class="p">(</span><span class="n">imX</span><span class="o">.</span><span class="n">shape</span><span class="p">[</span><span class="mi">0</span><span class="p">],</span>
<span class="n">np</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="n">match</span><span class="p">)),</span>
<span class="n">ha</span><span class="o">=</span><span class="s">'right'</span><span class="p">,</span> <span class="n">va</span><span class="o">=</span><span class="s">'top'</span><span class="p">,</span> <span class="n">transform</span><span class="o">=</span><span class="n">ax</span><span class="o">.</span><span class="n">transAxes</span><span class="p">)</span>
<span class="n">ax</span><span class="o">.</span><span class="n">set_xlim</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span> <span class="mf">0.2</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>
</div>
<div align="right"><p><a class="reference download internal" href="../../_downloads/plot_crossmatch.py"><tt class="xref download docutils literal"><span class="pre">[download</span> <span class="pre">source:</span> <span class="pre">plot_crossmatch.py]</span></tt></a></p>
</div></div>
</div>
</div>
</div>
<div class="clearer"></div>
</div>
</div>
<div class="footer">
<p style="text-align: center">This documentation is relative
to astroML version 0.2<br/>
© 2012, Jake Vanderplas.
Created using <a href="http://sphinx.pocoo.org/">Sphinx</a> 1.1.3. Design by <a href="http://webylimonada.com">Web y Limonada</a>.
<span style="padding-left: 5ex;">
<a href="../../_sources/examples/algorithms/plot_crossmatch.txt"
rel="nofollow">Show this page source</a>
</span></p>
</div>
</body>
</html>