经纬度坐标下求多边形的面积

最新推荐文章于 2024-09-13 16:04:01 发布

avcodec-decoder

最新推荐文章于 2024-09-13 16:04:01 发布

阅读量7.8k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： vs2008 文章标签：面积 boost

本文链接：https://blog.youkuaiyun.com/zhouhangjay/article/details/45822371

vs2008 专栏收录该内容

3 篇文章

订阅专栏

本文介绍了如何在经纬度坐标系下计算多边形面积，遇到不同投影方法不可行的问题，并最终采用boost库解决。通过实例演示了使用boost库求球面坐标面积的方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

经纬度坐标下求多边形的面积，mapinfo的值作为参考，试了几种投影都不行，原来应该计算球面坐标的面积。

最后采用了boost库的方法求面积。

对boost不熟悉的可以参考一下：

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>
namespace bg = boost::geometry;
#define MAPINFO_SEMIMAJOR   6370997.0
typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > DPoint;
typedef boost::geometry::model::polygon<DPoint> APolygon;
using boost::geometry::append;
double GetArea(vector<SGeoPoint>& vPt)
{
    APolygon sph_poly;
	for (size_t i = 0; i < vPt.size(); i++)
	{
		append(sph_poly,DPoint(vPt[i].x,vPt[i].y));
	}

	double area = bg::area(sph_poly);

	return fabs(area) *MAPINFO_SEMIMAJOR*MAPINFO_SEMIMAJOR;

	
}