Integrated Geodetic and Hydrographic Measurements of the Yacht Port for Nautical Charts and Dynamic Spatial Presentation

: The use of geodetic and hydrographic systems, methods and measuring devices to obtain geospatial data enables integrated modeling of 3D space and related objects. Factors such as the seaﬂoor relief of the marina and land infrastructure need to be considered to ensure safe stay of vessels in the marina. The article presents the results of bathymetric and geodetic measurements of the yacht port of the National Sailing Centre in Gda´nsk. They were made with the use of several measuring methods such as a hydrographic unmanned surface vessel, GNSS positioning, enabling geodetic inventory of hydrotechnical structures and port infrastructure, as well as laser scanning that resulted in three-dimensional presentation of land objects. As a result of research and analysis, a two-dimensional bathymetric chart of the yacht port was made. Thanks to the integration of geodetic and hydrographic data, three-dimensional visualization of the marina area was obtained. The charts can be used to ensure ship navigation safety when maneuvering in port and mooring at quays.


Introduction
An important aspect of the infrastructure safety of a yacht port is bathymetric information. It is part of a broader term, i.e. navigational information [1,2]. In order for this information to be useful, it needs to satisfy a number of basic conditions. Apart from the substantive (thematic) content, the information needs to have features that facilitate its verification and proper use, for example through an assessment of its reliability [3]. Moreover, the information should meet certain general requirements arising from users' needs. To meet these requirements, the information should be relevant, complete, accessible, timely, up-to-date, easy to use (to image), easy to store, and correctable. In order to ensure that it is up to date, periodic bathymetric plans for a water area are devised every two years for irregular use of water areas and marine structures [4,5].
Hydrographic surveys which provide a basis for the creation of a navigational chart comprise [6,7] bathymetric surveys of the sea bottom and onshore surveys [8,9] of quays and bitts with an accuracy specified in [10,11] depending on the sounding (water area) category.
As the port is approaching and the water area limitations are taken into account, bathymetric surveys are carried out using ever-smaller sounding vessels. Hydrographic ships and vessels carry out bathymetric surveys in open water areas, at great depths and in the 24/7 system. In the coastal area and in major harbour basins, surveys are carried out using hydrographic cutters and motor boats whose size prevents them from maneuvering within a marina, particularly where vessels are at berth. For this reason, depth measurements in small yacht harbours are taken on an irregular basis, which is why such marinas often have no up-to-date information of the seabed surface. In recent years, unmanned maritime vessels have been developing very rapidly. Unmanned surface vessels (USVs) are a variety of design solutions in the construction of the hull and the boat propulsion: single hull, double hull, with a screw or screwless propulsion with a small draft [12]. Due to the unprecedented functionalities, they found a wide application in many measurement applications in the second decade of the 21st century. One of the most important of these applications are hydrographic measurements: port basins, lakes, rivers and small reservoirs, where their aim is to measure the seafloor relief with the appropriate accuracy [13,14]. In addition, USVs are increasingly used for tasks related to supporting the navigation process [15], in underwater photogrammetry [16], or in geological works [17]. Owing to the significant increase in accuracy of GNSS (Global Navigation Satellite System) systems' position, unmanned hydrographic vessels no longer have to use expensive GNSS RTN (Real Time Network) systems (accuracy of 2-3 cm, p = 0.95) but can successfully use DGPS (Differential Global Positioning System) [18,19] and EGNOS (European Geostationary Navigation Overlay Service) [19] systems (accuracy of 2 m, p = 0.95) or even multi-system multi-GNSS receivers (accuracy of 3-5 m, p = 0.95) [20]. An equally important aspect of progress has been the emergence of miniature multibeam echosounders which have helped to significantly increase the efficiency of AUV (Autonomous Underwater Vehicle) hydrographic work [21].
The onshore part of hydrographic surveys includes geodetic inventory of quays and hydro-engineering devices. Nowadays, the most effective method of taking an inventory is the use of laser scanning. It ensures cloud point accuracy at a level of 2-5 cm [22][23][24]. Thanks to advanced error-reduction methods [25], cartographic imaging algorithms [26] and the effective processing of a large amount of data reduction [27], it is possible to precisely represent the 3D shape of objects from single objects [28] to that of an entire city [29]. An alternative to laser scanning is the method of quay inventory using GNSS networks with a measurement accuracy of 2-5 cm, p = 0.95 [30]; however, its limitation is the lack of the option of 3D modelling of objects.
The acquired geospatial data is usually presented in the form of two-dimensional cartographic products, such as geodetic maps, nautical and bathymetric charts. Digital data is also used to present in SENC (System of Electronic Navigational Chart) onboard vessels. The 3D capability can be utilised in laboratory as well as operational use. A chart, which depicts three-dimensional data, is the best vehicle to visualise real world [31,32]. Dynamic perspective projection [33][34][35][36] can be particularly useful when presenting the surroundings during the movement of a floating vessel [37].
Geodetic and hydrographic surveys were carried out in the yacht port of the National Sailing Centre at the Gdańsk University of Physical Education and Sport. It is located at the mouth of thé Smiała Wisła river ( Figure 1) naturally surrounded by a landscape park, several kilometres to the east of the TriCity agglomeration. The port has 51 berths at the quay and floating landing stage. On land, there are lecture and seminar rooms, a shed, a biological regeneration complex with a sauna, a navigation and meteorological office with a look-out tower, a sailmaker's workshop with sail sewing machines, and a boatbuilding and mechanical workshop. In 2015, the marina was reconstructed. In addition to infrastructure upgrading, the basin was deepened to a minimum depth of 4.5 m. Prior the reconstruction, the depths shallow by 50 cm.
In the research: scanning method for hydrography, multisensory bathymetry, 3D visualization for navigation, USV surveys in restricted area and surveys between vessels and yachts are advances in bathymetry. In the research: scanning method for hydrography, multisensory bathymetry, 3D visualization for navigation, USV surveys in restricted area and surveys between vessels and yachts are advances in bathymetry.

