我有一个 (x,y) 坐标列表，它们代表线条骨架。 该列表是直接从二值图像中获得的:
import numpy as np
list=np.where(img_skeleton>0)
现在列表中的点根据它们在图像中沿其中一个轴的位置进行排序。
我想对列表进行排序，使顺序代表沿线的平滑路径。 (目前这不是线弯曲回来的情况)。 随后，我想将样条曲线拟合到这些点。
已使用 arcPy here 描述并解决了类似的问题.有没有一种方便的方法可以使用 python、numpy、scipy、openCV(或其他库？)来实现这一点？
以下是示例图像。它会生成一个包含 59 个 (x,y) 坐标的列表。
当我将列表发送到 scipy 的样条拟合例程时，我遇到了一个问题，因为线上的点没有“排序”:
最佳答案
我提前为冗长的回答道歉:P(问题不是那么简单)。
让我们从改写问题开始。找到一条连接所有点的线，可以重新表述为图中的最短路径问题，其中(1)图节点是空间中的点，(2)每个节点都连接到它的 2 个最近的邻居，并且( 3) 最短路径通过每个节点只经过一次。最后一个约束非常重要(并且很难优化)。本质上，问题是找到长度为N的排列，其中排列是指每个节点的顺序(N是节点的总数)路径。
找到所有可能的排列并评估它们的成本太昂贵(如果我没记错的话，有 N! 个排列，这对于问题来说太大了)。下面我提出了一种方法，该方法找到 N 最佳排列(每个 N 点的最佳排列)，然后找到排列(从那些 N)，最大限度地减少错误/成本。
现在，让我们开始创建一个示例问题:
import matplotlib.pyplot as plt
import numpy as np
x = np.linspace(0, 2 * np.pi, 100)
y = np.sin(x)
plt.plot(x, y)
plt.show()
这里，未排序版本的点[x, y] 来模拟空间中的随机点连接成一条线:
idx = np.random.permutation(x.size)
x = x[idx]
y = y[idx]
plt.plot(x, y)
plt.show()
然后问题是对这些点进行排序以恢复其原始顺序，以便正确绘制线。
我们可以先重新排列 [N, 2] 数组中的点:
points = np.c_[x, y]
然后，我们可以从创建一个最近邻图开始，将每个节点连接到它的 2 个最近邻:
from sklearn.neighbors import NearestNeighbors
clf = NearestNeighbors(2).fit(points)
G = clf.kneighbors_graph()
G 是一个稀疏的 N x N 矩阵，其中每一行代表一个节点，列的非零元素是到这些点的欧几里得距离。
然后我们可以使用 networkx 从这个稀疏矩阵构造一个图:
import networkx as nx
T = nx.from_scipy_sparse_matrix(G)
然后，魔法开始了:我们可以使用 dfs_preorder_nodes 提取路径，这将在给定一个起始节点(如果没有给出，将选择 0 节点)的情况下，本质上创建一条通过所有节点(仅通过每个节点一次)的路径。
order = list(nx.dfs_preorder_nodes(T, 0))
xx = x[order]
yy = y[order]
plt.plot(xx, yy)
plt.show()
嗯，还不错，但我们可以注意到重建不是最优的。这是因为无序列表中的 0 点位于行的中间，也就是说它先朝一个方向前进，然后返回并朝另一个方向结束。
所以，为了得到最优的顺序，我们可以直接得到所有节点的最优顺序:
paths = [list(nx.dfs_preorder_nodes(T, i)) for i in range(len(points))]
现在我们有了从每个 N = 100 节点开始的最佳路径，我们可以丢弃它们并找到最小化连接之间距离的路径(优化问题):
mindist = np.inf
minidx = 0
for i in range(len(points)):
p = paths[i] # order of nodes
ordered = points[p] # ordered nodes
# find cost of that order by the sum of euclidean distances between points (i) and (i+1)
cost = (((ordered[:-1] - ordered[1:])**2).sum(1)).sum()
if cost < mindist:
mindist = cost
minidx = i
为每条最优路径对点进行排序，然后计算成本(通过计算所有点对 i 和 i+1 之间的欧式距离)。如果路径从 start 或 end 点开始，它将具有最小的成本，因为所有节点都是连续的。另一方面，如果路径从位于线中间的节点开始，则成本在某个点会非常高，因为它需要从线的末端(或起点)行进到初始位置探索另一个方向。使该成本最小化的路径是从最佳点开始的路径。
opt_order = paths[minidx]
现在，我们可以正确地重构顺序了:
xx = x[opt_order]
yy = y[opt_order]
plt.plot(xx, yy)
plt.show()
关于python - 对点进行排序以形成一条连续的线，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/37742358/
很好奇，就使用rubyonrails自动化单元测试而言，你们正在做什么？您是否创建了一个脚本来在cron中运行rake作业并将结果邮寄给您？git中的预提交Hook？只是手动调用？我完全理解测试，但想知道在错误发生之前捕获错误的最佳实践是什么。让我们理所当然地认为测试本身是完美无缺的，并且可以正常工作。下一步是什么以确保他们在正确的时间将可能有害的结果传达给您？ 最佳答案 不确定您到底想听什么，但是有几个级别的自动代码库控制:在处理某项功能时，您可以使用类似autotest的内容获得关于哪些有效，哪些无效的即时反馈。要确保您的提
关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点？我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等)，但我确实想创建一个输出文件。
在控制台中反复尝试之后，我想到了这种方法，可以按发生日期对类似activerecord的(Mongoid)对象进行分组。我不确定这是完成此任务的最佳方法，但它确实有效。有没有人有更好的建议，或者这是一个很好的方法？#eventsisanarrayofactiverecord-likeobjectsthatincludeatimeattributeevents.map{|event|#converteventsarrayintoanarrayofhasheswiththedayofthemonthandtheevent{:number=>event.time.day,:event=>ev
我试图获取一个长度在1到10之间的字符串，并输出将字符串分解为大小为1、2或3的连续子字符串的所有可能方式。例如:输入:123456将整数分割成单个字符，然后继续查找组合。该代码将返回以下所有数组。[1,2,3,4,5,6][12,3,4,5,6][1,23,4,5,6][1,2,34,5,6][1,2,3,45,6][1,2,3,4,56][12,34,5,6][12,3,45,6][12,3,4,56][1,23,45,6][1,2,34,56][1,23,4,56][12,34,56][123,4,5,6][1,234,5,6][1,2,345,6][1,2,3,456][123
我正在编写一个包含C扩展的gem。通常当我写一个gem时，我会遵循TDD的过程，我会写一个失败的规范，然后处理代码直到它通过，等等......在“ext/mygem/mygem.c”中我的C扩展和在gemspec的“扩展”中配置的有效extconf.rb，如何运行我的规范并仍然加载我的C扩展？当我更改C代码时，我需要采取哪些步骤来重新编译代码？这可能是个愚蠢的问题，但是从我的gem的开发源代码树中输入“bundleinstall”不会构建任何native扩展。当我手动运行rubyext/mygem/extconf.rb时，我确实得到了一个Makefile(在整个项目的根目录中)，然后当
这是一道面试题，我没有答对，但还是很好奇怎么解。你有N个人的大家庭，分别是1,2,3,...,N岁。你想给你的大家庭拍张照片。所有的家庭成员都排成一排。“我是家里的friend，建议家庭成员安排如下:”1岁的家庭成员坐在这一排的最左边。每两个坐在一起的家庭成员的年龄相差不得超过2岁。输入:整数N，1≤N≤55。输出:摄影师可以拍摄的照片数量。示例->输入:4，输出:4符合条件的数组:[1,2,3,4][1,2,4,3][1,3,2,4][1,3,4,2]另一个例子:输入:5输出:6符合条件的数组:[1,2,3,4,5][1,2,3,5,4][1,2,4,3,5][1,2,4,5,3][
我已经构建了一些serverspec代码来在多个主机上运行一组测试。问题是当任何测试失败时，测试会在当前主机停止。即使测试失败，我也希望它继续在所有主机上运行。Rakefile:namespace:specdotask:all=>hosts.map{|h|'spec:'+h.split('.')[0]}hosts.eachdo|host|begindesc"Runserverspecto#{host}"RSpec::Core::RakeTask.new(host)do|t|ENV['TARGET_HOST']=hostt.pattern="spec/cfengine3/*_spec.r
我们的git存储库中目前有一个Gemfile。但是，有一个gem我只在我的环境中本地使用(我的团队不使用它)。为了使用它，我必须将它添加到我们的Gemfile中，但每次我checkout到我们的master/dev主分支时，由于与跟踪的gemfile冲突，我必须删除它。我想要的是类似Gemfile.local的东西，它将继承从Gemfile导入的gems，但也允许在那里导入新的gems以供使用只有我的机器。此文件将在.gitignore中被忽略。这可能吗？ 最佳答案 设置BUNDLE_GEMFILE环境变量:BUNDLE_GEMFI
这似乎非常适得其反，因为太多的gem会在window上破裂。我一直在处理很多mysql和ruby-mysqlgem问题(gem本身发生段错误，一个名为UnixSocket的类显然在Windows机器上不能正常工作，等等)。我只是在浪费时间吗？我应该转向不同的脚本语言吗？ 最佳答案 我在Windows上使用Ruby的经验很少，但是当我开始使用Ruby时，我是在Windows上，我的总体印象是它不是Windows原生系统。因此，在主要使用Windows多年之后，开始使用Ruby促使我切换回原来的系统Unix，这次是Linux。Rub
这个问题在这里已经有了答案:关闭10年前。PossibleDuplicate:Pythonconditionalassignmentoperator对于这样一个简单的问题表示歉意，但是谷歌搜索||=并不是很有帮助；)Python中是否有与Ruby和Perl中的||=语句等效的语句？例如:foo="hey"foo||="what"#assignfooifit'sundefined#fooisstill"hey"bar||="yeah"#baris"yeah"另外，类似这样的东西的通用术语是什么？条件分配是我的第一个猜测，但Wikipediapage跟我想的不太一样。