Many radar-gauge merging methods have been developed to produce improved rainfall data by leveraging the advantages of gauge and radar observations. Two popular merging methods, Regression Kriging and Bayesian Regression Kriging were utilized and compared in this study to produce hourly rainfall data from gauge networks and multi-source radar datasets. The authors collected, processed, and modeled the gauge and radar rainfall data (Stage IV, MRMS and RTMA radar data) of the two extreme storm events (i.e., Hurricane Harvey in 2017 and Tropical Storm Imelda in 2019) occurring in the coastal area in Southeast Texas with devastating flooding. The analysis of the modeled data on consideration of statistical metrics, physical rationality, and computational expenses, implies that while both methods can effectively improve the radar rainfall data, the Regression Kriging model demonstrates its superior performance over that of the Bayesian Regression Kriging model since the latter is found to be prone to overfitting issues due to the clustered gauge distributions. Moreover, the spatial resolution of rainfall data is found to affect the merging results significantly, where the Bayesian Regression Kriging model works unskillfully when radar rainfall data with a coarser resolution is used. The study recommends the use of high-quality radar data with properly spatial-interpolated gauge data to improve the radar-gauge merging methods. The authors believe that the findings of the study are critical for assisting hazard mitigation and future design improvement.