EFDC从零开始|01制作EFDC所需网格文件(ArcGIS和Python)

最新推荐文章于 2025-02-12 15:11:42 发布
最新推荐文章于 2025-02-12 15:11:42 发布
阅读量3.8k 收藏 28
点赞数 10
分类专栏： 水文模型-EFDC 文章标签： python arcgis
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/smilesnailGIS/article/details/120688476
版权
该博客介绍了如何使用ArcGIS将矢量地图转换为栅格，并通过Python处理栅格文件以生成EFDC所需的cell.inp、dxdy.inp和lxly.inp文件。涉及的操作包括数据格式转换、特定值替换以及文件整理，为水动力学模型的预处理提供了详细步骤。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

一. EFDC网格文件的格式说明

（这部分以后有时间再补充）

二. ArcGIS操作

1.矢量地图转栅格

给shp文件增加一列属性，命名为value，值为5；

 

工具：Conversion Tools - To Raster - Polygon to Raster

将Cellsize设置为你需要的栅格大小，我设置了500m。

 before&after

 备注：我所使用的坐标系

投影坐标系：Krasovsky_1940_Transverse_Mercator

地理坐标系：GCS_Krasovsky_1940

2.栅格导出为TXT

工具：Conversion Tools - From Raster - Taster to ASCII

 打开taihu500.txt查看一下数据格式，局部如图所示。NoData的部分为-9999，与EFDC不同。

三. 制作EFDC所需文件（Python）

1. cell.inp

①手动删除前六行。提前记录一下，后面可能会用到。

 ②将5和-9999交界处（包括角）的-9999替换为9，其它-9999替换为9

代码如下：

f = open("taihu500.txt", "r")
f2 = open("cell-init.txt", "wt")

lines = f.read()
data = lines.split('\n')
row = len(data)
for i in range(1, row-2):
    line0 = data[i-1].split()
    line = data[i].split()
    line1 = data[i+1].split()
    col = len(line)
    for j in range(1, col-1):
        '''
        if line[j].strip('\n') == '-9999' and line[j+1].strip('\n') == '5':
            line[j] = '9'
        elif line[j].strip('\n') == '-9999' and line[j-1].strip('\n') == '5':
            line[j] = '9'
        elif line[j].strip('\n') == '-9999' and line0[j].strip('\n') == '5':
            line[j] = '9'
        elif line[j].strip('\n') == '-9999' and line1[j].strip('\n') == '5':
            line[j] = '9'
        elif line[j].strip('\n') == '-9999':
            line[j] = '0'
        '''
        if line[j].strip('\n') == '-9999':
            if line[j+1].strip('\n') == '5':
                line[j] = '9'
            if line[j-1].strip('\n') == '5':
                line[j] = '9'
            if line0[j].strip('\n') == '5':
                line[j] = '9'
            if line1[j].strip('\n') == '5':
                line[j] = '9'
            # corner
            if line0[j - 1].strip('\n') == '5':
                line[j] = '9'
            if line0[j + 1].strip('\n') == '5':
                line[j] = '9'
            if line1[j - 1].strip('\n') == '5':
                line[j] = '9'
            if line1[j + 1].strip('\n') == '5':
                line[j] = '9'
    for j in range(1, col-1):
        if line[j].strip('\n') == '-9999':
            line[j] = '0'
        f2.write(str(line[j]))
        f2.write(' ')
    f2.write('\n')

f2.close()
f.close()

 结果，局部如图所示：（顶部手动加了两行0缓冲，不为什么，想加）

 

 整体如图所示：注意不要有空行（像我图里这样最后有一行空行是错的，手动检查一下就行）

 

 

 

 ③整理成cell.inp格式

代码：

f = open("cell-init.txt", "r")
f2 = open("input_cell.txt", "wt")

f2.write('C Cell.inp\n')
f2.write(
    'C    0        1         2         3         4         5         6         7         8         9        10        11        12        13        14        15        16        17        18        19        20        21\n')
f2.write(
    'C    12345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012\n')
f2.write('C\n')

