Contents 

Mira Pro x64
Contents
New Features
Getting Application Help
About this Document
Tech Support
Getting Started
Glossary of Terms
Selecting Source Images
Setting Application Properties
Command Dialogs
Screen Capture
Mira's Special Folders
Changing the GUI Theme
Critical Concepts
Tip of the Day
Tutorials
Overview
Displaying an Image
Displaying an Image Set
Creating Plots from Images
Introduction to Image Registration
Making a Scatter Plot of Measurements
Using FITS Keywords to Analyze Image Data
Introduction to Aperture Photometry
Doing Time Series Photometry
Introduction to File Event Scripting
Cleaning Artifacts from an Image Set
Introduction to Source Extraction
Introduction to Astrometric Calibration
Examples
File Lists
Creating File Lists
Opening Files from File Lists
Images
The Boost Transfer Function Option
Image Rotation
Adding Labels to an Image
Editing the FITS Header
Image Registration
Plots
Plot Windows
Plotting Examples
Plotting an Averaged Line
Row Plots
Adding Series to a Line Profile Plot
Changing Plot Series Attributes
Plotting in World Coordinates
Comparison of Histogram Preferences
Contour Plotting Smoothing
3-D Plots
3-D Pixel Representations
3-D Z-Scaling
Measurements
Distance and Angle
Selecting Cells in a Grid Control
Tutorial: Making a Scatter Plot of Measurements
Fixing Header Problems in Photometric Data
Importing a Photometry Catalog
Image Basics
Image Types
Image Sets
Overview
Image Set Flags
Image Set Properties
Image Set Toolbar
Image Set Applications
Image Sets Menu
Duplicating an Image Set
File Lists
Overview
Creating file lists
Editing file lists
Opening files from file lists
Opening Images
Overview
File opening and saving commands
Using Drag & Drop
Opening image files
Creating an image catalog
Duplicating an Image
Duplicating an Image Set
Opening files from file Lists
Using File Opening Plug-ins
Selecting Source Images
Importing a Binary Image
Importing a Text Image
Saving Images
Saving a single image
Saving an image set
Saving all image set members
Saving an image set to a Folder
Saving an image set with a filename Suffix
Saving an image set to a file list
Exporting as a binary image
Exporting as a text image
Palettes
Overview
Image Palettes
Palette Pane
Palette Properties
Using the Palette Property Trackbars
Rubber Banding the Palette Graph
Importing Palettes
Transfer Functions
Overview
Choosing Transfer Function Properties
Using the Boost Transfer Function Option
Default Transfer Function Properties
Transfer Function Pane
Transfer Function Popup Menus
Image Cursor
Overview
Image Cursor Properties
Setting default and local properties
Image Window Properties
Image Plot Properties
Histogram Properties
Centroid Properties
FWHM Preferences
Radial Profile Properties
Palette Properties
Transfer Function Properties
General Properties
Selecting source images for commands
Overview
Image Display
Image Windows
Image Bar
Animation Bar
User Interface
User Interface Overview
Working with Toolbars
Initializing Toolbars and Docking Panes
Docking Panes
Command Dialogs
Command Toolbars
Main Message Pane
Image Cursor Toolbar
Image Keys Pane
Status Bar
Accelerator Keys
Profile Control
Status Bar
Save Window Position
GUI Theme
Changing the GUI theme
Context Menus
Image Context Menu
Image Context Menu for an image set
Plot Context Menu
3-D Plot Context Menu
Grid Context Menu
Dialog Styles
Command Dialogs
Dialog Buttons
Docking Panes
Overview of Docking Panes
Docking Panes Toolbar
Image Bar
Animation Bar
Image Keys Pane
Measurement Panes
Initializing Toolbars and Docking Panes
Property Profiles for Commands
Overview
Profile Control
Grid Controls
Overview
Report Windows
Grid Context Menu
Measurement Panes
Table Properties
Overview
Table Appearance Properties
Table Printing Properties
Menus
Tear-off menus
Main Context Menu
Default Menu Bar
Align Images Menu
Pictures of Menus
Creating new text documents
Default File Menu
Calibration Menu
Diagnostics menu
Edit Menu
File Menu
Measure Menu
Mira Windows
Image Windows
Plot Windows
Grid Controls
Image
Overview
Animation Pane
Image Bar
Image Set Toolbar
Image Window Keystroke Commands
Image Keys Pane
Adjusting the Image Cursor
Command Modes
Cursor Command Modes
Roam Mode
Cursor Mode
Plot
Overview
Plotting Images and Data
Creating a Plot from Table Data
Plotting Examples
Plot Bar
Plot Animation Bar
Plot Properties
Plot Series Properties
Plot Marker Properties
Image Plot Properties
Save Window Position
Chart
Chart Windows
Plotting Images and Data
Creating a Chart from Table Data
Chart Command Menu
Chart Axis Limits
Chart Properties Pane
Chart Series Properties
Chart Titles
Chart Legend Properties
Save Window Position
3-D Plot
Overview
3-D Animation Toolbar
3-D Rotation Toolbar
3-D Drawing Attributes
3-D Light Source Dialog
3-D Viewpoint Dialog
3-D Window Dialog
Save Window Position
File Open
Open Command
Opening & Saving Files
Opening Files from File Lists
Report
Overview
Grid Controls
Report Commands
Report Window Menu
Save Window Position
Grid Properties
Overview
Table Appearance Properties
Table Printing Properties
Text Editor
Overview
Creating a new text document
Text Editor Edit Menu
Text Editor View Menu
Save Window Position
Script Editor Window
Script Editor Toolbar
Script Editor Context Menu
Script Editor View Menu
Save Window Position
Plug-ins
File Opening Plug-ins
Image Processing Plug-ins
Toolbars
Working with Toolbars
Command Toolbars
Image Keys Pane
3-D Animation Toolbar
3-D Rotation Toolbar
Cursor Tools Toolbar
Image Cursor Toolbar
Image Plot Toolbar
Image Set Toolbar
Image Tools Toolbar
Image Bar
Calibration Tools Toolbar
Interactive Repair Toolbar
Main Toolbar
Measurements Toolbar
Working with FITS Format
Overview
FITS Format Definition
FITS keywords
FITS Header Editor
FITS Header Pane
FITS image orientation
List COMMENT Keywords
List HISTORY Keywords
List Statistics Keywords
List WCS keywords
MIRAPROC keyword
Editing the FITS header
Modifying keywords
Adding metadata at file open or file save
Creating an image list from keywords
Searching images for keyword values
Viewing the FITS File Header
Fixing Header Problems in Photometric Data
Image Information Editor (RGB images)
Coordinate Systems
Overview
Image Coordinate System
Pixel Coordinate Definition
Subpixel Coordinate Definition
World Coordinate System
Image Display Orientation
FITS Image Orientation
Image Coordinate Display
Coordinate Commands
Go to cursor coordinates
Go To Object
Astrometric Calibration
Set an arcsecond scale
Set an equatorial scale
Set a linear scale
Set independent linear scales
Set linear world coordinate units
Setting brightness units
Alignment Commands
High precision image registration
Align on Object
Align on Point
Align by WCS
Coordinate Readout
Set X-Axis Linear Format
Set Y-Axis Linear Format
Set Pixel Value Format
Select Image Properties
World Coordinate System
Overview
World Coordinate System Keywords
Deleting the WCS
Astrometric Calibration
Overview
Tutorial
Entering coordinate data
Changing reference point data
Calibration Preferences dialog
Plate Solution Strategy
Other preferences
Measuring Images
Overview
Measurement Panes
Image Measurements Toolbar
Image Cursor Toolbar
Image Keys Pane
Moving the image cursor to a coordinate
Setting Application Properties
Testing using synthetic images
Setting statistical estimator properties
Opening a table in Microsoft Excel
Tutorial: Working with Measurements in a Grid
Angle Measurements
Command overview
Marker properties
Aperture Photometry
Overview of aperture photometry
Apphot Pane for photometry results
Tutorial: Introduction to Aperture Photometry
Tutorial: Doing Time Series Photometry
How Magnitude and Errors are calculated
Editing target data and standard star data
Automatic detection of objects
Plotting a light curve
Kwee - van Woerden
Setting marker attributes
Preparing a report for the AAVSO
Setting the aperture size and shape
Setting photometry keywords
Fixing header problems in photometric data
Importing a photometry catalog
Calculate Julian Date
Definitions of measurements
Definitions of errors
Understanding signal to noise ratio
Aperture Photometry Toolbar
Aperture Photometry Properties
Aperture Photometry Properties dialog
Point Markers
Calibration Results
Image Information
Photometry Keywords
AAVSO Reporting
Other Properties
Measurement Output
Astrometric Calibration
Overview
Tutorial
Plate Solution Strategy
Centroid Measurements
Command overview
Centroid properties
Distance & Angle Measurements
Command overview
Measuring Distance and Angle
Example
Angle Measurement Definition
FWHM Measurements
Command overview
FWHM properties
PSF Fit Properties
Labeling Images
Overview
Label Properties
Examples
Line Measurements
Command overview
Marker properties
Plot Line Measurements
Command overview
Marker Properties
Plot Point Measurements
Plot Point Measurements
Marker Properties
Shape Properties
Point Measurements
Command overview
Marker properties
Region Measurements
Command overview
Marker properties
Region properties
Statistics Measurements
Command overview
Statistics Measurements
Statistics properties
List Statistics Keywords
Remove Statistics Keywords
Marker Properties
Overview
Drawing Properties
Centroid Properties
Region Properties
Copying & Pasting Markers
Copy Markers
Paste Markers
Scatter Plots
Scatter Plot command
Tutorial: Making a Scatter Plot of Measurements
Selecting Cells in a Grid Control
Point Sampling
Overview
Pixel Table Pane
Viewing and editing pixel values
Report Windows
Overview
Grid Controls
Grid Context Menu
Optimizing the column width
Opening a table in Excel
Table Properties
Table Printing Properties
Table Appearance Properties
Plotting Images and Data
Overview
Plot Bar Pane
Creating a Chart from Table Data
Creating a Plot from Table Data
Examples of plot types
Examples of row plots
Setting Application Properties
Horizontal Slice Plot
Vertical Slice Plot
Plot Windows
Overview
Creating a Plot from Table Data
Plot Bar
Plot Animation Bar
Plot Keys Pane
Plot Context Menu
Set X-Axis Format
Set Y-Axis Format
Plot Series Data
Plot Measurements
Go to Source Image
Plot Averaging Modes
Saving a plot as a text file
Save Window Position
Window Properties
Plot Properties
Plot Series Properties
Chart Windows
Chart Windows
Creating a Chart from Table Data
Chart Command Menu
Chart Axis Limits
Chart Properties Pane
Chart Series Properties
Chart Titles
Chart Legend Properties
Save Window Position
Working with Plot Series
Changing Plot Series Attributes
Overplotting & Animating
Go to Source Image
Copying & pasting plot series
Adding series to a line profile plot
Plot Default Marker Properties
Setting marker attributes before the plot
Default Series Properties
Plot Series Properties
Plot Series Data
Column and Row Profile Plots
Column profile plots
Making a column profile plot
Row profile plots
Making a row profile plot
Plotting an averaged line
Overview
Plot Averaging Modes
Changing the Plot Appearance
Changing Plot Attributes
The Plot Attributes dialog
Using a dialog shortcut
X Axis Properties
Y Axis Properties
Plot Labels
Frame Properties
Selecting a Coordinate System
Plot coordinate systems
Using World Coordinates
Comparison of coordinate systems
Line profile plots
Command
Creating a Line Profile plot
Adding series along different lines
Adding parallel series
Preferences
Setting Line Profile preferences
Centroid Preferences
Marker Properties
Radial profile plots
Making a radial profile plot
Radial Profile Properties
PSF Fit Properties
Histogram plots
Making a Histogram plot
Histogram Properties
Comparison of Histogram Preferences
Making a Bit Histogram plot
Plot Measurements
Overview
Plot Point Measurements
Plot Line Measurements
Statistics Measurements
Statistics Properties
Mark Statistic
Plot Marker Properties
Plot Notes Pane
Scatter plots
Scatter Plot command
Tutorial: Making a Scatter Plot of Measurements
Selecting Cells in a Grid Control
Pixel Series plots
Plotting a value through an image set
3-D plots
3-D Surface Plot command
3-D Plot Type Examples
3-D Z-Scaling Examples
Preferences
3-D Viewpoint
3-D Window Dialog
3-D Drawing Attributes
3-D Light Source
Surface Plot Z Limits
Contour plots
Contour Plot command
Contour Plot Properties
Contour Drawing Properties
Interactive Contour Plot command
Interactive Contour Properties
Contour Plot Smoothing
Copy & Paste
Example
Copying plot data to the clipboard
Pasting plots from the clipboard
Copying the window bitmap
Printing Plots
Overview
Properties
Overview
Printing Placement Properties
Printing Overlay Properties
Printing Heading Properties
Printing Text Properties
Calibrating Images
Selecting Source Images
Overview
Creating synthetic images
Trimming an image section
Adding metadata at File Open or Save
Modifying header keywords
Creating a Pixel Mask
Applying a Pixel Mask
Image Combining Strategy
Editing and Creating a Blemish Mask
Applying a Blemish Mask
File Event Scripting
Astrometric Calibration
Express Image Calibration
Express Image Calibration
Bias Correction Properties
Dark Correction Properties
Flat Correction Properties
Cosmetic Correction Properties
Creating a Master Bias
Create Master Bias
Image Combining Strategy
Creating a Master Dark
Create Master Dark
Image Combining Strategy
Create Master Flat
Create Master Flat
Flat Frame Selection Criteria
Image Combining Strategy
Create Illumination Flat
Create Illumination Flat
Flat Correction Properties
Image Combining Strategy
Create Pixel Flat
Create Pixel Flat
Flat Correction Properties
Image Combining Strategy
Repairing Artefacts and Cosmetic Defects
Overview of Methods
Cleaning artefacts from an image set
Interactively Repairing Defects
Evaluating repair techniques using synthetic images
Blemish Masks
Creating and editing a blemish mask
Applying a Blemish Mask
Pixel Masks
Creating a Pixel Mask
Editing a Pixel Mask
Applying a Pixel Mask
Removing Stars from Images
Overview
Changing the aperture properties
Viewing data for template samples
Properties
Overview
Point Properties
Other Properties
Processing Images
Overview
Command Dialogs
Selecting Source Images
Image Types
Image Region Commands
Testing calibration strategies using synthetic images
Image Keys Pane
File Event Scripting
Adding metadata at File Open or Save
Creating an image list from keywords
Searching images for keyword values
Cosmetic Repairs
Interactive Repair
Interactive Repair Properties
Removing cosmic rays from an image set
Blemish Masks
Creating and editing a blemish mask
Applying a blemish mask
Pixel Masks
Creating a pixel mask from an image
Editing a pixel mask
Applying a pixel mask
Removing Stars from Images
Overview
Changing the aperture properties
Viewing data for template samples
Properties
Overview
Markers
Background
Sample Data
Combining Images
Combine Image Set
Combine Files
Image combining methods
Setting normalization preferences
Region Statistics
Aligning Images
Overview
Aligning on a point
Aligning on an extended object
Aligning using the WCS calibration
High Precision Registration command
Overview
Tutorial: Introduction to Image Registration
Example using Point Tracking
Image Registration Properties
Image Registration Properties Dialog
Resampling
Background
Image Registration Residuals
Marker Properties
Resampling Properties
Tracking Properties
Image Region Statistics
Overview
Estimators
Editing & Copying
Overview
Duplicating to a new window
Emptying the Clipboard
Math Commands
Overview
Arithmetic Operations
Change Pixel Type
Combine Image Set
Fit Background
Image Arithmetic
Normalize Image
Region Statistics
Value Arithmetic
Using Plug-ins
Geometry Commands
Overview
Rotating Images
Rotation Menu
Affine Transformation
Expand Image
Imbed Image
Mirror
Reverse
Rotate
Rotate by angle (from Image Toolbar)
Scale
Shift
Transpose
Cropping to specified bounds
Cropping to the Image Cursor
Cropping by a rectangle
Align Horizontal
Align Vertical
Commands for RGB data
Working with RGB Images
Applying a palette to the pixels
Converting to RGB
Converting to Grayscale
Converting to 16 bit Luminance
Converting to 8 bit Luminance
Removing color noise
Extracting channels
Making an LRGB Image
Making an RGB Image
Commands for Intensity data
Extracting byte planes
Converting the data type
Converting to RGB
Undo
Using undo
Filtering Images
Overview
Removing Extreme Pixel Values
Cosmic Ray Filter
De-speckle Filter command
Clip High Values
Clip Low Values
Rank Filtering
Rank Filter
Median Filter
Minimum Filter
Maximum Filter
Smoothing Image Details
Rectangle Filter
Elliptical Filter
Binomial Filter
Gaussian Filter
Block Average Filter
Block Sum Filter
Sharpening Image Details
High Pass Filter
Unsharp Mask
Maximum Entropy command
Maximum Entropy example
Enhancing Edges
Gradient Filter
Rotational Gradient Filter
Rotational Gradient example
Laplacian Filter
Miscellaneous Kernel Filters
Custom Filter Kernel
3x3 Custom Filter command
Extracting Sources from Images
Using the MExtract Module
Definitions of Extracted Properties
Running the Pipeline
Tutorial
Preferences
Properties for the Source Extraction Tutorial
Source Extraction dialog
Procedure page
Background page
Detect page
Filter page
Match page
Difference page
Post Process page
Post Processing Options
Centroid Properties (Source Extraction)
FWHM Properties (Source Extraction)
PSF Fit Properties
Aperture Photometry Properties
Variability Properties
Image Diagnostics and Calculations
Overview
Create Synthetic Image
Bit Histogram Plot
Calculate Airmass
Calculate Gain
Calculate Image Scale
Calculate Julian Date
Calculate Readout Noise
Using the Pro Script Module
Working with scripts
Using the Script Editor Window
Using the Script Manager
Editing script keywords
Save Window Position
Related Topics
Script Editor Toolbar
Script Editor View Menu
Script Editor Context Menu
File Event Scripting
File Event Scripting
Tutorial: Introduction to File Event Scripting

