同じ図形の複数の省略記号をテキストファイルのデータを使ってプロットする

Question

私はPythonの初心者ですので、あまりにも些細なことがあればごめんなさい。

これは、テキストファイルの最初の2行の例です。

RA 12月major_axis minor_axisは、ファイルの各行は、1つの楕円をプロットするために必要とされる5つのパラメータを有する

149.20562 2.29594 0.00418 0.00310 83.40

をposition_angle。最初の2つの列は中央の列です。次の3つの列はそれぞれ長軸、短軸および位置角です。このファイルは、多数の行を持つ巨大なカタログの一部です。私はそれらの楕円をすべて1つの図にプロットしたい。

これは私が試みたものです。

import matplotlib.pyplot as plt 
import numpy as np 
import astropy.io.ascii as asciitable 
from matplotlib.patches import Ellipse 


path=/users/vishnu/Desktop/ 
fw=open(path + 'data_plot.txt', 'r') 
table = asciitable.read(path+ "data_plot.txt") 
ra_degrees=[table['ra']] 
dec_degrees=[table['dec']] 
major_axis_deg=[table['major_axis']] 
minor_axis_deg=[table['minor_axis']] 
position_angle_deg=[table['position_angle']] 



for ra, dec, w, h, angle in zip(ra_degrees, 
dec_degrees,major_axis_deg, minor_axis_deg, position_angle_deg): 
    ellipse = Ellipse(xy=(ra, dec), width=w, height=h, angle=angle) 
    ax.add_patch(ellipse) 
plt.axis('scaled') 
plt.show() 
fw.close() 

これはエラーログです。

runfile('/users/vishnu/.spyder2-py3/radio_sources.py', wdir='/users/vishnu/.spyder2-py3') 
Traceback (most recent call last): 

    File "<ipython-input-299-a0011c0326f5>", line 1, in <module> 
    runfile('/users/vishnu/.spyder2-py3/radio_sources.py', wdir='/users/vishnu/.spyder2-py3') 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/spyderlib/widgets/externalshell/sitecustomize.py", line 699, in runfile 
    execfile(filename, namespace) 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/spyderlib/widgets/externalshell/sitecustomize.py", line 88, in execfile 
    exec(compile(open(filename, 'rb').read(), filename, 'exec'), namespace) 

    File "https://stackoverflow.com/users/vishnu/.spyder2-py3/radio_sources.py", line 63, in <module> 
    ax.add_patch(ellipse) 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 1783, in add_patch 
    self._update_patch_limits(p) 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/matplotlib/axes/_base.py", line 1803, in _update_patch_limits 
    xys = patch.get_patch_transform().transform(vertices) 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/matplotlib/patches.py", line 1409, in get_patch_transform 
    self._recompute_transform() 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/matplotlib/patches.py", line 1398, in _recompute_transform 
    .scale(width * 0.5, height * 0.5) \ 

    File "https://stackoverflow.com/users/vishnu/anaconda3/lib/python3.5/site-packages/matplotlib/transforms.py", line 1965, in scale 
    np.float_) 

ValueError: setting an array element with a sequence. 

これを行うには、配列を作成する必要がなく、よりスマートな方法があれば教えてください。

Answer 1

PatchCollectionを使用する必要があるようです（this exampleを参照）。私はこの問題を推測している。この

plt.figure() 
ax = plt.gca() 
ellipses = [] 
plt.xlim([0,100]) 
plt.ylim([0,100]) 
for ra, dec, w ... zip(...): 
    ellipse = Ellipse(xy=(ra, dec), width=w, height=h, angle=angle) 
    ellipses.append(ellipse) 
p = PatchCollection(ellipses) 
ax.add_collection(p) 
plt.show() 

Answer 2

のようなものを試してみてくださいすることは、あなたが意志[「RA」] 1つのあまりにも多くのなどra_degreesにネストのレベルは、あなたがテーブルを印刷する場合、ある

を持っているということですおそらくそれは既に配列であることがわかります。角括弧で囲むと、[table ['ra']]は長さ1のリストになります。 zip関数は一度に各シーケンスから1つの項目を取ります。したがって、forループでは、最初の唯一の時間にraにtable ['ra']が割り当てられます。

別の問題は、ファイルを2回開いていることです。あなたはasciitable.readにファイル名を渡していますが、決して実際にはfwを使用していません。

以下のコードはうまくいく場合があります。

import numpy as np 
import matplotlib.pyplot as plt 
from matplotlib.patches import Ellipse 
import astropy.io.ascii as asciitable 

path=/users/vishnu/Desktop/ 
table = asciitable.read(path + "data_plot.txt") 

ra_degrees = table['ra'] 
dec_degrees = table['dec'] 
major_axis_deg = table['major_axis'] 
minor_axis_deg = table['minor_axis'] 
position_angle_deg = table['position_angle'] 

for ra, dec, w, h, angle in zip(ra_degrees, dec_degrees, 
     major_axis_deg, minor_axis_deg, position_angle_deg): 
    ellipse = Ellipse(xy=(ra, dec), width=w, height=h, angle=angle) 
    ax.add_patch(ellipse) 
plt.axis('scaled') 
plt.show() 

asciitableは、開いているファイルではなくファイル名で動作し、代わりに次を使用している場合：

with open(path + "data_plot.txt") as fw: 
    table = asciitable.read(fw) 

これは、ファイルがasciitable.readにエラーがあっても閉じていることを確認しますそれによって例外が発生します。上記のエラーを考えると、すでにデータを読み込んでいるようです。