def generate_points():
global triangles,cellsize
f = open('67250_5950_25.asc',mode='r')
ncols = int(f.readline().split()[1]);
nrows = int(f.readline().split()[1]);
xllcenter = float(f.readline().split()[1]);
yllcenter = float(f.readline().split()[1]);
cellsize = float(f.readline().split()[1]);
nodata_value = int(f.readline().split()[1]);
# print "Columns in data file: ",ncols
# print "Rows in data file: ",nrows
# print "Coordinate X-wise center (SWEREF99): ",xllcenter
# print "Coordinate Y-wise center (SWEREF99): ",yllcenter
# print "Cell size in meters: ",cellsize
# print "Value if no data available in point: ",nodata_value
row = 0
col = 0
index = 0
arr = np.zeros((ncols*nrows,3));
arr = arr.astype(np.float32)
for line in f:
valarr = line.split()
for val in valarr:
row = index/ncols
col = index%ncols
val = (float(val))
arr[index] = [row,col,val]
index += 1
#Delaunay triangulation of laser points
tri = Delaunay(np.delete(arr,2,1))
triangles = arr[tri.simplices]
triangles = np.vstack(triangles)
评论列表
文章目录