File: colourimage.shtml

package info (click to toggle)
munipack 0.6.2-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 33,104 kB
  • sloc: cpp: 29,677; sh: 4,909; f90: 2,872; makefile: 278; python: 140; xml: 72; awk: 12
file content (212 lines) | stat: -rw-r--r-- 5,965 bytes parent folder | download | duplicates (2)
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
<!DOCTYPE HTML>
<html lang="en">
<head>
<!-- #include virtual="/meta.shtml" -->
<title>Munipack ‒ Colour images tutorial</title>
</head>
<body>
<!-- #include virtual="/head.shtml" -->
<section>

<h1>Colour images</h1>

<p>
How to create of images in natural colours with standard photometric filters.
</p>

<h1>Dumbbell nebula in colours</h1>

<p>
  <a href="https://en.wikipedia.org/wiki/Dumbbell_Nebula">Dumbbell nebula</a>
  (M 27) is a planetary nebula in Saggita constellation. Dumbbell is
  one from favourite objects for small telescopes, although its colours are
  too dim for visual observers.
</p>


<h2>Sample Data</h2>

<p>
A sample data are available as
<a href="https://integral.physics.muni.cz/ftp/munipack/munipack-data-m27.tar.gz">munipack-data-m27.tar.gz</a>.
Save it to arbitrary directory (for example <samp>/tmp/</samp>) like:
</p>
<pre>
$ cd /tmp
$ tar zxf munipack-data-m27.tar.gz
$ cd munipack-data-m27/
</pre>
<p>
  The directory <samp>munipack-data-m27/</samp> contains exposures in
  the instrumental colour system (nearly Johnson) in BVR filters and
  corresponding dark frames. B filter exposures has duration 60 seconds.
  Both V and R filter exposures has 40 seconds. Flat-fields frames
  are not available.
</p>

<h2>Data Processing</h2>

<p>
  There is a complete list of Munipack's commands
  for composing of colour frame of Dumbbell nebula.
  This part prepares frames:
</p>

<h3>Corrections</h3>
  <p>
    Prepare correction frames and pre-correct object's frames.
  </p>
<pre>
$ munipack dark -o d60.fits d60_*.fits
$ munipack phcorr -dark d60.fits m27_*B.fits
$ munipack dark -o d40.fits d40_*.fits
$ munipack phcorr -dark d40.fits m27_*[VR].fits
</pre>

<h3>Stars detection and photometry</h3>
<pre>
$ munipack find m27_*.fits
$ munipack aphot m27_*.fits
</pre>
<p>To check detected stars and its preliminary photometry, run xmunipack:</p>
<pre>
$ xmunipack m27_01R.fits
</pre>

<h3>Catalogue</h3>
<p>
  Search the reference astrometry catalogue around Dumbbell's position
  (required for astrometry calibration) in Virtual Observatory:
</p>
<pre>
$ munipack cone -r 0.1 --magmin 10 --magmax 13 -- 299.87 22.71
</pre>

<h3>Astrometry</h3>
<p>Prepare the astrometric calibration of images</p>
<pre>
$ munipack astrometry m27_*.fits
</pre>

<h3>Stack images</h3>
<p>Specify centre of projection of the output image. A good choice is
a point near of centre of object; also the centre of a selected frame
can be recommended. We has chooses <i>α</i>=299.9° and <i>δ</i>=22.72°.
The specification of the common centre of projection (--rcen, --dcen)
for all frames is important for correct alignment of all particular frames.
</p>
<p>
We are merging all images in single filter to get a deeper exposure:
</p>
<pre>
$ munipack kombine -o m27_B.fits --rcen 299.9 --dcen 22.72 m27_*B.fits
$ munipack kombine -o m27_V.fits --rcen 299.9 --dcen 22.72 m27_*V.fits
$ munipack kombine -o m27_R.fits --rcen 299.9 --dcen 22.72 m27_*R.fits
</pre>

<h3>Photometry calibration</h3>
<p>
  To get proper colours, the final frames must be calibrated
  (this adds suffix <samp>_cal</samp> to original filenames):
</p>
<pre>
  $ munipack find m27_?.fits
  $ munipack aphot m27_?.fits
  $ for F in B V R; do
      munipack phcal --verbose  --photsys-ref Johnson --area 0.2 -c cone.fits \
                 -f $F --col-mag ${F}mag --col-magerr e_${F}mag m27_${F}.fits
    done
</pre>


<h2>Colouring</h2>

<p>
  Once data are kombined into a deep exposure frames, everything
  is prepared for colouring.
</p>

<p>
  A direct way to compose a frame in natural colours is to select
  a white star, the star with its colour index near zero:
  our suitable white star is at pixel coordinates (391,54):
</p>
<pre>
  $ munipack colouring --white-star 391,54 m27_B.fits m27_V.fits m27_R.fits
</pre>
<p>
  It leaves frame in <samp>colouring.fits</samp>. Be carefull, result
  does depends on proper choice of the white star. Option
  <samp>--white-radius</samp>, and <samp>--backs</samp>
  can be used for fine tuning.
</p>

<p>
  More simple way is offered by calibrated frames:
</p>
<pre>
  $ munipack colouring m27_B_cal.fits m27_V_cal.fits m27_R_cal.fits
</pre>
<p>
  The photometry calibration of frames does job for us,
  the calibrated frames has balanced the white, by definition.
</p>

<div class="twocolumn">
  <div class="column">
    <figure>
      <img class="figure" src="color-best.png" alt="M27.png" title="Dumbell in Colours">
      <figcaption>Dumbbell nebula presented in natural colours</figcaption>
    </figure>
  </div>
  <div class="column">
    <figure>
      <img class="figure" src="color-scotopic.png" alt="M27.png" title="Dumbell in Night">
      <figcaption>Dumbbell nebula shown in grayscales simulating perception
	by a visual observer adapted to	low-light conditions.</figcaption>
    </figure>
  </div>
</div>

<h2>Astrophysical background</h2>
<p>
  The output image nicely shows regions where radiation by forbidden lines
  dominates (green) and Hα regions heated by shock-waves (red).
</p>

<p>
  Also, try how this image would be visible by human vision at night,
  when colour vision receptors (rods) is not activated.
  The appearance will be similar to visual sights.
  Use xmunipack and Tune->Colour_menu, or try helper utility:
</p>
<pre>
$ fitspng -l 100,7e3 m27.fits
$ fitspng -l 100,7e3 -n 1e9,1e3 m27.fits
</pre>
<p>
  The image, by scotopic vision, is near to the pictures drawn
  by observers from the past, and looks more like a dumbbell.
</p>

<h2>See Also</h2>

<p>
Manuals:
<a href="man_colouring.html">Colouring</a>,
<a href="man_kombine.html">Kombine</a>,
<a href="man_astrometry.html">Astrometry</a>,
<a href="man_aphot.html">Aperture Photometry</a>,
<a href="man_phcorr.html">Photometry pre-corrections</a>.
</p>

<p>
Tutorial:
<a href="kombitut.html">Image Compositions</a> and
the shell script <a href="m27.sh">m27.sh</a>.
</p>

</section>
<!-- #include virtual="/foot.shtml" -->
</body>
</html>