请教一个关于坐标点生成算法

[chengxiao]

现在已知 4 个坐标经纬度围成了一个矩形，如何在这个已知矩形中生成 500 个等距的经纬度坐标点呢？

[xuanbg]

除非矩形的长宽比是固定的几个值，否则不可能存在解

[no1xsyzy]

等距是怎么个定义法？是矩形排列还是六边形排列？平面距离还是空间距离？是必须包含边界还是不需要？
你什么都不说别人怎么回答？楼上楼下除了给出 8 个答案的其他都是在瞎猜。

[bilibilifi]

经纬度的话涉及到不同国家和地区的参考系，是一个很麻烦的问题

[SeanChense]

两个点是线段，三个点是等边三角形，四个点是啥

[bilibilifi]

同一个参考系的话可以利用库就 ok，不同的话要计算投影

[necomancer]

不知道你说的等距是啥概念，如果想求的是均分球面，等价于求 \theta\in(0, \theta), \phi\in(0, \phi) 的 500 个单位向量中，\sum_ij 1/|r_i - r_j| 的最小值，是一个有区域限制的汤普森问题。

初始化生成 500 个单位向量，在 0<theta<经度, 0<phi<维度 满足你的矩形限制条件下，最小化能量函数 U(r1,r2,...r500)=\sum_{i=1}^500 1/|r_i-r_j|，mathematica 有比较简单的代码。要自己撸需要用梯度下降法，或者参考 scipy.optimize 里的 minimize 函数。

[necomancer]

如果是均分经度、维度，那么 1 楼说了，需要特定的长宽比。上面均分球面算法保证弧线、弦距离相等。

[necomancer]

如果你不需要那么精确的结果，也可以考虑直接用密度，也就是如果你的矩形所占球台面积是 10，那么密度就是 50，如果球面积是 200，那么总点数就是 200x50=10000，Mathematica 有个 SpherePoints 函数，可以生成很均匀的点，大量点精确度会变高而且速度很快，详见 Fejes Tóth 问题。速度很快。你可以生成好以后用矩形去框一下，密度法可能有一定的涨落，如果要求结果不用那么精确的话（ 500+-10 ），这个方法是最高效的。

[opengps]

这种需求是网格处理方面的算法吗？

[chengxiao]

@opengps @necomancer @SeanChense @bilibilifi @xuanbg
sorry 各位 题目已补充 麻烦大家帮我看一下

[necomancer]

球面均匀点的算法我在之前已经提供了，看起来你并不是绝对需要 500 个点，你可以考虑尝试 Mathematica 的 SpherePoints 函数，等距性很好。我没能理解你的 dbscan 是怎么做的，空洞是什么意思？从图 1 来看感觉所有点都是连在一起的。你是直接对 (经，纬，1 ）这样相当于对经纬度分类还是每个小区有不同的权重？ dbscan 基于密度，加上权重也许你的问题就能解决了。另外，你可以试试 hdbscan 算法。

[chengxiao]

@necomancer 权重是一样的 直接队经纬度进行聚类 图二上能看到上海最密的区域是个空洞 没有 1 个点，但从图 1 可以看出那里是最密集的区域

[chengxiao]

@chengxiao 感谢 信息量有点大 我先消化消化

[chengxiao]

@necomancer 感谢 信息量有点大 我先消化消化....激动到发错了

[Raymon111111]

geohash 块的中心就行