Mira Pro x64 User's Guide

Prev Page Next Page

Tutorial: Introduction to Image Registration


This tutorial shows you how to register a set of images. If the alignment involves correcting a simple field drift or a wobble, then you can align them using a single reference point, such as a star. However, in this example we will need the full power of the Image Registration command because one of the images has a serious parallelogram distortion.

Overview

In this exercise we will use 4 images that are distorted relative to each other. The pictures below show the images combined into a single image using the Mean method (a simple average). Each picture shows features from all 4 images together to illustrate the misalignment at various places around the field of view. Within each picture, each group of 4 points shows the geometric relationship between the 4 images. It is clear that what works to align the images at one location does not work at another location. No simple combination of shifting, scaling, and rotation can bring the images into alignment.

To register this set of images requires using a 6-coefficient Shear transformation, which requires using at least 3 reference points per image. But relax! Even though this sounds complicated, Mira handles most of the details behind the scenes. You simply click some points on one image, let Mira copy them to the other images, and click two other buttons to calculate and apply the corrections. The Shear method corrects for image shift, field rotation, non-perpendicular axes, and unequal directional image scale. As you will see at the end of this tutorial, that is sufficient to correct the images so that features overlap with positional uncertainties near the 0.01 pixel level.

Getting Started

To begin, we need to open 4 images of Comet Hyakutake that are supplied with the Mira installation. These were acquired using a CCD camera with a 105 mm lens, which makes them severely undersampled with an image scale of 48 arcseconds per pixel. As you can see, point sources (stars) are typically 1 or 2 pixels wide.