Materials
The project was carried out using measuring instruments and systems in three areas: a hydrographic system to investigate the shape of the sea bottom, a GNSS system receiver for measurements of the coastline, platforms and bitts, and a laser scanner for the measurement of a geospatial point cloud of the onshore infrastructure.
Hydrographic surveys were carried out using an unmanned hydrographic vessel which enabled the performance of surveys in both manual and autonomous mode. The study used two operating modes: autonomous in the open space of the marina along parallel profiles, and manual in areas difficult to access. The hydrographic equipment of the vessel included a SonarMite hydrographic singlebeam echosounder operating at a frequency of 200 kHz, and a Trimble R10 geodetic receiver operating in the VRS Net at a frequency of 10 Hz. On shallow water (a depth of up to 6 m) and on quiet water, an MRU (Motion Reference Unit) sensor was not used for the HPR (heave, pitch, roll) compensation.
Geodetic surveys of quays, platforms and bitts were carried out using a Topcon Hyper II receiver operating in the TPI NetPro network with NET RTCM 3.0 which ensures positioning with a horizontal accuracy better than 3 cm [9,38]. They were conducted in the PL-2000 system with the EVRF2007 geoid model determined with 0.0001 m vertical resolution.
The surveys of port infrastructure were carried out using a Trimble TX8 laser scanner without the photo-taking option. Hence, the obtained point clouds only had colours resulting from the calculated laser beam reflection intensity. Basic parameters of Trimble TX8 laser scanner are presented in Table 1.

Materials
The project was carried out using measuring instruments and systems in three areas: a hydrographic system to investigate the shape of the sea bottom, a GNSS system receiver for measurements of the coastline, platforms and bitts, and a laser scanner for the measurement of a geospatial point cloud of the onshore infrastructure.
Hydrographic surveys were carried out using an unmanned hydrographic vessel which enabled the performance of surveys in both manual and autonomous mode. The study used two operating modes: autonomous in the open space of the marina along parallel profiles, and manual in areas difficult to access. The hydrographic equipment of the vessel included a SonarMite hydrographic singlebeam echosounder operating at a frequency of 200 kHz, and a Trimble R10 geodetic receiver operating in the VRS Net at a frequency of 10 Hz. On shallow water (a depth of up to 6 m) and on quiet water, an MRU (Motion Reference Unit) sensor was not used for the HPR (heave, pitch, roll) compensation.
Geodetic surveys of quays, platforms and bitts were carried out using a Topcon Hyper II receiver operating in the TPI NetPro network with NET RTCM 3.0 which ensures positioning with a horizontal accuracy better than 3 cm [9,38]. They were conducted in the PL-2000 system with the EVRF2007 geoid model determined with 0.0001 m vertical resolution.
The surveys of port infrastructure were carried out using a Trimble TX8 laser scanner without the photo-taking option. Hence, the obtained point clouds only had colours resulting from the calculated laser beam reflection intensity. Basic parameters of Trimble TX8 laser scanner are presented in Table 1.

