[SciPy-User] helmert implementation (7 parameters - geometric transformation)

Pauli Virtanen pav+sp at iki.fi
Wed Mar 3 10:22:07 EST 2010


Wed, 03 Mar 2010 00:29:42 +0100, Massimo Di Stefano wrote:
> from here :
> 
> http://docs.scipy.org/doc/numpy/reference/generated/
> numpy.linalg.cholesky.html
> 
> seems it is applicable only on a square matrix. is that true ?

Cholesky decomposition works only on Hermitian matrices -- so being 
square is one requirement.

Also, if you are going to do least-squares fitting, use linalg.lstsq; it 
does the same thing as (A^T A)^-1 A^T b, but should be numerically more 
stable.

-- 
Pauli Virtanen




More information about the SciPy-User mailing list