surfer制作冲淤变化图.doc

冲淤变化图制图要点注意：在研究一个新的区域的冲淤变化时，最好首先去确定冲淤变化的研究范围，根据这个范围取点、白化等，这很重要1、原始数据的处理：水深点、背景底图目的是得到可以在surfer中处理的数据，主要包括水深数据和工程图（范围、建筑、标记，用作背景）原始数据：一般为各种cad工程图（如果有处理过的水深点数据就可以直接用），其中用到的主要数据是各时间段某个区域的水深数据处理方法：（1）检查cad图中水深点数据的属性（特性），有时有些数据点的属性中没有水深值，则必须提取出图中标记在点上的那个值，此时可用韩志远给的插件提取，得到TXT文件，如图1所示。（2）如果是用作背景的工程范围图，则可直接另存为dxf格式的数据，然后在surfer中作为背景图直接使用。图1 用插件从dwg数据中提取的水深点数据2、用水深数据生成网格源数据：水深点数据，dat、txt等格式处理方法：（1）打开surfer，点击Grid-Data，弹出后选择水深数据文件。（2）根据水深文件的具体情况及制图的需求对数据进行网格化，保存为.grd网格数据，如图2、图3所示。（3）对各个时期各个工程区域的水深数据都进行网格化注意：sufer计算冲淤变化时要求参与计算的两个grd文件有相同大小和数量的网格，即他们的xy坐标范围和节点数量（node）都相同，根据此设置图3中的参数。（可参考网格化文件夹中的参数截图）图2 选择水深文件后注意文件内点数据的排序图3 注意XYZ对应的数据列，根据点的密度选择合适的网格搜索半径，根据需要设置网格大小和XY的范围3、白化处理 blankSurfer对水深数据插值计算生成网格数据时是覆盖整个矩形范围的，无论每个网格区域内是否有数据点都会在该网格中进行插值计算。因此计算冲淤变化之前先要确定是否有需要白化的区域，如果有就要先制作bln白化文件源数据：网格化后的grd数据处理方法：打开网格化后的grd数据（以contour map的形式），可看到有效工程范围内有等深线，右键点击地图，弹出菜单中选择Digitize，此时可在地图上点选一个多边形，将无效区域围起来（可能需要多个多边形才能将无效区域包围，不过尽可能用一个多边形包围），Digitized 对话框会记录所有点的坐标，如图4所示。多边形绘制玩后，在Digitized 对话框中选择File-Save as，将其保存为bln文件，完成后还要再用记事本打开改文件，把第一个坐标行复制一份到最后一行，使多边形闭合，第一行第一个数字表示点数，加1，第二个数字为标志位，选择0，保存，如图5所示。如果有多个多边形，则对相应的bln文件修改后可将其内容合并到一个bln文件中。最后用Grid-Blank对grd数据进行白化处理图4 用数字化工具选择需要白化的多边形区域图5 bln文件格式4、用网格数据计算冲淤变化图源数据：白化后的grd数据（如不需要白化，则用网格化后的grd数据），一般分为不同时期不同区域处理方法：选择Grid-Math，弹出的对话框中先后选择同一区域不同时间的grd数据，用时间小的数据减去时间大的数据（如13年-14年），保存计算结果，如图6所示。图6 计算冲淤变化5、调整冲淤变化图在冲淤变化图（contour map）的属性管理器中设置其填充、等深线级别和颜色卡等等。颜色卡可保存为clr文件供其他图使用。调整颜色条color scale属性，使其位置、大小合适。用划线和文字的方式添加图例，指示所计算的冲淤图的2个时间点。6、添加2个时间的等深线以contour map的形式分别添加两个白化后的网格数据，设置其属性为不填充，并设置显示线的颜色，标记，间隔等等。如图7所示。图7 等深线的属性设置，不填充，没有color scale，线的颜色和level的级别7、添加背景图如图8所示。图8 为冲淤图添加工程背景图图9 如果背景图范围大于冲淤图，则会弹出该对话框，选择否，使范围以冲淤图为准8、添加冲淤厚度值（此处是我会的方法，应该有其他方法？）冲淤厚度值就是冲淤变化的数值，即把数字尽量均匀的添加到冲淤变化图上。首先要生成dat格式的厚度值数据文件，选择Grid-Convert，将冲淤图的grd文件转换为dat文件，此时dat文件中所有白化区域的网格会填充无效值，如图10所示。编程将这些含有无效值的行去掉，只保留有效值。图10 冲淤厚度值，白化区域均填充为无效值，需要过滤掉处理完后，类似添加背景图一样，为冲淤图添加一个post layer，在其属性管理器中设置属性，如图11、图12所示。图11 设置数据列、显示密度、字体大小；Frequency最重要，控制显示密度，多试几次选最佳效果图12 设置数据列和字体、格式等属性9、遮挡添加的背景图、冲淤厚度值等可能会超出有效区域，此时