lines = f.read()
data = lines.split('\n')
row = len(data)
print(row)
for i in range(0, row):
    # print(i)
    f2.write('%3d  ' % (row - i))
    line = data[i].split()
    col = len(line)
    for j in range(0, col):
        f2.write(str(line[j]))

    f2.write('  %d  ' % (row - i))
    f2.write('\n')

f2.write('C\n')
f2.write(
    'C    12345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012345678901234567890123456789012\n')
f2.write(
    'C    0        1         2         3         4         5         6         7         8         9        10        11        12        13        14        15        16        17        18        19        20        21\n')
f2.close()
f.close()

 结果如图所示：

 

2. dxdy.inp

DX和DY是Cellsize = 500；

DEPTH，BOTTOM ELEV等目前写成了固定的，其具体赋值方法有待探索。

代码：

f = open("cell-init.txt", "r")
f2 = open("input_dxdy.txt", "wt")

f2.write('C dxdy.inp file, in free format across columns\n')
f2.write('C\n')
f2.write('C    I      J      DX      DY      DEPTH      BOTTOM ELEV      ZROUTH      VEG TYPE\n')
f2.write('C\n')

lines = f.read()
data = lines.split('\n')
row = len(data)
for j in range(0, row):  # Y方向-行-文件第二列
    line = data[row-1-j].split()
    col = len(line)
    for i in range(0, col):  # X方向-列-文件第一列
        if line[i].strip('\n') == '5':
            f2.write('   %3d   %3d    500.0     500.0     5.0     -5.0     0.02     0\n' % (i+1, j+1))

f2.write('C\n')
f2.write('C  I         ARRAY INDEX IN X DIRECTION\n')
f2.write('C  J         ARRAY INDEX IN Y DIRECTION \n')
f2.write('C  DX        CELL DIMENSION IN X DIRECTION, METERS \n')
f2.write('C  DY        CELL DIMENSION IN Y DIRECTION, METERS\n')
f2.write('C  DEPTH     INITIAL WATER DEPTH, METERS \n')
f2.write('C  BOTTOM    ELEV  BOTTOM BED ELEVATION, METERS\n')
f2.write('C  ZROUGH    LOG LAW ROUGHNESS HEIGHT, ZO, METERS \n')
f2.write('C  VEG TYPE  VEGETATION TYPE CLASS, INTEGER VALUE\n')
f2.close()

 结果如下图所示：

 

3. dxdy.inp

注意XLNUTME和YLNUTME的计算（有待考证）；

CCUE等目前写成了固定的，其具体赋值方法有待探索。

代码：

f = open("cell-init.txt", "r")
f2 = open("input_lxly.txt", "wt")

f2.write('C lxly.inp file, in free format across columns\n')
f2.write('C\n')
f2.write('C    I   J    XLNUTME    YLNUTME   CCUE    CCVE    CCUN    CCVN    TMPVAL\n')
f2.write('C\n')

lines = f.read()
data = lines.split('\n')
row = len(data)
for j in range(0, row):  # Y方向-行-文件第二列
    line = data[row-1-j].split()
    col = len(line)
    for i in range(0, col):  # X方向-列-文件第一列
        if line[i].strip('\n') == '5':
            f2.write('   %3d   %3d   %7.1f   %7.1f     1.0     0.0     0.0     1.0     -0.5000000\n' % (i+1, j+1, (i+1)*500+250, (j+1)*500+250))
            #  有待讨论：上一句有错误，实际情况中不应该+250，还要结合坐标系原点

f2.write('C\n')    
f2.write('C  I         ARRAY INDEX IN X DIRECTION\n')
f2.write('C  J         ARRAY INDEX IN Y DIRECTION \n')
f2.write('C  XLNUTME   X CELL CENTER COORDINATE, LONGITUDE, METERS,OR KM \n')
f2.write('C  YLTUTMN   Y CELL CENTER COORDINATE, LONGITUDE, METERS,OR KM  \n')
f2.write('C  CCUE      ROTATION MATRIX COMPONENT \n')
f2.write('C  CCVE      ROTATION MATRIX COMPONENT \n')
f2.write('C  CCUN      ROTATION MATRIX COMPONENT \n')
f2.write('C  CCVN      ROTATION MATRIX COMPONENT \n')   
f2.close()

 四. 总结

