Abnormal tumor hemodynamics are a critical determinant of a tumor's microenvironment (TME), and profoundly affect drug delivery, therapeutic efficacy and the emergence of drug and radio-resistance. Since multiple hemodynamic variables can simultaneously exhibit transient and spatiotemporally heterogeneous behavior, there is an exigent need for analysis tools that employ multiple variables to characterize the anomalous hemodynamics within the TME. To address this, we developed a new toolkit called HemoSYS for quantifying the hemodynamic landscape within angiogenic microenvironments. It employs multivariable time-series data such as in vivo tumor blood flow (BF), blood volume (BV) and intravascular oxygen saturation (Hbsat) acquired concurrently using a wide-field multicontrast optical imaging system. The HemoSYS toolkit consists of propagation, clustering, coupling, perturbation and Fourier analysis modules. We demonstrate the utility of each module for characterizing the in vivo hemodynamic landscape of an orthotropic breast cancer model. With HemoSYS, we successfully described: (i) the propagation dynamics of acute hypoxia; (ii) the initiation and dissolution of distinct hemodynamic niches; (iii) tumor blood flow regulation via local vasomotion; (iv) the hemodynamic response to a systemic perturbation with carbogen gas; and (v) frequency domain analysis of hemodynamic heterogeneity in the TME. HemoSYS (freely downloadable via the internet) enables vascular phenotyping from multicontrast in vivo optical imaging data. Its modular design also enables characterization of non-tumor hemodynamics (e.g. brain), other preclinical disease models (e.g. stroke), vascular-targeted therapeutics, and hemodynamic data from other imaging modalities (e.g. MRI).