River water surface extent can be extracted from optical and radar satellite images; this is useful for estimating streamflow from space. The radiation characteristics of open water from the visible and microwave bands are different and provide independent information. In this study, for the purpose of improving streamflow estimation from space for data-sparse regions, a method that combines satellite optical and radar images data for streamflow estimation using a machine learning technique was proposed. The method was demonstratedthrough a case study in the river segment upstream of the Ganzi gauging station on the Yalong River, China. Utilizing the support vector regression (SVR) model, the feasibility of different combinations of water surface area derived from Sentinel-1 synthetic aperture radar images (AREA_SAR), modified normalized difference water index derived from Landsat 8 images (MNDWI), and reflectance ratios between NIR and SWIR channels derived from MODIS images (RNIR/RSWIR) for streamflow estimation were evaluated through three experiments. In Experiment I, three models using AREA_SAR (Model 1), MNDWI (Model 2), and a combination of AREA_SAR and MNDWI (Model 3) were built; the mean relative error (MRE) and mean absolute error (MAE) of streamflow estimates corresponding to the SVR model using both AREA_SAR and MNDWI (Model 3) were 0.19 and 31.6 m3/s for the testing dataset, respectively, and were lower than two models using AREA_SAR (Model 1) or MNDWI (Model 2) solely as inputs. In Experiment II, three models with AREA_SAR (Model 4), RNIR/RSWIR (Model 5), and a combination of AREA_SAR and RNIR/RSWIR (Model 6) as inputs were developed; the MRE and MAE for the model using AREA_SAR and RNIR/RSWIR (Model 6) were 0.25 and 56.5 m3/s, respectively, which outperformed the two models treating AREA_SAR (Model 4) or MNDWI (Model 5) as single types of inputs. In Experiment III, three models using AREA_SAR (Model 7), MNDWI, and RNIR/RSWIR (Model 8) and the combination of AREA_SAR, MNDWI and RNIR/RSWIR (Model 9) were built; combining all three types of satellite observations (Model 9) exhibited the highest accuracy, for which the MRE and MAE were 0.18 and 18.4 m3/s, respectively. The results of all three experiments demonstrated that integrating optical and microwave observations could improve the accuracy of streamflow estimates using a data-driven model; the proposed method has great potential for near-real-time estimations of flood magnitude or to reconstruct past variations in streamflow using historical satellite images in data-sparse regions.