<div class="abstract"><div class="abstract_item"><em>Modern techniques for the map analysis allow for the creation of full or partial geometric reconstruction of its content. The projection is described by the set of estimated constant values: transformed pole position, standard parallel latitude, longitude of the central meridian, and a constant parameter. Analogously the analyzed map is represented by its constant values: auxiliary sphere radius, origin shifts, and angle of rotation. Several new methods denoted as M6-M9 for the estimation of an unknown map projection and its parameters differing in the number of determined parameters, reliability, robustness, and convergence have been developed. However, their computational demands are similar. Instead of directly measuring the dissimilarity of two projections, the analyzed map in an unknown projection and the image of the sphere in the well-known (i.e., analyzed) projection are compared. Several distance functions for the similarity measurements based on the location as well as shape similarity approaches are proposed. An unconstrained global optimization problem poorly scaled, with large residuals, for the vector of unknown parameters is solved by the hybrid BFGS method. To avoid a slower convergence rate for small residual problems, it has the ability to switch between first- and second-order methods. Such an analysis is beneficial and interesting for historic, old, or current maps without information about the projection. Its importance is primarily referred to refinement of spatial georeference for the medium- and small-scale maps, analysis of the knowledge about the former world, analysis of the incorrectly/inaccurately drawn regions, and appropriate cataloging of maps. The proposed algorithms have been implemented in the new version of the <span style="font-family: monospace;">detectproj</span> software.</em></div></div>
Plotting the map projection graticule involving discontinuities based on combined sampling