Use the File > Open command or Ctrl+O to open the Open dialog, then select the 4 images as shown. Be sure that the Options box shows Open as Image Set so that the images open into the same window.

On the Open dialog, click [Open] and the images quickly appear in an Image Window as shown below.

As shown here, we have de-magnified the image and shrunk the window to keep the pictures small. You may choose to use a large window with a magnified image or to shrink the image or window. The display has no effect on the registration quality since the reference markers use coordinates computed from the image pixel data, not from the display.

Examining the Images

Before starting to register the images, it is a good idea first to examine the images you have opened. Perhaps some images are poor quality or you have other reasons not to align them. You need to do any image set adding, deleting, or re-ordering before starting Image Registration (or Photometry) since these adjustments are not allowed while the image registration package is active. To examine the images, move the red thumb on the trackbar in the Animation Bar at the bottom of the Image Window. If you need to change the images in the image set, see Changing the Image Set.

After making any necessary changes to the image set, we are ready to register the images. From the Image Window menu bar, select the Process > Alignment > Image Registration command. This opens the Image Registration Toolbar (see Image Registration ), as shown below (note: the Register command is also available as the left-most button on the Animation Bar). When you are finished with the image registration task, close the Registration Toolbar by re-using the Process > Alignment > Image Registration command that opened the toolbar. You can also close the Image Registration Toolbar by clicking the left-most button ("Register") on the Animation Bar.

