python算法. 在球面上取随机N个均匀的点(或者间距不小于某距离的点)
希望能在球上获得均匀分布, 或者 每两个点之间的间距不小于某个值的N个点的坐标.
点的数量不需要太大, 在100到200之间就够用了.
球的中心点就是坐标系原点.
有看到另外一个大牛写的.
https://www.oschina.net/code/...
但是传入100个点的时候, 相邻很近的点出现几率非常大. 导致在球面上的点上放东西的时候, 就叠在一起了.
求教, 有没有什么其他算法能实现.
球面上要实现均匀采样不难,用正态分布随机变量产生三维向量再单位化就可以了。
#include <iostream>
#include <fstream>
#include <random>
using namespace std;
int main()
{
std::default_random_engine gen;
std::normal_distribution<float> distrib(0.f, 1.f);
ofstream ofs("sphere.txt");
for (int i = 0; i < 1000; i++) {
float x = distrib(gen);
float y = distrib(gen);
float z = distrib(gen);
float r = sqrt(x*x + y*y + z*z);
ofs << x / r << ' ' << y / r << ' ' << z / r << endl;
}
return 0;
}
不过不知道满不满足相邻点之间的要求。如果要保证相邻点比较远,可以借鉴一下
jittering或者stratified sampling之类的思路。
Java版
import java.util.Random;
import java.io.*;
class SphericalSampling{
public static void main(String[] args){
Random rnd = new Random();
try{
PrintWriter writer = new PrintWriter("sphere.txt", "UTF-8");
for(int i = 0; i < 1000; i++){
double x = rnd.nextGaussian();
double y = rnd.nextGaussian();
double z = rnd.nextGaussian();
double r = Math.sqrt(x*x + y*y + z*z);
writer.println(x/r + " " + y/r + " " + z/r);
}
}catch (Exception e) {
e.printStackTrace(System.out);
}
}
}
另外,保存的sphere.txt可以用CloudCompare打开查看点云。
题主的意思是想让球面上的点间距尽量大,而均匀随机分布无法保证不出现距离任意小的两点,所以这个题与球面上的随机分布无关(标题太坑人)。
这个问题本质是求解球面点最小间距的最大值,各种优化工具包可以各显神通了。如果目标函数直接写成球面点最小间距的形式,函数稳定性很差,不容易求到最优解。这里将目标函数取为所有点间距平方的倒数和并求最小值:
$$\text{minimize:} \quad \sum_{i\lt{}j}\frac{1}{d^2(i,j)}$$
这样既突出了相邻点间距又保持函数相对平滑。
我用的是Mathematica提供的NMinimize
函数,点数比较多时需要很长计算。比如在我机器上算160个点需要四个小时。结果画图:
Spherical Codes
网上一搜才发现,原来这个问题还是蛮有来头的,叫做Tamme's problem,问题的解称为“spherical codes”。这里有一些计算好的结果。同时也知道,当点数比较多时寻找和证明最优解是很困难的。所以题主找到个还不错的次优解就可以啦。
球面均匀随机分布
啰嗦一句球面均匀分布吧。前面大神给出的神奇算法我百思不得其解,为啥用正态分布?后来从单位化上窥见了端倪:单位化其实是体分布到球面的投影。因为正态分布是球对称的,因此它投影到球面上就一定是均匀的了。也就是说,真正重要的是分布的球对称性,具体形式无所谓。比如圆内的面积均匀分布投影可以得到圆上的均匀分布:
玩蛇网文章,转载请注明出处和文章网址:https://www.iplaypy.com/wenda/wd13731.html
相关文章 Recommend
- • 2019年3月最新消息: Python 3.4.10 现已推出
- • [上海]招Python量化系统开发工程师
- • 优集品网络科技有限公司招Python中/高级工程师
- • 爱因互动科技发展有限公司招募Python开发攻城狮
- • mozio招聘Python/Django工程师
- • Kavout金融科技公司招Python研发工程师
- • Python数组逆向输出,编程练习题实例四十
- • Python数组插入排序,编程练习题实例三十九
- • Python矩阵for循环应用,编程练习题实例三十八
- • Python操作Redis数据库方面的问题
- • 请python高手帮我看看这段python代码中函数setter的
- • Python什么方法可以快速将两个队列变成字典
必知PYTHON教程 Must Know PYTHON Tutorials
- • python 解释器
- • python idle
- • python dir函数
- • python 数据类型
- • python type函数
- • python 字符串
- • python 整型数字
- • python 列表
- • python 元组
- • python 字典
- • python 集合
- • python 变量
- • python print
- • python 函数
- • python 类定义
- • python import
- • python help
- • python open
- • python 异常处理
- • python 注释
- • python continue
- • python pass
- • python return
- • python global
- • python assert
- • python if语句
- • python break
- • python for循环
- • python while循环
- • python else/elif
- • lambda匿名函数
必知PYTHON模块 Must Know PYTHON Modules
- • os 模块
- • sys 模块
- • re 正则表达式
- • time 日期时间
- • pickle 持久化
- • random 随机
- • csv 模块
- • logging 日志
- • socket网络通信
- • json模块
- • urlparse 解析URL
- • urllib 模块
- • urllib2 模块
- • robotparser 解析
- • Cookie 模块
- • smtplib 邮件
- • Base64 编码
- • xmlrpclib客户端
- • string 文本
- • Queue 线程安全
- • math数学计算
- • linecache缓存
- • threading多线程
- • sqlite3数据库
- • gzip压缩解压
最新内容 NEWS
- • Python程序员解决棘手问题的常用库
- • 求助关于restfull api接口几个问题
- • qiniu pythonsdk提示ImportError错误求解
- • 问一个关于Hadoop Python中读写文件统计分析
- • 求问str()同__str__原理上有什么不同,分别在
- • 大神帮忙看下20行的python代码,文件io和数
- • python 爬虫爬wiki 报错 [Errno 65] No route to
- • python续点上传问题None bad token...
- • python3环境下文本中超链接出错,要如何修
- • Python环境保存操作思路问题求助
图文精华 RECOMMEND
-
Python程序员解决棘手问题的常用库
-
求问str()同__str__原理上有什么不同
-
scrapy框架里面用link extractor怎么能
-
python {}.fromkeys创建字典append添加操
-
python3 类型Type str doesn't support th
-
python里面为什么系统的时区是东八
热点文章 HOT
- 学习Python有什么好的书籍推荐?
- Python匿名函数 Lambda表达式作用
- Python与Java、C、Ruby、PHP等编程语言有什么
- Python 正则中文网页字符串提取问题
- 如何为实时性应用存取经纬度?django my
- 想用python做个客户端,在二维码登录这个地
- 有让IDE可识别Python函数参数类型的方法吗
- Python字符串转换成列表正则疑问