我與校準軟件PEST,這只是一個示例站點文件類型兼容內的公用事業工作。我想用此實用程序來校準我的模型,但我的模型輸出格式非常不同(CSV)。 SSF是製表符分隔的,僅包含4個關鍵數據列:示例站點ID,日期(MM/DD/YYYY),時間和Streamflow。 SSF文件的轉換CSV文件的水文模型
例子:
134 01/01/1980 00:00:00 34
134 01/02/1980 00:00:00 30
134 01/03/1980 00:00:00 28
134 01/04/1980 00:00:00 38
從我的水文模型的輸出被寫入到包含日期(年月日),仿真流量(QSIM),觀察流量(Qobs)CSV文件,氣溫,降水,實際蒸發量,潛在蒸發量和雪水當量。模型輸出的
例子:
134
Date, Qsim, Qobs, Temp, Precip, AET, PET, SWE
19800101, 34, 31, 11, 21, 3.4, 4.0, 0
19800102, 30, 30, 11, 15, 3.0, 4.4, 0
19800103, 28, 25, 12, 0, 3.1, 4.0, 0
19800104, 38, 45, 8, 30, 0.5, 3.8, 0
所以我決定,爲了使該轉換我會:
- 添加包含網站ID列0(行0在CSV文件)CSV文件
- 刪除行0和1
- CSV文件
- 添加AC刪除行2-7 olumn的時間(在這種情況下,我只需要在「00:00:00」的第2列的每一行)
- 格式化日期從年月日到MM/DD/YYYY
- 寫文件的製表符分隔
我絕對有蟒蛇一個初學者,這是我的代碼已經能夠到目前爲止寫:
import csv
HBVout = csv.reader(open('C:\\ENVpest\\Output\\Results.csv', 'rb'))
HBVout.next()
newSSF = csv.writer(open('SSF1.txt', 'wb+'), delimiter='\t')
for cline in HBVout:
new_line = [val for col, val in enumerate(cline) if col not in (2,3,4,5,6,7)]
newSSF.writerow(new_line)
我能跳過第一行,刪除不必要的列,並重寫爲製表符分隔的文件。我非常感謝將任何列添加到文件並重新設置日期的幫助!我已經搜索了SO和其他網站上的其他問題,但還沒有找到任何運氣。
謝謝!是的,「刪除第2-7行」是一個錯字,意思是「刪除第2-7列」。 雖然我在運行代碼時遇到了語法錯誤,關於「open()as inf,open()as outf」)之間的逗號以及我不確定爲什麼這是語法錯誤。你有什麼可以解決它的想法嗎? 除此之外,我真的很喜歡你爲此寫的劇本,它非常簡單,並以非常直接的方式完成了所有的事情。我絕對需要更多地關注DictReader,它似乎是重組CSV文件的一種非常強大的方式。 – Geogrammer
Python 2.7之前不允許多個上下文表達式。解決方法是將一個嵌套在另一個內部,這將使得還需要縮進另外兩個級別中的所有語句。 – martineau
是的,使用嵌套的「work ... as」語句,但我現在收到「KeyError:Date」。對於DictReader()幾乎沒有經驗,我剛開始研究一下,但還沒有找到明確的解決方案來解決這個關鍵錯誤。你有什麼建議嗎? – Geogrammer