Choosing Image Registration Settings

Image registration involves many complex mathematical operations that can be applied in different ways for different situations. This means that many different settings are involved in getting the best results. Typically, you will use the same settings so you will not need to change them each time you do registration. However, you should have a sense of what adjustments and settings are available and how they affect the results. Since this is your first time doing image registration, we will describe these settings in detail, even though you will most often just allow them to operate silently, behind the "front end" of the simple toolbar user interface.

On the Image Registration Toolbar, click to open the Marker Properties dialog. This is a standard dialog for setting "marker Properties" used by numerous Mira commands. This dialog contains the two pages shown below.

The pictures above show the Drawing and Centroid pages that control how registration markers are located and displayed on the images. The Drawing page shows that the markers are red and 1.5 points (about 2 screen pixels) wide. The Centroid page settings are more important because they are somewhat sensitive to the size of the alignment points on the images. These settings may change for different image sets. In this case, the images are severely undersampled and the stars are only 1 to 2 pixels wide. We use a small centroid Sample Radius of 3 pixels. The tracking Search Radius should be 20 to allow for 20 pixel jumps between images. You should seldom use a value larger than around 40. In addition, these are positive images with brighter things having higher signal, so the centroiding should be done on Positive objects, not Negative ones. You should be using the centroid positions near where you click the mouse on the image. You do not have to use centroiding; instead Mira can use the exact coordinate where you click the mouse. But the Centroid First Point box must be checked if you intend to use centroid positions. When you are finished changing the settings, click [OK] or [Apply] to save your changes.

