Search Astrofriend's homepage:
 All pages shall now have been validated Info Cookies (Kakor) / GDPR

• Svenska
 Follow Astrofriend Travel News 20230818 Travel Travel to Turku (Åbo) in Finland, 2023 20230722 Historia, Stockholm Cykeltur till Brunnsviken 20230311 Travel Travel itinerary for a round trip in Southern Germany, 2023 20221028 Travel Croatia Five coast cities in Northern Croatia 20220821 Travel Germany By train in Northern Germany 20220730 Travel Sweden A visit to the island Gotland in Seden, 2022 2021201 Travel Fuerteventura, Spain 2021 20211028 Travel: Crete, Greek 2021 20210815 Travel: Cruise Sweden: From Stockholm to Härnösand and Höga Kusten 20200902 Travel in Sweden: The village Röttle in Sweden, 2007 20200723 Sveriges Historia: Befästningslinjen Södra Fronten Full list - - All News

Tutorial: AstroImageJ Align and Drizzle

Content:

Note:
I take no responsibility or liability for what are written here, you use the information at your own risk!

4, Inverse matrix:

Inverse matrix:

In part one I talked about the matrix T that transform our sub images to fit the reference image. But how to find out what T is?

We need to know how to calculate the matrix T and its internal parameters.

X' = T X.

With help of the invers matrix X-1 we can get T alone on the right side, in matrix manipulation we normally can't change order of the matrices. It looks a bit strange if you are not used with it.

X-1 is of the construction that X-1 X give the unit matrix I (I=1). X must be a square matrix!

X' X-1 = T X X-1 = T I = T

T = X' X-1

Now we get one more problem, to calculate X-1.

Later when we have find T values for each sub image we can put in each sub images coordinates one pixel by one to the transformation T matrix and get the new coordinates in the reference image. In more advanced programs like Matlab we put in the whole image, in more general software like AstroImageJ we have too work with loops.

Inverse matrix, a deeper look inside:

This was the easy case when we have square matrix and an unique solution.

Generalized inverse:

Many times we don't have square matrices, it could be that we have more reference stars than our translation matrix need, they are called to be over-determined. Or the reference points maybe depends on each other, maybe we by mistake have chosen the same star twice.

We can solved that too, but this time it's a mean value more than exact. This is the most common situation, the real life.

Take the example of a warp matrix from the first page:

Z t = X'

Lets say that we have taken 50 reference stars, much more than needed for our translation matrix t, in the example with second degree it's only needs 12 equations, not 50.

You cant just use the inverse of Z like (left side): Z-1 Z t = Z-1 Z t = (Z-1 Z) t = I t = t.

The generalized inverse of Z is the Z-1left matrix or it's components (ZT Z)-1 ZT. ZT is the transpose of the matrix Z. There is an Z-1right too.

Then we have that: Z-1left Z t= Z-1left X' : t = (ZT Z)-1 ZT X'.

In Matlab this is very simple, looks like this: t = X' \ Z.

Transpose of a matrix, in depth:

Generalized inverse, in depth:

 Next part: V, Determinant