We used intrinsic Second Harmonic Generation (SHG) by fibrillar collagen to visualize the three-dimensional architecture of collagen fibrosis at the micrometer scale using laser scanning nonlinear microscopy. We showed that SHG signals are highly specific to fibrillar collagen and provide a sensitive probe of the micrometer-scale structural organization of collagen in tissues. Moreover, recording simultaneously other nonlinear optical signals in a multimodal setup, we visualized the tissue morphology using Two-Photon Excited Fluorescence (2PEF) signals from endogenous chromophores such as NADH or elastin. We then compared different methods to determine accurate indexes of collagen fibrosis using nonlinear microscopy, given that most collagen fibrils are smaller than the microscope resolution and that second harmonic generation is a coherent process. In order to define a robust method to process our three-dimensional images, we either calculated the fraction of the images occupied by a significant SHG signal, or averaged SHG signal intensities. We showed that these scores provide an estimation of the extension of renal and pulmonary fibrosis in murine models, and that they clearly sort out the fibrotic mice.