Next we need to check the method used to align the images. Click the button on the Registration Toolbar to open the Image Registration Properties dialog. This dialog has 2 pages, as shown below.

Settings in the Image Registration Properties dialog control the way image pixels are aligned and how the pixel values are adjusted to fall at new coordinates in the output images.

  • Settings on the Resampling page relate to "resampling", which is the mathematical prescription for how the pixel values are shifted and adjusted to align the images.

  • Settings on the Background page control the "background" (or empty image) into which the resampled images are saved. For example, you can choose to make the registered images larger than the original ones or have all the aligned images offset as a group.

For the purposes of this tutorial, make all settings like those shown above. After making changes, click [OK] or [Apply] to save the calibration and close the dialog. We are now ready to mark registration points.

Marking Reference Points

Enter marking mode by clicking the button on the Registration Toolbar and use the Animation Bar to go to the first image. Now click on each of the 6 stars shown in the picture below. As you move onto each point, notice how the magnified view on the Image Bar shows you a detailed view, despite the lack of detail in the de-magnified image. You can adjust the tracker magnification by double-clicking on it to cycle between 1, 3, 5, and other factors of the original image.

Seldom do you ever need more than 5 to 6 reference points on a CCD image. So long as the markers are well separated and placed all around the image, the alignment should be excellent. But you must choose "good" points. The points we have marked were chosen using these guidelines:

  • bright, but not saturated.

  • fairly round.

  • well isolated from nearby stars.

  • distributed around the image, rather than clumped in one area.

