The intertidal mudflats that are exposed subaerially between the high and low tides have a complex morphology comprising ebb and flood tidal deltas, spits, shoals, and tidal sandbanks intervened with a network of meandering and branching tidal channels and inlets. The intertidal zones are highly transient and understanding their morphodynamics is critical to the ecology and economy of the region. Monitoring intertidal mudflat dynamics is imperative for sustainable development activities along coastal regions. The topography of intertidal mudflat is estimated for the eastern coast of the Gulf of Khambhat (north-western coast of India) using the waterline detection method from a series of satellite images. We have used Modified Normalized Difference Water Index and image threshold techniques to delineate waterlines from satellite images of high to low tidal conditions at different tidal heights. The waterlines are considered as elevation contours corresponding to tidal height during satellite observation, and a series of elevation contours delineated from satellite images of various tidal heights are then used to calculate the mudflat topography. The estimated topography is validated with the Differential Global Positioning System (DGPS), the Sound and navigation ranging (Sonar), and ICESat-2 ATL03 dataset observations, which results in R2 values of 0.97, 0.89, and 0.9, respectively, as well as RMSE values of 0.39 m, 0.33 m, and 0.49 m, respectively. The mudflat dynamics are studied by comparing the topographical changes estimated between the satellite images of 1998-99 and 2019. The major morphological changes are the landward shift of the low waterline, lateral migration of tidal bars, narrowing of the river mouth, as well as anthropogenic influences such as mudflat reclamation and impedance of long-shore sediment transport. The region with deposition of sediment above Mean Sea Level (MSL) has flattened the slope, making the mudflat stable and promoting mangroves expansion, while erosion has steepened the mudflats with loss of mangrove cover. These mudflats exhibit highly complex morphological changes mainly induced by sediment transport, two-way tidal currents, and fluvial fluxes. The estimated topography shall find its application in modeling coastal processes, identifying vulnerable coastal stretches, and site suitability for mangrove plantations.