質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
OpenCV

OpenCV(オープンソースコンピュータービジョン)は、1999年にインテルが開発・公開したオープンソースのコンピュータビジョン向けのクロスプラットフォームライブラリです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Q&A

解決済

1回答

728閲覧

Python フラクタルツリーを扇型にする方法が知りたい

mia_flower

総合スコア66

OpenCV

OpenCV(オープンソースコンピュータービジョン)は、1999年にインテルが開発・公開したオープンソースのコンピュータビジョン向けのクロスプラットフォームライブラリです。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

0グッド

1クリップ

投稿2022/04/30 09:27

編集2022/04/30 09:28

前提

PythonでOpenCVを用いてフラクタルツリーを作っています。

発生している問題・エラーメッセージ

以下のコードを書いたところ、添付画像のような図形が描画されました。
枝を広げた扇型を目指していたのですが、枝の結び目のy座標が揃ってしまっています。

扇型のツリーにするにはどのようにしたらよろしいのでしょうか。

イメージ説明

該当のソースコード

import numpy as np import cv2 img = np.zeros((1200,1200,3))#height, width l = 0 Deltal = 0 Tree = [] PreTree = [] Prel = 1 def Main(number): global Prel global l global Deltal global Tree global PreTree if number == 1: PreEndX = 600 PreEndY = 300 PreRadius = 300 PreTheta = 10 R = Branch(PreEndX, PreEndY, PreRadius, PreTheta, 1) #右方向へのびる枝 r = Branch(PreEndX, PreEndY, PreRadius, PreTheta, -1) #左 Tree.append(R) Tree.append(r) print(Tree) print('::::::::::::') else: for i in range(1, Deltal + 1): #増えた配列の長さの分、for文を実行 PreEndX= PreTree[Prel + i-1:Prel + i][0][2] PreEndY= PreTree[Prel + i-1:Prel + i][0][3] PreRadius=PreTree[Prel + i-1:Prel + i][0][4] PreTheta= PreTree[Prel + i-1:Prel + i][0][5] R = Branch(PreEndX, PreEndY, PreRadius, PreTheta, 1) r = Branch(PreEndX, PreEndY, PreRadius, PreTheta, -1) Tree.append(R) Tree.append(r) Prel=l l= len(Tree) Deltal=l-Prel for i in range(1,l+1): cv2.line(img, (int(Tree[i-1:i][0][0]), int(Tree[i-1:i][0][1])), (int(Tree[i-1:i][0][2]), int(Tree[i-1:i][0][3])),(255,255,255),1) print((int(Tree[i-1:i][0][0]), int(Tree[i-1:i][0][1])), (int(Tree[i-1:i][0][2]), int(Tree[i-1:i][0][3]))) PreTree=Tree #枝の座標を計算する関数 def Branch(ThisEndX, ThisEndY, ThisRadius, ThisTheta, ThisDirection): ThisBeginX=ThisEndX ThisBeginY=ThisEndY ThisRadius = ThisRadius * 0.67 ThisTheta = ThisTheta +10 if ThisTheta < 0: ThisTheta += 360; if ThisTheta >= 360: ThisTheta -= 360; ThisThetaRad = np.radians(ThisTheta) ThisEndX = ThisBeginX + np.cos(ThisThetaRad) * ThisRadius * ThisDirection ThisEndY = ThisBeginY + np.sin(ThisThetaRad) * ThisRadius print('=',ThisBeginX, ThisBeginY, ThisEndX, ThisEndY, ThisRadius, ThisTheta) return (ThisBeginX, ThisBeginY, ThisEndX, ThisEndY, ThisRadius, ThisTheta) for number in range(1, 8): Main(number) print("-------", number, "回実行されました。") #1〜7まで出力される cv2.line(img, (600, 0), (600, 1200), (255, 255, 255)) cv2.imshow('forest of objects', img) cv2.waitKey(0)

試したこと

y座標を求める式に誤りがあるのは、と考え
見直したのですが、原因がわかりませんでした。

補足情報(FW/ツールのバージョンなど)

Mac OS
python 3

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

自己解決

ThisEndX = ThisBeginX + np.cos(ThisThetaRad) * ThisRadius
にThisDirectionを掛け算していたことが誤りでした。
無事に目的を達成することができました。

投稿2022/04/30 13:03

mia_flower

総合スコア66

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問