Meteoinfo画散点密度图，读取数据出现问题

呼小喵喵

  参照王老师之前发的脚本，对读取数据那块进行了修改，结果还是错误。求指导脚本如下：

#Create normal distribution point data fn = 'f:/lightning-0514/0514IC.txt' tf = open(fn) x1 = [] y1 = [] for aline in tf:   if aline.isspace():      continue      datalist = aline.split()       x = float(datalist[1])       y = float(datalist[0])       x1.append(x)       y1.append(y)   x = array(x1)   y = array(y1)      #Calculate point density x_g = linspace(x.min(), x.max(), 100) y_g = linspace(y.min(), y.max(), 100) z2,x2,y2,z = griddata((x,y), None, xi=(x_g, y_g), method='inside_count') #Sort the points by density, so that the densest points are plotted last idx = argsort(z) x = x[idx] y = y[idx] z = z[idx] #Plot scatter point colored by density ls = scatter(x, y, s=4, c=z, edgecolor=None, cnum=20) colorbar(ls) title('Scatter density plot example')

提示错误如下：
File "<iostream>", line 10
    x = float(datalist[1])
    ^
SyntaxError: mismatched input '' expecting DEDENT


数据就两列：
0514IC.txt (189.58 KB, 下载次数: 20)

2017-5-24 16:58 上传

点击文件名下载附件

呼小喵喵

#Create normal distribution point data
fn = 'f:/lightning-0514/0514IC.txt'
tf = open(fn)
x1 = []
y1 = []
for aline in tf:
  if aline.isspace():
     continue
     datalist = aline.split()
      x = float(datalist[1])
      y = float(datalist[0])
      x1.append(x)
      y1.append(y)
  x = array(x1)
  y = array(y1)  
  
#Calculate point density
x_g = linspace(x.min(), x.max(), 100)
y_g = linspace(y.min(), y.max(), 100)
z2,x2,y2,z = griddata((x,y), None, xi=(x_g, y_g), method='inside_count')

#Sort the points by density, so that the densest points are plotted last
idx = argsort(z)
x = x[idx]
y = y[idx]
z = z[idx]

#Plot scatter point colored by density
ls = scatter(x, y, s=4, c=z, edgecolor=None, cnum=20)
colorbar(ls)
title('Scatter density plot example')

MeteoInfo

应该是代码缩进的问题。

呼小喵喵

MeteoInfo 发表于 2017-5-24 19:24
应该是代码缩进的问题。

好的，谢谢王老师，我自己研究研究~~

呼小喵喵

MeteoInfo 发表于 2017-5-24 19:24
应该是代码缩进的问题。

王老师，我解决问题了，谢谢老师。~
#Create normal distribution point data
fn = 'D:/temp/0514IC.txt'
tf = open(fn)
tf.readline()
aline = tf.readline()
x1 = []
y1 = []
for aline in tf:
    datalist = aline.split()
    y = float(datalist[0])
    x = float(datalist[1])
    x1.append(x)
    y1.append(y)
x = array(x1)
y = array(y1)  

#Calculate point density
x_g = linspace(x.min(), x.max(), 100)
y_g = linspace(y.min(), y.max(), 100)
z2,x2,y2,z = griddata((x,y), None, xi=(x_g, y_g), method='inside_count')

#Sort the points by density, so that the densest points are plotted last
idx = argsort(z)
x = x[idx]
y = y[idx]
z = z[idx]

#Plot scatter point colored by density
ls = scatter(x, y, s=4, c=z, edgecolor=None, cnum=20)
colorbar(ls)
title('Scatter density plot example')

木船渔人

Traceback (most recent call last):
  File "<iostream>", line 20, in <module>
TypeError: 'NoneType' object is not iterable

为什么报错了呢？