Hydrography
Where there is enough time to take measurements freely, and where there is no limitation resulting from a contract or a task assigned, hydrometeorological conditions should be used as a guideline. As regards hydrographic surveys, the implementation is significantly affected by waving, water level and acoustic properties of water. Although waving hinders the control of a sounding unit which, in this case was a small unmanned vessel, it also affects the accuracy of determining the seabed shape. AUV rolls cause the deflection of the positioning system antenna to one side, and of the echosounder beam to another when measuring the sloping distance. The second factor, i.e., the water level, is offset in postprocessing based on water gauge observations accessible on an Internet to the nearest station(s). Where a few stations are used, the water level is interpolated based on the position in relation to them. Changes in the speed of sound in the water are long-term (seasonal) and short-term (daily) in nature. The smallest changes in the speed of sound in water are observed in the winter and, from a hydrographic perspective, this is the best period.
When carrying out a sounding using a singlebeam echosounder, the profiles are planned in a parallel system, perpendicular to the contour lines. On open water, when bathymetry (contour lines system) is known, the location of profiles is obvious. Regulations [4,5], which describe in detail the distances between profiles and their lengths, are applied in port basins, at the time of inspections of quays and bathymetric surveys in the areas where vessels are stationed. When no risk to the stability of a marine structure has been identified, or where no such obligation has been found to have been imposed by a specialised construction supervision authority, parallel profiles are positioned perpendicular to the quay at a distance of 10 m over a length of three depths, however no less than 20 m. In addition, a survey needs to be carried out in the immediate vicinity of a hydraulic structure, at a distance of 1 m away from it, another at a distance of 2 m from the previous one, and the next ones every 5 m. This is achieved by taking measurements on profiles parallel to the quay.
On this basis, it was undertaken to plan profiles in an area unoccupied by yachts by using the photographs taken by a UAV (Unmanned Aerial Vehicle) at the time of reconnaissance. Given the movement of vessels in the marina, the dynamics of the situation need to be taken into consideration. Each photograph is only up-to-date at the time when it is taken. Therefore, the best period for carrying out surveys is outside of the sailing season, when a proportion of the yachts have not yet been launched. Figure 2 presents an example of a sounding profile plan for the automatic mode.

Hydrography
Where there is enough time to take measurements freely, and where there is no limitation resulting from a contract or a task assigned, hydrometeorological conditions should be used as a guideline. As regards hydrographic surveys, the implementation is significantly affected by waving, water level and acoustic properties of water. Although waving hinders the control of a sounding unit which, in this case was a small unmanned vessel, it also affects the accuracy of determining the seabed shape. AUV rolls cause the deflection of the positioning system antenna to one side, and of the echosounder beam to another when measuring the sloping distance. The second factor, i.e., the water level, is offset in postprocessing based on water gauge observations accessible on an Internet to the nearest station(s). Where a few stations are used, the water level is interpolated based on the position in relation to them. Changes in the speed of sound in the water are long-term (seasonal) and shortterm (daily) in nature. The smallest changes in the speed of sound in water are observed in the winter and, from a hydrographic perspective, this is the best period.
When carrying out a sounding using a singlebeam echosounder, the profiles are planned in a parallel system, perpendicular to the contour lines. On open water, when bathymetry (contour lines system) is known, the location of profiles is obvious. Regulations [4,5], which describe in detail the distances between profiles and their lengths, are applied in port basins, at the time of inspections of quays and bathymetric surveys in the areas where vessels are stationed. When no risk to the stability of a marine structure has been identified, or where no such obligation has been found to have been imposed by a specialised construction supervision authority, parallel profiles are positioned perpendicular to the quay at a distance of 10 m over a length of three depths, however no less than 20 m. In addition, a survey needs to be carried out in the immediate vicinity of a hydraulic structure, at a distance of 1 m away from it, another at a distance of 2 m from the previous one, and the next ones every 5 m. This is achieved by taking measurements on profiles parallel to the quay.
On this basis, it was undertaken to plan profiles in an area unoccupied by yachts by using the photographs taken by a UAV (Unmanned Aerial Vehicle) at the time of reconnaissance. Given the movement of vessels in the marina, the dynamics of the situation need to be taken into consideration. Each photograph is only up-to-date at the time when it is taken. Therefore, the best period for carrying out surveys is outside of the sailing season, when a proportion of the yachts have not yet been launched. Figure 2 presents an example of a sounding profile plan for the automatic mode. Bathymetric surveys aimed at obtaining geospatial data for the seabed within the yacht port area were carried out twice. In the first stage, in April 2019, basic surveys were carried out in the port, in view of the beginning of the sailing season, as many positions were not occupied by vessels. They were carried out in the port due to unfavourable navigation and hydrographic conditions outside the internal port zone. The second phase of surveys covering the external zone of the NSC (National Sailing Centre) yacht port and the entire external zone of the adjacent objects from the southern Bathymetric surveys aimed at obtaining geospatial data for the seabed within the yacht port area were carried out twice. In the first stage, in April 2019, basic surveys were carried out in the port, in view of the beginning of the sailing season, as many positions were not occupied by vessels. They were carried out in the port due to unfavourable navigation and hydrographic conditions outside the internal port zone. The second phase of surveys covering the external zone of the NSC (National Sailing Centre) yacht port and the entire external zone of the adjacent objects from the southern direction were carried out in the summer. The depth measurements, with account taken of the transducer draught, were related to the MSL (Mean Sea Level). Figure 3 shows a sounding vessel and its trajectory in the first stage of surveys.
Geosciences 2020, 10, x FOR PEER 5 of 14 direction were carried out in the summer. The depth measurements, with account taken of the transducer draught, were related to the MSL (Mean Sea Level). Figure 3 shows a sounding vessel and its trajectory in the first stage of surveys. The singlebeam echosounder was calibrated using the bar check method, which involves adjusting the average speed of sound in water at known depths to echosounder indications. The basic method for measuring the speed of sound is the use of a CTD/STD meter or an SVP probe. An alternative is a measurement of water temperature and the determination of the speed of sound based on an empirical equation when the salinity is known. Some echosounders determine the speed of sound in water based on the temperature by means of a sensor located in the transducer.

