1. A method for mapping color gamut, comprising: defining at least one hue leaf associated with a source device with respect to at least two arbitrary color gamuts utilizing a vector math function in order to thereafter determine a most saturated point with respect to said at least one hue leaf, wherein said at least one hue leaf is mapped with a color gamut of a target space based on a continuous, one-one and onto function from an edge of said at least two color gamuts, wherein a polygon is selected for computing a color transformation based on a location of a selected point, if said selected point is in a safe area, above an upper polygon, or below a lower polygon; estimating said safe area with respect to an intersection point by approximating said most saturated point in order to thereafter continuously sub-divide an upper hull and a lower hull associated with at least one hue leaf into at least one section by constructing at least one vector, wherein said safe area is estimated with respect to an intersection point of hue leaf line segments by approximating a most saturated point in said at least one hue leaf, wherein said approximating said most saturated point comprises utilizing a vector function comprising p=aV p1 +bV p2 +V p0 , wherein p comprises a plane, wherein scalar point a comprises vector V p1 constrained to a Lab space from point (0, 0, 0) to point (100, 0, 0), wherein scalar point b comprises vector V p2 constrained to a Lab space from point (0, 0, 0) to a point to be mapped (point M), and wherein V p0 comprises a displacement vector; and determining an appropriate polygon section from a selected point in said safe area to compute a vector relationship with respect to said at least two arbitrary color gamuts in order to map said at least two arbitrary color gamuts based on said continuous, one-one and onto function, thereby creating an invertible transformation.