vispy_pointcloud.py 文件源码

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)