The coastal zone is an area that includes the sea coast and adjacent parts of the land and sea, where the mutual interaction of these environments is clearly marked. Hence, the modelling of the land and seabed parts of the coastal zone is crucial and necessary in order to determine the dynamic changes taking place in this area. The accurate determination of the terrain in the coastal zone is now possible thanks to the use of Unmanned Aerial Vehicles (UAVs) and Unmanned Surface Vehicles (USVs). The aim of this article is to present land and seabed surface modelling in the coastal zone using UAV/USV-based data integration. Bathymetric and photogrammetric measurements were carried out on the waterbody adjacent to a public beach in Gdynia (Poland) in 2022 using the DJI Phantom 4 Real Time Kinematic (RTK) UAV and the AutoDron USV. As a result of geospatial data integration, topo-bathymetric models in the coastal zone were developed using the following terrain-modelling methods: Inverse Distance to a Power (IDP), kriging, Modified Shepard’s Method (MSM) and Natural Neighbour Interpolation (NNI). Then, the accuracies of the selected models obtained using the different interpolation methods, taking into account the division into land and seabed parts, were analysed. Research has shown that the most accurate method for modelling both the land and seabed surfaces of the coastal zone is the kriging (linear model) method. The differences between the interpolated and measurement values of the R95 measurement are 0.032 m for the land part and 0.034 m for the seabed part. It should also be noted that the data interpolated by the kriging (linear model) method showed a very good fit to the measurement data recorded by the UAVs and USVs.