在阵列中找到差异

Question

我在Python中有一个阵列，该阵列由电线之间的米（链条）距离（链条）组成。这是

d=([0, 98, 198, 300, 400, 500, 600, 700, 800, 900, 997, 1102, 1187, 1282, 1382, 1480, 1570, 1670, 1775, 1885, 1980, 2083, 2178])

我想编写一个脚本以计算每1000m或接近1000m但不超过1000m的差异（部分长度）的示例。

这是一个示例

• pection1 = 997-0

• pection2 = 1980-997

• pection3 = 2178-1980

结果应该在阵列中，

Section = ([997, 983, 198])

我是Python的初学者。我还在学习。我没有附加任何代码，因为我对从哪里开始不了解。

Answer 1

因此，这将稍久，因为您是一个初学者，我将尝试分解我在这里所做的事情，而不仅仅是粘贴代码。如果您有任何疑问，请随时回答和询问，但是无论如何，我们就可以了！因此，听起来您想开始使用您测量的所有值来创建一个列表。在Python中，您只需使用[]（并且您将使用（）创建一个元组）。因此，从您的列表开始...

d=[0, 98, 198, 300, 400, 500, 600, 700, 800, 900, 997, 1102,1187, 1282, 1382, 1480, 1570, 1670, 1775, 1885, 1980, 2083, 2178]

接下来，我们将需要找到差异。对于初学者来说，让我们创建一个变量，该变量将跟踪我们使用的最后一个值（将要减去），在这种情况下，该变量将以零为零。接下来，我们还将创建我们的第二个列表，我们将在其中存储差异。

current = 0
section = []

现在是通过它们循环的最棘手的部分。因此，在Python中，您可以遍历列表，但是要检查我们目前正在使用的项目，并确保它是最接近1000的项目，而无需越过，我们也需要访问下一个项目。因此，我们将通过一组数字进行迭代，而不是直接迭代列表，这些数字将是从第一个开始的列表中项目的地址。

for i in range(len(d) - 1):

此语句在这里说我们是从将I设置为零开始，并且将通过代码循环我们接下来输入，直到达到列表D负1的长度。因此，现在检查当前项目是否最接近1000不走过。

for i in range(len(d) - 1):
    if ((d[i] % 1000) == 0) or ((d[i] % 1000) > (d[i + 1] % 1000)):
        section.append(d[i] - current)
        current = d[i]

因此，现在本节说，如果D [i] mod 1000等于零（地址i的D项目是1000），或者如果正在研究的当前项目Mod 1000大于下一个（这表明下一个大于大于1000的倍数）然后，我们将一个元素添加到等于d [i]的差异列表中 - 然后将电流更改为等于d [i]。

这将负责大部分工作，现在我们只需要完成最后一行即可。

section.append(d[len(d) - 1] - current)

这条线采用了最后一个元素，并始终添加它和我们当前价值之间的区别，就像您似乎想要的那样。因此，最后我们的整个代码是...

d=[0, 98, 198, 300, 400, 500, 600, 700, 800, 900, 997, 1102,1187, 1282, 1382, 1480, 1570, 1670, 1775, 1885, 1980, 2083, 2178]
current = 0
section = []
for i in range(len(d) - 1):
    if ((d[i] % 1000) == 0) or ((d[i] % 1000) > (d[i + 1] % 1000)):
        section.append(d[i] - current)
        current = d[i]

section.append(d[len(d) - 1] - current)

希望这可以帮助您解决问题并学习一些东西。有一个很棒的！

Answer 2

既然您需要一个 numpy 解决方案，那么让我们取消那些 for 循环：

d=np.array([0, 98, 198, 300, 400, 500, 600, 700, 800, 900, 997, 1102, 1187, 1282, 1382, 1480, 1570, 1670, 1775, 1885, 1980, 2083, 2178])

def func(d, k = 1000):
    insert = (np.arange(d.max() // k) + 1) * k                   #intermediate points
    vals = np.r_[d[0], d[np.searchsorted(d, insert) - 1], d[-1]] #values below the intermediate points
    return np.diff(vals)                                         #difference between the intermediate points

func(d)
Out[]: array([997, 983, 198])

Answer 3

即使你不知道如何有效地解决问题，尝试写一些东西，它可能只是如何获得 1000m 间隔，甚至是找到最接近值的繁琐脚本。这是学习如何编码的唯一方法，你必须练习，当你遇到困难时就去询问。从阅读别人的代码中获得的价值是有限的。

为了解决您的问题，我计算了与给定的特定列表相关的间隔。然后，我创建一个列表来存储列表中与间隔最接近的值，并使用 for 循环计算间隔。我使用列表理解来计算相关部分的差异。

    intervals = np.arange(min(d),max(d)+1000,1000)
    closest_values = [min(d)]
    for i in intervals:
        diff = abs((d-i)/i)
        closest_values.append(d[np.argmin(diff)])
    section = [closest_values[i+1] - closest_values[i] for i,_ in enumerate(closest_values[:-1]]

使用您的输入列表，将其转换为 numpy 数组后，给出了以下结果：

[997, 983, 198]

Answer 4

另一个选择：这可能是一个很好的机会，可以看一下 > groupby（） 来自标准库的 模块。

尝试

from itertools import groupby

d = [0, 98, 198, 300, 400, 500, 600, 700, 800, 900, 997, 1102, 1187, 1282, 1382, 1480, 1570, 1670, 1775, 1885, 1980, 2083, 2178]

for key, group in groupby(d, key=lambda n: n // 1000):
    print(f"{key = }, {list(group) = }")

，您将获得

key = 0, list(group) = [0, 98, 198, 300, 400, 500, 600, 700, 800, 900, 997]
key = 1, list(group) = [1102, 1187, 1282, 1382, 1480, 1570, 1670, 1775, 1885, 1980]
key = 2, list(group) = [2083, 2178]

groupby chops d根据键>键函数的结果，这是1000 。

考虑到这一点，您可以尝试

result, last = [], 0
for _, group in groupby(d, key=lambda n: n // 1000):
    *_, current = group
    result.append(current - last)
    last = current

获得

[997, 983, 198]