推 QQ1512:感恩 已解決^^ 07/06 16:59
※ 引述《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