mesh文件生成

第三章MESH文件生成3.1MESH文件生成步骤第一步在MIKEZero中选择MeshGenerator。第二步选择投影带选择投影带主要用来计算地球自转引起的科氏力，当计算区域较小可以不考虑科氏力，此时可选择NON-UTM。国内常用的北京54及西安80均已经内置，用户直接选择即可。第三步设置模拟区域根据实际计算范围，设置工作范围。依次打开Options→Workspace打开相应界面进行设置。这一步需要指定计算区域左下角和右上角坐标。第四步导入背景图片（1）导入背景图片之前，需要将背景图片左下角和右上角在工作范围内的坐标。这一步可以用Tool中设置（如左图，MIKE2009之前版本没有此功能），也可以直接在图片所在文件夹中放置一个与图片同名的属性文件，用文本编辑器编辑里面的数据。更为简单的方法是在（2）步中导入图片后，点击旁边的edit来设置。打开顺序：MIKEZero→MikeZeroToolbox→FileConverter（2）完成以上设置后，通过Options→ImportGraphicLayers.点击下图箭头处的图标，新建一个背景图片。（3）导入后，点击工具栏中缩小按钮，可以显示图片。只需要设置x0,y0,x1,y1的值即可，gifw文件需要与背景图片同名。导入新建当不需要显示图片时，可在第（2）步中界面，选择OverlayManager（见下图）,去掉图片前的对号，即可隐藏图片。第五步描绘边界边界可以在底图上描绘，也可以直接导入数据确定边界。（1）在导入底图后，可以用工具栏中Drawacr工具进行描绘边界。边界描绘完成后，可能由于底图的原因，无法看到描绘的边界，此时可以按照第四步（3）中方法隐藏图片。（2）通过selectacr工具（）选择已经绘制的边界，右击选择RedistributedVertices,选择适当的距离对Vertices进行均匀化处理。第六步定义边界（1）将需要定义为开办界的线段两端的vertice(图中显示为红色的点)转换为node（图中显示为蓝色的点）。用工具选择vertice，右击选择Vertices→Nodes（2）转换后，用selectarc选择该边界，右击选择property，将表中数字改为大于0的数字，即可将该线段定义为开办界。第七步导入地形数据依次打开Data→ImportScatterData打开相应界面，点击Add可导入后缀为.xyz的地形文件。地形文件中数据依次为X，Y，Z坐标，可事先制作好。第八步生成网格（1）在导入地形数据后，依次点击Mesh→GenerateMesh打开网格设置界面，可以设置网格最大面积等相应参数。参数设置完后，点击Generate可生成网格。（2）点击Mesh→SmoothMesh可对网格进行平滑处理。第九步地形插值点击Mesh→Interpolate，对地形进行插值，保证每一个网格都有地形数据。第十步导出mesh点击Mesh→Exportmesh，输出mesh文件，供模型文件使用。注意文件名及路径中尽量不要有中文。3.2常用数据格式（1）边界文件格式3547.994484.421004022.694162.301004314.883541.011003987.782647.471003527.051886.561002577.661432.811001970.331572.431001698.081991.271001656.202717.281001572.433484.171001614.314071.551001872.604314.881002207.684484.421002744.204567.191003352.534581.151003547.994484.42000上述数据共分为5列，第一列为x坐标，第二列为y坐标，第三列为连接属性，第四列为z坐标，第五列为节点类型，其中前三列为必须数据，后两列可以省略。第三列一般为1，在连接处为0。（2）地形数据格式401374.82302624554.151014.4180401482.11102624627.833017.8190401369.05202624551.003014.4150401374.63002624554.045014.0690401378.25602624556.023014.3890401386.73202624560.645010.1910401491.09202624617.560017.7460401378.56202624556.190011.9320401456.77302624447.292011.8750401521.46402624484.087011.4660401442.01902624439.127014.4550401448.86102624442.913012.1810401550.78902624499.312012.5140401357.28902624544.588014.3330401542.93202624494.966011.0470401544.27102624494.3.706011.7430上述数据第一列为x坐标，第二列为y坐标，第三列为z坐标。在做糙率场时，第三列为（x,y）处的糙率。3.3局部加密在生成网格时，可以对某一部分区域网格进行局部加密，具体方法如下：（1）使用Drawarc工具将拟加密区域，描绘出来，形成封闭区域。（2）用工具栏中InsertPloygons工具，在拟加密区域中点绘一个点。拟加密现有边（3）用工具栏中SelectPloygons工具，选择在（2）中点绘的点，右击选择properties,选择applytriangularmesh，并勾选uselocalmaximumarea，指定加密区域最大划分面积。注意：若选择Excludefrommesh，则该区域作为固边界，不打网格。在水域中心有岛，桥墩等情况时，常采用此方法将其排除在计算区域之外。（4）设定好后，在进行全局打网格时，该区域将会被加密。若选择Excludefrommesh，则网格是：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-33度带中央子午线计算公式：当地经度/3=N；中央子午线L=3XN我国的经度范围西起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坐标系带区表示方法Beijing19543DegreeGKCM75E.prj表示三度分带法的北京54坐标系，中央经线在东75度的分带坐标，横坐标前不加带号；Beijing19543DegreeGKZone24.3.prj表示三度分带法的北京54坐标系，中央经线在东75度的分带坐标，横坐标前加带号；Beijing1954GKZone14.2.prj表示六度分带法的北京54坐标系，分带号为13，横坐标前加带号；Beijing1954GKZone13N.pr表示六度分带法的北京54坐标系，分带号为13，横坐标前不加带号。