mysql lat long_看看lat / long是否使用mysql落在多边形内

我已经创建了下表

CREATE TABLE geom (g GEOMETRY);

并插入了很多行,例如:

INSERT INTO geom (g)

VALUES(PolygonFromText('POLYGON((

9.190586853 45.464518970,

9.190602686 45.463993916,

9.191572471 45.464001929,

9.191613325 45.463884676,

9.192136130 45.463880767,

9.192111509 45.464095594,

9.192427961 45.464117804,

9.192417811 45.464112862,

9.192509035 45.464225851,

9.192493139 45.464371079,

9.192448471 45.464439002,

9.192387444 45.464477861,

9.192051402 45.464483037,

9.192012814 45.464643592,

9.191640825 45.464647090,

9.191622331 45.464506215,

9.190586853 45.464518970))')

);

现在我想搜索所有的数据,并返回一个lat / long我与任何多边形下降的条目.

这怎么可以使用mysql?或者是否知道任何链接会指向正确的方向?

这段代码的目的是将给定的经纬度划分到以 `accuracy` 为大小的栅格中,然后计算出栅格的中心点经纬度。以下是一些代码优化的建议: 1. 可以将常用的参数值保存在变量中,以便多次使用。例如,可以将栅格划分范围和栅格大小保存在变量中,然后在计算经纬度增量和栅格编号时使用这些变量。 2. 可以使用更具描述性的变量名,以便更容易理解代码的功能。例如,可以将 `lon_col` 改为 `longitude_index`,将 `hb_lon` 改为 `center_longitude`。 3. 可以使用更简洁的语法来计算经纬度增量和栅格编号。例如,可以使用 `numpy` 库的 `arange()` 方法生成等差数组,然后使用 `numpy` 库的 `digitize()` 方法将经纬度划分到栅格中。 4. 可以使用文档字符串来解释代码的功能和参数,以便更易读和维护。 下面是优化后的代码示例: ``` import math import numpy as np # 设置常用参数 min_longitude = 113.75194 max_longitude = 114.624187 min_latitude = 22.447837 max_latitude = 22.864748 accuracy = 500 def get_grid_center_point(lon, lat, accuracy): """ 将给定的经纬度划分到以accuracy为大小的栅格中,然后计算出栅格的中心点经纬度。 :param lon: 经度 :param lat: 纬度 :param accuracy: 栅格大小(单位m) :return: 栅格中心点的经纬度 """ # 计算栅格的经纬度增加量大小 delta_lon = accuracy * 360 / (2 * math.pi * 6371004 * math.cos((min_latitude + max_latitude) * math.pi / 360)) delta_lat = accuracy * 360 / (2 * math.pi * 6371004) # 计算栅格的经纬度编号 longitude_index = np.arange(min_longitude, max_longitude, delta_lon) latitude_index = np.arange(min_latitude, max_latitude, delta_lat) lon_col = np.digitize(lon, longitude_index) - 1 lat_col = np.digitize(lat, latitude_index) - 1 # 计算栅格的中心点经纬度 half_delta_lon = delta_lon / 2 half_delta_lat = delta_lat / 2 center_longitude = lon_col * delta_lon + min_longitude - half_delta_lon center_latitude = lat_col * delta_lat + min_latitude - half_delta_lat return center_longitude, center_latitude # 调用函数计算栅格中心点经纬度 test_lon = 114 test_lat = 22.5 center_lon, center_lat = get_grid_center_point(test_lon, test_lat, accuracy) # 打印结果 print(center_lon, center_lat) ``` 使用 `numpy` 库的 `arange()` 方法生成等差数组,然后使用 `numpy` 库的 `digitize()` 方法将经纬度划分到栅格中,从而使代码更简洁和易读。同时,添加了文档字符串来解释函数的功能和参数。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值