Matrix2Image := proc(  A :: Matrix )
   local P, m, n, i, j;

   m, n := LinearAlgebra:-Dimension(A);
   P := Array(1..m,1..n,datatype=float[8],storage=rectangular,order=C_order):
   for i from 1 to m do
      for j from 1 to n do
         P[i,j] := A[i,j]
      end do
   end do;

   return P
end proc;