最后生成了3个网格文件

修改文件名和后缀

 

 

EFDC视频教程.rar
11-09
efdc教学视频，有二十几个，喜欢的赶紧下载
EFDC软件操作步骤
04-30
EFDC软件操作具体步骤详解，从小白到上手的最快途径，帮助初学者快速掌握EFDC建模
EEMS 12.2(EFDC+ Explorer 12.2.0 and Grid+ 1.2)
最新发布
qq_62821232的博客
02-12 521
在开发过程中，我们改进了我们的建模系统，提高了处理速度、更多的分析工具、新的子模型和错误修复。在 EE 之前，如果要使用任何版本的 EFDC 创建流体动力学模型，则需要一到两年的时间来创建、编辑、运行和校准模型。现在，借助 EE 的输入工具、计算效率、一致性检查、改进的处理速度和整体可用性改进，您可以在短短一到两个月内获得一个可行的复杂模型！一种用于流体动力学模型的网格生成工具，使构建复杂的二维曲线正交网格的过程变得快速、直观且足够强大，即使是最具挑战性的水体也是如此。1、快速构建复杂的网格。
EFDC水模型 初学者入门 及软件下载学习指导
weixin_34008933的博客
06-29 2845
EFDC模型有三个不同的版本，它们分别代表了计算机不同的操作菜单和功能。其中的两个版本功能齐全，包括了EFDC 所有的产品版本。另外一个是简化版，涵盖了水动力学（包括盐度和温度的影响），毒物，沉积物的3-time level 计算方法。这一简化版称为EFDC水力学版，主要用于研究地区IV，两个完全版继续由Tera Tech, Inc（中国最大的国外企业级软件供应商）和Dyn...
EFDC建模操作流程
04-30
EFDC软件操作具体步骤详解，从小白到上手的最快途径，帮助初学者快速掌握EFDC建模
环境流体动力学模型(EFDC_Explore)中文手册
05-08
环境流体动力学模型(EFDC_Explore)中文手册
EFDC_Explorer8.1 中文使用手册.pdf
05-29
EFDC_Explorer8.1 中文使用手册 用于使用efdc模型的人使用和参考，方便英文不好的人学习使用。
浓度文件转换程序.rar_EFDC_plt文件_浓度文件转换程序
09-20
总的来说，这个“浓度文件转换程序”是一个专门设计来解决环境科学和工程领域中特定数据格式问题的工具，它帮助用户将EFDC模型的输出转换为TechPlot可以方便处理的格式，从而促进更高效的数据分析和可视化工作流程。...
EFDC8.3免费版
02-11
这个文件是Windows操作系统下的可执行文件，用户可以通过运行此文件来安装EFDC模型。在安装过程中，系统通常会提示用户选择安装路径、接受许可协议，并可能需要配置一些基本的设置，如模拟的水体类型、输出选项等。 ...
EFDC控制输入文件
11-17
`Isresti`用于在模拟开始时从`restart.inp`读取初始条件，`Isresto`则决定在何时输出`restart.out`文件，可以按时间周期或在模拟结束时输出。 `ISLOG`参数控制日志文件`EFDC.LOG`的生成，当设置为1时，会记录模型...
EFDC_Explore（动力学模型模拟）中文使用手册.rar
08-25
EFDC_Explore（动力学模型模拟）中文使用手册，EFDC_Explorer是微软基于Windows™的预处理程序和后处理程序开发的三维环境流体动力学程序，最初由John Hamrick（1996）年设计形成。EFDC由美国环境保护署资助（EPA）。EFDC_Explorer是一个基于微软视窗Windows、由Dynamic Solutions-International, LLC 公司为环境流体动力学软件EFDC开发的前处理和后处理计算机图形界面软件。
EFDC建模操作流程注意事项
04-30
EFDC软件操作具体步骤详解，从小白到上手的最快途径，帮助初学者快速掌握EFDC建模
著名的水质模拟软件，EFDC.本人亲测安装成功。
11-15
网上下了好多EFDC的程序，不是试用版就是不能用。这个是可以用的，但是仅限于一维。二维已经商业化了。
EFDC（The Environmental Fluid Dynamics Code）模型源码
07-09
EFDC（The Environmental Fluid Dynamics Code）模型是由威廉玛丽大学维吉尼亚海洋科学研究所（VIMS，Virginia Institute of Marine Science at the College of William and Mary）的John Hamrick等人开发的三维地表水水质数学模型，可实现河流、湖泊、水库、湿地系统、河口和海洋等水体的水动力学和水质模拟，是一个多参数有限差分模型。经过近20年的发展和完善，目前该模型已在大学，政府机关和环境咨询公司等组织中被广泛使用，并成功用于美国和欧洲其他国家100多个水体区域的研究，在我国已被应用于云南滇池水质模拟，重庆两江汇流水动力模拟、密云水库营养物模拟等以及内蒙古乌梁素海地区水体富营养化模拟等。[1] 该模型系统包括水动力、泥沙、有毒物质、水质、底质、风浪等模块，模拟计算过程中首先完成流场计算，获得三维流速场的时空分布特征，在此基础上计算泥沙迁移、冲淤作用，进而模拟受粘性泥沙吸附影响的各水质变量动态变化过程。为更好的拟合研究区地形条件，模型在水平方向除可采用传统的 直角坐标外还可在水平向使用正交曲线坐标，垂直方向采用σ坐标。 EFDC水动力学模块可计算如下内容：流速，示踪剂，温度，盐度，近岸羽流和漂流。水动力学模型输出变量可直接与水质，底泥迁移和毒性物质等模块耦合，作为物质运移的驱动条件。同时EFDC也提供了与WASP等软件的接口，输出可供水质模拟使用的.HYD文件。EFDC泥沙模块可进行多组分泥沙的模拟，根据在水体里面的迁移特征把泥沙分为悬移质和推移质；悬移质根据粒径大小分为粘性泥沙和非粘性泥沙，进而还可细分为若干组。可根据物理或经验模型模拟泥沙的沉降、沉积、冲刷及再悬浮等过程。EFDC有毒污染物模块可以模拟各类型污染物在水体中的迁移转化过程，该模块需要研究者针对特定有毒污染物提供具体反应过程设定反应系数。EFDC的水质模块，主要模拟水体中以藻类生长为中心的各变量间相互关系。而底质模块模拟沉积物与水体之间的物质交换过程。
EFDC软件操作步骤指南
gitblog_06691的博客
10-31 729
EFDC软件操作步骤指南 【下载地址】EFDC软件操作步骤指南分享 EFDC软件操作步骤指南欢迎来到EFDC（Environmental Fluid Dynamics Code）软件操作详细步骤指南 项目地址: https://g...
EFDC从零开始|00前言（一些废话）
smilesnailGIS的博客
10-12 1178
看了很多资料，还是很凌乱 买不起商业版模型软件，只能靠开源 找到了很多代码，没一个能用起来 装软件，学Fortran 好不容易用起来了，运行时间贼长 道阻且长 道阻且跻 我自己也比较菜吧，代码只会写一些简单的 记忆力还超级差，过段时间不看就忘了之前是怎么做的了 反反复复在原地打转 所以把这个过程记录下来 一是好记性不如烂笔头 二是也许会对别人有一些帮助 三是，如果有大神看到我在苦苦挣扎，如果能指点一下我将不胜感激 ...
EFDC建模方法及在地表水环境评价、水源地划分、排污口论证中的应用
PhyliciaFelicia的博客
05-11 694
EFDC模型、水质模型、水动力、地表水环评、地下水、水质、水环境、MIKE、Delft3D、Qual2K、水源地划分、排污口论证
评论 6
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值