file_MIKE21 mesh文件制作教程

第三章 MESH 文件生成 3.1 MESH 文件生成步骤 第一步 在 MIKE Zero 中选择 Mesh Generator。 第二步 选择投影带 选择投影带主要用来计算地球自转引起的科氏力，当计算区域较小可以不考 虑科氏力，此时可选择 NON-UTM。国内常用的北京 54 及西安 80 均已经内置， 用户直接选择即可。 第三步 设置模拟区域 根据实际计算范围，设置工作范围。依次打开 Options→ Workspace 打开相 应界面进行设置。这一步需要指定计算区域左下角和右上角坐标。 第四步 导入背景图片 （1）导入背景图片之前，需要将背景图片左下角和右上角在工作范围内的 坐标。这一步可以用 Tool 中设置（如左图，MIKE2009 之前版本没有此功能）， 也可以直接在图片所在文件夹中放置一个与图片同名的属性文件，用文本编辑器 编辑里面的数据。更为简单的方法是在（2）步中导入图片后，点击旁边的 edit 来设置。 打开顺序： MIKE Zero→Mike Zero Toolbox→ File Converter （2）完成以上设置后，通过 Options→ Import Graphic Layers .点击下图箭头 处的图标，新建一个背景图片。 （3）导入后，点击工具栏中缩小按钮 ，可以显示图片。 只需要设置 x0,y0,x1,y1的值即可， gifw 文件需要与背景图片同名。 导入 新建 当不需要显示图片时，可在第（2）步中界面，选择 Overlay Manager（见下 图）,去掉图片前的对号，即可隐藏图片。 第五步 描绘边界 边界可以在底图上描绘，也可以直接导入数据确定边界。 （1）在导入底图后，可以用工具栏中 Draw acr 工具 进行描绘边界。边 界描绘完成后，可能由于底图的原因，无法看到描绘的边界，此时可以按照第四 步（3）中方法隐藏图片。 （2）通过 select acr 工具（ ）选择已经绘制的边界，右击选择 Redistributed Vertices,选择适当的距离对 Vertices 进行均匀化处理。 第六步 定义边界 （1） 将需要定义为开办界的线段两端的 vertice(图中显示为红色 的点)转换为 node（图中显示为蓝色的点）。用工具 选择 vertice，右 击选择 Vertices→Nodes （2） 转换后，用 select arc 选择该边界，右击选择 property，将 中数字改为大于 0 的数字，即可将该线段定义为开办界。 第七步 导入地形数据 依次打开 Data→Import Scatter Data 打开相应界面，点击 Add 可导入后缀 为.xyz 的地形文件。地形文件中数据依次为 X，Y，Z 坐标，可事先制作好。 第八步 生成网格 （1）在导入地形数据后，依次点击 Mesh→Generate Mesh 打开网格设置界 面，可以设置网格最大面积等相应参数。参数设置完后，点击 Generate 可生成 网格。 （2）点击 Mesh→Smooth Mesh 可对网格进行平滑处理。 第九步 地形插值 点击 Mesh→Interpolate，对地形进行插值，保证每一个网格都有地形数据。 第十步 导出 mesh 点击 Mesh→Export mesh，输出 mesh 文件，供模型文件使用。注意文件名 及路径中尽量不要有中文。 3.2 常用数据格式 （1）边界文件格式 3547.99 4484.42 1 0 0 4022.69 4162.30 1 0 0 4314.88 3541.01 1 0 0 3987.78 2647.47 1 0 0 3527.05 1886.56 1 0 0 2577.66 1432.81 1 0 0 1970.33 1572.43 1 0 0 1698.08 1991.27 1 0 0 1656.20 2717.28 1 0 0 1572.43 3484.17 1 0 0 1614.31 4071.55 1 0 0 1872.60 4314.88 1 0 0 2207.68 4484.42 1 0 0 2744.20 4567.19 1 0 0 3352.53 4581.15 1 0 0 3547.99 4484.42 0 0 0 上述数据共分为 5 列，第一列为 x 坐标，第二列为 y 坐标，第三列为连接属 性，第四列为 z 坐标，第五列为节点类型，其中前三列为必须数据，后两列可以 省略。第三列一般为 1，在连接处为 0。 （2）地形数据格式 401374.8230 2624554.1510 14.4180 401482.1110 2624627.8330 17.8190 401369.0520 2624551.0030 14.4150 401374.6300 2624554.0450 14.0690 401378.2560 2624556.0230 14.3890 401386.7320 2624560.6450 10.1910 401491.0920 2624617.5600 17.7460 401378.5620 2624556.1900 11.9320 401456.7730 2624447.2920 11.8750 401521.4640 2624484.0870 11.4660 401442.0190 2624439.1270 14.4550 401448.8610 2624442.9130 12.1810 401550.7890 2624499.3120 12.5140 401357.2890 2624544.5880 14.3330 401542.9320 2624494.9660 11.0470 401544.2710 2624494.3.7060 11.7430 上述数据第一列为 x 坐标，第二列为 y 坐标，第三列为 z 坐标。在做糙率场 时，第三列为（x,y）处的糙率。 3.3 局部加密 在生成网格时，可以对某一部分区域网格进行局部加密，具体方法如下： （1） 使用 Draw arc 工具将拟加密区域，描绘出来，形成封闭区域。 （2） 用工具栏中 Insert Ploygons 工具，在拟加密区域中点绘一个 点。 拟加密 现有边 （3） 用工具栏中 Select Ploygons 工具，选择在（2）中点绘的点， 右击选择 properties,选择 apply triangular mesh，并勾选 use local maximum area，指定加密区域最大划分面积。 注意：若选择 Exclude from mesh，则该区域作为固边界，不打网格。在水 域中心有岛，桥墩等情况时，常采用此方法将其排除在计算区域之外。 （4） 设定好后，在进行全局打网格时，该区域将会被加密。 若选择 Exclude from mesh，则网格是： 3.4 北京 54 坐标投影选择 北京坐标系下，投影带的选择方法 1∶2.5 万及 1∶5 万的地形图采用 6 度分带投影，即经差为 6 度，从零度子 午线开始，自西向东每个经差 6 度为一投影带，全球共分 60 个带，用 1，2，3， 4，5，……表示．即东经 0～6 度为第一带，其中央经线的经度为东经 3 度，东 经 6～12 度为第二带，其中央经线的经度为 9 度…… 1∶1 万的地形图采用 3 度分带，从东经 1.5 度的经线开始，每隔 3 度为一带， 用 1，2，3，……表示，全球共划分 120 个投影带，即东经 1.5～ 4.5 度为第 1 带，其中央经线的经度为东经 3 度，东经 4.5～7.5 度为第 2 带，其中央经线的经 度为东经 6 度． 如何计算当地的中央子午线？ 当地中央子午线决定于当地的直角坐标系统，首先确定您的直角坐标系统是 3 度带还是 6 度带投影公式推算: 6 度带中央子午线计算公式：当地经度/6=N；中央子午线 L=6 * N （带号） 当没有除尽，N 有余数时， 中央子午线 L=6*N - 3 3 度带中央子午线计算公式： 当地经度/3=N；中央子午线 L=3 X N 我国的经度范围西起 73°东至 135°，可分成 六度带十一个（13 号带————23 号带），各带中央经线依次为（75°、 81°、……123°、129°、135°）； 三度带二十二 个（24 号带————45 号带）。各带中央经线依次为（72°、 75°、……132°、135°）； 六度带可用于中小比例尺（如 1：250000）测图，三度带可用于大比例尺（如 1：10000）测图，城建坐标多采用三度带的高斯投影 如何判断投影坐标是 3 度带坐标还是 6 度带坐标 如(4231898,21655933)其中 21 即为带号，同样所定义的东伪偏移值也需要加 上带号，如 21 带的东伪偏移值为 21500000 米。 假如你的工作区经度在 120 度 至 126 度范围，则该坐标系为 6 度带坐标系，该带的中央经度为 123 度。 如(2949320,36353822)其中 36 即为带号，已知该地点位于贵阳市附近，而从 地图上我们看到贵阳大概的经度是东经 108 度左右，因此可以 36*3=108，所以 该坐标系为 3 度带坐标系，该带的中央经度为 108 度。而不可能为 6 度带： 36*6=216。 MIKE 中北京 54 坐标系带区表示方法 Beijing 1954 3 Degree GK CM 75E.prj 表示三度分带法的北京 54 坐标系，中 央经线在东 75 度的分带坐标，横坐标前不加带号； Beijing 1954 3 Degree GK Zone 24.3.prj 表示三度分带法的北京 54 坐标系， 中央经线在东 75 度的分带坐标，横坐标前加带号； Beijing 1954 GK Zone 14.2.prj 表示六度分带法的北京 54 坐标系，分带号为 13，横坐标前加带号； Beijing 1954 GK Zone 13N.pr 表示六度分带法的北京 54 坐标系，分带号为 13，横坐标前不加带号。