import bpy
import math
# mesh 數(shù)組（點(diǎn)、面、邊）
verts = []
faces = []
edges = []
#3D supershape 參數(shù)
m = 14.23
a = -0.06
b = 2.78
n1 = 0.5
n2 = -.48
n3 = 1.5
scale = 3
Unum = 50
Vnum = 50
Uinc = math.pi / (Unum/2)
Vinc = (math.pi/2)/(Vnum/2)
# 套用公式，填充頂點(diǎn)數(shù)組
theta = -math.pi
for i in range (0, Unum + 1):
   phi = -math.pi/2
   r1 = 1/(((abs(math.cos(m*theta/4)/a))**n2+(abs(math.sin(m*theta/4)/b))**n3)**n1)
   for j in range(0,Vnum + 1):
       r2 = 1/(((abs(math.cos(m*phi/4)/a))**n2+(abs(math.sin(m*phi/4)/b))**n3)**n1)
       x = scale * (r1 * math.cos(theta) * r2 * math.cos(phi))
       y = scale * (r1 * math.sin(theta) * r2 * math.cos(phi))
       z = scale * (r2 * math.sin(phi))
       vert = (x,y,z) 
       verts.append(vert)
       #增加 phi
       phi = phi + Vinc
   #增加 theta
   theta = theta + Uinc
# ------------------------------------------------------------------------------- 
# 填充面數(shù)組
count = 0
for i in range (0, (Vnum + 1) *(Unum)):
    if count < Vnum:
        A = i
        B = i+1
        C = (i+(Vnum+1))+1
        D = (i+(Vnum+1))
        face = (A,B,C,D)
        faces.append(face)
        count = count + 1
    else:
        count = 0
# 創(chuàng)建 mesh 和 object
mymesh = bpy.data.meshes.new("supershape")
myobject = bpy.data.objects.new("supershape",mymesh)
# 設(shè)置 mesh 的 location
myobject.location = bpy.context.scene.cursor.location # *
bpy.context.scene.collection.objects.link(myobject) # *
# 從 python 數(shù)據(jù)創(chuàng)建 mesh
mymesh.from_pydata(verts,edges,faces)
mymesh.update(calc_edges=True)
# 設(shè)置 object 為編輯模式
bpy.context.view_layer.objects.active = myobject # *
bpy.ops.object.mode_set(mode='EDIT')
# 移除重復(fù)的頂點(diǎn)
bpy.ops.mesh.remove_doubles() 
# 重新計(jì)算法線
bpy.ops.mesh.normals_make_consistent(inside=False)
bpy.ops.object.mode_set(mode='OBJECT')
# 新建細(xì)分修改器（subdivide modifier）
myobject.modifiers.new("subd", type='SUBSURF')
myobject.modifiers['subd'].levels = 3
# 平滑 mesh
mypolys = mymesh.polygons
for p in mypolys:
    p.use_smooth = True