(答案满意加到100分)涉及到 GIS 经纬度坐标的几个问题,帮帮忙解决一下,偶不会已经好几天没有弄出来了

Mister 2013-03-20 04:56:56
首先说明一下
1)发计算求两点距离的朋友,就麻烦别发了,我有,网上到处也有。
2)我需要的是代码实现,不要发计算公式,偶看不懂。

问题一:
已知线段(折线) A 的经纬度数组,和某个目标点 B 的经纬度,求 B 到 A 的最短距离。

问题二:
已知多边形 A 的经纬度数组,和某个线段 B 的经纬度集合,判断 B 和 A 之间的关系(1.线段在多边形之内? 2.线段不在多边形之内? 3.线段与多边形相交?)

问题三:
已知多边形 A 的经纬度数组,和某个目标点 B 的经纬度,求 B 到 A 的最短距离。



注1:计算必须考虑传入参数是 GIS 经纬度坐标,而不是普通的平面坐标。
注2:需要实现源代码能是 C# 最好,Java、VB、C++ 也行。
...全文
284 11 打赏 收藏 转发到动态 举报
写回复
用AI写文章
11 条回复
切换为时间正序
请发表友善的回复…
发表回复
Mister 2013-04-03
  • 打赏
  • 举报
回复
自顶一下,希望达人能看见。
Mister 2013-03-25
  • 打赏
  • 举报
回复
还有 liuchaolin 的代码,是计算什么东东?说明一下啊。
Mister 2013-03-25
  • 打赏
  • 举报
回复
引用 5 楼 unicorn_dsx 的回复:
gis开发,除了知道经纬度,你需要有对应的地图才可以计算最短距离吧. 因为两点只能有障碍物.例如,建筑物,河流,高山等等. 有地图才能知道,有那些通道,才会有最短距离. 而且早些年,arcgis 超图,灵图的开发平台,都有实现最短路径的方法.你用的是什么开发平台.到他们公司论坛问一下.
你说的无需考虑,我没有说要实现逻辑上的最短距离。
md5e 2013-03-22
  • 打赏
  • 举报
回复
引用 7 楼 Mister 的回复:
引用 6 楼 liuchaolin 的回复: C# code?1234567891011121314 protected double distHaversine(double lat1, double lon1, double lat2, double lon2) { Int16 R = 6371; double dLat = toRad(lat2 - ……

    protected string destPoint(double lat, double lon, double brng, double d) 
    {
        Int16 R = 6371; // earth's mean radius in km
        double lat1 = toRad(lat);
        double lon1 = toRad(lon);
        brng = toRad(brng);

        double lat2 = Math.Asin(Math.Sin(lat1) * Math.Cos(d / R) +
                              Math.Cos(lat1) * Math.Sin(d / R) * Math.Cos(brng));
        double lon2 = lon1 + Math.Atan2(Math.Sin(brng) * Math.Sin(d / R) * Math.Cos(lat1),
                                     Math.Cos(d / R) - Math.Sin(lat1) * Math.Sin(lat2));
        lon2 = (lon2 + Math.PI) % (2 * Math.PI) - Math.PI;  // normalise to -180...+180

        //if (isNaN(lat2) || isNaN(lon2)) return null;
        return "<lat>" + toDeg(lat2) + "</lat><lon>" + toDeg(lon2) + "</lon>";
    }
    protected double toRad(double Rad)
    {
        return Rad * Math.PI / 180;
    }

    protected double toDeg(double Deg)
    {  // convert radians to degrees (signed)
        return Deg * 180 / Math.PI;
    }
Mister 2013-03-22
  • 打赏
  • 举报
回复
引用 6 楼 liuchaolin 的回复:
C# code?1234567891011121314 protected double distHaversine(double lat1, double lon1, double lat2, double lon2) { Int16 R = 6371; double dLat = toRad(lat2 - lat1) ; ……
这段代码是计算什么的?toRad 方法定义呢?
a569673493 2013-03-21
  • 打赏
  • 举报
回复
学GIS的飘过,接分
md5e 2013-03-21
  • 打赏
  • 举报
回复

    protected double distHaversine(double lat1, double lon1, double lat2, double lon2)
    {
        Int16 R = 6371;
        double dLat = toRad(lat2 - lat1) ;
        double dLon = toRad(lon2 - lon1) ;
        lat1 = toRad(lat1);
        lat2 = toRad(lat2);
        double a = Math.Sin(dLat / 2) * Math.Sin(dLat / 2) +
        Math.Cos(lat1) * Math.Cos(lat2) *
        Math.Sin(dLon / 2) * Math.Sin(dLon / 2);
        double c = 2 * Math.Atan2(Math.Sqrt(a), Math.Sqrt(1 - a));
        double d = R * c;
        return d;
    }
