Thanks for J.M's suggestions, with the help of FindIntersectionsGraphics`Mesh`FindIntersections
newSolution[mat1_, mat2_] := Module[{graph, pts}, graph = ParametricPlot[ {mat1.{Sin[θ], Cos[θ]]Cos[θ], 1}, mat2.{Sin[θ], Cos[θ], 1}}, {θ, 0, 2 Pi}]; pts = FindIntersections[First[graph]]; ParametricPlot[ ,Epilog -> {mat1.Point[{Sin[θ], Cos[θ].1, 1.2},]}]; pts = mat2.{Sin[θ],Graphics`Mesh`FindIntersections[First[graph]]; Cos[θ], 1}},graph {θ,/. 0, 2 Pi}, (Epilog -> _) -> Epilog -> {Red, PointSize[Medium], Point[pts]}] ]