等经纬 等距离 格点数据

cmy

等经纬   等距离  格点数据
micaps 第4类等经纬数据可以直接转换为surfer格式的格点文件，  等距离的micaps 第4类文件如何转成surfer格式的格点文件？

holz

等间距数据可以不经过网格化就直接转成 grd 文件，请看帮助，有具体例子。

cmy

请问 帮助  及  具体例子  在什么地方？

holz

Surfer 的 ASCII 格式的网格文件规范

格式如下：

1. 
   
2. id                        第一行，只能是 DSAA
   
3. nx ny                第二行，只能是整数，nx 是 X 方向的网格列数；ny 是 Y 方向的网格行数
   
4. xlo xhi                第三行，xlo 是 X 的最小值；xhi 是 X 的最大值
   
5. ylo yhi         第四行。ylo 是 Y 的最大值；yhi 是 Y 的最大值
   
6. zlo zhi                第五行。zlo 是 Z 的最小值；zhi 是 Z 的最大值。
   
7. grid row 1        第六行。从这一行开始，一直到第 ny+5 行，都是网格 Z 值。此行对应 Y=ylo。
   
8. grid row 2        第七行。每一行有 nx 个 Z 值，从左到右对应的 X 值是 xlo 到 xhi。
   
9. grid row 3        第八行。如果某个节点值被白化，其 Z 值应该是 1.71041e38。
   
10. ……
    
11. grid row ny        第 ny+5 行，也是最后一行，此行对应 Y=yhi。
    

复制代码

例子如下：

假设有一个 10 列 9 行的网格， X 的范围 0 到 9， Y 的范围 0 到 7， Z 的范围 25 到 97.19。文件全部内容如下，可以将其保存为 demo.grd 试试。

1. 
   
2. DSAA
   
3. 10 9
   
4. 0.0 9.0
   
5. 0.0 7.0
   
6. 25.00 97.19
   
7. 91.03 77.21 60.55 46.67 52.73 64.05 41.19 54.99 44.30 25.00
   
8. 96.04 81.10 62.38 48.74 57.50 63.27 48.67 60.81 51.78 33.63
   
9. 92.10 85.05 65.09 53.01 64.44 65.64 52.53 66.54 59.29 41.33
   
10. 94.04 85.63 65.56 55.32 73.18 70.88 55.35 76.27 67.20 45.78
    
11. 97.19 82.00 64.21 61.97 82.99 80.34 58.55 86.28 75.02 48.75
    
12. 91.36 78.73 64.05 65.60 82.58 81.37 61.16 89.09 81.36 54.87
    
13. 86.31 77.58 67.71 68.50 73.37 74.84 65.35 95.55 85.92 55.76
    
14. 80.88 75.56 74.35 72.47 66.93 75.49 86.39 92.10 84.41 55.00
    
15. 74.77 66.02 70.29 75.16 60.56 65.56 85.07 89.81 74.53 51.69
    

复制代码

cmy

关键是micaps第4类文件（等距离），  数据文件头没有提供四个点的平面坐标，   提供的是经纬度值。
而且把这个点的经纬度值转换为平面坐标以后，不是严格的长方形或正方形，  有一定的旋转角度。

cmy

请高人出手。针对等距离 格点数据，必须按离散点重新插值吗？

flowers940815

请问楼主解决了这个问题吗？

chengxf

Holz 大侠的答复是完美的。
补充一点：
1、等距离的格点数据中，是没有其他坐标的，只有总距离长度除以格点数以后得到的序列数。X、Y轴都一样，从1到Xn，以及1到Yn。这就是它的坐标。
2、等经纬度不等于等距离，因为地球是球体。同样的经度值，在不同的纬度，具有不同的距离(长度)。所以你得到的“不是严格的长方形或正方形”，你感觉“有一定的旋转角度”，这很正确。为了消除你这样的感觉，只有将你的格点数据限制在一个合理的范围内，使得最终绘制的图形不要畸变太严重。
3、所谓“合理”的范围，需要根据你的数据精度取舍，精度越高，范围就应该取得越小。一般在一个城市大小的范围内可以近似地看作矩形，一个省的范围就有点不适合了，特别是南北向特别长的省份(江苏、宁夏、新疆、西藏、河南等等，大部分的省份都不合适，海南省可以)。
这些是参考意见。