#  簡単なチェック
	#  線形状が平らであるかのチェックは省略。
	#  平らでない閉じた線形状の場合、エラーにはならないが、得られる結果は意味を持たない
def check_object(ob) :
	if ob.type != 4 :	#  線形状でない
		return False
	elif not ob.closed :	#  閉じていない
		return False
	elif ob.number_of_control_points < 2 :	#  2つ以上のアンカーポイントがない
		return False
	else :
		ob2 = ob.parent
		return not (ob2.type == 2 and ob2.part_type == 1)	#  自由曲面内の線形状は対象外



# bezier の4つの制御点座標を取得
def get_bezier_data(ob) :
	n = ob.total_number_of_control_points
	ap = []
	inh = []
	outh = []
	for i in range(n) :
		cp = ob.control_point(i)
		ap.append(cp.position)
		inh.append(cp.in_handle)
		outh.append(cp.out_handle)
	p = []
	for i in range(n - 1) :
		bz = []
		bz.append(ap[i])
		bz.append(outh[i])
		bz.append(inh[i + 1])
		bz.append(ap[i + 1])
		p.append(bz)
	bz = []
	bz.append(ap[n - 1])
	bz.append(outh[n - 1])
	bz.append(inh[0])
	bz.append(ap[0])
	p.append(bz)
	
	return p


# ベクトル長さ
def get_len_vec(v) :
	return math.sqrt(v[0]**2 + v[1]**2 + v[2]**2)


#  ベクトルの差
def get_sub_vec(p2, p1) :
	return [p2[0] - p1[0], p2[1] - p1[1], p2[2] - p1[2]]
	

#  ベクトルの和
def get_add_vec(p2, p1) :
	return [p2[0] + p1[0], p2[1] + p1[1], p2[2] + p1[2]]
	
	
#  ベクトルの外積
def get_cross_vec(v1, v2) :
	return[v1[1]*v2[2] - v1[2]*v2[1], v1[2]*v2[0] - v1[0]*v2[2], v1[0]*v2[1] - v1[1]*v2[0]]
	

# p1 から p2 への単位ベクトル
def get_unit_vec(p2, p1) :
	v = get_sub_vec(p2, p1)
	r = get_len_vec(v)
	if abs(r) > 0.000001 :
		return [v[0]/r, v[1]/r, v[2]/r]
	else :
		return None
	

#  二つの単位ベクトルの交差角
def get_angle_of_unit_vec(v1, v2) :
	a = v1[0]*v2[0] + v1[1]*v2[1] + v1[2]*v2[2]
	a = max(-1, a)
	a = min(1, a)
	return math.acos(a)

	
#  bezier 区間の全曲率の近似値を簡易的に取得
def get_total_curvature(bz) :
	v1 = get_unit_vec(bz[1], bz[0])
	v2 = get_unit_vec(bz[2], bz[1])
	v3 = get_unit_vec(bz[3], bz[2])	
	b1 = (v1 != None)
	b2 = (v2 != None)
	b3 = (v3 != None)
	
	if b1 and b2 and b3 :
		return get_angle_of_unit_vec(v1, v2) + get_angle_of_unit_vec(v2, v3)
	elif b1 and b2 :
		return get_angle_of_unit_vec(v1, v2)
	elif b2 and b3 :
		return get_angle_of_unit_vec(v2, v3)
	elif b1 and b3 :
		return get_angle_of_unit_vec(v1, v3)
	else :
		return 0
		

# bezier の階差を求める
def get_bezier_difference(bz) :
	d11 = get_sub_vec(bz[1], bz[0])
	d12 = get_sub_vec(bz[2], bz[1])
	d13 = get_sub_vec(bz[3], bz[2])
	d21 = get_sub_vec(d12, d11)
	d22 = get_sub_vec(d13, d12)
	d31 = get_sub_vec(d22, d21)
	return [bz[0], d11, d21, d31]


#  bezier parameter t における座標を返す
def get_point_on_line(b, t) :
	t2 = t*t
	t3 = t*t*t
	x = b[3][0]*t3 + 3*b[2][0]*t2 + 3*b[1][0]*t + b[0][0]
	y = b[3][1]*t3 + 3*b[2][1]*t2 + 3*b[1][1]*t + b[0][1]
	z = b[3][2]*t3 + 3*b[2][2]*t2 + 3*b[1][2]*t + b[0][2]
	return [x, y, z]
	
	
def main() :
	#  簡単なチェック
	ob = scene.active_shape()
	if not check_object(ob) :
		xshade.beep()
		return

	#  線形状の座標を取得
	scene.inhibit_update()
	scene.copy_object(None, None, None, None)
	ob2 = scene.active_shape()
	scene.create_part()
	ob3 = scene.active_shape()
	ob3.cancel_transformation()
	
	ob2.place_child(1)
	p = get_bezier_data(ob2)	# bezier の4つの制御点座標を取得
	c = ob2.center_point		#  線形状の中心座標
	ob3.remove()
	scene.allow_update()
	scene.force_update()
	
	nb = len(p)
	areaV = [0, 0, 0]
	for i in range(nb) :		#  各 bezier 区間について
		bz = p[i]
		theta = get_total_curvature(bz)		#  曲線の曲がり具合を簡易的に求める
		epsilon = 0.05/180*math.pi			# epsilon を小さくすれば精度は上がる
		n = max(1, int(math.ceil(theta/epsilon)))	#  当該 bezier 区間の分割数

		b = get_bezier_difference(bz)			#  bezier の階差表示を求める
		
		v2 = get_sub_vec(bz[0], c)
		v1 = []
		for i in range(n) :
			v1 = v2
			p2 = get_point_on_line(b, float(i + 1)/n)	#  bezier parameter t = float(i + 1)/n) における座標を求める
			v2 = get_sub_vec(p2, c)
			a = get_cross_vec(v1, v2)
			areaV = get_add_vec(areaV, a)
	
	area = get_len_vec(areaV)/2
	k = scene.native_to_user_unit
	area2 = area*k*k
		
	print ''
	print str(area) + ' ( mm^2 )'		
	print str(area2) + ' ( user unit ^2 )'	


	

import math
scene = xshade.scene()
main()