Sea ice models have become essential components of weather, climate, and ocean models. A realistic representation of sea ice affects the reliability of process representation, environmental forecast, and climate projections. Realistic simulations of sea ice kinematics require the consideration of both large-scale and finescale geomorphological structures such as linear kinematic features (LKF). We propose a multiscale directional analysis (MDA) that diagnoses the spatial characteristics of LKFs. The MDA is different from previous analyses in that it (i) does not detect LKFs as objects, (ii) takes into account the width of LKFs, and (iii) estimates scale-dependent orientation and intersection angles. The MDA is applied to pairs of deformation fields derived from satellite remote sensing data and from a numerical model simulation with a horizontal grid spacing of ~4.5 km. The orientation and intersection angles of LKFs agree with the observations and confirm the visual impression that the intersection angles tend to be smaller in the satellite data compared to the model data. The MDA distributions can be used to compare satellite data and numerical model fields using conventional metrics such as a Euclidean distance, the Bhattacharyya coefficient, or the Earth mover’s distance. The latter is found to be the most meaningful metric to compare distributions of LKF orientations and intersection angles. The MDA proposed here provides a tool to diagnose if modified sea ice rheologies lead to more realistic simulations of LKFs.