傅里叶的忧伤 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 求差集那个函数有问题,有无佬可以解决
  • ¥15 【提问】基于Invest的水源涵养
  • ¥20 微信网友居然可以通过vx号找到我绑的手机号
  • ¥15 寻一个支付宝扫码远程授权登录的软件助手app
  • ¥15 解riccati方程组
  • ¥15 display:none;样式在嵌套结构中的已设置了display样式的元素上不起作用?
  • ¥15 使用rabbitMQ 消息队列作为url源进行多线程爬取时,总有几个url没有处理的问题。
  • ¥15 Ubuntu在安装序列比对软件STAR时出现报错如何解决
  • ¥50 树莓派安卓APK系统签名
  • ¥65 汇编语言除法溢出问题