編集履歴

質問編集履歴

各関数の補足やソースコード作成時の参考サイトの記入を行いました。また、その際、浮動小数点をまるめる関数を不要とした為削除しました。

2020/09/27 13:02

投稿

スコア16

title CHANGED Viewed

	@@ -1,1 +1,1 @@
1	- 任意の~~データ~~の曲率半径を算出したい。
1	+ 指定のグラフの曲率半径を算出したい。

body CHANGED Viewed

@@ -1,7 +1,13 @@
 ### 前提・実現したいこと
-下記のデータの曲率半径を算出したいです。
+下記のグラフの曲率半径をすべて算出したいです。
 ![イメージ説明](7e99694c16978d154238406299cbc6d8.png)
+最終的にグラフの曲率半径最小位置を求めたいと考えている為、
+グラフの端からは端まですべての曲率半径を出力したいです。
+該当ソースコード作成時に参考にしたサイトは下記の通りです。
+[Three Points Circle (Electronic Station) - 3点を通る円の方程式を求める](http://py3.hateblo.jp/entry/2014/02/23/172305)
 ### 補足
 原因がx,yの配列がnumpyである事だと思うのですが、
 numpyではない通常(?)の配列への変換方法がわからず対処出来ていません。
@@ -72,7 +78,6 @@
 ```Python
 import math
 import numpy as np
-import matplotlib.pyplot as plt
 x =np.array([[930.49552541],
      [930.50241252],
@@ -278,33 +283,16 @@
      [0.00428208]])
+#(x1,y1)(x2,y2)(x3,y3)の3点を通る円の中心と半径を取得
 def get_circle_center_and_radius(x1, y1, x2, y2, x3, y3):
     r = []
     d = 2 * ((y1 - y3) * (x1 - x2) - (y1 - y2) * (x1 - x3))
     x = ((y1 - y3) * (y1 ** 2 - y2 ** 2 + x1 ** 2 - x2 ** 2) - (y1 - y2) * (y1 ** 2 - y3 ** 2 + x1 ** 2 - x3 ** 2)) / d
     y = ((x1 - x3) * (x1 ** 2 - x2 ** 2 + y1 ** 2 - y2 ** 2) - (x1 - x2) * (x1 ** 2 - x3 ** 2 + y1 ** 2 - y3 ** 2)) / -d
-    r = math.sqrt((x - x1) ** 2 + (y - y1) ** 2)
+    r = math.sqrt((x - x1) ** 2 + (y - y1) ** 2)　#ここで求めたrを配列Rに格納していきたい
-    print(r)
     return (x, y), r
-def remove_suffix(s, suffix):
-    return s[:-len(suffix)] if s.endswith(suffix) else s
-def make_equation(x, y, r):
+#上記で算出した円の中心座標と半径より円の方程式の作成
-    def format_float(f):
-        result = str(round(f, 2))
-        result = remove_suffix(result, '.00')
-        result = remove_suffix(result, '.0')
-        return result
-    def make_part(name, value):
-        num = format_float(abs(value))
-        sign = '-' if value > 0 else '+'
-        return name if num == '0' else '({0}{1}{2})'.format(name, sign, num)
-    return "{}^2+{}^2={}^2".format(make_part('x', x), make_part('y', y), format_float(r))
 def checkio(data):
     (x1, y1), (x2, y2), (x3, y3) = eval(data)
     (x0, y0), r = get_circle_center_and_radius(x1, y1, x2, y2, x3, y3)
@@ -315,20 +303,20 @@
 row_num = i
 R = []
 while i <= 100 :
-    X = x[row_num]
+    X = x[row_num] #1つ目の座標
     Y = y[row_num]
     i = i + 1
     row_num = i
-    X2 = x[row_num]
+    X2 = x[row_num] #2つ目の座標
     Y2 = y[row_num]
     i = i + 1
     row_num = i
-    X3 = x[row_num]
+    X3 = x[row_num] #3つ目の座標
     Y3 = y[row_num]
-    checkio("(X,Y),(X2,Y2),(X3,Y3)")
+    checkio("(X,Y),(X2,Y2),(X3,Y3)") #上記3点から作られる円の方程式を算出
     R = r
 print(R)