GNSS Surveying
As part of the measurement campaign, measurements of 680 points were taken: 160-quays, 310platforms and berth moorings, and 310-bitts and electricity intake points. A Topcon HyPer II receiver with a Sokkia SHC-25 controller operating in the TPI NETpro network with the NET RTCM 3.0 with a virtual reference station were used. This ensured obtaining the accuracy of the determination of position coordinates: horizontal of 0.5-1.8 cm, and vertical of 0.6-2.1 cm.
Overall, 160 characteristic points of the permanent quay were used for comparison with the Electronic Navigational Chart (Figure 4) contained in the PL5WISLA.000 cell along with six (the latest in January of 2018) updates taking into account the port reconstruction and bathymetry. In general, geodetic surveys of the main structure of platforms (red colour) coincide with the ENC imaging. A small difference (approx. 13 cm) was found on the components marked in Figure 5 in grey colour and presented in Figure 6 (on the left). Out of 265 points, 208 supplement the 57 points The singlebeam echosounder was calibrated using the bar check method, which involves adjusting the average speed of sound in water at known depths to echosounder indications. The basic method for measuring the speed of sound is the use of a CTD/STD meter or an SVP probe. An alternative is a measurement of water temperature and the determination of the speed of sound based on an empirical equation when the salinity is known. Some echosounders determine the speed of sound in water based on the temperature by means of a sensor located in the transducer.

GNSS Surveying
As part of the measurement campaign, measurements of 680 points were taken: 160-quays, 310-platforms and berth moorings, and 310-bitts and electricity intake points. A Topcon HyPer II receiver with a Sokkia SHC-25 controller operating in the TPI NETpro network with the NET RTCM 3.0 with a virtual reference station were used. This ensured obtaining the accuracy of the determination of position coordinates: horizontal of 0.5-1.8 cm, and vertical of 0.6-2.1 cm.
Overall, 160 characteristic points of the permanent quay were used for comparison with the Electronic Navigational Chart (Figure 4) contained in the PL5WISLA.000 cell along with six (the latest in January of 2018) updates taking into account the port reconstruction and bathymetry.
Geosciences 2020, 10, x FOR PEER 5 of 14 direction were carried out in the summer. The depth measurements, with account taken of the transducer draught, were related to the MSL (Mean Sea Level). Figure 3 shows a sounding vessel and its trajectory in the first stage of surveys. The singlebeam echosounder was calibrated using the bar check method, which involves adjusting the average speed of sound in water at known depths to echosounder indications. The basic method for measuring the speed of sound is the use of a CTD/STD meter or an SVP probe. An alternative is a measurement of water temperature and the determination of the speed of sound based on an empirical equation when the salinity is known. Some echosounders determine the speed of sound in water based on the temperature by means of a sensor located in the transducer.

GNSS Surveying
As part of the measurement campaign, measurements of 680 points were taken: 160-quays, 310platforms and berth moorings, and 310-bitts and electricity intake points. A Topcon HyPer II receiver with a Sokkia SHC-25 controller operating in the TPI NETpro network with the NET RTCM 3.0 with a virtual reference station were used. This ensured obtaining the accuracy of the determination of position coordinates: horizontal of 0.5-1.8 cm, and vertical of 0.6-2.1 cm.
Overall, 160 characteristic points of the permanent quay were used for comparison with the Electronic Navigational Chart (Figure 4) contained in the PL5WISLA.000 cell along with six (the latest in January of 2018) updates taking into account the port reconstruction and bathymetry. In general, geodetic surveys of the main structure of platforms (red colour) coincide with the ENC imaging. A small difference (approx. 13 cm) was found on the components marked in Figure 5 in grey colour and presented in Figure 6 (on the left). Out of 265 points, 208 supplement the 57 points In general, geodetic surveys of the main structure of platforms (red colour) coincide with the ENC imaging. A small difference (approx. 13 cm) was found on the components marked in Figure 5 in grey colour and presented in Figure 6 (on the left). Out of 265 points, 208 supplement the 57 points that describe the mooring equipment contained in the ENC under the acronym of MORFAC (Mooring/Warping facility). An ENC cell contains only bitts located on the permanent quay.
Geosciences 2020, 10, x FOR PEER 6 of 14 that describe the mooring equipment contained in the ENC under the acronym of MORFAC (Mooring/Warping facility). An ENC cell contains only bitts located on the permanent quay. Figure 6 shows a view of the platform along with bitts (on the left), a bitt with a water take-off point (in the centre), and an electricity intake point (on the right).

