i.pan.sharpen contains python code to do an inverse pca. It uses some of the 
new r.mapcalc capabilities to do parallel processing on multi-core processors. 
So it's pretty fast.

Michael
____________________
C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University

voice:  480-965-6262 (SHESC), 480-727-9746 (CSDC)
fax:          480-965-7671 (SHESC),  480-727-0709 (CSDC)
www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu











On Mar 26, 2013, at 4:24 AM, 
<[email protected]<mailto:[email protected]>>
 wrote:

From: Martin Landa <[email protected]<mailto:[email protected]>>
Subject: [GRASS-dev] calculating inverse matrix
Date: March 26, 2013 4:20:08 AM MST
To: GRASS developers list 
<[email protected]<mailto:[email protected]>>


Hi all,

I found several methods in GRASS code which compute inverse matrix

1) v.generalize

http://trac.osgeo.org/grass/browser/grass/trunk/vector/v.generalize/matrix.c#L134

2) i.orto.photo

http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.ortho.photo/lib/m_inverse.c#L16

3) GMath lib (depends on lapack/blas)

http://trac.osgeo.org/grass/browser/grass/trunk/lib/gmath/la.c#L566

It would be nice to consolidate it including data structures (modules
noted above define their own data structures for matrix). In other
worlds to incorporate v.generalize/i.orto.photo code into GMath lib.
Any pointers here?

Thanks, Martin

--
Martin Landa <landa.martin gmail.com<http://gmail.com/>> * 
http://geo.fsv.cvut.cz/~landa


_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev

Reply via email to