Polarization dependence second harmonic generation (P-SHG) microscopy is gaining increase popularity for in situ quantification of fibrillar protein architectures. In this report, we combine P-SHG microscopy, new linear least square (LLS) fitting and modeling to determine and convert the complex second-order non-linear optical anisotropy parameter ρ of several collagen rich tissues into a simple geometric organization of collagen fibrils. Modeling integrates a priori knowledge of polyhelical organization of collagen molecule polymers forming fibrils and bundles of fibrils as well as Poisson photonic shot noise of the detection system. The results, which accurately predict the known sub-microscopic hierarchical organization of collagen fibrils in several tissues, suggest that they can be subdivided into three classes according to their microscopic and macroscopic hierarchical organization of collagen fibrils. They also show, for the first time to our knowledge, intrahepatic spatial discrimination between genuine fibrotic and non-fibrotic vessels. CCl4-treated livers are characterized by an increase in the percentage of fibrotic vessels and their remodeling involves peri-portal compaction and alignment of collagen fibrils that should contribute to portal hypertension. This integrated P-SHG image analysis method is a powerful tool that should open new avenue for the determination of pathophysiological and chemo-mechanical cues impacting collagen fibrils organization.