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

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

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

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

Q&A

2回答

692閲覧

2次元配列から3次元配列へ

halu0301

総合スコア8

NumPy

NumPyはPythonのプログラミング言語の科学的と数学的なコンピューティングに関する拡張モジュールです。

Python 3.x

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

0グッド

0クリップ

投稿2018/08/14 14:31

編集2022/01/12 10:55

DICOM画像を扱っています。3次元配列から2次元配列へ分割。その後、2次元配列の画像に処理を加えて上書きし、3次元配列へ戻したいのですが。分割したものを元にもどせない。
python3.x
コード
print(dataset.pixel_array)
print(dataset.pixel_array.shape)

ds = dataset.pixel_array
ds = np.split(ds,380,axis=0)

以下のような配列です。
マトリックス1024*1024の画像が380枚
[[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]
...
[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]]



[[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]
...
[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]
[0 0 0 ... 0 0 0]]]
(380, 1024, 1024)

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

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

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

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

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

hayataka2049

2018/08/14 14:43

どんな方法で分割していますか?
halu0301

2018/08/14 14:56

返信ありがとうございます。質問内容を更新させて頂きましたが、伝わりましたでしょうか?
guest

回答2

0

380枚それぞれに処理したいとわかっているのであればこんな感じでよいのでは。

python

1import numpy as np 2 3 4ds = dataset.pixel_array 5result = np.zeros(ds.shape) 6 7for idx in range(ds.shape[0]): 8 print(ds[idx:idx+1].shape) 9 # (1, 1024, 1024) 10 11 # 何かの処理 12 new_array = ds[idx:idx+1] 13 result[idx:idx+1] = new_array 14

投稿2018/08/15 02:18

wakame

総合スコア1170

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

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

halu0301

2018/08/17 23:52

ありがとうございました。解決できました。
guest

0

(380, 1024, 1024)を380個の2次元配列にして、それぞれ独立に処理したあと元に戻したいということですね。

一番簡単な方法は、380個の2次元配列が要素になっているリストにした上で

python

1lst = 380個の2次元配列が要素になっているリスト 2np.array(lst)

こんな感じでいけますね。

他にnp.newaxisを使う方法もありますが、今回のケースだと記述が面倒くさいだけであまりメリットはないかもしれません。

以下簡単な実例。

python

1>>> import numpy as np 2>>> a = np.array([[[1,2],[3,4]],[[5,6],[7,8]]]) 3>>> a 4array([[[1, 2], 5 [3, 4]], 6 7 [[5, 6], 8 [7, 8]]]) 9>>> a[0] 10array([[1, 2], 11 [3, 4]]) 12>>> a[1] 13array([[5, 6], 14 [7, 8]]) 15>>> np.array([a[0], a[1]]) # リストに入れてnumpy配列に変換 16array([[[1, 2], 17 [3, 4]], 18 19 [[5, 6], 20 [7, 8]]]) 21>>> np.vstack([a[0][np.newaxis,...], a[1][np.newaxis,...]]) # newaxisを使う場合(こちらで書く場合、内包表記やループ等を適宜活用してください) 22array([[[1, 2], 23 [3, 4]], 24 25 [[5, 6], 26 [7, 8]]])

投稿2018/08/14 15:04

hayataka2049

総合スコア30933

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだベストアンサーが選ばれていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.50%

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

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

質問する

関連した質問