surfer制作冲淤演变图

软件准备：surfer，arcgis（或者SMS，AutoCAD）

数据准备：①不同年份的水深散点数据，若水深散点的投影坐标不一致，可在arcgis中转为一致（可参考这个链接：ArcGIS平面坐标互转（从122度15分转到123度）），高程不一致也需要转到一致，②岸线shp文件，可用做白化文件；

操作步骤（step by step）:

1.打开surfer，gird模块，点Gird Data,在右边框里选择水深文件，左侧插值方法就选默认的克里金插值（Kriging），然后一路下一步(next)，直到Output选项。

2.在输出网格时：①注意不同年份的网格设置要一致，才能做后面的冲淤分析；②在这个页面可以导入白化文件，白化文件可以是提前制作好的shp文件，或者是根据数据范围在surfer中勾出来的多边形文件(.bln)，选择NoData Inside项，即在多边形内的数据被白化（删除），该步骤也可以单独操作，在Grid->Assign Nodata中选择制作好的Grid和白化多边形文件(shp,bln格式的多边形文件)，同样选择Nodata Inside选项，注意Output Gird时的文件名的命名。

3.生成的grid文件会自动导入到surfer左侧的目录栏（contents）中，此时可以对gird显示的等值线进行修改。Color scale 打勾可以显示颜色轴，Fill contours可以根据选择的fill colors对等值线进行填充。

4.使用Grids模块中的Calculate工具 Math进行冲淤计算。具体操作，将第2步生成的gird文件，添加到对话框中，导入后进行计算，导出计算后的网格。（注：若水深点为负值，则大年份减小年份，若为正值则相反）

5.冲淤图，按照自己的要求对冲淤进行修改。

6.断面提取：根据断面线和2年的Girdwe文件可以做不同年份的断面的高程变化。具体步骤：在Map Tools模块，Add to map中添加冲淤区域的断面线，然后再Grids模块中Calculate框里使用Slice来提取断面上的点数据，在提取的时候需要输入gird数据和断面线数据，保存数据选择第2个.dat格式。然后用记事本（或UE）打开，可以看到四列数据，分别是X，Y坐标，水深值，距离起始点的距离。（注：若断面超出gird范围，出现的水深值会不合理，处理的时候要把这些点删除）