December 21, 2014 12:36 · Dec 21, 2014
diff --git a/lidar2.py b/lidar2.py
@@ -0,0 +1,70 @@
+import bpy  
+import time
+import csv
+
+A = time.time()
+
+
+dfile = r"C:\Users\dealga\Desktop\Archive\SU8606_DSM_1M.asc"
+
+getval = lambda i: int(next(i).split()[1])
+
+with open(dfile) as ofile:
+    ncols = getval(ofile)
+    nrows = getval(ofile)
+    xllcorner = getval(ofile)
+    yllcorner = getval(ofile)
+    cellsize = getval(ofile)
+    NODATA_value = getval(ofile)
+
+    print(ncols, nrows, xllcorner, yllcorner, cellsize, NODATA_value)
+
+    # this will read the rest
+    verts = []
+    add_vert = verts.append
+    asc_reader = csv.reader(ofile, delimiter=' ')
+
+    # ni = nrows #? +1 -1
+    # nj = ncols #? +1 -1
+    ni = 1000
+    nj = 1000
+
+    for i, row in enumerate(asc_reader):
+        if i >= ni:
+            break
+        for j in range(int(ncols)):
+            if j >= nj:
+                break
+            z = (float(row[j]) / 60)
+            x = j * 0.01   #   cell x width
+            y = i * 0.01   #   cell y width
+            add_vert((x,y,z))
+
+    print('done')
+    print('last vertex:', verts[-1])
+
+B = time.time()
+
+total_time = B-A
+print('total_time:', total_time)
+
+faces = []
+add_face = faces.append
+
+# generate_edges, i = verts y, j = verts x
+total_range = ((ni-1) * (nj))
+indices = []
+for i in range(total_range):
+    if not ((i+1) % nj == 0):
+        add_face([i, i+nj, i+nj+1, i+1])
+
+# do your own error handling
+mesh_data = bpy.data.meshes.new("LIDAR_mesh_data3")  
+mesh_data.from_pydata(verts, [], faces)  
+mesh_data.update()
+
+LIDAR_object = bpy.data.objects.new("LIDAR_Object3", mesh_data)  
+
+scene = bpy.context.scene    
+scene.objects.link(LIDAR_object)    
+LIDAR_object.select = True
No results found