2012-11-12 67 views
2

我使用python(特別是numpy的,SciPy的,和matplotlib)提取數據以從多個文件讀取數據。具體而言,我試圖提取Z值9個不同的點雲數據集(.ply格式:X,Y,Z,R,G,B)。從多個文件

與提取的Z(或高度值),我試圖讓高度的相對頻率直方圖,平均相對頻率,標準偏差,並從9個文件每個倉的置信區間。

我有一個工作代碼這樣做了奇異的文件,但我有困難想它的規模與9個文件的工作。我可以簡單地用許多冗餘代碼蠻力,但我試圖比這更有效一些。

我有一些代碼的效果:

filelist = os.listdir('path to 9 ply files') 
    for files in filelist: 
    CAM_XYZ = loadtxt(files,skiprows = 14, usecols = (0,1,2)) 
    CAMt = transpose(CAM_XYZ) 
    CAM_Z = CAMt[2] 

我想提取每個文件的Z值/列,然後將它們附加到另一個陣列。理想情況下,這個最終陣列將有9列,每個文件的Z值。我只是堅持如何實際做到這一點,我的當前代碼的設置方式,它只是覆蓋值。

任何建議/地方尋找將不勝感激。

回答

1

任何原因,你不能只是追加的所有值到一個列表?

filelist = os.listdir('path to 9 ply files') 
z_values = [] 
for files in filelist: 
    CAM_XYZ = loadtxt(files,skiprows = 14, usecols = (0,1,2)) 
    CAMt = transpose(CAM_XYZ) 
    z_values.append(CAMt[2]) 

這似乎是最簡單的解決方案。