傅里叶的忧伤 2019-01-11 16:55 采纳率: 25%
浏览 434
已采纳

python怎么计算殴氏距离(经纬坐标),元数据有59W行,数据为CSV格式(见图),求代码,感谢!

图片说明

            起初用的FOR循环,但数据量比较大内存不够运算(代码如下),想寻求一些简单的算法。
print('--------for循环计算殴氏距离试验---------')
number=len(X)
print(number)
print('-----for循环------')
knn={}
for i in range(number):
    for j in range(number):
        knn[i,j]=np.sqrt((X['LONGITUDE'][i]-X['LONGITUDE'][j])**2+(X['LATITUDE'][i]-X['LATITUDE'][j])**2)
print(knn)
  • 写回答

1条回答 默认 最新

  • ndcy 2021-06-26 16:15
    关注
    # -*- coding:utf-8 -*-
    """
    author:lqr
    date:20210626 10:00
    theme:Dimension reduction
    参考:https://blog.csdn.net/A793539835/article/details/116265977
    """
    import pandas as pd
    import numpy as np
    import csv
    import math
    import copy
    f = open(r'mydata.csv')
    reader = csv.reader(f)
    m = []
    for row in reader:
        m.append(row)
    od = []  # outputdata
    for i in range(1, 32401):   # (1, 53)
        subOD = []  # sub outputdata
        for j in range(1, 121):  # (1, 53)
            deltX = float(m[i][1])-float(m[j][1])
            deltY = float(m[i][2])-float(m[j][2])
            len = math.sqrt((deltX**2)+(deltY**2))   # a**2 ==>求a的平方    a**3 ==>求a的立方
            subOD.append(round(len, 3))
        od.append(subOD)
    print(od)
    
    csvFile = open('distance.csv', 'w+', newline="")
    writer = csv.writer(csvFile)
    writer.writerows(od)
    csvFile.close()
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

问题事件

  • 已采纳回答 1月9日

悬赏问题

  • ¥15 CARSIM前车变道设置
  • ¥50 三种调度算法报错 有实例
  • ¥15 关于#python#的问题,请各位专家解答!
  • ¥200 询问:python实现大地主题正反算的程序设计,有偿
  • ¥15 smptlib使用465端口发送邮件失败
  • ¥200 总是报错,能帮助用python实现程序实现高斯正反算吗?有偿
  • ¥15 对于squad数据集的基于bert模型的微调
  • ¥15 为什么我运行这个网络会出现以下报错?CRNN神经网络
  • ¥20 steam下载游戏占用内存
  • ¥15 CST保存项目时失败