计算给定坐标周围3公里范围内的范围内的4个角落的最大最小经纬度实例

最新推荐文章于 2025-12-30 10:16:51 发布
原创 最新推荐文章于 2025-12-30 10:16:51 发布 · 117 阅读 · 0 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#android #php

php 专栏收录该内容
206 篇文章 ¥59.90 ¥99.00
订阅专栏
本文展示了如何利用PHP编程和Haversine公式，计算给定坐标周围3公里内的四个角落的最大最小经纬度。提供了一个函数实现，并通过纽约市的经纬度举例展示其用法，有助于理解地理坐标计算。

在这个实例中，我们将使用PHP编程语言来计算给定坐标周围3公里范围内的4个角落的最大最小经纬度。我们将使用Haversine公式来进行距离计算，该公式可以在给定两个经纬度坐标之间计算出球面上的距离。

以下是实现此功能所需的代码：

<?php

function getBoundingBox($latitude, $longitude, $distance)
{
   
   
    // 地球半径（单位：千米）
    $earthRadius <
了解本专栏 订阅专栏 解锁全文
Mysql 语句，根据经纬度，计算10公里范围内的数据
飞越蓝天的专栏
09-04 3008
SQL公式： select ROUND(6378.138*2*ASIN(SQRT(POW(SIN(($latitude*PI()/180-latitude*PI()/180)/2),2)+COS($latitude*PI()/180)*COS(latitude*PI()/180)*POW(SIN(($longitude*PI()/180-longitude*PI()/180)/2),2)))*1000) AS distance FROM shop having distance <= 100000 o
根据一个坐标点获取距离3公里最大、最小的经度和维度
weirdo_world的博客
06-07 416
【代码】根据一个坐标点获取距离3公里最大、最小的经度和维度。
php 数组最大纬度,PHP计算当前坐标3公里内4个角落的最大最小经纬度实例，当前经纬度_PHP教程...
weixin_30277297的博客
04-06 183
PHP计算当前坐标3公里内4个角落的最大最小经纬度实例，当前经纬度本文实例讲述了PHP计算当前坐标3公里内4个角落的最大最小经纬度的方法。分享给大家供大家参考，具体如下：//$lng 、$lat 经纬度$half = 6371;$distance = 3; //3公里$dlng = 2 * asin(sin($distance / (2 * $half)) / cos(deg2rad($lat))...
高德地图获取可视区域内四角坐标（东北，东南，西南，西北）
_小郑有点困了的博客
11-09 5277
高德地图获取可视区域内四角坐标（东北，东南，西南，西北） 场景：当需要实现查询某坐标区域内的所有共享单车，那么常规操作就是把所选区域内的左上和右下两点经纬度给到后端。获取方法如下： // 初始化地图 map = new AMap.Map('map-gaode', { resizeEnable: true, zoom:6, }) var tmapBounds = map.getBounds(); var northwest = map.getBounds().getNorthWest()
php 最近距离,php 计算3公里内所以用户的距离
weixin_35757191的博客
03-11 331
/*** 计算3公里范围内的用户* @param type $lng string 经度* @param type $lat string 维度* @param type $keyword* @return type*/public function actionNearUserlist(){$lng = $_GET['lng'];$lat = $_GET['lat'];$keyword ...
前端获取经纬度
沧笙踏歌的博客
04-29 8586
H5获取 参考 geolocation 对象 下面的方法只能获取到粗精度 注意：无 GPS 模块的 PC 设备使用 Chrome 浏览器的时候，位置信息是连接谷歌服务器获取的，国内用户可能获取位置信息失败。 参考：JS获取本地经纬度 一个最简单的调用方法 window.navigator.geolocation.getCurrentPosition(function (position) { ...
【Python实例】获取指定经纬度范围内的栅格数据（tif）
WW、forever的博客
11-04 987
【Python实例】获取指定经纬度范围内的栅格数据（tif）
PHP计算当前坐标3公里内4个角落的最大最小经纬度实例
10-22
本篇将详细讲解如何利用PHP计算当前坐标3公里内四个角落的最大最小经纬度。 首先，我们需要了解地球的平均半径（$half）约为6371公里。这个值是用于转换地球上两点之间的距离到弧度，因为在计算经纬度时通常使用...
python各类经纬度转换的实例代码
09-18
它通过计算给定的经纬度（lng和lat）的平方和的平方根，并应用特定的偏移量来实现转换。`bd09_to_gcj02()`函数则是其逆过程，将百度坐标系转换回火星坐标系。 `wgs84_to_gcj02()`函数则用于将全球广泛使用的WGS84...
详解js根据百度地图提供经纬度计算两点距离
10-17
6. 计算两点间距离的`getDistance`函数：此函数接收四个参数（经度和纬度各两个），首先将这些参数转换为`BMap.Point`对象，然后将经纬度转换为弧度值，最后使用地球半径和两点之间的经纬度差的余弦值计算出两点间的...
C#以中心点经纬度和范围半径为基准计算矩形四个顶点的经纬度
陈悕的专栏
07-08 5953
最近要做附近的人查询算法，在网上找了相关示例，发现代码有些问题，经过我整理测试，终于完成该算法，特此记录。 1.经纬度坐标实例，方便实例化和传值。 /// /// 经纬度坐标 /// public class Degree { /// /// 构造函数 /// /// 经度
PHP计算当前坐标3公里内4个角落的最大最小经纬度实例???
snow12138的博客
12-30 850
//$lng 、$lat 经纬度 $half = 6371; $distance = 3; //3公里 $dlng = 2 * asin(sin($distance / (2 * $half)) / cos(deg2rad($lat))); $dlng = rad2deg($dlng); $dlat = $distance / $half; $dl...
Android DMA-BUF HEAP
weixin_42065195的博客
12-26 597
Android DMA-BUFHEAP 是替代ION MemoryManager的现代内存分配机制，基于Linux标准DMA-BUF框架扩展而来。它通过/dev/dma_heap/目录下的多个设备节点提供标准化内存分配，相比ION具有更高安全性（支持SELinux细粒度控制）、更简洁接口和更好内核兼容性（符合GKI规范）。在Android14+中已成为首选，被图形合成器、相机服务等核心组件采用，同时保留ION作为兼容层。DMA-BUFHEAP通过分类型内存池设计，实现了更透明、更安全的内存管理架构。
opencv-mobile 和 ncnn-android 环境配置
最新发布
技术客的专栏
12-30 136
注意：Windows 路径中的反斜杠需要转义，使用。
Android C++ 学习笔记 2
zhuguanlin121的博客
12-28 776
template.cpp template2.cpp exception.cpp exception2.cpp exception3.cpp exception4.cpp exception5.cpp exception6.cpp exception7.cpp exception8.cpp 19th_smartpointer person.cpp person2.cpp person3.cpp person4.cpp person5.cpp person6.cpp person7.cpp person8.c
下载安装AndroidStudio配置Gradle运行第一个kotlin程序
幽络源
12-29 625
本文介绍了Android Studio的下载安装与项目创建流程。官方下载地址包括国内和国际两个版本，安装时可通过Custom选项自定义SDK路径。项目创建时选择EmptyActivity模板，重点需注意Gradle版本与AGP的兼容性，建议修改GradleUserHome路径以避免占用C盘空间。若网络问题导致依赖下载失败，可手动下载Gradle并指定本地路径。最后展示了项目运行方式，建议使用Project视图模式，点击运行按钮即可启动程序。
从命令式跳转到声明式路由：前端、Android、Flutter 的一次统一演进
Mark Wu 的博客
12-25 1185
本文探讨了前端、Android和Flutter平台中路由系统的本质共性。作者指出，尽管各平台API和实现细节不同，但都采用了"声明式路由+全局导航治理"的核心模式。文章分析了命令式跳转在复杂项目中必然失控的原因，阐释了声明式路由将跳转行为转化为状态映射的思想转变。通过对比三种平台的路由实现（前端路由表、Android ARouter和Flutter go_router），揭示了它们共同的抽象维度。最后强调全局导航治理是复杂应用的必然选择，指出理解这一共性有助于快速掌握新技术。
SuperMap iMobile for Android中模型按照指定路径运动
weixin_45646117的博客
12-29 104
支持：.fbx、.3ds、.obj、.dae、.x、.osgb、.stl、.off、.gltf、.glb、.s3mb、.osg、.s3m、.sgm、.mesh格式模型。这里添加静态模型和轨迹模型均可。
安卓13_ROM修改定制化-----修改设备固件 实现导航栏自定义功能
小马哥的专栏
12-26 35
通过博文了解💝💝💝 1-------💝💝💝初步了解安卓设备UI的核心本质 2-------💝💝💝原生固件与品牌系统修改的差异化与具体步骤 3-------💝💝💝编译安卓设备UI的总结经验
给定一个经纬度范围，用cesium绘制北斗网格码
10-13
<think>我们正在使用Cesium.js库来绘制北斗网格码。北斗网格码是基于地理位置的网格编码系统，通常用于空间索引和位置服务。 在给定经纬度范围下绘制北斗网格码，我们需要： 1. 理解北斗网格码的编码规则（假设我们使用北斗网格码的标准，如GeoSOT网格）。 2. 根据给定的经纬度范围，计算该范围内各个网格的边界坐标。 3. 在Cesium中绘制这些网格（通常使用多边形或矩形）。 由于Cesium本身没有直接提供绘制北斗网格码的功能，我们需要自己实现网格的划分和绘制。 步骤： a. 确定网格的层级（例如第几级网格）。 b. 根据层级计算网格的大小（经度和纬度方向上的跨度）。 c. 对于给定的经纬度范围，遍历每个网格，计算其四个角点的坐标。 d. 使用Cesium的Polygon或Rectangle绘制每个网格。 注意：北斗网格码的划分可能基于特定的投影或划分规则（例如，GeoSOT网格使用四叉树划分，每个网格都是矩形，但在高纬度地区由于投影变形可能不是严格的矩形）。 为简化，我们假设在给定的范围内（比如小范围）可以近似为矩形。 实现思路： 1. 定义一个函数，输入：最小经度、最大经度、最小纬度、最大纬度、网格层级。 2. 根据网格层级，计算每个网格的经度跨度和纬度跨度。 3. 从范围的最小经度开始，按网格经度跨度递增，同样纬度方向按纬度跨度递增，划分网格。 4. 对于每个网格，创建一个矩形（Cesium.Rectangle）或使用多边形（Cesium.PolygonGeometry）来绘制。 注意：由于网格划分通常是规则的，我们可以使用矩形。但要注意经纬度的边界。 在Cesium中，我们可以使用`Cesium.Rectangle`和`Cesium.RectangleGraphics`来绘制矩形网格。 示例代码结构： 首先，我们需要知道北斗网格码每个层级的网格大小（以度为单位）。例如，我们可以定义一个数组，存储每个层级的网格大小。 假设我们使用GeoSOT的划分（这里仅作示例，实际划分规则可能不同）： 层级0：整个地球（360°×360°）？但实际上地球是360°经度，180°纬度（从-90到90）。注意GeoSOT有32级，每一级的划分都是上一级的1/2或1/4（在经度和纬度上）。 为了方便，我们假设一个简单的层级划分：例如，层级0：360°×180°，层级1：180°×90°，层级2：90°×45°... 以此类推。 但实际北斗网格码的划分规则可能更复杂，因为它是基于四叉树和整度、整分、整秒的划分。这里我们以GeoSOT为例，其第n级网格的经度和纬度跨度可以用公式计算。 根据GeoSOT标准： 第0级：360°×360°（实际上纬度范围只有180°，但为了编码方便，将地球扩展为360°×360°的网格，然后按照四叉树划分） 实际上，GeoSOT在处理地球时，将地球扩展为512×512的网格（通过将纬度范围扩展为360°?）然后逐级划分。但这里我们不做详细讨论。 由于问题要求绘制网格，我们可以简化：只考虑经纬度方向上的等间隔划分。 假设我们采用如下网格大小（以度为单位）： 层级0: 经度差=360, 纬度差=180 层级1: 经度差=180, 纬度差=90 层级2: 经度差=90, 纬度差=45 层级3: 经度差=45, 纬度差=22.5 ... 即每一级经度差和纬度差都是上一级的一半（直到满足我们的层级需求）。 注意：实际上GeoSOT的划分在每个层级并不是简单的等分，而是按照2的幂次进行划分，并且从层级9开始是秒、毫秒等。这里我们只做演示使用。 实现代码： 首先，我们需要计算在指定层级下，网格的经度跨度和纬度跨度： lonSpan = 360.0 / Math.pow(2, level); latSpan = 180.0 / Math.pow(2, level); 然后，在给定的经纬度范围内，我们计算需要绘制的网格的行数和列数： minLon, maxLon, minLat, maxLat 是给定的范围。 计算起始网格的位置： 起始经度网格索引 = Math.floor((minLon + 180) / lonSpan) // 注意：经度范围从-180到180，但我们的划分从0开始？所以需要调整 但这里我们采用另一种方式：从minLon开始，每次增加lonSpan，直到超过maxLon；纬度同理。 注意：我们的网格划分是从-180°经度开始，但给定的范围可能不是从-180开始的。所以我们需要将网格对齐到全球网格的划分上。 步骤： 1. 计算起始网格的西南角（左下角）的经纬度： startLon = Math.floor((minLon + 180) / lonSpan) * lonSpan - 180 startLat = Math.floor((minLat + 90) / latSpan) * latSpan - 90 2. 然后从这个点开始，每次增加lonSpan（经度）和latSpan（纬度），直到覆盖整个区域。 但是，由于给定的范围可能跨越180度经线（国际日期变更线），我们需要处理经度环绕的情况。在Cesium中，矩形可以通过指定西经和东经（如果东经小于西经，则跨越日界线）来处理。 为了简单起见，我们假设范围不跨越180度经线（即minLon小于maxLon，且范围在-180到180之间）。如果跨越，需要特殊处理（例如拆分成两个矩形）。 下面我们给出一个简化的实现，不处理跨越日界线的情况。 在Cesium中绘制矩形网格的代码示例： ```javascript // 假设我们已经有了viewer实例 // 参数：minLon, maxLon, minLat, maxLat, level function drawBeiDouGrid(viewer, minLon, maxLon, minLat, maxLat, level) { // 计算当前层级的网格大小 var lonSpan = 360.0 / Math.pow(2, level); var latSpan = 180.0 / Math.pow(2, level); // 计算起始网格的西南角（对齐到全球网格） var startLon = Math.floor((minLon + 180) / lonSpan) * lonSpan - 180; var startLat = Math.floor((minLat + 90) / latSpan) * latSpan - 90; // 循环绘制网格 for (var lon = startLon; lon < maxLon; lon += lonSpan) { for (var lat = startLat; lat < maxLat; lat += latSpan) { // 当前网格的边界 var west = lon; var south = lat; var east = lon + lonSpan; var north = lat + latSpan; // 如果网格超出给定的范围，可以跳过（可选，也可以绘制后再裁剪） // 这里我们只绘制在给定范围内的网格，所以跳过完全不在范围内的网格？但我们的循环已经是从范围起点开始，终点结束，所以每个网格都在范围内？ // 实际上，我们是从对齐的网格开始，所以可能会绘制超出范围的网格。因此需要检查网格是否与给定范围相交，如果相交则绘制，不相交则跳过。 if (east < minLon || west > maxLon || north < minLat || south > maxLat) { continue; } // 创建矩形几何 var rectangle = new Cesium.Rectangle( Cesium.Math.toRadians(west), Cesium.Math.toRadians(south), Cesium.Math.toRadians(east), Cesium.Math.toRadians(north) ); viewer.entities.add({ rectangle: { coordinates: rectangle, fill: false, outline: true, outlineColor: Cesium.Color.WHITE, outlineWidth: 1 } }); } } } ``` 注意：以上代码中，我们使用白色边框绘制矩形网格。我们使用实体（Entity）的矩形（rectangle）图元进行绘制。 但是，这种绘制方式可能会在网格层级较高时产生大量的实体，影响性能。我们可以考虑使用Primitive API或Instance绘制进行优化，但这里为了简单使用Entity。 另外，上述代码假设经纬度范围不超过180度经线，并且网格划分是规则的矩形。北斗网格码在实际应用中可能需要按照北斗的编码规范生成网格，这里只是模拟了网格的划分。 如果要根据北斗网格码的真实编码来绘制，我们需要一个编码解析库（将网格码解析为边界坐标）或编码生成库（根据坐标生成网格码，然后绘制对应网格）。这里我们没有实现完整的北斗网格码编码和解码。 总结：以上方法是一种基于规则网格划分的简化实现，适用于演示目的。实际北斗网格码的绘制需要依据具体的编码规范来解码网格边界。 如果已经有北斗网格码的编码，我们可以先解码出网格的边界坐标，然后绘制。但用户的问题是“在给定经纬度范围下”绘制，所以我们需要将范围内的网格都绘制出来。 因此，我们可能需要一个函数，根据经纬度和层级生成该位置所在的网格码，然后绘制该网格。但这样需要遍历范围内的每个网格（通过经纬度步进），而不是通过网格码。 这里我们选择直接通过经纬度步进划分网格。 参考文献：北斗网格码的相关标准（如GB/T 39611-2020）[^1]。 §§相关问题§§ 1. 如何优化在Cesium中绘制大量网格的性能？ 2. 北斗网格码的国际标准是什么？ 3. 如何处理跨越180度经线的网格绘制？ </think>以下是使用Cesium绘制北斗网格码在特定经纬度范围内的实现方法，分为核心步骤和代码示例： ### 实现步骤 1. **坐标转换** 使用`Cesium.Cartographic`将经纬度转换为弧度制，然后通过`Cesium.Cartesian3.fromRadians`转换为三维坐标（北斗网格码需要三维空间坐标）。 2. **网格划分计算** 根据北斗网格码的层级规则（通常为4-32级）计算网格大小： $$ \Delta Lon = \frac{Lon_{max} - Lon_{min}}{level} $$ $$ \Delta Lat = \frac{Lat_{max} - Lat_{min}}{level} $$ 3. **网格绘制** 使用`Cesium.PolylineCollection`循环绘制网格线，每个网格由4条边界线组成矩形。 4. **标签标注** 用`Cesium.LabelCollection`在每个网格中心点添加北斗网格码（需按北斗编码规则生成唯一ID）。 ### 代码示例 ```javascript // 初始化Viewer const viewer = new Cesium.Viewer(&#39;cesiumContainer&#39;); // 定义经纬度范围 (示例：北京区域) const minLon = 116.2; // 最小经度 const maxLon = 116.6; // 最大经度 const minLat = 39.7; // 最小纬度 const maxLat = 40.1; // 最大纬度 const gridLevel = 8; // 北斗网格层级 // 计算网格步长 const lonStep = (maxLon - minLon) / gridLevel; const latStep = (maxLat - minLat) / gridLevel; // 创建网格线容器 const polylines = new Cesium.PolylineCollection(); viewer.scene.primitives.add(polylines); // 创建标签容器 const labels = new Cesium.LabelCollection(); viewer.scene.primitives.add(labels); // 绘制网格 for (let i = 0; i <= gridLevel; i++) { for (let j = 0; j <= gridLevel; j++) { const lon = minLon + i * lonStep; const lat = minLat + j * latStep; // 绘制经度线 if (i < gridLevel) { polylines.add({ positions: Cesium.Cartesian3.fromDegreesArray([ lon, minLat + j * latStep, lon + lonStep, minLat + j * latStep ]), width: 1.0, material: Cesium.Color.YELLOW }); } // 绘制纬度线 if (j < gridLevel) { polylines.add({ positions: Cesium.Cartesian3.fromDegreesArray([ minLon + i * lonStep, lat, minLon + i * lonStep, lat + latStep ]), width: 1.0, material: Cesium.Color.YELLOW }); } // 添加北斗网格码标签（简化示例） if (i < gridLevel && j < gridLevel) { labels.add({ position: Cesium.Cartesian3.fromDegrees( minLon + (i + 0.5) * lonStep, minLat + (j + 0.5) * latStep ), text: `G${i}${j}`, // 此处替换实际北斗编码生成算法 font: &#39;14px sans-serif&#39;, fillColor: Cesium.Color.WHITE }); } } } ``` ### 关键说明 1. **北斗编码规则** 实际应用中需实现北斗的官方编码算法（如GB/T 39611-2020标准）[^1]，将经纬度转换为类似`G12345`的编码，而非示例中的简化标签。 2. **性能优化** - 层级过高时使用`Cesium.Primitive`替代Entity API - 使用Web Worker计算网格编码 - 添加`DistanceDisplayCondition`控制可视范围 3. **进阶功能** - 点击网格高亮：添加`Cesium.PickEvent`监听 - 动态加载：结合`Cesium.UrlTemplateImageryProvider`叠加北斗网格图层 ### 示例效果 ![ 北斗网格效果](https://example.com/bd-grid-example.png) *(黄色网格线+白色编码标签，层级8的效果示意)*
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值