Laser Scanning
Laser scanning, frequently used in open spaces, is also determined by atmospheric conditions. In particular, the devices are adversely affected by negative temperatures. Situations in which the difference in ambient temperature and the temperature inside the device are significantly different should be avoided. Specific values are determined by the particular scanner model. However, solar rays may cause interference with the receiver in terms of the returning laser beam, and thus errors when taking a measurement of distance. When taking measurements, one should also avoid rainfall as rain drops can act as lenses, which may cause laser beam deflection. Dust and water particles found in the air increase the risk of beam dispersion, and thus the emitted signal strength is reduced [39].
The three-dimensional presentation of infrastructural components of the marina of the National Sailing Centre at the Gdańsk University of Physical Education and Sport was a significant supplement to bathymetric data. This task was performed using terrestrial laser scanning technology. The instrument that was used (a Trimble TX8 scanner) is equipped with a module that enables taking Geosciences 2020, 10, x FOR PEER 6 of 14 that describe the mooring equipment contained in the ENC under the acronym of MORFAC (Mooring/Warping facility). An ENC cell contains only bitts located on the permanent quay. Figure 6 shows a view of the platform along with bitts (on the left), a bitt with a water take-off point (in the centre), and an electricity intake point (on the right).

Laser Scanning
Laser scanning, frequently used in open spaces, is also determined by atmospheric conditions. In particular, the devices are adversely affected by negative temperatures. Situations in which the difference in ambient temperature and the temperature inside the device are significantly different should be avoided. Specific values are determined by the particular scanner model. However, solar rays may cause interference with the receiver in terms of the returning laser beam, and thus errors when taking a measurement of distance. When taking measurements, one should also avoid rainfall as rain drops can act as lenses, which may cause laser beam deflection. Dust and water particles found in the air increase the risk of beam dispersion, and thus the emitted signal strength is reduced [39].
The three-dimensional presentation of infrastructural components of the marina of the National Sailing Centre at the Gdańsk University of Physical Education and Sport was a significant supplement to bathymetric data. This task was performed using terrestrial laser scanning technology. The instrument that was used (a Trimble TX8 scanner) is equipped with a module that enables taking  Figure 6 shows a view of the platform along with bitts (on the left), a bitt with a water take-off point (in the centre), and an electricity intake point (on the right).

Laser Scanning
Laser scanning, frequently used in open spaces, is also determined by atmospheric conditions. In particular, the devices are adversely affected by negative temperatures. Situations in which the difference in ambient temperature and the temperature inside the device are significantly different should be avoided. Specific values are determined by the particular scanner model. However, solar rays may cause interference with the receiver in terms of the returning laser beam, and thus errors when taking a measurement of distance. When taking measurements, one should also avoid rainfall as rain drops can act as lenses, which may cause laser beam deflection. Dust and water particles found in the air increase the risk of beam dispersion, and thus the emitted signal strength is reduced [39].
The three-dimensional presentation of infrastructural components of the marina of the National Sailing Centre at the Gdańsk University of Physical Education and Sport was a significant supplement to bathymetric data. This task was performed using terrestrial laser scanning technology. The instrument that was used (a Trimble TX8 scanner) is equipped with a module that enables taking photographs of the surroundings and then assigning a colour to particular surveyed points. Within the marine area, 16 scanner measurement stands were established, at which a distance between the target and the device as well as horizontal and vertical angles were measured for each of the cloud points.
Based on the above geodetic observations, the coordinates of particular points in the local coordinate system of the stand were calculated. It is worth noting that each measuring stand had a separate coordinate system with a starting point in the emission-receiving centre of the laser scanner. The angular orientation of the coordinate system resulted from the way (direction) of mounting the instrument on a tripod. In order to create a uniform spatial set for the entire marina in one (global) coordinate system, transformation parameters for particular point clouds with local coordinated systems were determined. The effect of the registration (matching) process is a geometrically and logically consistent set of spatial data that present the surveyed object. In the registration, the matching points were the coordinates of spherical centres of tags (styrofoam balls) placed prior to the laser measurement at locations ensuring the best possible visibility from the neighbouring measurement stands. The registration error for 16 (part of them presented in Figure 7) point clouds did not exceed 1 cm. The resulting set is a digital representation of the marina (Figure 8) that can be displayed, viewed and processed in software dedicated (e.g., CloudCompare) to this purpose.
Geosciences 2020, 10, x FOR PEER 7 of 14 photographs of the surroundings and then assigning a colour to particular surveyed points. Within the marine area, 16 scanner measurement stands were established, at which a distance between the target and the device as well as horizontal and vertical angles were measured for each of the cloud points. Based on the above geodetic observations, the coordinates of particular points in the local coordinate system of the stand were calculated. It is worth noting that each measuring stand had a separate coordinate system with a starting point in the emission-receiving centre of the laser scanner. The angular orientation of the coordinate system resulted from the way (direction) of mounting the instrument on a tripod. In order to create a uniform spatial set for the entire marina in one (global) coordinate system, transformation parameters for particular point clouds with local coordinated systems were determined. The effect of the registration (matching) process is a geometrically and logically consistent set of spatial data that present the surveyed object. In the registration, the matching points were the coordinates of spherical centres of tags (styrofoam balls) placed prior to the laser measurement at locations ensuring the best possible visibility from the neighbouring measurement stands. The registration error for 16 (part of them presented in Figure 7) point clouds did not exceed 1 cm. The resulting set is a digital representation of the marina (Figure 8) that can be displayed, viewed and processed in software dedicated (e.g., CloudCompare) to this purpose.