一个点一个点去比较吧 Response.Write("<br>" +distHaversine(24.316029, 109.366053, 40.7118740691058, -74.0005159378052));
xixihaha_2011_098 2013-03-21
  • 打赏
  • 举报
回复
gis开发,除了知道经纬度,你需要有对应的地图才可以计算最短距离吧. 因为两点只能有障碍物.例如,建筑物,河流,高山等等. 有地图才能知道,有那些通道,才会有最短距离. 而且早些年,arcgis 超图,灵图的开发平台,都有实现最短路径的方法.你用的是什么开发平台.到他们公司论坛问一下.
木头小菜 2013-03-21
  • 打赏
  • 举报
回复
貌似arcgis 有这个功能
yanasdf789 2013-03-20
  • 打赏
  • 举报
回复
arcgis 二次开发 组件 自己调用 自己写算法 要琢磨一阵子了
言多必失 2013-03-20
  • 打赏
  • 举报
回复
还是帮你顶下吧。
 OpenGL-自主高性能三维GIS平台架构与实现/第二季:实现三维GIS球体+ 高程数据章节名称DEM基础1DEM基础知识1.介绍基本的DEM知识2.什么是DEM,作用是什么2DEM数据1.如何获取/ 传统测量/激光扫描/无人机测量/ 点云数据/ 倾斜摄影2.如何使用/局部小规模(栅格数据,图片/tif),3. 组织方式4. 根据使用目的不同,介绍多种优化方法3DEM图层的实现原理14DEM数据结构定义struct  V3U3N4顶点数据的生成和计算WGS84投影计算5wgs84 投影球体被切成一个个小圆弧,一共60个投影带,别为01,02.........60WGS的最新版本为WGS 84(也称作WGS 1984、EPSG:4326),1984年定义、最后修订于2004年。接口定义坐标转换Wgs84 数据加载6瓦片编号计算生成算法1. 经纬度到大地坐标的转换2.大地坐标经纬度坐标转换3. 根据经纬度获取瓦片编号框架重构7智能指针重构框架1. 基类定义(所有的类继承自基类),基类派生自 std::enbale_shared_from_this2. 实现智能指针的动态转换接口3. 实现向下转换4. 已有的类实现全部使用智能指针重构5. 任务系统(多线程加载任务)8引入图层(Layer)1. 介绍图层的概念以及重要性2. 图层类实现3. 修改框架(使用图层的方式重构框架)9Layer-bug排查(绘制过程中出现错位,发)1. 框架重构后遇到问题(绘制结果错误)2. 瓦片索引方式发生变化,多线程中引起内存问题3. 修改索引方式,解决绘制发错误问题10引入数据源(TileSource)1. 数据源的作用与设计目的2. 当前存在的问题,数据调度中存在问题3. 数据源(TileSource)类实现11数据格式管理(FormatMgr)1. 数据格式管理(FormatMgr) 提出的目的,需要解决问题2. CELLFormat基类接口抽象3. 实现几个标准格式类4. 修改框架流程,使用FormatMgr重构流程5. 扩展支持,后续支持任务格式数据加入系统12Task(任务)优化1. 任务中低耦合数据结构,目的是让Task更加的通用2. 修改任务读取代码与任务处理代码,完善处理流程DEM高程13DEM-数字高程定义1. 什么是数字化高程数据2. 当下GIS系统中有哪些常见的高程格式3. 课程体体系中使用的哪种格式4. 高程类定义以及实现,并加入到FormatMgr 管理系统中14高程瓦片数据读取1. 介绍GIS系统相关的工具(在数据转换)数据生成方面可以解决大量时间2. 自定义高程瓦片格式说明3. 自定义高程格式文件解析,并以智能对象的方式引入到系统中4. 完善框架代码,适配高程数据15高程瓦片文件的读取1. 实现基本的读取算法2. 增加格式化组件,并加入到系统中3. 配置高程图层以及高程数据源,并加载数据,验证数据正确性16瓦片数据结构重构1.顶点生成2.UV坐标计算3.面数据生成17DEM重构绘制流程1. 修改绘制数据结构,去除无用字段2. 增加Mesh类,实现光栅数据转换成三角面数据,计算UV数据,提炼接口3. 修改系统调度,实现顶点数据,UV数据,以及面数据的生成与更新4. 按需更新数据,而不是每一帧更新18DEM-数据精度问题(CPU)1. 因为瓦片数据使用大地坐标作为系统输入,造成瓦片坐标很大,单浮点数据精度不够2. 使用局部坐标的方式解决单浮点精度问题3. 调整相机参数,解决投影矩阵数据计算深度精度问题4. 修改绘制shader 实现对瓦片数据的绘制19DEM-数据精度问题(LogDepth)1. 使用对数深度(log depth )算法在GPU中 计算解决单浮点经纬计算问题2. 修改shader ,增加对(logDepth)算法支持3. 修改C++端代码,实现对shader数据的输入20DEM-数据结构优化1.当下使用CPU端数据通过接口的方式传递给GPU,速度慢2. 使用Instance 方式降低Vertex Buffer 的大小,优化渲染系统21DEM-GPU缓冲区优化1. 使用Vertex Buffer Object / Index Buffer Object  / Instance  方式优化渲染系统2. 修改绘制接口,使用DrawElementsInstanceBaseInstance方式提升系统性能内存池与对象池22瓦片生成优化/对象池1. 相机移动过程中会频繁的建立与释放瓦片,对CPU有较大的消耗2. 引入内存池,避免频繁的内存申请与释放,降低CPU时间3. 改造智能指针对象,对象释放通知到内存管理,回收对象内存23改造任务系统支持对象池1. 任务系统是一个公用模块,被多个模块使用,避免频繁的内存操作,引起的内存碎片2. 实现对象池,并应用到任务模块法线计算24法线计算1. 修改现有顶点结构,增加法线支持2. 修改shader,增加法线顶点输入,使用平行光光照模型3. 修改绘制流程,支持光照计算,使用探照灯作为光源输入25顶点法线计算/共享法线计算1. 增加数据结构保存顶点数据被多个面共享的次数2. 计算面法线,并累加到顶点法线中3. 根据顶点被面共享的次数做平均法线计算4. 修改流程,按需更新法线数据26法线数据压缩1. 法线数据使用3 * float 数据存储,大大的增加了系统的数据2. 实现算法,将3 * float 数据压缩成4字节数据3. 改造绘制代码,支持压缩数据输入27GPU中计算产生法线数据(去掉CPU中计算)1. 引擎支持 Geometry Shader 阶段2. 编写 Geometry Shader,实现法线计算系统功能优化28重构CPU拾取流程1. 当下的拾取流程,只支撑二维数据拾取,无法准群的拾取三维数据2. Terrain中增加拾取接口,输入射线,输出拾取到顶点数据29绘制拾取结果1. 增加一个绘制点的方法,实现绘制代码2. 修改shader,增加logdepth3. 调试代码,花费了很多时间排查错误,最总排查到是因为uniform参数笔误写错造成。30任务系统完善,避免任务队列无线膨胀1. 任务系统中,没有限制队列的大小,生产者的能力远大于消费者的能力,造成任务队列膨胀2. 处理办法,限制生产者的生产能力,而不是限制任务队列大小(这种方式会造成业务逻辑异常复杂)3. 使用sleep休眠方式(这种方式是严重错误的)31如何避免瓦片数据抖动1. 产生瓦片抖动的原因 ? 裂算法与回退算法中间没有过度2. 引入过度流程,避免内存抖动,参数因子是一个重要的数据,需要谨慎使用3. 有必要结合瓦片自身数据动态计算参数因子32瓦片数据管理-fepk文件格式支持-全球数据加载1. 支持fepk文件格式,增加fepk读取组件,适配fepk文件2. fepk管理数据方式:一般情况选择全球前10级别作为基础级别,因数据量不大(1G)左右,后续以8级作为基础级别,全球19级别数据被划为 2^8 * 2^7(512 * 256)个块。每个块中包含了256 * 256 张小瓦片33fepk高程数据读取 34高程裂处理当瓦片没有高程数据,那么子节点以及其他后代节点该如何共享父节点的数据35lesson-734-高程瓦片裂处理(2)-算法实现高程数据裂算法实现实现对高程数据的切,并对特殊数据进行处理36高程瓦片裂处理(3)-问题排查 37高程瓦片裂处理(4)-(后代节点更新问题)当一个瓦片高程数据更新后,他的儿子节点,孙子节点...该如何处理?38瓦片视锥裁剪错误高程数据更新后,没有技术计算瓦片包围盒信息,造成包围盒错误,进而引视锥计算错误39http支持1.引入三方库 Libcurl2.http类封装,支持http读取数据40fepk.server使用 生成三维地球41改造四叉树-统一使用经纬度输入42地形网络生成算法重构 43引入球体坐标系 44使用球体坐标改造瓦片 45多图层(加载标签数据) 课时截图:镜头拉近后,显示细节数据加载矢量SHP国界线数据:加载矢量三维白膜数据截图高程数据加载点云数据 加载倾斜摄影数据 

110,538

社区成员

发帖
与我相关
我的任务
社区描述
.NET技术 C#
社区管理员
  • C#
  • Web++
  • by_封爱
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告

让您成为最强悍的C#开发者

试试用AI创作助手写篇文章吧