Tortuosity of corneal nerve fibers acquired by in vivo Confocal Microscopy (IVCM) are closely correlated to numerous diseases. While tortuosity assessment has conventionally been conducted through labor-intensive manual evaluation, this warrants an automated and objective tortuosity assessment of curvilinear structures. This paper proposes a method that extracts the image-level features for corneal nerve tortuosity grading. For an IVCM image, all corneal nerve fibers are first segmented and then, their tortuosity are calculated by morphological measures. The ordered weighted averaging (OWA) approach, and the k-Nearest-Neighbor guided dependent ordered weighted averaging (kNNDOWA) approach are proposed to aggregate the tortuosity values and form a set of extracted features. This is followed by running the Wrapper method, a supervised feature selection, with an aim to identify the most informative attributes for tortuosity grading. Validated on a public and an in-house benchmark data sets, experimental results demonstrate superiority of the proposed method over the conventional averaging and length-weighted averaging methods with performance gain in accuracy (15.44% and 14.34%, respectively). The simultaneous use of multiple aggregation operators could extract the image-level features that lead to more stable and robust results compared with that using average and length-weighted average. The OWA method could facilitate the explanation of derived aggregation behavior through stress functions. The kNNDOWA method could mitigate the effects of outliers in the image-level feature extraction.