Dynamic Perspective Projection
The perspective projection belongs to the class of geometrical planar projections. The projection is made on the plane, however, linear projection beams are used. Visual effect of the single-point photographs of the surroundings and then assigning a colour to particular surveyed points. Within the marine area, 16 scanner measurement stands were established, at which a distance between the target and the device as well as horizontal and vertical angles were measured for each of the cloud points. Based on the above geodetic observations, the coordinates of particular points in the local coordinate system of the stand were calculated. It is worth noting that each measuring stand had a separate coordinate system with a starting point in the emission-receiving centre of the laser scanner. The angular orientation of the coordinate system resulted from the way (direction) of mounting the instrument on a tripod. In order to create a uniform spatial set for the entire marina in one (global) coordinate system, transformation parameters for particular point clouds with local coordinated systems were determined. The effect of the registration (matching) process is a geometrically and logically consistent set of spatial data that present the surveyed object. In the registration, the matching points were the coordinates of spherical centres of tags (styrofoam balls) placed prior to the laser measurement at locations ensuring the best possible visibility from the neighbouring measurement stands. The registration error for 16 (part of them presented in Figure 7) point clouds did not exceed 1 cm. The resulting set is a digital representation of the marina (Figure 8) that can be displayed, viewed and processed in software dedicated (e.g., CloudCompare) to this purpose.

Dynamic Perspective Projection
The perspective projection belongs to the class of geometrical planar projections. The projection is made on the plane, however, linear projection beams are used. Visual effect of the single-point

Dynamic Perspective Projection
The perspective projection belongs to the class of geometrical planar projections. The projection is made on the plane, however, linear projection beams are used. Visual effect of the single-point perspective projection is similar to the effect occurring in the human optic system. In the perspective projection the real image of an observed 3D object originates as a result of projection of each point by projection beams radiated from the projection centre and crossing the projection plane.
Assuming that a chart should be prepared for observer's position (for vessel's position), it is suggested that the projection plane Π and reference system Π are given with respect to ' in the following way [33][34][35][36]: The reference system = {O, (e 1 , e 2 , e 3 )} (called the geocentrical) is connected with the Earth in this way, that its versors e 1 , e 2 , e 3 determine the axes X, Y, Z of this system. The reference system ' = {O', (e 1 ', e 2 ', e 3 ')} (called the observer's) is connected with the position of the vessel and its beginning and the base is determined with respect to the reference system . It follows that the projection is made in the observer reference system '. However, before projection, the base (e 1 ', e 2 ', e 3 ') transformation and the centre O of the geocentrical reference system in relation to the base (e 1 ', e 2 ', e 3 ') and the centre O' of the observer reference system ', should be made.
The transformation, which consists of translation of the point O to the point O' and rotation of the base (e 1 , e 2 , e 3 ) with respect to the base (e 1 ', e 2 ', e 3 '), is done using the following transformation matrix [33][34][35][36]: where e 1 = x 1 , y 1 , z 1 , e 2 = x 2 , y 2 , z 2 , e 3 = x 3 , y 3 , z 3 , O = (x 0 , y 0 , z 0 ) When location of vectors e 1 ', e 2 ', e 3 ' with respect to the geocentrical reference system of the observer reference system ' in the location space E 3 are functions of the time variable t (t∈τ), then the system ' is in movement • with respect to the reference system in the time range τ.
There is assumed, that the system ' is in rigid (Euclidean) movement with respect to the system in the time range τ. • unambiguously describes movement of the location space E 3 with respect to the reference system [33][34][35][36]: for t∈τ. This method was used for dynamic visualisation of the surroundings built on the basis of the geospatial data included in the ENC [40,41]. In this model it was assumed that: the land surface is flat, all buildings are of equal height, and the sea bottom model is built on the basis of objects including the geospatial data: LNDARE (land area), DEPCNT (depth contour) and SOUNDG (sounding) [40,41]. Spatial presentation of the NSC area is shown in Figure 9.

Hydrographic and Geodetic Measurements
A bathymetric chart, as one of the reporting documents from bathymetric surveys, was prepared in accordance with the rules adopted in hydrography and presented in Figure 10. The basic components include: scale, frame, coastline, contour lines and (depth) soundings.
The map scale applied in most bathymetric studies is usually 1:500 or 1:1000. Where surveys are carried out in accordance with [4,5], the map scale is 1:100, 1:500 or 1:1000. It is selected to the printout format, in order to maintain its cartometric properties. The frame contains coordinates in which the study was carried out. This is usually in the UTM system, although this study was prepared in the PL-2000 system because of integration of the hydrographic and geodetic data for land management. The grid of crosses is enriched with lines of the grid of geographical coordinates.

