As a unique and important ecosystem, tidal flats provide a variety of ecosystem functions and services. Mapping tidal flats is essential for the protection and management of coastal ecosystems. However, large-scale tidal flats mapping still faces challenges due to the tidal variation and spectral similarity between tidal flats and inland wetlands. Previous methods rely on the coastlines or maximum seawater extent to exclude inland areas, which is limited by its inability to effectively differentiate tidal flats from spectrally similar inland wetlands. To address these issues, we proposed a new tidal-flat mapping method by integrating Sentinel-2 time series imagery with tide height (TH) data from ground-based tide stations on Google Earth Engine. We first generated images at the lowest and highest tidal stages, established the statistical relationship between the Normalized Difference Water Index (NDWI) of each pixel and TH, and then concatenated them into a Random Forest classifier for further classification. The statistical relationship between NDWI and TH amplified the difference between tidal flats and inland wetlands, thus significantly reducing the influence of spectral similarity. This method could produce a high-precision tidal flats map with an overall accuracy of 97.30% in the coastal zone of the Chinese mainland. By quantitatively comparing with the previous tidal flat maps, we found that the strategies of tidal-level information simulation and inland area exclusion were the two main reasons producing the differences among the maps. The proposed method does not rely on space constraints to exclude inland wetlands and can capture more estuarine tidal flats, so it can be used as a reliable means to monitor the tidal flats in large-scale areas.