123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229 |
- #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
- {
- 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;
- }
- 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);
- }
- }
- 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)
- {
-
-
-
- if (pj_gridlist_find_all(gridname, grids.gridinfo, gridindexes))
- return true;
- std::size_t orig_size = grids.gridinfo.size();
-
- typename StreamPolicy::stream_type is;
- stream_policy.open(is, gridname);
- if (! pj_gridinfo_init(gridname, is, grids.gridinfo))
- {
- return false;
- }
-
- pj_gridlist_add_seq_inc(gridindexes, orig_size, grids.gridinfo.size());
- return true;
- }
- 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)
- {
-
-
-
- {
- typename SharedGrids::read_locked lck_grids(grids);
- if (pj_gridlist_find_all(gridname, lck_grids.gridinfo, gridindexes))
- return true;
- }
-
- typename StreamPolicy::stream_type is;
- stream_policy.open(is, gridname);
- pj_gridinfo new_grids;
- if (! pj_gridinfo_init(gridname, is, new_grids))
- {
- return false;
- }
-
- std::size_t orig_size = 0;
- std::size_t new_size = 0;
- {
- typename SharedGrids::write_locked lck_grids(grids);
-
-
- 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;
- }
- 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)
- {
-
- 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);
- }
- }
- }}}
- #endif
|