资源简介
本人写的高斯空气扩散模型,对于了解arcgis 编写Python 模型的有一定的帮助。
代码片段和文件信息
# ---------------------------------------------------------------------------
# test.py
# Created on: 星期五 一月 07 2011 01:14:21 下午
# 高斯空气扩散模型
# Author:黄兆明
# ---------------------------------------------------------------------------
import win32com.clientsys string os arcgisscriptingmathfileinput
#空气污染类型对应类型,PType:污染物类型;I:污染指数
def GetType(PType I):
if PType == “TSP“ :
if I == 50 :
return 0.120
elif I == 100 :
return 0.300
elif I == 200 :
return 0.500
elif I == 300 :
return 0.625
elif I == 400 :
return 0.875
else :
return 1.000
if PType == “SO2“ :
if I == 50 :
return 0.050
elif I == 100 :
return 0.150
elif I == 200 :
return 0.250
elif I == 300 :
return 1.600
elif I == 400 :
return 2.100
else :
return 2.620
if PType == “NOx“ :
if I == 50 :
return 0.050
elif I == 100 :
return 0.100
elif I == 200 :
return 0.150
elif I == 300 :
return 0.565
elif I == 400 :
return 0.750
else :
return 0.940
#获取稳定度对应Qy,x为距离,Q大气稳定度
def GetQy(xQ) :
if Q==“A“ :
if x<=1000:
a=0.901074
r=0.425809
else :
a=0.850934
r=0.602052
elif Q==“B“ :
if x<=1000:
a=0.914370
r=0.281846
else :
a=0.865014
r=0.396353
elif Q==“B-C“ :
if x<=1000:
a=0.919325
r=0.229500
else :
a=0.875086
r=0.314238
elif Q==“C“ :
if x<=1000:
a=0.924279
r=0.177154
else :
a=0.885157
r=0.232123
elif Q==“C-D“ :
if x<=1000:
a=0.926849
r=0.143940
else :
a=0.886940
r=0.189396
elif Q==“D“ :
if x<=1000:
a=0.929418
r=0.110726
else :
a=0.888723
r=0.146669
elif Q==“D-E“ :
if x<=1000:
a=0.925118
r=0.0985631
else :
a=0.892794
r=0.124308
elif Q==“E“ :
if x<=1000:
a=0.920818
r=0.0864001
else :
a=0.896864
r=0.101947
else :
if x<=1000:
a=0.929418
r=0.0553634
else :
a=0.888723
r=0.0733348
return r*pow(xa)
#获取稳定度对应Qz,x为距离,Q大气稳定度
def GetQz(xQ) :
if Q==“A“ :
if x<=300:
a=1.12154
r=0.0799904
elif x>300 and x<=500:
a=1.51360
r=0.00854771
else :
a=2.10881
r=0.000211545
elif Q==“B“ :
if x<=500:
a=0.964435
r=0.127190
else :
a=0.109356
r=0.057025
elif Q==“B-C“ :
if x<=500:
a=0.941015
r=0.114682
else :
a=1.00770
r=0.0757182
elif Q==“C“ :
a=0.917595
r=0.106803
elif Q==“C-D“ :
if x<=2000:
a=0.838628
r=0.126152
elif x>2000 and x<=10000 :
a=0.756410
r=0.235667
else :
a=0.815575
r=0.136659
elif Q==“D“ :
if x<=1000:
a=0.826212
r=0.104634
elif x>1000 and x<=10000 :
a=0.632023
r=0.400167
else :
a=0.55536
r=0.810763
elif Q==“D-E“ :
if x<=2000:
a=0.776864
r=0.111771
elif x>2000 and x<=10000 :
a=0.572347
r=0.5289922
else :
a=0.499149
r
相关资源
- python3+Tkinter+GUI界面+pyserial+串口+串口
- Python.CLF
- python实现2048游戏
- 大学排名爬虫
- annoy-1.15.2-cp36-cp36m-win_amd64.whl
- python矩阵乘法小程序
- python 音乐播放器
- Python分词系统jieba代码有注释,无报错
- Python生成3D图,饼图,合图,散点图,
- 树莓派目标跟踪代码python
- 简单的Python蚁群算法
- Python爬虫教程千万别错过
- get_wuyou.zip
- python入门到精通某智付费课程超级完
- 用python代码把背景颜色过滤掉
- python所有源代码
- PyHook3-1.6.1-cp36-cp36m-win32.whl
- L2正则化python实现案例()
- tensorflow手写数字识别python源码案例
- 批量excel转shp(面)python
- python ds evidence theory code
- Python串口调试助手
- 推荐系统物质扩散代码python
- 基于python的行人重识别代码
- 基于Python语言的SAD算法进行双目立体
- Python 聊天室 客户端和服务端 聊天
- python随机森林应用
- 肺实质分割python代码
- 基于python2.7的LeNet5源代码实现
- caffe_log绘制accuracy和loss曲线python3
评论
共有 条评论