看板 GIS 關於我們 聯絡資訊
※ 引述《QQ1512 (Anderson)》之銘言: : GIS轉出來的DTM : 前六行會出現像 : ncols 238 : nrows 194 : xllcorner 260122.500000000 : yllcorner 2655950 : cellsize 10 : NODATA_value -9999 : 這樣的訊息 : 我要用FORTRAN讀取 前六行不讀 從第七行開始讀取的矩陣 : 請問要怎麼寫程式 ? 因為每個檔案都會有這前六行 所以每讀一個檔,就先下6次read把這些值讀進來 ncol代表有幾欄 nrows代表有幾列 xllcorner左下角的x坐標值 yllcorner左下角的y坐標值 cellsize 每個點的間隔是十公尺 NODATA_value -9999代表該點無資料 再依據ncols與nrows值用迴圈方式把每個點讀進來, 看你的用途,是要存矩陣或是什麼都可以 : 還有假設我有1000筆這樣的DTM資料(矩陣數一樣) : 要怎麼使用"迴圈"的方法讀取 : 總不能open(1,file='1.txt')一直寫到open(1000,file='1000.txt') ~"~ : 然後read(1,*)a1 一直寫到 read(1000,*)a1000 吧 : 小弟不才 又再此地方發問 關於以上兩個小問題 : 知道的大大希望能幫忙解惑 感激不盡 >"< 如果有多個檔,做法很多,如果檔名是流水號,就用迴圈方式一個一個讀 對fortran也不熟,希望對你有幫助,有錯請指正!! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 123.205.124.70
QQ1512:感恩 已解決^^ 07/06 16:59