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

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

新規登録して質問してみよう
ただいま回答率
85.47%
Python 3.x

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

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

8083閲覧

python 3次元のscipy.interpolateについて

NR4200

総合スコア41

Python 3.x

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

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

2クリップ

投稿2020/09/11 03:44

前提・実現したいこと

質問を見てくださりありがとうございます.

現在私は,3次元のグラフ上に,点をプロットして,それを直線ではなく,曲線でつなげようとしています.(調べた感じは,Bスプライン補間)

イメージとして,図1のようなものから,図2のようなものを作ろうと考えています.
イメージ説明
図1
イメージ説明
図2

以下のサイトを見本にして作ったのですがうまくいきません.
https://qiita.com/Ken227/items/aee6c82ec6bab92e6abf

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

File "q.py", line 43, in <module> tck, u= interpolate.splprep(data, k=5) File "C:\Users\N\AppData\Local\Programs\Python\Python38\lib\site-packages\scipy\interpolate\fitpack.py", line 156, in splprep res = _impl.splprep(x, w, u, ub, ue, k, task, s, t, full_output, nest, per, File "C:\Users\N\AppData\Local\Programs\Python\Python38\lib\site-packages\scipy\interpolate\_fitpack_impl.py", line 220, in splprep idim, m = x.shape ValueError: too many values to unpack (expected 2)

該当のソースコード

python

1from scipy import interpolate 2from mpl_toolkits.mplot3d import Axes3D 3import numpy as np 4import matplotlib.pyplot as plt 5 6 7x = [0, 0, 1, 0] 8y = [0, -1, 0, -1] 9z = [0, -1, -2, -3] 10 11data = [[x],[y],[z]] 12 13tck, u= interpolate.splprep(data, k=5) 14 15new = interpolate.splev(np.linspace(0,1,500), tck, der=0) 16 17 18fig = plt.figure() 19ax = Axes3D(fig) 20ax.plot(data[0], data[1], data[2]) 21ax.plot(new[0], new[1], new[2], label='fit', lw =2, c='red') 22ax.legend() 23plt.show()

試したこと

2次元の場合にはうまくいくので,interpolate自体が3次元に対応していないということ?
3次元の場合は違うモジュールを使う必要がある?

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

ここにより詳細な情報を記載してください。

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

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

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

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

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

guest

回答1

0

ベストアンサー

完成形

###変更点

data = [[x],[y],[z]] -->data = [x,y,z]`

  • TypeError: m > k must hold対策でkを小さくしました。

tck, u= interpolate.splprep(data, k=5)
--> tck, u= interpolate.splprep(data, k=3)

Python3

1from scipy import interpolate 2from mpl_toolkits.mplot3d import Axes3D 3import numpy as np 4import matplotlib.pyplot as plt 5 6 7x = [0, 0, 1, 0] 8y = [0, -1, 0, -1] 9z = [0, -1, -2, -3] 10 11# data = [[x],[y],[z]] 12data = [x,y,z] 13 14# tck, u= interpolate.splprep(data, k=5) 15tck, u= interpolate.splprep(data, k=3) 16 17new = interpolate.splev(np.linspace(0,1,500), tck, der=0) 18 19 20fig = plt.figure() 21ax = Axes3D(fig) 22ax.plot(data[0], data[1], data[2], label='raw', lw =2, c='blue') 23ax.plot(new[0], new[1], new[2], label='fit', lw =2, c='red') 24ax.legend() 25plt.show() 26

投稿2020/09/21 01:33

退会済みユーザー

退会済みユーザー

総合スコア0

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.47%

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

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

質問する

関連した質問