※ 引述《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