2024年6月9日发(作者：)

ascii 码的grd文件及插值问题

在第七部分给出的grd文件是ascii码的直接grd文件格式，无需进行grid过程，这

样 做简单，准确，可以在工作和学习过程中提供很多方便。下面做一详细介绍：

ascii码的直接grd文件格式包括5个文件头信息(关于网格的限制)，然后是z的值，

值得注 意的是，这种文件格式必须用空格来区分两个数据。其一般格式为:

DSAA (注意,必须大写且从第一行第一列开始,即DSAA前面不能有空格)

nx ny (x 方向网格点的个数nx, y方向网格点的个数ny)

xlo xhi (x 方向网格点的最小值和最大值)

ylo yhi (y 方向网格点的最小值和最大值)

zlo zhi (z 方向网格点的最小值和最大值)

grid row 1 (z 的值在点(x1,y1), (x2,y1),....(xn,y1))

grid row 2 (z 的值在点(x1,y2), (x2,y2),....(xn,y2))

grid row ........ grid row n (z 的值在点(x1,yn), (x2,yn),....(xn,yn)) 注意,此处x1

y1

下面是一个具体的例子:

DSAA

10 10

0.0 9.0

11.0 20.0

25.00 97.19

91.03 77.21 60.55 46.67 52.73 64.05 41.19 54.99 44.30 25.00 96.04 81.10 62.38

48.74 57.50 63.27 48.67 60.81 51.78 33.63 92.10 85.05 65.09 53.01 64.44 65.64

52.53 66.54 59.29 41.33 94.04 85.63 65.56 55.32 73.18 70.88 55.35 76.27 67.20

45.78 97.19 82.00 64.21 61.97 82.99 80.34 58.55 86.28 75.02 48.75 91.36 78.73

64.05 65.60 82.58 81.37 61.16 89.09 81.36 54.87 86.31 77.58 67.71 68.50 73.37

74.84 65.35 95.55 85.92 55.76 80.88 75.56 74.35 72.47 66.93 75.49 86.39 92.10

84.41 55.00 74.77 66.02 70.29 75.16 60.56 65.56 85.07 89.81 74.53 51.69 70.00

54.19 62.27 74.51 55.95 55.42 71.21 74.63 63.14 44.99

在这个例子中，共有100个z的值，其中x为10个，y为10个，x的最小最大值是

0.0和 9.0, y的最小最大值是11.0和20.0，z的最小最大值是25.00和97.19，最后是100

个z的 值， 每行10个值，共10行。 当然，我们也没有必要这样排列，只有顺序不变，

照样可以 每行只有一个值，即类似下面的排法: .... (省略了文件头信息) 91.03 77.21

60.55 ... 74.63 63.14 44.99 这里我给出一个FORTRAN语言生成ascii码grd文件的子程

序,仅供参考: