CN109033681B - 一种边坡计算模型导入时的自动修复方法 - Google Patents

一种边坡计算模型导入时的自动修复方法 Download PDF

Info

Publication number
CN109033681B
CN109033681B CN201810916739.5A CN201810916739A CN109033681B CN 109033681 B CN109033681 B CN 109033681B CN 201810916739 A CN201810916739 A CN 201810916739A CN 109033681 B CN109033681 B CN 109033681B
Authority
CN
China
Prior art keywords
line
chain
point
line segment
intersection
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810916739.5A
Other languages
English (en)
Other versions
CN109033681A (zh
Inventor
柳侃
柳文强
黄树根
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Geological Engineering Surveying Institute Of Fujian Province
Original Assignee
Geological Engineering Surveying Institute Of Fujian Province
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Geological Engineering Surveying Institute Of Fujian Province filed Critical Geological Engineering Surveying Institute Of Fujian Province
Priority to CN201810916739.5A priority Critical patent/CN109033681B/zh
Publication of CN109033681A publication Critical patent/CN109033681A/zh
Application granted granted Critical
Publication of CN109033681B publication Critical patent/CN109033681B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/203—Drawing of straight lines or curves

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种边坡计算模型导入时的自动修复方法，包括图形线条自动修复步骤：步骤1、导入边坡计算模型；步骤2、计算图形的最大高度和宽度，设置带标尺的边框，下边框的长度为图形的宽度，左边框和右边框上的顶端纵坐标值大于图形的最大高度值；步骤3、手动识别出滑动面和水位线，然后对图形中的其余线条进行修正；步骤4、自动识别坡面线并连接，生成连续坡面线；步骤5、自动识别岩土线，并对不连续的岩土线进行延长或缩短，生成连续岩土线；步骤6、完成图形线条的修复。本发明自动修复方法能够自动纠错图形线条，用户可先进行滑动面、水位线定义，并自动识别坡面线和岩土线；能够自动构建封闭区域，形成土层模块，自动完成填土。

Description

一种边坡计算模型导入时的自动修复方法
技术领域
本发明涉及一种边坡计算模型导入时的自动修复方法，属于地质岩土软件领域。
背景技术
目前，国内一些岩土边坡(滑坡)计算软件，用户在导入边坡计算模型时，由于线型问题或者人为建模(一般都是通过AUTOCAD软件建模)时的失误，例如线条断开、重叠等，导致岩土边坡计算软件不能辨识计算模型，用户须重新检查、修改原有CAD模型后再次导入；其次，倘若用户根据个人经验或者后期补充勘察认为某段地层线须做适当调整，亦只能重新CAD建模而不能直接在原有的计算界面上直接进行调整；对于自动搜索最危险滑动面时，由于现有岩土边坡计算软件未对计算结果进行分类填充形成封闭区间云图，因此，用户无法对滑动面进行快速准确搜索。
请参阅图1，现有的边坡计算模型导入过程，仅判断线条是否为线段、圆弧、且相互连接，如果不是，要求用户再CAD软件中修正后再导入，如果是，判断图形是否分层且独立成区域，不是，则进行人工手动分层并手动填充区域，是，则，则可进行边坡安全系数计算、支档结构计算、滑坡推理计算、普通抗滑桩计算、锚索抗滑桩、挡土墙计算。边坡计算模型的应用都必须以导入后的图形为基础，且该图形必须由线段组成且线段连续，并形成独立的封闭区间云图。
现有技术中线条的修复大都必须通过CAD软件手动修正后才能再导入。图形线条未形成独立封闭区域，则采用人工手动分层填充，处理速度慢，耗时较长。
发明内容
为了解决上述技术问题，本发明提供一种边坡计算模型导入时的自动修复方法，其能够对导入的图形线条进行自动修复，自动生成连续坡面线和连续岩土线，同时，还能进一步生成封闭区域，完成自动填土。
本发明的技术方案如下：
一种边坡计算模型导入时的自动修复方法，包括图形线条自动修复步骤：
步骤1、导入边坡计算模型；
步骤2、计算图形的最大高度和宽度，设置带标尺的边框，下边框的长度为图形的宽度，左边框和右边框上的顶端纵坐标值大于图形的最大高度值；
步骤3、手动识别出滑动面和水位线，然后对图形中的其余线条进行修正；
步骤4、自动识别坡面线并连接，生成连续坡面线；
步骤5、自动识别岩土线，并对不连续的岩土线进行延长或缩短，生成连续岩土线；
步骤6、完成图形线条的修复。
更优地，所述步骤3中，对其余线条的修正为：对于线段，保持不变；对于折线，将折点处的线条打断，使折线转变为多个线段依次相连，从而将滑动面和水位线之外的所有线条均转变为线段，将得到的线段集合在一起，生成第一线段集合。
更优地，所述步骤4中，生成连续坡面的过程为：
4-1、逐条判断第一线段集合中的线段是否重叠，若重合，则保留重叠线段中的一条，删除多余的重叠线段；
4-2、对第一线段集合中的各线段逐条做中点辅助线，若线段的中点辅助线和第一线段集合中的其他线段无交点，将该线段定义为坡面线；若有交点，说明该线段不是坡面线；判断完毕将该中点辅助线删除；
4-3、将被定义为坡面线的线段按照各线段的横坐标从小到大进行排序，依次连接，若相邻两线段不连续，则自动将后一条坡面线与前一条坡面线连接，直至将全部坡面线连接完毕，形成连续坡面。
更优地，所述步骤4-2为：第一线段集合中的各线段逐条执行如下操作：获取线段的中点坐标值(x，y)，以(x，y+K)作起点，以(x，y0)做终点，画中点辅助线，若该中点辅助线与其他线段无交点，则将该线段定义为坡面线，所述y0为左边框或右边框上标尺的最大值，所述K≥0.001。
更优地，所述步骤5的执行过程如下：
5-1、将所述第一线段集合中除坡面线外的其余线段均定义为岩土线，并生成岩土线集合，从而完成岩土线的自动识别；
5-2、岩土线集合中的各线段逐条执行如下操作：
5-2-1、判断岩土线两端点是否均与其他岩土线有相交，若是，说明该岩土线不存在独立端点，不需要处理，选择下一条岩土线重复执行步骤5-2，若否，执行步骤5-2-2；所述独立端点为与其他岩土线不连接且不相交的端点；
5-2-2、判断该岩土线上是否有交点S1，所述交点S1是该岩土线与其他岩土线的交点，若岩土线存在交点S1，执行步骤5-2-3，若岩土线上无交点S1且两端点均为独立端点，则将其删除，选择下一条岩土线重复执行步骤5-2，若该岩土线无交点S1，且仅有一端点为独立端点，则从该独立端点出发，反向延长该岩土线直至与边框或者坡面线相交，生成交点A，将交点A作为该岩土线的新的端点，取消原来的独立端点；然后选择下一条岩土线重复执行步骤5-2；
5-2-3、对该岩土线上的独立端点逐一执行如下操作：计算独立端点与该岩土线上最近交点S1的距离L1，然后以该独立端点为起点做该岩土线的反向延长辅助线，直至与边框或者坡面线相交，判断该反向延长辅助线与其他岩土线是否相交，若无，则将与边框或坡面线的交点命名为交点S2，若有，则将反向延长线上与独立端点最近的交点命名为交点S2，计算该独立端点与该交点S2的距离L2，若距离L1≥距离L2，将该岩土线延长到交点S2，即交点S2作为该岩土线的端点，取消原独立端点，若距离L1＜距离L2，将该岩土线缩短到交点S1，即交点S1作为该岩土线的端点，取消原独立端点；然后删除反向延长辅助线，选择下一条岩土线重复执行步骤5-2；
5-3、各条岩土线均经过步骤5-2的操作，使任意一岩土线的两端点必然与其他岩土线的端点或边框或坡面线连接，生成多条连续岩土线。
更优地，所述自动修复方法还包括自动填土步骤：
步骤7、生成点链辅助线；
步骤8、创建初点链；
步骤9、构建每一个初始点链对应的整体性点链；
步骤10、创建由整体性点链形成的封闭区域，生成土层模块。
更优地，所述步骤7为：
7-1、自动记录所有坡面线、岩土线以及边框线段，并将统计到的线段集合在一起，生成第二线段集合；
7-2、从第二线段集合中依次导入其中一条线段Li，i为自然数,i的初始值为1，每循环一次，i＝i+1；
7-3、依次导入第二线段集合中其余线段的一条，每循环一次，导入一条其余线段；
7-4、判断线段Li是否与导入的其余线段，是，记录交点S3，执行步骤7-5，否，执行步骤7-5；
7-5、判断第二线段集合中的其余线段是否导入完毕，是，执行步骤7-6，否，返回步骤7-3；
7-6、将线段Li的端点和交点S3按照横坐标排序，并做点链辅助线，相连两点之间做一条点链辅助线，将所有点链辅助线记录在点链辅助线集合中；
7-7、判断线段Li是否全部导入完毕，若否，执行步骤7-2，若是，执行步骤8。
更优地，所述步骤8为：
8-1、统计出所有点链辅助线，并用集合A记录每一条点链辅助线的端点坐标；
8-2、从集合A中查找与左边框的横坐标大小相同的端点坐标，并将该端点坐标标记为端点坐标Bj，然后记录在集合B中；j为自然数，j的初始值为1，每标记一次，j＝j+1；全部查找完毕，执行步骤8-3；
8-3、依次导入集合B中的端点坐标Bj，每循环一次，j＝j+1；
8-4、依次输入点链辅助线集合中的点链辅助线，每循环一次导入一条点链辅助线；
8-5、判断导入的点链辅助线的端点坐标是否为当前导入的坐标Bj，若是，执行步骤8-6，若否，执行步骤8-7；
8-6、创建点链Ck，其记录坐标Bj和该点链辅助线的另一个端点坐标，每创建一个点链Ck，k＝k+1；该点链Ck为初点链；
8-7、判断点链辅助线集合中点链辅助线是否导入完毕，是，执行步骤8-8，否，执行步骤8-4；
8-8、判断集合B中的坐标点Bj是否导入完毕，是，执行步骤9，否，返回步骤8-3。
更优地，所述步骤9为：
9-1、将所有点链Ck导入点链集合C，从点链集合中选取一个点链Ck，并将被选取的点链Ck的最后一个坐标点记录为E；每循环一次，k＝k+1；
9-2、依次导入点链辅助线集合中的点链辅助线，每循环一次导入一条点链辅助线；
9-3、判断点链辅助线的端点坐标是否为E，若是，创建点链Pkm，其记录点链Ck上所有坐标，以及该点链辅助线的另一端点坐标，然后执行步骤9-4；若否，执行步骤9-4；m为自然数，每创建一次点链Pkm，m＝m+1；
9-4、判断点链辅助线集合中的各点链辅助线是否均导入完毕，若否，返回步骤9-2，若是，删除点链Ck，执行步骤9-5；
9-5、判断点链Pkm的最后一个坐标的横坐标是否等于右边框的横坐标，若是，执行步骤9-6，若否，将点链Pkm重命名为点链Ck，此步骤中k值大于初点链的总数量，避免出现相同k值编号的点链，然后返回步骤9-1；
9-6、判断点链集合C中的各点链是否均导入完毕，若是，则点链Pkm为整体性点链，执行步骤10，若否，返回步骤9-1。
更优地，所述步骤10为：
10-1、将所有点链Pkm统一重命名为点链Dx；将所述点链Dx集合在一起，生成待排序点链集合M；
10-2、依次导入点链Dx，且将该点链Dx的第一个坐标点标记为坐标点F；每循环导入一次，x＝x+1；
10-3、依次导入除步骤10-2中的点链Dx外的其余的点链Dx，并将导入的该条其余点链Dx的第一个坐标点标记为坐标点G，每循环一次，导入一条其余点链Dx；
10-4、判断坐标点F的纵坐标是否大于坐标点G的纵坐标：
若是，创建点链Hy，使点链Hy＝点链Dx，第一次创建点链Hy时，y＝1，以后每创建一次，y＝y+1；然后执行步骤10-5；
若否，判断坐标点F的纵坐标是否等于坐标点G的纵坐标：若是，则依次分别将点链Dx和该条其余点链Dx的下一个坐标点标记为坐标点F和坐标点G，该下一个坐标点的横坐标相同，再次执行步骤10-4；若否，创建点链Hy，使点链Hy＝该条其余点链Dx，然后执行步骤10-5；
10-5、判断待排序点链集合M中的其余点链Dx是否均已逐一导入完毕：
若是，从待排序点链集合M中删除与点链Hy等值的点链Dx，或删除与点链Hy等值的其余点链Dx，然后执行步骤10-6；
若否，返回步骤10-3；
10-6、判断待排序点链集合M的点链Dx是否均已执行步骤10-2：
若是，则完成排序，由点链Hy和点链Hy+1前后两组点链分别，构成多边形，形成封闭区域，生成土层模块，自动填土步骤完成；
若否，返回步骤10-2。
本发明具有如下有益效果：
1、本发明自动修复方法能够自动纠错图形线条，用户可先进行滑动面、水位线定义，并自动识别坡面线和岩土线；
2、能够自动构建封闭区域，形成土层模块，自动完成填土；
3、导入dxf计算模型方便快捷，减少因线型问题或者人为操作失误产生的重复建模工作量；
4、边坡模型导入后，经过自动修复，使得图形能够用于边坡安全系数计算、支档结构计算、滑坡推理计算、普通抗滑桩计算、锚索抗滑桩、挡土墙计算。
附图说明
图1为现有技术中边坡计算模型导入流程；
图2为本发明自动修复方法的总流程图；
图3为本发明中图形线条自动修复流程图；
图4为本发明中线条修正流程图；
图5为本发明线条修正前后的图形线条对比图；
图6为本发明中生成连续坡面线的流程图；
图7为本发明中做中点辅助线的流程示意图；
图8为本发明中点辅助线的图形示意图；
图9为本发明中生成连续岩土线的流程图；
图10为本发明岩土线延长或缩短的图形示意图；
图11为本发明中生成点链辅助线的流程图；
图12为本发明生成点链辅助线的图形示意图；
图13为本发明中创建初点链的流程图；
图14为本发明中构建整体性点链的流程图；
图15为本发明中创建封闭区域的流程图。
具体实施方式
下面结合附图和具体实施例来对本发明进行详细的说明。
实施例一
请参阅图2和图3，一种边坡计算模型导入时的自动修复方法，包括图形线条自动修复步骤：
步骤1、导入边坡计算模型；
步骤2、计算图形的最大高度和宽度，设置带标尺的边框，下边框的长度为图形的宽度，左边框和右边框上的顶端纵坐标值大于图形的最大高度值；
步骤3、手动识别出滑动面和水位线，然后对图形中的其余线条进行修正；其余线条的修正过程为：请参阅图4和图5，对于线段，保持不变；对于折线，将折点处的线条打断，使折线转变为多个线段依次相连，从而将滑动面和水位线之外的所有线条均转变为线段，(图5中，虚线为水位线、双点划线为滑动面)将得到的线段集合在一起，生成第一线段集合；本发明所述线段包括直线线段和弧线线段，直线线段的两端点之间的连线为直线，弧线线段的两端点之间的连线为弧线；
请参阅图6，步骤4、自动识别坡面线并连接，生成连续坡面线；
4-1、逐条判断第一线段集合中的线段是否重叠，若重合，则保留重叠线段中的一条，删除多余的重叠线段；
4-2、对第一线段集合中的各线段逐条做中点辅助线，若线段的中点辅助线和第一线段集合中的其他线段无交点，将该线段定义为坡面线；若有交点，说明该线段不是坡面线；判断完毕将该中点辅助线删除；
请参阅图7和图8，所述中点辅助线的实施方式为：
第一线段集合中的各线段逐条执行如下操作：获取线段的中点坐标值(x，y)，以(x，y+K)作起点，以(x，y0)做终点，画中点辅助线，若该中点辅助线与其他线段(即第一线段集合中的其他线段)无交点，则将该线段定义为坡面线，所述y0为左边框或右边框上标尺的最大值，若左边框上标尺最大值≥右边框上标尺最大值，则y0为左边框上标尺的最大值，反之，y0为右边框上标尺的最大值；所述K≥0.001；图8中，K＝0.001，图中左边的重点辅助线与其他线段没有交点，则对应的线段为坡面线，图中右边的重点辅助线与其他线段产生交点，因此，该对应的线段非坡面线；
4-3、将被定义为坡面线的线段按照各线段的横坐标从小到大进行排序，依次连接，若相邻两线段不连续，则自动将后一条坡面线与前一条坡面线连接，直至将全部坡面线连接完毕，形成连续坡面；
请参阅图9，步骤5、自动识别岩土线，并对不连续的岩土线进行延长或缩短，生成连续岩土线；
5-1、将所述第一线段集合中除坡面线外的其余线段均定义为岩土线，并生成岩土线集合，从而完成岩土线的自动识别；
5-2、岩土线集合中的各线段逐条执行如下操作：
5-2-1、判断岩土线两端点是否均与其他岩土线有相交，若是，说明该岩土线不存在独立端点，不需要处理，选择下一条岩土线重复执行步骤5-2，若否，执行步骤5-2-2；所述独立端点为与其他岩土线不连接且不相交的端点；
5-2-2、判断该岩土线上是否有交点S1，所述交点S1是该岩土线与其他岩土线的交点：
若该岩土线存在交点S1(交点S1可能是1个，也可能是多个)，执行步骤5-2-3，
若该岩土线无交点S1且两端点均为独立端点，则将其删除，选择下一条岩土线重复执行步骤5-2，
若该岩土线无交点S1且仅有一端点为独立端点，则从该独立端点出发，反向延长该线段直至与边框或者坡面线相交，判断该反向延长辅助线与其他岩土线是否相交，若无，则将与边框或坡面线的交点命名为交点S2，若有，则将反向延长线上与独立端点最近的交点命名为交点S2，将交点A作为该线段的新的端点，取消原来的独立端点；选择下一条岩土线重复执行步骤5-2；
5-2-3、对该岩土线上的独立端点逐一执行如下操作：计算独立端点与该岩土线上最近交点S1的距离L1，然后以该独立端点为起点做该岩土线的反向延长辅助线，直至与边框或者坡面线相交，生成交点S2，计算该独立端点与该交点S2的距离L2，若距离L1≥距离L2，将该岩土线延长到交点S2，即交点S2作为该岩土线的端点，取消原独立端点，若距离L1＜距离L2，将该岩土线缩短到交点S1，即交点S1作为该岩土线的端点，取消原独立端点；然后删除反向延长辅助线；选择下一条岩土线重复执行步骤5-2；请参阅图10，线段C1C2上，端点C1为独立端点，端点C1与该线段上距离最近的交点S1之间的距离为L1，从端点C1反向延长后，仅与左边框产生交点S2，端点C1与交点S2之间的距离为L2，由于L1＞L2，因此将该岩土线延长至交点S2，将交点S2作为新的端点，取消原独立端点C1；再看线段C3C4，其中端点C3为独立端点，端点C3与该线段上距离最近的交点S1之间的距离为L1，从端点C3反向延长后，首先与其他线段产生交点S2，因此，将端点C3与交点S2之间的距离记为L2，由于L1＜L2，因此将该岩土线缩短至交点S1，将交点S1作为新的端点，取消原独立端点C3；
5-3、各条岩土线均经过步骤5-2的操作，使任意一岩土线的两端点必然与其他岩土线的端点或边框或坡面线连接，生成多条连续岩土线。
步骤6、完成图形线条的修复。
实施例二
请参阅图2，本实施例二在完成图形线条修复的基础上，完成自动填土。
所述自动修复方法的自动填土步骤，包括：
请参阅图11，步骤7、生成点链辅助线：
7-1、自动记录所有坡面线、岩土线以及边框线段，并将统计到的线段集合在一起，生成第二线段集合；
7-2、从第二线段集合中依次导入其中一条线段Li，i为自然数，i的初始值为1，每循环一次，i＝i+1；
7-3、依次导入第二线段集合中其余线段的一条，每循环一次，导入一条其余线段；所述其余线段是指除步骤7-2导入的线路Li之外的线段；
7-4、判断线段Li是否与导入的其余线段相交，是，记录交点S3，执行步骤7-5，否，执行步骤7-5；
7-5、判断第二线段集合中的其余线段是否导入完毕，是，执行步骤7-6，否，返回步骤7-3；
7-6、将线段Li的端点和交点S3按照横坐标排序，并做点链辅助线，相连两点之间做一条点链辅助线，将所有点链辅助线记录在点链辅助线集合中；
7-7、判断线段Li是否全部导入完毕，若否，执行步骤7-2，若是，执行步骤8；
请参阅图12(为避免虚线混淆，图12中将滑动面改为实线，且未示出水位线)，以先导入线段L1为例，多次循环判断后，仅有一条线段与其存在交点S3，因此，生成对应的点链辅助线为T1和T2，第二次导入线段L2，多次循环判断后，仅有一条线段与其存在交点S3，因此，生成对应的点链辅助线T3和T4；
请参阅图13，步骤8、创建初点链：
8-1、统计出所有点链辅助线，并用集合A记录每一条点链辅助线的端点坐标；
8-2、从集合A中查找横坐标与左边框的横坐标大小相同的端点坐标，并将该端点坐标标记为端点坐标Bj，然后记录在集合B中；j为自然数，j的初始值为1，每标记一次，j＝j+1；全部查找完毕，执行步骤8-3；
8-3、依次导入集合B中的端点坐标Bj，可以从端点坐标B1开始导入，每循环一次，j＝j+1；
8-4、依次输入点链辅助线集合中的点链辅助线，每循环一次导入一条点链辅助线；
8-5、判断导入的点链辅助线的端点坐标是否为当前导入的坐标Bj，若是，执行步骤8-6，若否，执行步骤8-7；
8-6、创建点链Ck，其记录坐标Bj和该点链辅助线的另一个端点坐标，每创建一个点链Ck，k＝k+1；该点链Ck为初点链；
8-7、判断点链辅助线集合中点链辅助线是否导入完毕，是，执行步骤8-8，否，执行步骤8-4；
8-8、判断集合B中的坐标点Bj是否导入完毕，是，执行步骤9，否，返回步骤8-3。
请参阅图14，步骤9、构建每一个初始点链对应的整体性点链；
9-1、将所有点链Ck导入点链集合C，从点链集合中选取一个点链Ck，并将被选取的点链Ck的最后一个坐标点记录为E；每循环一次，k＝k+1；
9-2、依次导入点链辅助线集合中的点链辅助线，每循环一次导入一条点链辅助线；
9-3、判断点链辅助线的端点坐标是否为E，若是，创建点链Pkm，其记录点链Ck上所有坐标，以及该点链辅助线的另一端点坐标，然后执行步骤9-4；若否，执行步骤9-4；m为自然数，每创建一次点链Pkm，m＝m+1；
9-4、判断点链辅助线集合中的各点链辅助线是否均导入完毕，若否，返回步骤9-2，若是，删除点链Ck，执行步骤9-5；
9-5、判断点链Pkm的最后一个坐标的横坐标是否等于右边框的横坐标，若是，执行步骤9-6，若否，将点链Pkm重命名为点链Ck，此步骤中k值大于初点链的总数量，避免出现相同k值编号的点链，然后返回步骤9-1；
9-6、判断点链集合C中的各点链是否均导入完毕，若是，则点链Pkm为整体性点链，执行步骤10，若否，返回步骤9-1；
请参阅图15，步骤10、创建由整体性点链形成的封闭区域，生成土层模块：
10-1、将点链Pkm统一重命名为点链Dx；将所述点链Dx集合在一起，生成待排序点链集合M；即将所有点链Pkm进行排序编号，X的最大值等于点链Pkm；例如，将点链P12、P13、P21、P22进行重命名，对应得到点链D1、D2、D3、D4；
10-2、依次导入点链Dx，且将该点链Dx的第一个坐标点标记为坐标点F；每循环导入一次，x＝x+1；
10-3、依次导入除步骤10-2中的点链Dx外的其余的点链Dx，并将导入的该条其余点链Dx的第一个坐标点标记为坐标点G，每循环一次，导入一条其余点链Dx；
10-4、判断坐标点F的纵坐标是否大于坐标点G的纵坐标：
若是，创建点链Hy，使点链Hy＝点链Dx，第一次创建点链Hy时，y＝1，以后每创建一次，y＝y+1；然后执行步骤10-5；
若否，判断坐标点F的纵坐标是否等于坐标点G的纵坐标：若是，则依次分别将点链Dx和该条其余点链Dx的下一个坐标点标记为坐标点F和坐标点G，该下一个坐标点的横坐标相同，再次执行步骤10-4；若否，创建点链Hy，使点链Hy＝该条其余点链Dx，然后执行步骤10-5；
10-5、判断待排序点链集合M中的其余点链Dx是否均已逐一导入完毕：
若是，从待排序点链集合M中删除与点链Hy等值的点链Dx，或删除与点链Hy等值的其余点链Dx，然后执行步骤10-6；
若否，返回步骤10-3；
10-6、判断待排序点链集合M的点链Dx是否均已执行步骤10-2：
若是，则完成排序，由点链Hy和点链Hy+1前后两组点链分别，构成多边形，形成封闭区域，生成土层模块，自动填土步骤完成；
若否，返回步骤10-2。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

Claims (10)

1.一种边坡计算模型导入时的自动修复方法，其特征在于，包括图形线条自动修复步骤：
步骤1、导入边坡计算模型；
步骤2、计算图形的最大高度和宽度，设置带标尺的边框，下边框的长度为图形的宽度，左边框和右边框上的顶端纵坐标值大于图形的最大高度值；
步骤3、手动识别出滑动面和水位线，然后对图形中的其余线条进行修正；
步骤4、自动识别坡面线并连接，生成连续坡面线；
步骤5、自动识别岩土线，并对不连续的岩土线进行延长或缩短，生成连续岩土线；
步骤6、完成图形线条的修复。
2.根据权利要求1所述的一种边坡计算模型导入时的自动修复方法，其特征在于：所述步骤3中，对其余线条的修正为：对于线段，保持不变；对于折线，将折点处的线条打断，使折线转变为多个线段依次相连，从而将滑动面和水位线之外的所有线条均转变为线段，将得到的线段集合在一起，生成第一线段集合。
3.根据权利要求2所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤4中，生成连续坡面的过程为：
4-1、逐条判断第一线段集合中的线段是否重叠，若重合，则保留重叠线段中的一条，删除多余的重叠线段；
4-2、对第一线段集合中的各线段逐条做中点辅助线，若线段的中点辅助线和第一线段集合中的其他线段无交点，将该线段定义为坡面线；若有交点，说明该线段不是坡面线；判断完毕将该中点辅助线删除；
4-3、将被定义为坡面线的线段按照各线段的横坐标从小到大进行排序，依次连接，若相邻两线段不连续，则自动将后一条坡面线与前一条坡面线连接，直至将全部坡面线连接完毕，形成连续坡面。
4.根据权利要求3所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤4-2为：第一线段集合中的各线段逐条执行如下操作：获取线段的中点坐标值(x，y)，以(x，y+K)作起点，以(x，y0)做终点，画中点辅助线，若该中点辅助线与其他线段无交点，则将该线段定义为坡面线，所述y0为左边框或右边框上标尺的最大值，所述K≥0.001。
5.根据权利要求3所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤5的执行过程如下：
5-1、将所述第一线段集合中除坡面线外的其余线段均定义为岩土线，并生成岩土线集合，从而完成岩土线的自动识别；
5-2、岩土线集合中的各线段逐条执行如下操作：
5-2-1、判断岩土线两端点是否均与其他岩土线有相交，若是，说明该岩土线不存在独立端点，不需要处理，选择下一条岩土线重复执行步骤5-2，若否，执行步骤5-2-2；所述独立端点为与其他岩土线不连接且不相交的端点；
5-2-2、判断该岩土线上是否有交点S1，所述交点S1是该岩土线与其他岩土线的交点，若岩土线上存在交点S1，执行步骤5-2-3，若岩土线上无交点S1且两端点均为独立端点，则将其删除，选择下一条岩土线重复执行步骤5-2；若该岩土线无交点S1，且仅有一端点为独立端点，则从该独立端点出发，反向延长该岩土线直至与边框或者坡面线相交，生成交点A，将交点A作为该岩土线的新的端点，取消原来的独立端点，选择下一条岩土线重复执行步骤5-2；
5-2-3、对该岩土线上的独立端点逐一执行如下操作：计算独立端点与该岩土线上最近交点S1的距离L1，然后以该独立端点为起点做该岩土线的反向延长辅助线，直至与边框或者坡面线相交，判断该反向延长辅助线与其他岩土线是否相交，若否，则将与边框或坡面线的交点命名为交点S2，若是，则将反向延长线上与独立端点最近的交点命名为交点S2，计算该独立端点与该交点S2的距离L2，若距离L1≥距离L2，将该岩土线延长到交点S2，即交点S2作为该岩土线的端点，取消原独立端点，若距离L1＜距离L2，将该岩土线缩短到交点S1，即交点S1作为该岩土线的端点，取消原独立端点；然后删除反向延长辅助线，选择下一条岩土线重复执行步骤5-2；
5-3、各条岩土线均经过步骤5-2的操作，使任意一岩土线的两端点必然与其他岩土线的端点或边框或坡面线连接，生成多条连续岩土线。
6.根据权利要求2所述的一种边坡计算模型导入时的自动修复方法，其特征在于，自动修复方法还包括自动填土步骤：
步骤7、生成点链辅助线；
步骤8、创建初点链；
步骤9、构建每一条初始点链对应的整体性点链；
步骤10、创建由整体性点链形成的封闭区域，生成土层模块。
7.根据权利要求6所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤7为：
7-1、自动记录所有坡面线、岩土线以及边框线段，并将统计到的线段集合在一起，生成第二线段集合；
7-2、从第二线段集合中依次导入其中一条线段Li，i为自然数,i的初始值为1，每循环一次，i＝i+1；
7-3、依次导入第二线段集合中其余线段的一条，每循环一次，导入一条其余线段；
7-4、判断线段Li是否与导入的其余线段相交，是，记录交点S3，执行步骤7-5，否，执行步骤7-5；
7-5、判断第二线段集合中的其余线段是否导入完毕，是，执行步骤7-6，否，返回步骤7-3；
7-6、将线段Li的端点和交点S3按照横坐标排序，并做点链辅助线，相连两点之间做一条点链辅助线，将所有点链辅助线记录在点链辅助线集合中；
7-7、判断线段Li是否全部导入完毕，若否，执行步骤7-2，若是，执行步骤8。
8.根据权利要求7所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤8为：
8-1、统计出所有点链辅助线，并用集合A记录每一条点链辅助线的端点坐标；
8-2、从集合A中查找横坐标与左边框的横坐标大小相同的端点坐标，并将该端点坐标标记为端点坐标Bj，然后记录在集合B中；j为自然数，j的初始值为1，每标记一次，j＝j+1；全部查找完毕，执行步骤8-3；
8-3、依次导入集合B中的坐标点Bj，每循环一次，j＝j+1；
8-4、依次输入点链辅助线集合中的点链辅助线，每循环一次导入一条点链辅助线；
8-5、判断导入的点链辅助线的端点坐标是否为当前导入的坐标Bj，若是，执行步骤8-6，若否，执行步骤8-7；
8-6、创建点链Ck，其记录坐标Bj和该点链辅助线的另一个端点坐标，k的初始值为1，每创建一个点链Ck，k＝k+1；该点链Ck为初点链；
8-7、判断点链辅助线集合中点链辅助线是否导入完毕，是，执行步骤8-8，否，执行步骤8-4；
8-8、判断集合B中的坐标点Bj是否导入完毕，是，执行步骤9，否，返回步骤8-3。
9.根据权利要求8所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤9为：
9-1、将所有点链Ck导入点链集合C，从点链集合中选取一个点链Ck，并将被选取的点链Ck的最后一个坐标点记录为E；每循环一次，k＝k+1；
9-2、依次导入点链辅助线集合中的点链辅助线，每循环一次导入一条点链辅助线；
9-3、判断点链辅助线的端点坐标是否为E，若是，创建点链Pkm，其记录点链Ck上所有坐标，以及该点链辅助线的另一端点坐标，然后执行步骤9-4；若否，执行步骤9-4；m为自然数，每创建一次点链Pkm，m＝m+1；
9-4、判断点链辅助线集合中的各点链辅助线是否均导入完毕，若否，返回步骤9-2，若是，删除点链Ck，执行步骤9-5；
9-5、判断点链Pkm的最后一个坐标的横坐标是否等于右边框的横坐标，若是，执行步骤9-6，若否，将点链Pkm重命名为点链Ck，此步骤中k值大于初点链的总数量，避免出现相同k值编号的点链，然后返回步骤9-1；
9-6、判断点链集合C中的各点链是否均导入完毕，若是，则点链Pkm为整体性点链，执行步骤10，若否，返回步骤9-1。
10.根据权利要求9所述的一种边坡计算模型导入时的自动修复方法，其特征在于，所述步骤10为：
10-1、将所有点链Pkm统一重命名为点链Dx；将所述点链Dx集合在一起，生成待排序点链集合M；
10-2、依次导入点链Dx，且将该点链Dx的第一个坐标点标记为坐标点F；每循环导入一次，x＝x+1；
10-3、依次导入除步骤10-2中的点链Dx外的其余的点链Dx，并将导入的该条其余点链Dx的第一个坐标点标记为坐标点G，每循环一次，导入一条其余点链Dx；
10-4、判断坐标点F的纵坐标是否大于坐标点G的纵坐标：
若是，创建点链Hy，使点链Hy＝点链Dx，第一次创建点链Hy时，y＝1，以后每创建一次，y＝y+1；然后执行步骤10-5；
若否，判断坐标点F的纵坐标是否等于坐标点G的纵坐标：若是，则依次分别将点链Dx和该条其余点链Dx的下一个坐标点标记为坐标点F和坐标点G，该下一个坐标点的横坐标相同，再次执行步骤10-4；若否，创建点链Hy，使点链Hy＝该条其余点链Dx，然后执行步骤10-5；
10-5、判断待排序点链集合M中的其余点链Dx是否均已逐一导入完毕：
若是，从待排序点链集合M中删除与点链Hy等值的点链Dx，或删除与点链Hy等值的其余点链Dx，然后执行步骤10-6；
若否，返回步骤10-3；
10-6、判断待排序点链集合M的点链Dx是否均已执行步骤10-2：
若是，则完成排序，由点链Hy和点链Hy+1前后两组点链两两分别构成多边形，形成封闭区域，生成土层模块，自动填土步骤完成；
若否，返回步骤10-2。
CN201810916739.5A 2018-08-13 2018-08-13 一种边坡计算模型导入时的自动修复方法 Active CN109033681B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810916739.5A CN109033681B (zh) 2018-08-13 2018-08-13 一种边坡计算模型导入时的自动修复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810916739.5A CN109033681B (zh) 2018-08-13 2018-08-13 一种边坡计算模型导入时的自动修复方法

Publications (2)

Publication Number Publication Date
CN109033681A CN109033681A (zh) 2018-12-18
CN109033681B true CN109033681B (zh) 2019-07-05

Family

ID=64633853

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810916739.5A Active CN109033681B (zh) 2018-08-13 2018-08-13 一种边坡计算模型导入时的自动修复方法

Country Status (1)

Country Link
CN (1) CN109033681B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111914330B (zh) * 2020-08-07 2022-06-17 合肥市市政设计研究总院有限公司 基于图解试算法的土岩结合边坡稳定分析方法
CN111951256A (zh) * 2020-08-20 2020-11-17 南京中车浦镇城轨车辆有限责任公司 一种切割下料图形缺陷自动检测修复的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413000A (zh) * 2013-08-19 2013-11-27 浙江浙大万维科技有限公司 自动绘制斜坡线方法和装置
CN106650127A (zh) * 2016-12-28 2017-05-10 北京华大九天软件有限公司 一种自动将多段线转换成封闭图形的方法
CN108345751A (zh) * 2018-02-13 2018-07-31 福建省地质工程勘察院 用AutoCAD插件进行岩土工程边坡建模方法及设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20100078433A (ko) * 2008-12-30 2010-07-08 주식회사 동부하이텍 사다리꼴 금속배선의 시뮬레이션을 위한 모델링 구조

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413000A (zh) * 2013-08-19 2013-11-27 浙江浙大万维科技有限公司 自动绘制斜坡线方法和装置
CN106650127A (zh) * 2016-12-28 2017-05-10 北京华大九天软件有限公司 一种自动将多段线转换成封闭图形的方法
CN108345751A (zh) * 2018-02-13 2018-07-31 福建省地质工程勘察院 用AutoCAD插件进行岩土工程边坡建模方法及设备

Also Published As

Publication number Publication date
CN109033681A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109033681B (zh) 一种边坡计算模型导入时的自动修复方法
CN104868465B (zh) 基于模糊机会约束的电力系统网架重构优化方法
CN107330030B (zh) 一种面向海量路网数据压缩存储的层次网络构建方法
CN103886023B (zh) Excel数据表的存储、提取方法及系统
CN106228183A (zh) 一种半监督学习分类方法与装置
CN109635414A (zh) 一种风力发电机组风轮叶片有限元建模方法
CN113223116B (zh) 一种基于AutoCAD的滑边坡地质模型、构建方法及评价方法
CN106648608B (zh) 一种灵活方便的流程图构建方法
CN104574449B (zh) 基于dem的投影面积计算方法
CN101604419A (zh) 一种建设项目拟选址自动提取的方法
CN104820733A (zh) 一种高速列车需求元模型建立方法和装置
CN106548278A (zh) 一种高标准基本农田建设项目区优选自动化计算方法
CN106844979A (zh) 一种城市虚拟公园交互式组合设计系统及方法
CN108256218A (zh) 一种基于实测腰线数据的地下连通巷道精细建模方法
CN105138711B (zh) 一种图元索引、检测方法及装置
CN102339473B (zh) 一种基于遗传算法的动漫模型设计方法
CN101398943A (zh) 基于双层规则化的渠道地面网格模型表示方法
CN103886129B (zh) 将测井数据离散到储层网格模型的方法和装置
CN106446397A (zh) 一种基于状态目标性的启发式航天器任务规划方法
CN108629036A (zh) 一种道路缩编方法及装置
Satoh et al. Interactive system for arranging issues based on PROLEG in civil litigation
CN101894063B (zh) 一种用于微处理器功能验证的测试程序生成方法及装置
CN102411646A (zh) 配网数据的处理方法及装置
CN108052743A (zh) 一种阶梯接近中心度确定方法及系统
KR101063827B1 (ko) 한국토지정보시스템 연속지적도와 수치지형도의 기하학적 지도 변환을 위한 반자동화된 공액점 쌍 추출방법

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant