123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229 |
- // Boost.Geometry
- // This file is manually converted from PROJ4
- // This file was modified by Oracle on 2018, 2019.
- // Modifications copyright (c) 2018-2019, Oracle and/or its affiliates.
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- // Use, modification and distribution is subject to the Boost Software License,
- // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- // This file is converted from PROJ4, http://trac.osgeo.org/proj
- // PROJ4 is originally written by Gerald Evenden (then of the USGS)
- // PROJ4 is maintained by Frank Warmerdam
- // This file was converted to Geometry Library by Adam Wulkiewicz
- // Original copyright notice:
- // Author: Frank Warmerdam, [email protected]
- // Copyright (c) 2000, Frank Warmerdam
- // Permission is hereby granted, free of charge, to any person obtaining a
- // copy of this software and associated documentation files (the "Software"),
- // to deal in the Software without restriction, including without limitation
- // the rights to use, copy, modify, merge, publish, distribute, sublicense,
- // and/or sell copies of the Software, and to permit persons to whom the
- // Software is furnished to do so, subject to the following conditions:
- // The above copyright notice and this permission notice shall be included
- // in all copies or substantial portions of the Software.
- // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
- // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
- // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
- // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
- // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
- // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
- // DEALINGS IN THE SOFTWARE.
- #ifndef BOOST_GEOMETRY_SRS_PROJECTIONS_IMPL_PJ_GRIDLIST_HPP
- #define BOOST_GEOMETRY_SRS_PROJECTIONS_IMPL_PJ_GRIDLIST_HPP
- #include <boost/geometry/srs/projections/exception.hpp>
- #include <boost/geometry/srs/projections/grids.hpp>
- #include <boost/geometry/srs/projections/impl/pj_gridinfo.hpp>
- #include <boost/geometry/srs/projections/impl/pj_strerrno.hpp>
- #include <boost/geometry/srs/projections/par_data.hpp>
- namespace boost { namespace geometry { namespace projections
- {
- namespace detail
- {
- /************************************************************************/
- /* pj_gridlist_merge_grid() */
- /* */
- /* Find/load the named gridfile and merge it into the */
- /* last_nadgrids_list. */
- /************************************************************************/
- // Originally one function, here divided into several functions
- // with overloads for various types of grids and stream policies
- inline bool pj_gridlist_find_all(std::string const& gridname,
- pj_gridinfo const& grids,
- std::vector<std::size_t> & gridindexes)
- {
- bool result = false;
- for (std::size_t i = 0 ; i < grids.size() ; ++i)
- {
- if (grids[i].gridname == gridname)
- {
- result = true;
- gridindexes.push_back(i);
- }
- }
- return result;
- }
- // Fill container with sequentially increasing numbers
- inline void pj_gridlist_add_seq_inc(std::vector<std::size_t> & gridindexes,
- std::size_t first, std::size_t last)
- {
- gridindexes.reserve(gridindexes.size() + (last - first));
- for ( ; first < last ; ++first)
- {
- gridindexes.push_back(first);
- }
- }
- // Generic stream policy and standard grids
- template <typename StreamPolicy, typename Grids>
- inline bool pj_gridlist_merge_gridfile(std::string const& gridname,
- StreamPolicy const& stream_policy,
- Grids & grids,
- std::vector<std::size_t> & gridindexes,
- grids_tag)
- {
- // Try to find in the existing list of loaded grids. Add all
- // matching grids as with NTv2 we can get many grids from one
- // file (one shared gridname).
- if (pj_gridlist_find_all(gridname, grids.gridinfo, gridindexes))
- return true;
- std::size_t orig_size = grids.gridinfo.size();
- // Try to load the named grid.
- typename StreamPolicy::stream_type is;
- stream_policy.open(is, gridname);
- if (! pj_gridinfo_init(gridname, is, grids.gridinfo))
- {
- return false;
- }
- // Add the grid now that it is loaded.
- pj_gridlist_add_seq_inc(gridindexes, orig_size, grids.gridinfo.size());
- return true;
- }
- // Generic stream policy and shared grids
- template <typename StreamPolicy, typename SharedGrids>
- inline bool pj_gridlist_merge_gridfile(std::string const& gridname,
- StreamPolicy const& stream_policy,
- SharedGrids & grids,
- std::vector<std::size_t> & gridindexes,
- shared_grids_tag)
- {
- // Try to find in the existing list of loaded grids. Add all
- // matching grids as with NTv2 we can get many grids from one
- // file (one shared gridname).
- {
- typename SharedGrids::read_locked lck_grids(grids);
- if (pj_gridlist_find_all(gridname, lck_grids.gridinfo, gridindexes))
- return true;
- }
- // Try to load the named grid.
- typename StreamPolicy::stream_type is;
- stream_policy.open(is, gridname);
- pj_gridinfo new_grids;
- if (! pj_gridinfo_init(gridname, is, new_grids))
- {
- return false;
- }
- // Add the grid now that it is loaded.
- std::size_t orig_size = 0;
- std::size_t new_size = 0;
- {
- typename SharedGrids::write_locked lck_grids(grids);
- // Try to find in the existing list of loaded grids again
- // in case other thread already added it.
- if (pj_gridlist_find_all(gridname, lck_grids.gridinfo, gridindexes))
- return true;
- orig_size = lck_grids.gridinfo.size();
- new_size = orig_size + new_grids.size();
- lck_grids.gridinfo.resize(new_size);
- for (std::size_t i = 0 ; i < new_grids.size() ; ++ i)
- new_grids[i].swap(lck_grids.gridinfo[i + orig_size]);
- }
- pj_gridlist_add_seq_inc(gridindexes, orig_size, new_size);
- return true;
- }
- /************************************************************************/
- /* pj_gridlist_from_nadgrids() */
- /* */
- /* This functions loads the list of grids corresponding to a */
- /* particular nadgrids string into a list, and returns it. The */
- /* list is kept around till a request is made with a different */
- /* string in order to cut down on the string parsing cost, and */
- /* the cost of building the list of tables each time. */
- /************************************************************************/
- template <typename StreamPolicy, typename Grids>
- inline void pj_gridlist_from_nadgrids(srs::detail::nadgrids const& nadgrids,
- StreamPolicy const& stream_policy,
- Grids & grids,
- std::vector<std::size_t> & gridindexes)
- {
- // Loop processing names out of nadgrids one at a time.
- for (srs::detail::nadgrids::const_iterator it = nadgrids.begin() ;
- it != nadgrids.end() ; ++it)
- {
- bool required = (*it)[0] != '@';
- std::string name(it->begin() + (required ? 0 : 1), it->end());
- if ( ! pj_gridlist_merge_gridfile(name, stream_policy, grids, gridindexes,
- typename Grids::tag())
- && required )
- {
- BOOST_THROW_EXCEPTION( projection_exception(error_failed_to_load_grid) );
- }
- }
- }
- template <typename Par, typename ProjectionGrids>
- inline void pj_gridlist_from_nadgrids(Par const& defn, ProjectionGrids & proj_grids)
- {
- pj_gridlist_from_nadgrids(defn.nadgrids,
- proj_grids.grids_storage().stream_policy,
- proj_grids.grids_storage().hgrids,
- proj_grids.hindexes);
- }
- } // namespace detail
- }}} // namespace boost::geometry::projections
- #endif // BOOST_GEOMETRY_SRS_PROJECTIONS_IMPL_PJ_GRIDLIST_HPP
|