At any time when marking points, you can change your mind about using a point. For example, you can delete a point and re-mark it elsewhere, or just delete it. To delete a point, press the [Shift] key and click on it. In Marking Mode, the [Shift] key is a toggle in marking mode: up=mark, down=delete), or you can click the button on the Registration Toolbar to enter Delete mode. Using the [Shift] key is usually easiest.

With 6 points marked on the first image, we now need to transfer them to the other images. There are 2 ways to do this:

  • Change to the next image and click on all the same objects, repeating this for each image in the stack. This requires knowing or remembering which ones you chose, and perhaps flipping back and forth between images. However, it will be necessary if the images are misaligned in rotation.

  • If each successive image is within 10 to 20 pixels of alignment with the previous one, Mira can track the points through the remaining images. Tracking computes a new centroid position in each image based on its location in the previous image.

The Tracking method is best if the images are not far out of alignment. To track the objects marked in the current image, click the button on the Registration Toolbar. There are now reference points in all images, but how well matched are they? The only way to answer this is to look through the image set. Using the Animation Bar, step through each image and check if it found the correct point in all images. Actually, in this crowded field of undersampled stars, you will find that some of the markers locked onto the wrong star which is nearby the correct one. This is easily fixed! In Marking mode, click on the correct point to mark it, then press [Shift] and click on the misplaced marker to delete it. Repeat this for any incorrectly placed markers that you find. When you have scanned through all the markers on all images, you are ready to compute the mathematical equations that will be used to shift the pixels around and bring the images into alignment.

