Perspective factorization algorithm
Assume that Li are known, then PM is known.
Use the singular value decomposition
PM=US VT
In the noise-free case
S=diag(s1,s2,s3,s4,0, … ,0)
and a reconstruction can be obtained by setting:
P=the first four columns of US.
M=the first four rows of V.