Hydrographic and Geodetic Measurements
A bathymetric chart, as one of the reporting documents from bathymetric surveys, was prepared in accordance with the rules adopted in hydrography and presented in Figure 10. The basic components include: scale, frame, coastline, contour lines and (depth) soundings.

Hydrographic and Geodetic Measurements
A bathymetric chart, as one of the reporting documents from bathymetric surveys, was prepared in accordance with the rules adopted in hydrography and presented in Figure 10. The basic components include: scale, frame, coastline, contour lines and (depth) soundings.
The map scale applied in most bathymetric studies is usually 1:500 or 1:1000. Where surveys are carried out in accordance with [4,5], the map scale is 1:100, 1:500 or 1:1000. It is selected to the printout format, in order to maintain its cartometric properties. The frame contains coordinates in which the study was carried out. This is usually in the UTM system, although this study was prepared in the PL-2000 system because of integration of the hydrographic and geodetic data for land management. The grid of crosses is enriched with lines of the grid of geographical coordinates.  The map scale applied in most bathymetric studies is usually 1:500 or 1:1000. Where surveys are carried out in accordance with [4,5], the map scale is 1:100, 1:500 or 1:1000. It is selected to the printout format, in order to maintain its cartometric properties. The frame contains coordinates in which the study was carried out. This is usually in the UTM system, although this study was prepared in the PL-2000 system because of integration of the hydrographic and geodetic data for land management. The grid of crosses is enriched with lines of the grid of geographical coordinates.
Among the methods for creating the Digital SeaBottom Model [42][43][44][45], the IDW (Inverse Distance Weighted) model with separation every 0.5 m was applied to determine the contour lines. The size of the uniform (rectangular) grid cell was set at 0.5 m × 0.5 m. A matrix of 689 × 781 cells covering the NSC area and the adjacent water regions was created in this way ( Figure 10). DSBM-Digital SeaBottom Model can also be used for its spatial presentation. The (depth) soundings were plotted with a resolution of 0.1 m in accordance with [46] for a shallow water region with depths of up to 31 m. The study was carried out in the ArcGIS 10.5 environment using the ENC PL5WISLA cell.

Laser Scaning
The three-dimensional presentation of infrastructural components of the marina of the National Sailing Centre at the Gdańsk University of Physical Education and Sport was a significant supplement to bathymetric data. This task was performed using the terrestrial laser scanning technology. An image of infrastructure and the platform from TLS cloud are presented in Figure 11. Among the methods for creating the Digital SeaBottom Model [42][43][44][45], the IDW (Inverse Distance Weighted) model with separation every 0.5 m was applied to determine the contour lines. The size of the uniform (rectangular) grid cell was set at 0.5 m  0.5 m. A matrix of 689  781 cells covering the NSC area and the adjacent water regions was created in this way ( Figure 10). DSBM -Digital SeaBottom Model can also be used for its spatial presentation. The (depth) soundings were plotted with a resolution of 0.1 m in accordance with [46] for a shallow water region with depths of up to 31 m. The study was carried out in the ArcGIS 10.5 environment using the ENC PL5WISLA cell.

Laser Scaning
The three-dimensional presentation of infrastructural components of the marina of the National Sailing Centre at the Gdańsk University of Physical Education and Sport was a significant supplement to bathymetric data. This task was performed using the terrestrial laser scanning technology. An image of infrastructure and the platform from TLS cloud are presented in Figure 11. Where work is carried out on spatial information sets originating from various surveying instruments, it is recommended that data harmonisation (understood as the adjustment of coordinate systems and reference systems) should be carried out to reliably conclude and compare both elements [47]. In the case under consideration, it could be assumed that both spatial references are identical.
Using dynamic perspective projection on the basis of following parameters: vessel's position, course, height of the observer, and inclination angle, the image seen by the navigator and that created on the basis of geospatial data from TLS cloud are similar ( Figure 12).  Where work is carried out on spatial information sets originating from various surveying instruments, it is recommended that data harmonisation (understood as the adjustment of coordinate systems and reference systems) should be carried out to reliably conclude and compare both elements [47]. In the case under consideration, it could be assumed that both spatial references are identical.
Using dynamic perspective projection on the basis of following parameters: vessel's position, course, height of the observer, and inclination angle, the image seen by the navigator and that created on the basis of geospatial data from TLS cloud are similar ( Figure 12). Among the methods for creating the Digital SeaBottom Model [42][43][44][45], the IDW (Inverse Distance Weighted) model with separation every 0.5 m was applied to determine the contour lines. The size of the uniform (rectangular) grid cell was set at 0.5 m  0.5 m. A matrix of 689  781 cells covering the NSC area and the adjacent water regions was created in this way ( Figure 10). DSBM -Digital SeaBottom Model can also be used for its spatial presentation. The (depth) soundings were plotted with a resolution of 0.1 m in accordance with [46] for a shallow water region with depths of up to 31 m. The study was carried out in the ArcGIS 10.5 environment using the ENC PL5WISLA cell.