Computing the Registration Equations

To compute the registration equations, click on the Registration Toolbar. This matches all the patterns of markers to figure out which marker goes with which one in all other images and then calculated the coefficients of the equations that match them. This produces a listing of results in a Text Editor window:

The "Reference Image" chosen on the Resampling page is considered to be perfect and the other images are matched to it (by default, the reference image is the first image). For each of the other images Mira reports "x terms" and "y terms" plus RMS errors. For example, the x terms are x0, dCx, and dCy and the RMS error in the x axis is labeled x rms. Since the registration is optimal but not exact, there will be a small difference between the position of same detail in a pair of images. In this case, the RMS errors are typically around 0.02 to 0.03 pixels. That's extremely good! If any of the values were unusually large, it would probably indicate an erroneously marked point. You could go back through the images and find the bad point(s), delete and re-mark them, and click again to repeat the calculation. If you find that this does not help, use the button to open the Resampling page and click the [Residuals] button. The table of residuals shows the difference in calculated position and actual position for each registered point. A bad point has a residual that is very high or very low. Identify the bad points, close the dialog, and delete those points from the images. Then re-calculate.

Tip

You can copy, save, or print the Residuals table using commands in its menu. To open its menu, right click on the table. For example, you might copy the table, close the dialog, and then paste the table into a Mira Text Editor.

Now that the transformation equations are known, we perform the alignment. Click to do this and you are done with image registration. Use the Animation Bar to scan through the image set and view the results. If everything is well aligned (and it should be because the RMS errors were small), then you may wish to Save or Combine the images (or do whatever you wanted to do with registered images). If you don't like the result, Undo the alignment and try again.

Assessing the Results

Let's look at how well aligned the images have become. Stars in the original images were tiny but mostly round. Therefore, the image of a star should also be round if all the images are added together. If not, then one or more of the images is not well aligned. Use the Process > Combine command in the menu (or on the Animation Bar) to combine the images using the Mean method. The picture below shows a region magnified 16x. This shows one very bright star and a number of faint ones. Notice that each of the combined stars is extremely round, which means excellent alignment.

As a final step, let us assess the quality of the alignment in a quantitative way. Measure the centroid position for a nicely isolated star that was not one of the registration points. The positions should be very similar. In this example, use the star at coordinate (661,603). Move the Image Cursor to near this point and click the Centroid button to compute the centroid positions in all images. The result looks like this:

The X and Y centroid positions are intensity weighted coordinates of high precision. How do they compare? Actually, they are in good agreement for all 4 images. The typical difference, or scatter among them, gives an estimate of the exactness of the alignment. You can use your hand calculator to compute the mean and error of the mean position for each of the X and Y coordinates. We will quote the error of the mean, which is the uncertainty in the mean value; it is the standard deviation divided by the square root of the number of values in the sample (4 images gives a sample size of 4). Here are the results:

  • Mean X = 660.982 +- 0.0069

  • Mean Y = 602.995 +- 0.0053

The x and y coordinates agree to better than 1/100th of a pixel. That is good enough to write a book about!

Summary

In this tutorial we have taken a grand tour of image registration. We have examined all the factors that can affect the quality of the outcome. Typically, registration can be done much faster than in this case. If the images are of good quality and not far out of alignment, most of the detailed examination shown here can be bypassed and the results will be excellent.

Related Topics

Contents

Tutorials

Image Registration

Saving Image Sets

Combine Image Set

Image Registration Properties


Mira Pro x64 User's Guide, Copyright Ⓒ 2023 Mirametrics, Inc. All Rights Reserved.