実現したいこと
「該当のソースコード」はこのページを元にしたものですが、私の質問内容の具体例として挙げただけで、深い意味はありません。7、8行目のデータはこちらのCase BのCase 5: low particle density and small particleです。
「該当のソースコード」を実行すると、Outputx.txtが出力されて、その中身は以下のようになります。
[[ 1. 0. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1.
1. 1. 1. 1. 1. 0. 1. 1. 1. 1. 1. 0. 0. 0.
0. 0. 0. 0.]
[ 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. -2.
1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 0. 1. 1. 1.
0. 0. 0. 0.]
[ 1. 1. 1. 1. 1. 1. 1. 2. 1. 2. 2. 2. 2. 1.
1. 1. 1. 1. -15. 1. 1. 1. 1. 1. 0. 0. 1. 0.
0. 0. 0. 0.]
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
[ -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -15. -1.
-1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. 0. 0. -1.
-1. -1. -1. -1.]
[ -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1.
-1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1. -1.
-1. -1. -1. -1.]]
このデータをエクセルに貼り付けたいのですが、そのときに2つ問題があります。
① [[、[、]、]]を除きたい
② 改行されてしまうので、そうならないようにしたい
欲しい出力ファイルのイメージは
1 0 1 1 1 1・・・ 0 0 0 0 0 0
1 1 1 1 1 1・・・ 1 1 0 0 0 0
1 1 1 1 1 1・・・ 1 0 0 0 0 0
・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・
・・・・・・・・・・・・・・・
-1. -1. -1. -1. -1. ・・・ 0. -1. -1. -1. -1. -1.
-1. -1. -1. -1. -1. ・・・ -1. -1. -1. -1. -1. -1.
です。
上記のような出力にする方法を教えてください。
該当のソースコード
python
1import matplotlib.pyplot as plt 2import numpy as np 3np.set_printoptions(threshold=np.inf)4from skimage.io import imread 5from scipy.signal import correlate 6 7path1 = r'.\B005_1.tif'8path2 = r'.\B005_2.tif'9 10a = imread(path1)11b = imread(path2)12 13fig, axs = plt.subplots(1, 2, figsize=(9, 4))14axs[0].imshow(a, cmap=plt.cm.gray)15axs[1].imshow(b, cmap=plt.cm.gray)16plt.show()17 18win_size = 3219 20def vel_field(curr_frame, next_frame, win_size):21 ys = np.arange(0, curr_frame.shape[0], win_size)22 xs = np.arange(0, curr_frame.shape[1], win_size)23 dys = np.zeros((len(ys), len(xs)))24 dxs = np.zeros((len(ys), len(xs)))25 26 for iy, y in enumerate(ys):27 for ix, x in enumerate(xs):28 int_win = curr_frame[y : y + win_size, x : x + win_size]29 search_win = next_frame[y : y + win_size, x : x + win_size]30 cross_corr = correlate(31 search_win - search_win.mean(), int_win - int_win.mean(), method="fft"32 )33 dys[iy, ix], dxs[iy, ix] = ( 34 np.unravel_index(np.argmax(cross_corr), cross_corr.shape)35 - np.array([win_size, win_size])36 + 1 37 )38 39 ys = ys + win_size / 240 xs = xs + win_size / 241 42 return xs, ys, dxs, dys 43 44xs, ys, dxs, dys = vel_field(a, b, win_size)45norm_drs = np.sqrt(dxs ** 2 + dys ** 2)46 47x = xs 48y = ys[::-1]49x, y = np.meshgrid(x, y)50 51file1 = open('Outputx.txt', 'w')52file1.write(str(dxs))
試したこと
①に関しては、エクセルに貼り付けた後、データ→区切り位置で1回1回対応すればできるのですが、
もっとシンプルな方法があるのではないかと思ったので、質問させていただきました。
②については、ネットで調べてもよく分かりませんでした。
0 コメント