回答編集履歴
1
コードの追加
test
CHANGED
|
@@ -1 +1,23 @@
|
|
|
1
1
|
解決しました。cmapを定義すればいいのですね
|
|
2
|
+
|
|
3
|
+
```Python
|
|
4
|
+
import astropy.io.fits as iofits
|
|
5
|
+
import matplotlib.pyplot as plt
|
|
6
|
+
import numpy as np
|
|
7
|
+
|
|
8
|
+
list = iofits.open("a.fit")
|
|
9
|
+
hdu = list[0]
|
|
10
|
+
data = np.array([hdu.data])
|
|
11
|
+
header = hdu.header
|
|
12
|
+
re = data.reshape(1*301, 301) #1×301×301の3次元配列を、301×301の2次元配列に変換
|
|
13
|
+
|
|
14
|
+
cmap = plt.cm.nipy_spectral #cmapを定義
|
|
15
|
+
# cmap.set_under('white')
|
|
16
|
+
cmap.set_over('white')
|
|
17
|
+
|
|
18
|
+
plt.imshow(re, cmap=cmap, vmin=5e+12, vmax=2e+14)
|
|
19
|
+
# plt.imshow(re, origin='lower', cmap='nipy_spectral', vmin=1e+13, vmax=1e+14)
|
|
20
|
+
plt.colorbar()
|
|
21
|
+
plt.savefig('a.png', dpi=200)
|
|
22
|
+
|
|
23
|
+
```
|