在地质区划图上绘制伪震中

import numpy as np
import matplotlib.pyplot as plt
#matplotlib inline
from matplotlib import image
from matplotlib import pyplot as plt
#单纯画图
def loadDataSet(fileName):  
    dataMat = []              # 初始化一个空列表，文件的最后一个字段是类别标签
    fr = open(fileName)       # 读取文件
    for line in fr.readlines():   # 循环遍历文件所有行
        curLine = line.strip().split(' ')   # 切割每一行的数据
        fltLine = list(map(float, curLine))  # 映射所有的元素为 float（浮点数）类型
        dataMat.append(fltLine)   # 将数据追加到dataMat
    return dataMat       # 返回dataMat


datMat1 = np.mat(loadDataSet('./nolabel_expand.txt'))
data = image.imread('地图.png')
plt.figure(1)
plt.imshow(data,extent=(0, data.shape[1], 0, data.shape[0]))

x = list((datMat1)[:,0])
print(x)
y = list((datMat1)[:,1])
plt.scatter(x, y,c='green', s=10,)
plt.show()

 然后用聚类的方法求实际最大可能的震中位置即可