import shapefile import Part shp = shapefile.Reader("/home/yorik/3x3 grid AKZO filled.shp") shapes = [] for i,s in enumerate(shp.shapes()): pts = [FreeCAD.Vector(p[0],p[1],0) for p in s.points] pl = Part.makePolygon(pts) rec = shp.records()[i] elev = rec['ELEV'] shape.translate(App.Vector(0,0,elev)) shapes.append(shape) cp = Part.makeCompound(shapes) Part.show(cp)