图形/多段线内缩外扩思路

luuuuuuuuuuuuuuuuuu

已于 2023-05-11 22:06:05 修改

阅读量2.5k

点赞数 1

分类专栏：兴趣使然文章标签： c++ c语言算法

于 2022-09-29 21:26:58 首次发布

本文链接：https://blog.csdn.net/yoyo18520/article/details/127114050

版权

兴趣使然专栏收录该内容

18 篇文章 0 订阅

订阅专栏

博主分享了在处理多段线内缩外扩时遇到的问题，即现有算法库如Clipper和boost::geometry::buffer会改变点位顺序，不适用于需要保持原始点位信息的场景。为解决这一问题，提出了使用直线段平移的方法，保留顶点顺序，适用于带有圆弧的多段线。虽然方法较为简单，但能有效满足需求。博客还提及了几何理论和直线处理的相关知识，并链接了一个可能的解决方案实现。

摘要由CSDN通过智能技术生成

图形/多段线内缩外扩思路

前言

我在网上找了很多关于多边形内缩外扩的资料，也测试了一些算法库，如Clipper，boost::geometry::buffer，其实都可以得到不错的结果。

`Clipper`

例程

#include "clipper.hpp"   
... 
using namespace ClipperLib; 
  
int main() 
{ 
  Path subj; 
  Paths solution; 
  subj <<  
    IntPoint(348,257) << IntPoint(364,148) << IntPoint(362,148) <<  
    IntPoint(326,241) << IntPoint(295,219) << IntPoint(258,88) <<  
    IntPoint(440,129) << IntPoint(370,196) << IntPoint(372,275); 
  ClipperOffset co; 
  co.AddPath(subj, jtRound, etClosedPolygon); 
  co.Execute(solution, -7.0); 
    
  //draw solution ... 
  DrawPolygons(solution, 0x4000FF00, 0xFF009900); 
}

效果

`boost::geometry::buffer`

例程

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/geometries.hpp>


int main()
{
    typedef double coordinate_type;
    typedef boost::geometry::model::d2::point_xy<coordinate_type> point;
    typedef boost::geometry::model::polygon<point> polygon;

    // Declare strategies
    const double buffer_distance = 1.0;
    const int points_per_circle = 36;
    boost::geometry::strategy::buffer::distance_symmetric<coordinate_type> distance_strategy(buffer_distance);
    boost::geometry::strategy::buffer::join_round join_strategy(points_per_circle);
    boost::geometry::strategy::buffer::end_round end_strategy(points_per_circle);
    boost::geometry::strategy::buffer::point_circle circle_strategy(points_per_circle);
    boost::geometry::strategy::buffer::side_straight side_strategy;

    // Declare output
    boost::geometry::model::multi_polygon<polygon> result;

    // Declare/fill a linestring
    boost::geometry::model::linestring<point> ls;
    boost::geometry::read_wkt("LINESTRING(0 0,4 5,7 4,10 6)", ls);

    // Create the buffer of a linestring
    boost::geometry::buffer(ls, result,
                distance_strategy, side_strategy,
                join_strategy, end_strategy, circle_strategy);


    // Declare/fill a multi point
    boost::geometry::model::multi_point<point> mp;
    boost::geometry::read_wkt("MULTIPOINT((3 3),(4 4),(6 2))", mp);

    // Create the buffer of a multi point
    boost::geometry::buffer(mp, result,
                distance_strategy, side_strategy,
                join_strategy, end_strategy, circle_strategy);


    // Declare/fill a multi_polygon
    boost::geometry::model::multi_polygon<polygon> mpol;
    boost::geometry::read_wkt("MULTIPOLYGON(((0 1,2 5,5 3,0 1)),((1 1,5 2,5 0,1 1)))", mpol);

    // Create the buffer of a multi polygon
    boost::geometry::buffer(mpol, result,
                distance_strategy, side_strategy,
                join_strategy, end_strategy, circle_strategy);


    return 0;
}

效果

遇到的问题

我现在得到的图形是经过反向解析文件再生成的，每个点位都是包含了一些信息，原始的点位个数和顺序是不可以轻易被改变的，而且图形也不一定是封闭的，其组成是可能带圆弧的多段线，后面再生成的文件还得和原始图形的点位顺序对应上才行。而以上两种处理库，是会改变点传入点位的个数和顺序的，这样的结果并不能正确生成理想的文件，就是这个问题让我很是头疼。

后来尝试找些几何理论的资料来解决，而网上的不少资料都是通过向量的方式来求夹角的，在出现内缩到一定程度，会有自交的情况，需要自己跳过垂直平行重复的点位，显然这也不是我想要的结果。