An inverse problem of simultaneous estimation of radiation phase function, albedo and optical thickness