空间管理您的位置: MSCBSC | 移动通信网 » 中华儿女(QQ494317616) » 日志

交点

上一篇 / 下一篇 2008-06-16 10:19:01 / 个人分类：网上的一些源程序

查看( 14 ) / 评论( 0 ) / 评分( 0 / 0 )

include "mapbasic.def"
Declare sub Main

sub main
Dim obj1,obj2,obj3 as object
Dim iRow1,iRow2 as integer
Dim iNext as integer
Dim x1,y1,x2,y2,x,y as float
Dim objTemp as object
Dim iNode,iSec as integer
Dim i as integer
Dim len1,len2 as float
Dim min as float

min=30 ´最小线段长度，因为我发觉MapInfo中计算得到的两Line交点与真正截断的点之间存在误差
set distance units "m"
select * from 国省道 into STable

fetch first from STable ´first

while not eot(STable)
obj1=STable.obj
iRow1=STable.rowid

select * from 国省道 into temp where obj intersects obj1 and rowid<>iRow1

if TableInfo(temp,TAB_INFO_NROWS)>0 then
fetch first from temp
while not eot(temp)
obj2=temp.obj
iRow2=temp.rowid
obj3=IntersectNodes (obj1, obj2, INCL_CROSSINGS)
iNode=objectinfo(obj3,OBJ_INFO_NPNTS)
if iNode>0 and objectinfo(obj3,OBJ_INFO_TYPE)=OBJ_TYPE_PLINE then
´相交点
x=ObjectNodeX(obj3,1,1)
y=ObjectNodeY(obj3,1,1)

if objectinfo(obj2,OBJ_INFO_TYPE)=OBJ_TYPE_LINE then
´相交线起点
x1=ObjectGeography(obj2,OBJ_GEO_LINEBEGX)
y1=ObjectGeography(obj2,OBJ_GEO_LINEBEGY)

´相交线终点
x2=ObjectGeography(obj2,OBJ_GEO_LINEENDX)
y2=ObjectGeography(obj2,OBJ_GEO_LINEENDY)

else
if objectinfo(obj2,OBJ_INFO_TYPE)=OBJ_TYPE_PLINE and objectinfo(obj2,OBJ_INFO_NPOLYGONS)=1 then
iSec=objectinfo(obj2,OBJ_INFO_NPOLYGONS)
iNode=objectinfo(obj2,OBJ_INFO_NPOLYGONS+iSec)

´相交线起点
x1=ObjectNodeX(obj2,1,1)
y1=ObjectNodeY(obj2,1,1)

´相交线终点
x2=ObjectNodeX(obj2,iSec,iNode)
y2=ObjectNodeY(obj2,iSec,iNode)
else
goto NextStep
end if
end if

len1=distance(x,y,x1,y1,"m")
len2=distance(x,y,x2,y2,"m")
if (len1<=50*len2 and len2<=50*len1 ) and (len1>min and len2>min) then ´防止线段被分割地太小
select * from temp into temp1 where rowid=iRow2
set target on
objTemp=createCircle(x1,y1,len1)
insert into circle(obj) values (objTemp)
select * from circle into temp1
Objects split into target data ID=ID
set target off
delete from circle
end if

end if

NextStep:
fetch next from temp
wend
end if

fetch rec iRow1 from STable
fetch next from STable

wend
end sub

’＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
1。以上程序中，国省道.Tab是包含Line、PLine的线对象表，Circle.tab是临时表，用来创建截断所需的Circle区域对象，两个表必须是打开的。程序运行的结果就是将国省道.tab表中的所有相交线对象按相交点分割。
2。程序的思路按顺序访问表，每次取出一个对象obj1，查找与其相交的所有对象，顺序处理这些对象，比如obj2，找出Obj1和Obj2的交点Node；以Obj2的起点或终点为圆心，这个点与交点Node间的距离为半径作圆Circle，将Obj2设为目标，选择Circle，分割Obj2.
3。注意：在一次大循环中，Obj1并不分割，分割的只是与它相交的所有对象，Obj1本身要在以后的循环中由其它对象对其进行分割。
4。程序目前只能处理Line对象和只有一个Section的Pline对象，对于多个Section的PLine对象，望大侠们能给出方案。
5。请大家帮我复核一下这个现象：MapInfo中计算得到的两Line交点与真正截断发生的点之间存在误差，也就是根据Obj1将Obj2分隔成Obj3和obj4之后，Obj3或bj4与Obj1还是相交的，这一点放大之后看得很清楚。这是系统的误差还是计算的误差呢？

导入论坛引用链接收藏分享给好友推荐到圈子管理举报

TAG: 交点 integerDim sub floatDim objectDim

查看全部评论

用户菜单

我的栏目

标题搜索

日历

数据统计

访问量: 25715
日志数: 123
图片数: 13
文件数: 46
书签数: 1
建立时间: 2007-02-01
更新时间: 2008-10-25

RSS订阅