(* Záření v jednom směru *)
W=KroneckerProduct[IdentityMatrix[3],Transpose[{{1},{1},{1}}]]

(* Záření ve směru kolmém na první *)X= KroneckerProduct[{{1,1,1}},{{1,0,0},{0,1,0},{0,0,1}}]

(* Matice reprezentující průchod záření objektem *)
A=Join[W,X]

A//MatrixForm

(* Kočička před tomografem *)
kocicka = Image[]

(* Kočička natlačená do tomografu *)
kocicka = Image[ImageResize[kocicka,{300,300}]]

(* Rozdělení kočičky na mřížku odpovídající citlivosti tomografu. *)
kocickaParts = ImagePartition[ColorConvert[#,"Grayscale"]&@kocicka,100,100]

(* "Sumace" dat pro jednotlivé buňky mřížky *)
vysledek = Flatten[Map[ImageData[ImageResize[#,{1,1}],"Byte"]&,Flatten@kocickaParts]]

vysledek//MatrixForm

B=PseudoInverse[A]

(* Výsledky kočičky po tomografu *)
A.vysledek