mapinfo画扇形

mapinfo 中绘制基站扇形的算法讨论和工具共享 资料名称 mapinfo 中绘制基站扇形的算法讨论和工具共享cbsc 移动通信论坛拥有 30 万通信专业人员 超过 50 万份 GSM 3G 等通信技术资料 是国内领先专注于通信技术和通信人生活的社区 R t db w X C 9m0n l p0I i4 Y 方法方法 3 3 采用大地主题解算算法来实现 这种算法需要很多专业知识背景 我 就是猛学了一段时间的测量专业的教材才理解的 大地主题解算有很多算法 我采用了白赛尔算法实现 具体算法我已经放到了上帝之眼上 大家有兴趣可 以自行看看 googleearth 基站扇区绘制工具 YZL 就是按照这个算法弄的 方法方法 4 4 mapinfo 中 我采用 mapbasic 编程 将经纬度转换成直角坐标系 然后再在直接坐标系下 SETE COORDSYS nonearth 语句 画线和画圆弧 然后 再切换到经纬度坐标系下 这个方案比较简单 但是从 earth 坐标系到 nonearth 坐标系的转化 暂时没有 找到好的方法 而且误差也很大 我认为 不可取 方法方法 5 5 mapinfo 中 我采用 mapbasic 编程 想了另外一种方法 充分利用 MAPINFO 自身功能 完全抛开直角坐标到经纬度的转换 实现起来很方便而且 效果很好 国内领先的通信技术论坛 R f m3X k F 下面简单介绍一下方法 5 主要思路 采用类似于我们高中几何中常用到的辅助线的方法 直接再经纬度 坐标系下处理 MSCBSC 移动通信论坛 c7D R ZA8 Y g4d 这种做法缺点也有 由于 mapinfo 的处理精度问题 所以半径小于 10m 或者 大于 500km 均还是有一定误差 但是这已经不影响使用了 MSCBSC 移动通信论坛T 主要问题是 mapinfo 中 Create Arc 语句没有基于两个端点和圆心来绘制的方 式 而是采用了四个点的坐标方式 mscbsc 移动通信论坛拥有 30 万通信专业人员 超过 50 万份 GSM 3G 等通信技术资料 是国内领先专注于通信技术和通信人生活的社区 7i z3d4G0Q4J5G6f 采用辅助圆的方式取得需要绘制圆弧的四个点的坐标 国内领先的通信技术论坛 X7 2r7O E s e1N 通过这四个点绘制弧形 然后取得该弧形的两个端点坐标 通过端点坐标和圆心绘制直线 E 4J2K y7D q f4H 最后合并直线和圆弧就得到了标准的扇形 效果不错移动通信 通信工程师的家园 通信人才 求职招聘 网络优化 通信工程 出差住宿 通信企业黑名单 W Q zP4 9i miny ObjectGeography circle obj obj geo miny MSCBSC 移动通信论坛 6o6v6 2r e i8l maxx ObjectGeography circle obj obj geo maxx 移动通信 通信工程师的家园 通信人才 求职招聘 网络优化 通信工程 出差住宿 通信企业黑名单 5c1M V m7S S6G9s maxy ObjectGeography circle obj obj geo maxy 这四行取得绘制圆 弧的关键点坐标 if 方向角 270 and 半功率角360 then 此处的方向角是经 过处理后的方向角 Create Arc into Variable arc obj 移动通信 通信工程师的家园 通信人才 求职招聘 网络优化 通信工程 出差住宿 通信企业黑名单 q s7D5y9r Z minx miny maxx maxy start angle end angle MSCBSC 移动通信论坛 H4g V N0w T x y xx1 yy1 T s j v2s Create Line into Variable line2 obj mscbsc 移动通信论坛拥有 30 万通信专业人员 超过 50 万份 GSM 3G 等通信技术资料 是国内领先专注于通信技术和通信人生活的社区 4x5o l2O5I F x y xx2 yy2 这里是绘制两条直线 obj1 ConvertToRegion Combine arc obj Combine line1 obj line2 obj 合并移动通信 通信 工程师的家园 通信人才 求职招聘 网络优化 通信工程 出差住宿 通信企业黑名单 I q J M8E u2Q D e W else obj1 circle obj 全向站的处理 直接画圆 end if mycreateregion obj1 返回扇形 end functionMSCBSC 移动通信论坛 I t6b t G0X1s 2t g V6N Z c U 求圆弧端点的程序 mscbsc 移动通信论坛拥有 30 万通信专业人员 超过 50 万份 GSM 3G 等通信技术资料 是国内领先专注于通信技术和通信人生活的社区 e w A Z h1 b0p3k getArcBEGX 0 end if end function function getArcBEGY byval obj as object as float if objectinfo obj OBJ INFO TYPE 1 C s s7 3q 3g6B Q8p R b dim d as object 国内领先的通信技术论坛 4x w V j M 1W C7V S k end 3 cb D7p8f o E p1s getArcENDX ObjectNodeX d 1 objectinfo d OBJ INFO NPNTS 0 k G B u4j getArcENDX 0 end if end function function getArcENDY byval obj as object as float if objectinfo obj OBJ INFO TYPE 1 then dim d as object 国内领先的通信技术论坛 v f4Q1b q V d ConvertToPline obj z4o F 0C6T0H W9k r5h N getArcENDY ObjectNodeY d 1 objectinfo d OBJ INFO NPNTS mscbsc 移动通信论坛拥有 30 万通信专业人员 超过 50 万份 GSM 3G 等通信技术资料 是国内领先专注于通信技术和通信人生活的社区 l O 6P Z 4v Z1P q Q a E 3v u getArcENDY 0移动通信 通信工程师的家园 通信人才 求职招聘 网络优化 通信工程 出差住宿 通信企业黑名单 q Z0Q V9l i j9u F F end if end function 我基于方法 5 做 了一个简单的 mapinfo 基站扇区图形生成工具 大家试试看 大家使用超级用户密码登录就可以了 密码为 leizi y 这个工具实现了点 线 面图 形的自定义等等功能 但没有做 excel TXT 导入 mapinfo 这样基本的 功能 我认为这个 mapinfo 自身已经能够足够处理了 欢迎大家讨论 leizi y qq 31933208 大地主题解算大地主题解算 正算正算 代码 根据经纬度和方向角以及代码 根据经纬度和方向角以及 距离计算另外一点坐标距离计算另外一点坐标 新建模块 拷贝下面的大地主题 正算 代码 调用方法示例 起点经度 116 235 度 终点纬度 37 435 度 方向角 50 度 长度 500 米 终点经纬度 经度 纬度 Computation 37 435 116 235 50 500 Option Explicit Const pi 3 1415926535898 Private a b c alpha e e2 w V As Double Private B1 L1 B2 L2 As Double Private s As Double Private A1 A2 As Double Private Sub getellipseparameter a 6378245 b 6356752 3142 c a 2 b alpha a b a e Sqr a 2 b 2 a e2 Sqr a 2 b 2 b End Sub Private Function computerw w Sqr 1 e 2 Sin B1 2 V w a b End Function Function Computation STARTLAT STARTLONG ANGLE1 DISTANCE As Double As String 正算 Dim sinu1 cosu1 sinA0 cotq1 sin2q1 cos2q1 cos2A0 As Double Dim k2 q0 sin2q1q0 cos2q1q0 As Double Dim q As Double Dim theta As Double Dim aa BB cc EE22 AAlpha BBeta As Double Dim sinu2 lamuda As Double Dim e1 As Double Dim W1 As Double B1 STARTLAT L1 STARTLONG A1 ANGLE1 s DISTANCE Call getellipseparameter If B1 0 Then If A1 90 Then A2 270 B2 0 L2 L1 s a 180 pi End If If A1 270 Then A2 90 B2 0 L2 L1 s a 180 pi End If Exit Function End If B1 rad B1 L1 rad L1 A1 rad A1 Call computerw e1 e W1 w sinu1 Sin B1 Sqr 1 e1 e1 W1 cosu1 Cos B1 W1 sinA0 cosu1 Sin A1 cotq1 cosu1 Cos A1 sin2q1 2 cotq1 cotq1 2 1 cos2q1 cotq1 2 1 cotq1 2 1 cos2A0 1 sinA0 2 e2 Sqr a 2 b 2 b k2 e2 e2 cos2A0 aa b 1 k2 4 3 k2 k2 64 5 k2 k2 k2 256 BB b k2 8 k2 k2 32 15 k2 k2 k2 1024 cc b k2 k2 128 3 k2 k2 k2 512 e2 e1 e1 AAlpha e2 2 e2 e2 8 e2 e2 e2 16 e2 e2 16 e2 e2 e2 16 cos2A0 3 e2 e2 e2 128 cos2A0 cos2A0 BBeta e2 e2 32 e2 e2 e2 32 cos2A0 e2 e2 e2 64 cos2A0 cos2A0 q0 s BB cc cos2q1 sin2q1 aa sin2q1q0 sin2q1 Cos 2 q0 cos2q1 Sin 2 q0 cos2q1q0 cos2q1 Cos 2 q0 sin2q1 Sin 2 q0 q q0 BB 5 cc cos2q1q0 sin2q1q0 aa theta AAlpha q BBeta sin2q1q0 sin2q1 sinA0 theta AAlpha q BBeta sin2q1q0 sin2q1 sinA0 sinu2 sinu1 Cos q cosu1 Cos A1 Sin q B2 Atn sinu2 Sqr 1 e1 e1 Sqr 1 sinu2 sinu2 180 pi lamuda Atn Sin A1 Sin q cosu1 Cos q sinu1 Sin q Cos A1 180 pi If Sin A1 0 Then If Sin A1 Sin q cosu1 Cos q sinu1 Sin q Cos A1 0 Then lamuda Abs lamuda Else lamuda 180 Abs lamuda End If Else If Sin A1 Sin q cosu1 Cos q sinu1 Sin q Cos A1 0 Then lamuda Abs lamuda 180 Else lamuda Abs lamuda End If End If L2 L1 180 pi lamuda theta 180 pi A2 Atn cosu1 Sin A1 cosu1 Cos q Cos A1 sinu1 Sin q 180 pi If Sin A1 0 Then If cosu1 Sin A1 cosu1 Cos q Cos A1 sinu1 Sin q 0 Then A2 180 Abs A2