We evaluated the performance of nine machine learning regression algorithms and their ensembles for sub-pixel estimation of impervious areas coverages from Landsat imagery. The accuracy of imperviousness mapping in individual time points was assessed based on RMSE, MAE and R 2 . These measures were also used for the assessment of imperviousness change intensity estimations. The applicability for detection of relevant changes in impervious areas coverages at sub-pixel level was evaluated using overall accuracy, F-measure and ROC Area Under Curve. The results proved that Cubist algorithm may be advised for Landsat-based mapping of imperviousness for single dates. Stochastic gradient boosting of regression trees (GBM) may be also considered for this purpose. However, Random Forest algorithm is endorsed for both imperviousness change detection and mapping of its intensity. In all applications the heterogeneous model ensembles performed at least as well as the best individual models or better. They may be recommended for improving the quality of sub-pixel imperviousness and imperviousness change mapping. The study revealed also limitations of the investigated methodology for detection of subtle changes of imperviousness inside the pixel. None of the tested approaches was able to reliably classify changed and non-changed pixels if the relevant change threshold was set as one or three percent. Also for five percent change threshold most of algorithms did not ensure that the accuracy of change map is higher than the accuracy of random classifier. For the threshold of relevant change set as ten percent all approaches performed satisfactory.
One of the fundamental problems of modern geodesy is precise de fi nition of the gravitational fi eld and its changes in time. This is essential in positioning and navigation, geophysics, geodynamics, oceanography and other sciences related to the climate and Earth’s environment. One of the major sources of gravity data is satellite altimetry that provides gravity data with almost 75% surface of the Earth. Satellite altimetry also provides data to study local, regional and global geophysical processes, the geoid model in the areas of oceans and seas. This technique can be successfully used to study the ocean mean dynamic topography. The results of the investigations and possible products of altimetry will provide a good material for the GGOS (Global Geodetic Observing System) and institutions of IAS (International Altimetry Service). This paper presents the achievements in satellite altimetry in all the above disciplines obtained in the last years. First very shorly basic concept of satellite altimetry is given. In order to obtain the highest accuracy on range measurements over the ocean improved of altimetry waveforms performed on the ground is described. Next, signi fi cant improvements of sea and ocean gravity anomalies models developed presently is shown. Study of sea level and its extremes examined, around European and Australian coasts using tide gauges data and satellite altimetry measurements were described. Then investigations of the phenomenon of the ocean tides, calibration of altimeters, studies of rivers and ice-sheets in the last years are given.
Land consolidation procedures are an attempt to comprehensively change the existing spatial structure of land in rural areas. This treatment also brings many other social and economic benefits, contributing to the development of consolidated areas. Land consolidation in mountain areas differs in many respects from those implemented in areas with more favorable conditions for the functioning of agriculture. The unfavorable values of land fragmentation indices, terrain conditions and lower than the average soil quality affect both the dominant forms of agricultural activity and the limited opportunities to improve the distribution of plots in space, parameters of shape, and the area as a result of land consolidation. For this reason, the effectiveness of land consolidation in mountain areas can be achieved by improving the quality of transportation network and the accessibility of the plots, arranging ownership issues and improving the quality of cadastral documentation. This article presents the evaluation of the measures of effectiveness of land consolidation realized in mountain areas on the example of Łetownia Village in the Małopolska Province, located in the southern part of Poland. Selected village is an area with unfavorable conditions for the functioning of agriculture and high values of land fragmentation indices.
The cognitive aim of this study is to point to the optimum number of local government units and the optimum boundaries of spatial units in Poland with the assumption of minimizing the cumulated theoretical travel time to all settlement units in the country. The methodological aim, in turn, is to present the use of the ArcGIS location-allocation tool for the purposes of delimitation processes as exemplified by administrative boundaries in Poland. The rationale for the implementation of this study is that number and the boundaries of units of all levels of Poland’s current territorial division are far from optimum in the light of minimization of accumulated theoretical travel time to all settlement units in the country. It may be concluded that it would be justifiable to increase the number of voivodships from the current number of 16 to 18. Besides it would be necessary to introduce modifications in relation to units with regional functions. In contrast, the number of districts and communes should be reduced. A continuation of this research may go in the direction of including analysis of public transport network in the research, creating in this way a multimodal set of network data. This would illustrate, apart from the potential itself resulting from the infrastructure, also the actually existing connections.
The cadastral data, including land parcels, are the basic reference data for presenting various objects collected in spatial databases. Easy access to up-to-date records is a very important matter for the individuals and institutions using spatial data infrastructure. The primary objective of the study was to check the current accessibility of cadastral data as well as to verify how current and complete they are. The author started researching this topic in 2007, i.e. from the moment the Team for National Spatial Data Infrastructure developed documentation concerning the standard of publishing cadastral data with the use of the WMS. Since ten years, the author was monitoring the status of cadastral data publishing in various districts as well as participated in data publishing in many districts. In 2017, when only half of the districts published WMS services from cadastral data, the questions arise: why is it so and how to change this unfavourable status? As a result of the tests performed, it was found that the status of publishing cadastral data is still far from perfect. The quality of the offered web services varies and, unfortunately, many services offer poor performance; moreover, there are plenty services that do not operate at all.
In recent years, Global Navigation Satellite Systems (GNSS) have gained great importance in terms of the benefits it provides such as precise geodetic point positioning, determining crustal deformations, navigation, vehicle monitoring systems and meteorological applications etc. As in Turkey, for this purpose, each country has set up its own GNSS station networks like Turkish National Permanent RTK Network analyzed precise station coordinates and velocities together with the International GNSS Service , Turkish National Fundamental GPS Network and Turkish National Permanent GNSS Network (TNPGN) stations not only are utilized as precise positioning but also GNSS meteorology studies so total number of stations are increased. This work is related to the reactivated of the TRAB IGS station which was established in Karadeniz Technical University, Department of Geomatics Engineering. Within the COST ES1206 Action (GNSS4SWEC) KTU analysis center was established and Trop-NET system developed by Geodetic Observatory Pecny (GOP, RIGTC) in order to troposphere monitoring. The project titled “Using Regional GNSS Networks to Strengthen Severe Weather Prediction” was accepted to the scientific and technological research council of Turkey (TUBITAK). With this project, we will design 2 new constructed GNSS reference station network. Using observation data of network, we will compare water vapor distribution derived by GNSS Meteorology and GNSS Tomography. At this time, KTU AC was accepted as E-GVAP Analysis Centre in December 2016. KTU reference station is aimed to be a member of the EUREF network with these studies.
The dedicated gravity satellite missions, in particular the GRACE (Gravity Recovery and Climate Experiment) mission launched in 2002, provide unique data for studying temporal variations of mass distribution in the Earth’s system, and thereby, the geometry and the gravity field changes of the Earth. The main objective of this contribution is to estimate physical height (e.g. the orthometric/normal height) changes over Central Europe using GRACE satellite mission data as well as to analyse them and model over the selected study area. Physical height changes were estimated from temporal variations of height anomalies and vertical displacements of the Earth surface being determined over the investigated area. The release 5 (RL05) GRACE-based global geopotential models as well as load Love numbers from the Preliminary Reference Earth Model (PREM) were used as input data. Analysis of the estimated physical height changes and their modelling were performed using two methods: the seasonal decomposition method and the PCA/ EOF (Principal Component Analysis/Empirical Orthogonal Function) method and the differences obtained were discussed. The main findings reveal that physical height changes over the selected study area reach up to 22.8 mm. The obtained physical height changes can be modelled with an accuracy of 1.4 mm using the seasonal decomposition method.
MP estimation is a method which concerns estimating of the location parameters when the probabilistic models of observations differ from the normal distributions in the kurtosis or asymmetry. The system of Pearson’s distributions is the probabilistic basis for the method. So far, such a method was applied and analyzed mostly for leptokurtic or mesokurtic distributions (Pearson’s distributions of types IV or VII), which predominate practical cases. The analyses of geodetic or astronomical observations show that we may also deal with sets which have moderate asymmetry or small negative excess kurtosis. Asymmetry might result from the influence of many small systematic errors, which were not eliminated during preprocessing of data. The excess kurtosis can be related with bigger or smaller (in relations to the Hagen hypothesis) frequency of occurrence of the elementary errors which are close to zero. Considering that fact, this paper focuses on the estimation with application of the Pearson platykurtic distributions of types I or II. The paper presents the solution of the corresponding optimization problem and its basic properties. Although platykurtic distributions are rare in practice, it was an interesting issue to find out what results can be provided by MP estimation in the case of such observation distributions. The numerical tests which are presented in the paper are rather limited; however, they allow us to draw some general conclusions.
The objective of the paper was the hydrological analysis, in terms of categorizing main watercourses (based on coupled catchments) and marking areas covered by potential impact of the occurrence and activities of the European beaver Castor fi ber . At the analysed area – the Forest District Głogów Małopolski there is a population of about 200 beavers in that Forest District. Damage inflicted by beavers was detected on 33.0 ha of the Forest District, while in the area of 13.9 ha the damage was small (below 10%). The monitoring of the beavers’ behaviour and the analysis of their influence on hydrology of the area became an important element of using geoinformationtools in the management of forest areas. ArcHydro ArcGIS Esri module was applied, as an integrated set of tools for hydrographical analysis and modelling. Further steps of the procedure are hydrologic analyses such as: marking river networks on the DTM, filling holes, making maps of the flow direction, making the map of the accumulation flow, defining and segmentation of streams, marking elementary basins, marking coupled basins, making dams in the places, where beavers occur and localization of the area with a visible impact of damming. The result of the study includes maps prepared for the Forest District: the map of main rivers and their basins, categories of watercourses and compartments particularly threatened by beaver’s foraging.
The TerraSAR-X add-on for Digital Elevation Measurement ( TanDEM-X) mission launched in 2010 is another programme – after the Shuttle Radar Topography Mission (SRTM) in 2000 – that uses space-borne radar interferometry to build a global digital surface model. This article presents the accuracy assessment of the TanDEM-X intermediate Digital Elevation Model (IDEM) provided by the German Aerospace Center (DLR) under the project “Accuracy assessment of a Digital Elevation Model based on TanDEM-X data” for the southwestern territory of Poland. The study area included: open terrain, urban terrain and forested terrain. Based on a set of 17,498 reference points acquired by airborne laser scanning, the mean errors of average heights and standard deviations were calculated for areas with a terrain slope below 2 degrees, between 2 and 6 degrees and above 6 degrees. The absolute accuracy of the IDEM data for the analysed area, expressed as a root mean square error (Total RMSE), was 0.77 m.
The processing of cartographic data demands human involvement. Up-to-date algorithms try to automate a part of this process. The goal is to obtain a digital model, or additional information about shape and topology of input geometric objects. A topological skeleton is one of the most important tools in the branch of science called shape analysis. It represents topological and geometrical characteristics of input data. Its plot depends on using algorithms such as medial axis, skeletonization, erosion, thinning, area collapse and many others. Area collapse, also known as dimension change, replaces input data with lower-dimensional geometric objects like, for example, a polygon with a polygonal chain, a line segment with a point. The goal of this paper is to introduce a new algorithm for the automatic calculation of polygonal chains representing a 2D polygon. The output is entirely contained within the area of the input polygon, and it has a linear plot without branches. The computational process is automatic and repeatable. The requirements of input data are discussed. The author analyzes results based on the method of computing ends of output polygonal chains. Additional methods to improve results are explored. The algorithm was tested on real-world cartographic data received from BDOT/GESUT databases, and on point clouds from laser scanning. An implementation for computing hatching of embankment is described.
Unmanned aerial vehicles are increasingly being used in close range photogrammetry. Real-time observation of the Earth’s surface and the photogrammetric images obtained are used as material for surveying and environmental inventory. The following study was conducted on a small area (approximately 1 ha). In such cases, the classical method of topographic mapping is not accurate enough. The geodetic method of topographic surveying, on the other hand, is an overly precise measurement technique for the purpose of inventorying the natural environment components. The author of the following study has proposed using the unmanned aerial vehicle technology and tying in the obtained images to the control point network established with the aid of GNSS technology. Georeferencing the acquired images and using them to create a photogrammetric model of the studied area enabled the researcher to perform calculations, which yielded a total root mean square error below 9 cm. The performed comparison of the real lengths of the vectors connecting the control points and their lengths calculated on the basis of the photogrammetric model made it possible to fully confirm the RMSE calculated and prove the usefulness of the UAV technology in observing terrain components for the purpose of environmental inventory. Such environmental components include, among others, elements of road infrastructure, green areas, but also changes in the location of moving pedestrians and vehicles, as well as other changes in the natural environment that are not registered on classical base maps or topographic maps.
Time series of weekly and daily solutions for coordinates of permanent GNSS stations may indicate local deformations in Earth’s crust or local seasonal changes in the atmosphere and hydrosphere. The errors of the determined changes are relatively large, frequently at the level of the signal. Satellite radar interferometry and especially Persistent Scatterer Interferometry (PSI) is a method of a very high accuracy. Its weakness is a relative nature of measurements as well as accumulation of errors which may occur in the case of PSI processing of large areas. It is thus beneficial to confront the results of PSI measurements with those from other techniques, such as GNSS and precise levelling. PSI and GNSS results were jointly processed recreating the history of surface deformation of the area of Warsaw metropolitan with the use of radar images from Envisat and Cosmo- SkyMed satellites. GNSS data from Borowa Gora and Jozefoslaw observatories as well as from WAT1 and CBKA permanent GNSS stations were used to validate the obtained results. Observations from 2000–2015 were processed with the Bernese v.5.0 software. Relative height changes between the GNSS stations were determined from GNSS data and relative height changes between the persistent scatterers located on the objects with GNSS stations were determined from the interferometric results. The consistency of results of the two methods was 3 to 4 times better than the theoretical accuracy of each. The joint use of both methods allows to extract a very small height change below the level of measurement error.
The base map provides basic information about land to individuals, companies, developers, design engineers, organizations, and government agencies. Its contents include spatial location data for control network points, buildings, land lots, infrastructure facilities, and topographic features. As the primary map of the country, it must be developed in accordance with specific laws and regulations and be continuously updated. The base map is a data source used for the development and updating of derivative maps and other large scale cartographic materials such as thematic or topographic maps. Thanks to the advancement of science and technology, the quality of land surveys carried out by means of terrestrial laser scanning (TLS) matches that of traditional surveying methods in many respects. This paper discusses the potential application of output data from laser scanners (point clouds) to the development and updating of cartographic materials, taking Poland’s base map as an example. A few research sites were chosen to present the method and the process of conducting a TLS land survey: a fragment of a residential area, a street, the surroundings of buildings, and an undeveloped area. The entire map that was drawn as a result of the survey was checked by comparing it to a map obtained from PODGiK (pol. Powiatowy Ośrodek Dokumentacji Geodezyjnej i Kartograficznej – Regional Centre for Geodetic and Cartographic Records) and by conducting a field inspection. An accuracy and quality analysis of the conducted fieldwork and deskwork yielded very good results, which provide solid grounds for predicating that cartographic materials based on a TLS point cloud are a reliable source of information about land. The contents of the map that had been created with the use of the obtained point cloud were very accurately located in space (x, y, z). The conducted accuracy analysis and the inspection of the performed works showed that high quality is characteristic of TLS surveys. The accuracy of determining the location of the various map contents has been estimated at 0.02-0.03 m. The map was developed in conformity with the applicable laws and regulations as well as with best practice requirements.
In the following paper, geovisualisation will be applied to one spatial phenomenon and understood as a process of creating complementary visualisations: static two-dimensional, surface three-dimensional, and interactive. The central challenge that the researchers faced was to find a method of presenting the phenomenon in a multi- faceted way. The main objective of the four-stage study was to show the capacity of the contemporary software for presenting geographical space from various perspectives while maintaining the standards of cartographic presentation and making sure that the form remains attractive for the user. The correctness, effectiveness, and usefulness of the proposed approach was analysed on the basis of a geovisualisation of natural aggregate extraction in the Gniezno district in the years 2005–2015. For each of the three visualisations, the researchers planned a different range of information, different forms of graphic and cartographic presentation, different use and function, but as far as possible the same accessible databases and the same free technologies. On the basis of the final publication, the researchers pointed out the advantages of the proposed work flow and the correctness of the detailed flowchart.
A parcel is the most important object of real estate cadastre. Its primary spatial attribute are boundaries, determining the extent of property rights. Capturing the data on boundaries should be performed in the way ensuring sufficiently high accuracy and reliability. In recent years, as part of the project “ZSIN – Construction of Integrated Real Estate Information System – Stage I”, in the territories of the participating districts, actions were taken aimed at the modernization of the register of land and buildings. In many cases, this process was carried out basing on photogrammetric materials. Applicable regulations allow such a possibility. This paper, basing on the documentation from the National Geodetic and Cartographic Documentation Center and on the authors’ own surveys attempts to assess the applicability of the photogrammetric method to capture data on the boundaries of cadastral parcels. The scope of the research, most importantly, included the problem of accuracy with which it was possible to determine the position of a boundary point using photogrammetric surveys carried out on the terrain model created from processed aerial photographs. The article demonstrates the manner of recording this information in the cadastral database, as well as the resulting legal consequences. Moreover, the level of reliability of the entered values of the selected attributes of boundary points was assessed.