Laser Scaning
The three-dimensional presentation of infrastructural components of the marina of the National Sailing Centre at the Gdańsk University of Physical Education and Sport was a significant supplement to bathymetric data. This task was performed using the terrestrial laser scanning technology. An image of infrastructure and the platform from TLS cloud are presented in Figure 11. Where work is carried out on spatial information sets originating from various surveying instruments, it is recommended that data harmonisation (understood as the adjustment of coordinate systems and reference systems) should be carried out to reliably conclude and compare both elements [47]. In the case under consideration, it could be assumed that both spatial references are identical.
Using dynamic perspective projection on the basis of following parameters: vessel's position, course, height of the observer, and inclination angle, the image seen by the navigator and that created on the basis of geospatial data from TLS cloud are similar ( Figure 12).  Synthetic presentation of the geospatial data from bathymetric sounding and laser scanning in perspective projection is shown in Figure 13. In addition to static objects, such as berths, piers and buildings, it includes mobile objects, such as: yachts and boats, cranes and cars, captured during laser scanning.
Synthetic presentation of the geospatial data from bathymetric sounding and laser scanning in perspective projection is shown in Figure 13. In addition to static objects, such as berths, piers and buildings, it includes mobile objects, such as: yachts and boats, cranes and cars, captured during laser scanning. Figure 13. Synthesis of geospatial data from bathymetric sounding and laser scanning.

Conclusions
Bathymetric surveys of ports and yacht marinas are extremely difficult due to the limited maneuverability between berths and moored vessels. Autonomous measuring units that carry measuring equipment, i.e., a positioning system receiver and a depth measuring device, are perfect for this task. Precise positioning in dynamics is already possible thanks to the use of geodetic system satellite receivers that enable positioning in real time at a high frequency. A depth measurement is usually taken using a singlebeam echosounder due to its small size and weight. This helps the user plan a system of surveying profiles while taking account of the estimated seabed shape and the detail level of the sounding.
Geospatial data of the upper hemisphere which present an image of quays, platforms and harbour infrastructure provide more detail to the image of the yacht harbour. They can be obtained by laser scanning or aerial photogrammetry which can currently be done with a UAV (Unmanned Aerial Vehicle). The first method is more labour-intensive and time-consuming due to the need for carrying out surveys from multiple stands.
The article presents the synthesis of bathymetric surveys carried out using a singlebeam echosounder mounted on an autonomous vessel, with measurements taken with a laser scanner. Thanks to the use of a highly accurate hydrographic echosounder and a GNSS receiver with RTK corrections, echosounder calibration, and the consideration of changes in water level, an autonomous surveying unit enabled the determination of the seabed topography with an accuracy of single centimeters. The information is supplemented by a point cloud obtained by the laser scanning method with a millimeter accuracy of the coordinates being determined. Thanks to harmonisation, the coordinates were presented in a uniform rectangular coordinate system and the depths were referred to the average sea (North Sea) level determined for the automatic tide gauge in Amsterdam. Integrated GNSS, laser and hydrographic data can be used in visualization systems, such as ECS, ECDIS and navigational simulator for dynamic, spatial presentation. On the basis of the observer's position, the image generated and presented on the display is similar to the image seen by the navigator.

Conclusions
Bathymetric surveys of ports and yacht marinas are extremely difficult due to the limited maneuverability between berths and moored vessels. Autonomous measuring units that carry measuring equipment, i.e., a positioning system receiver and a depth measuring device, are perfect for this task. Precise positioning in dynamics is already possible thanks to the use of geodetic system satellite receivers that enable positioning in real time at a high frequency. A depth measurement is usually taken using a singlebeam echosounder due to its small size and weight. This helps the user plan a system of surveying profiles while taking account of the estimated seabed shape and the detail level of the sounding.
Geospatial data of the upper hemisphere which present an image of quays, platforms and harbour infrastructure provide more detail to the image of the yacht harbour. They can be obtained by laser scanning or aerial photogrammetry which can currently be done with a UAV (Unmanned Aerial Vehicle). The first method is more labour-intensive and time-consuming due to the need for carrying out surveys from multiple stands.
The article presents the synthesis of bathymetric surveys carried out using a singlebeam echosounder mounted on an autonomous vessel, with measurements taken with a laser scanner. Thanks to the use of a highly accurate hydrographic echosounder and a GNSS receiver with RTK corrections, echosounder calibration, and the consideration of changes in water level, an autonomous surveying unit enabled the determination of the seabed topography with an accuracy of single centimeters. The information is supplemented by a point cloud obtained by the laser scanning method with a millimeter accuracy of the coordinates being determined. Thanks to harmonisation, the coordinates were presented in a uniform rectangular coordinate system and the depths were referred to the average sea (North Sea) level determined for the automatic tide gauge in Amsterdam. Integrated GNSS, laser and hydrographic data can be used in visualization systems, such as ECS, ECDIS and navigational simulator for dynamic, spatial presentation. On the basis of the observer's position, the image generated and presented on the display is similar to the image seen by the navigator.