Computing the Differential Characteristics of Isointensity Surfaces