Real-world river networks are complex, depending on landscape features such as elevation, aspect, and lithology. Moreover, the fractal nature of river networks means that they are approximately scale-free (Tarboton et al., 1988) without a well-defined spatial resolution at which they should be represented (Davies & Bell, 2009;Yamazaki et al., 2009). As a result, hydrologic models often use a conceptual representation of river networks. To date, there are mainly three methods for representing river networks at different scales, each with its own advantages and disadvantages. However, limitations remain when representing river networks due to the resolution mismatch and their interactions with other hydrologic features (e.g., ocean).Two methods are useful at the watershed/regional scale, the first of which is the vector-based river networks analysis method. This method uses vector datasets to represent the river networks and their topological relationships (Lin et al., 2018;Mizukami et al., 2016). Vector datasets are often provided by public agencies, for example, the United States Geological Survey (USGS), or are digitized from satellite image processing, for example, vectorization-based river channel extraction. Various graph theory algorithms are then used to perform quality control and network analysis based on the vector river networks (Lindsay et al., 2019). Lindsay et al. reviewed several potential issues in existing vector river networks datasets: (a) the vertex coordinates of river flowlines may not exactly overlap with the actual locations due to digitization error and floating-point rounding errors. (b) A river flowline's starting and ending vertices may be reversed during spatial analysis, resulting in an opposite flow direction. (c) vector datasets obtained from different sources generally use different spatial references and cannot be used or combined directly. Even without these issues, vector datasets still require other pre-processing before use. For example, braided rivers are not universally supported in hydrologic models as multiple flow