CN109375285A - 一种应用于污染分析的玫瑰图制作方法 - Google Patents
一种应用于污染分析的玫瑰图制作方法 Download PDFInfo
- Publication number
- CN109375285A CN109375285A CN201810949653.2A CN201810949653A CN109375285A CN 109375285 A CN109375285 A CN 109375285A CN 201810949653 A CN201810949653 A CN 201810949653A CN 109375285 A CN109375285 A CN 109375285A
- Authority
- CN
- China
- Prior art keywords
- data
- compass azimuth
- rose
- time point
- wind
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Medicinal Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Physics & Mathematics (AREA)
- Food Science & Technology (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
Abstract
一种应用于污染分析的玫瑰图制作方法,包括以下步骤:选取预定地点的数据获取时间段,在所述时间段内选定若干测试时间点;获取每一时间点的风向数据、风速数据以及污染物浓度数据;以预定地点为中心,制作罗盘方位图;以罗盘方位图的中心为圆心,制作风速线性差值图;依据每一时间点的风向数据在罗盘方位图中生成风频廓线;依据所述时间点的风速数据和风向数据将每一时间点对应显示到所述罗盘方位图中;将污染物浓度映射成颜色深度;将与每一时间点的污染物浓度映射的颜色叠加到所述罗盘方位图中。由于在罗盘方位图中集成了污染物浓度以及风速数据,从而形象的表示了污染物传播方向与风速之间的关系。
Description
技术领域
本发明涉及气象分析领域,特别涉及一种应用于污染分析的玫瑰图制作方法。
背景技术
近年来,随着经济社会的持续快速发展,工业化和城市化的发展加剧,能源消耗迅速增加,空气污染日益严重,空气污染物排放总量居高不下。作为影响污染物传播的风一直是人们研究污染传播的主要手段。
风频玫瑰图(简称风玫瑰图)也叫风向频率玫瑰图,它是根据某一地区多年平均统计的各个风向和风速的百分数值,并按一定比例绘制,一般多用8个或16个罗盘方位表示,由于形状酷似玫瑰花朵而得名。
由于风频玫瑰图能够形象的表示风的风向频率,所以现有技术中,研究人员经常依据风频玫瑰图分析污染物的传播。但是其在研究过程中,由于风频玫瑰图并不能够直接反应与污染物之间的关系,所以需要研究人员对污染物的传播情况另做研究,用以配合污染物传播的预警。
发明内容
本发明的目的是提供一种应用于污染分析的玫瑰图制作方法,其优点是该方法能够制作出形象表现污染物传播的玫瑰图。
本发明的上述技术目的是通过以下技术方案得以实现的:
一种应用于污染分析的玫瑰图制作方法,包括以下步骤:
步骤1:选取预定地点的数据获取时间段,在所述时间段内选定若干测试时间点;
步骤2:获取每一时间点的风向数据、风速数据以及污染物浓度数据;
步骤3:以预定地点为中心,制作罗盘方位图;
步骤4:以罗盘方位图的中心为圆心,制作风速线性差值图;
步骤5:依据每一时间点的风向数据在罗盘方位图中生成风频廓线;
步骤6:依据所述时间点的风速数据和风向数据将每一时间点对应显示到所述罗盘方位图中;
步骤7:将污染物浓度映射成颜色深度;
步骤8:将与每一时间点的污染物浓度映射的颜色叠加到所述罗盘方位图中。
通过采用上述技术方案,风频廓线的生成表示为了在选定时间段内的主要风向,将污染物浓度映射成颜色深度后叠加到罗盘方位图中,实现了污染物浓度与风向风频数据的混合显示,由于在罗盘方位图中不但显示了风向风频数据,还集成了污染物浓度以及风速数据,从而形象的表示了污染物传播方向与风速之间的关系,在研究人员研究风对污染物传播的影响过程中,会提供形象而有效的数据支撑。当经过该方法制作的罗盘方位图表示未来一一段时间内,某地区会出现强污染天气时,可通过研究人员进行提前预警。
优选的,还包括步骤9:对罗盘方位图中的颜色做线性插值变换。
通过采用上述技术方案,由于单纯的通过选定的时间点的污染浓度和风速进行表示,会需要很多的时间点的数据支撑,不但增加了数据的获取难度,同时也会增加***的数据处理难度。通过采用线性插值变换,可在获取的时间点数据较少的情况下,形象的显示出罗盘方位图中各个点所代表的风向与浓度,从而方便研究人员研究罗盘方位图。
优选的,所述线性插值变换设置为双线性插值法。
优选的,所述罗盘方位图为八罗盘方位图。
优选的,所述罗盘方位图为十六罗盘方位图。
优选的,还包括有表示污染物浓度与颜色深度映射关系的浓度映射图。
通过采用上述技术方案,由于将颜色与污染物浓度一一映射,通过颜色的显示即可实现对污染物浓度分布的形象表示,但是在实际研究过程中,若研究人员并不清楚颜色的深度与污染物浓度之间的映射关系,就需要与技术人员沟通,或是不同的罗盘方位图中污染物浓度与颜色深度之间的映射关系不同,则会影响到研究人员对污染物浓度的判断,通过增加浓度映射图,可方便研究人员能够清楚的判断出罗盘方位图中的浓度分布。
综上所述,本发明具有以下有益效果:由于在罗盘方位图中不但显示了风向风频数据,还集成了污染物浓度以及风速数据,从而形象的表示了污染物传播方向与风速之间的关系,在研究人员研究风对污染物传播的影响过程中,会提供形象而有效的数据支撑。
附图说明
图1是应用于污染分析的玫瑰图制作方法的流程图;
图2是污染玫瑰图的示意图。
图中,1、罗盘方位图;2、风速线性差值图;3、风频廓线;4、浓度映射图。
具体实施方式
以下结合附图对本发明作进一步详细说明,其中相同的零部件用相同的附图标记表示。需要说明的是,下面描述中使用的词语“前”、“后”、“左”、“右”、“上”和“下”、“底面”和“顶面”指的是附图中的方向,词语“内”和“外”分别指的是朝向或远离特定部件几何中心的方向。
一种应用于污染分析的玫瑰图制作方法,如图1和图2所示,包括以下步骤:
步骤1:选定时间点和预定地点,选取预定地点的数据获取时间段,在所述时间段内选定若干测试时间点。
步骤2:获取时间点的气象数据和污染数据,气象数据包括风速数据和风速数据,污染数据包括污染物浓度数据。
步骤2-1:设置传感器,在选定地点设置污染物传感器,污染物传感器至少包括PM2.5传感器和二氧化硫传感器,涉及到空气污染物浓度检测的传感器均在本实施例的说明范围内,在此不做唯一限定。
步骤2-2:获取监测数据,通过数据中心采集各污染物传感器的检测数据,并上传至云数据库。
步骤2-3:获取气象数据,通过数据中心获取NOAA服务器中的gdas气象数据,并从gdas气象数据中抽取相应时间点的风向数据,将风向数据上传至云数据库。
步骤3:以预定地点为中心,制作罗盘方位图1,罗盘方位图1优选为八罗盘方位图1,但是对于方位表示更加清晰的十六罗盘方位图1也在本实施例的说明范围内。数据中心连接有若干研究终端,研究终端优选为设置在电脑或APP上的气象分析软件,在研究终端制作罗盘方位图1。
步骤4:以罗盘方位图1的中心为圆心,制作风速线性差值图2。
步骤5:生成风频廓线3,依据每一时间点的风向数据在罗盘方位图1中生成风频廓线3。
步骤6:罗盘布点,依据时间点的风速数据和风向数据将每一时间点对应显示到所述罗盘方位图1中。
步骤7:将污染物浓度映射成颜色深度。
步骤8:将与每一时间点的污染物浓度映射的颜色叠加到罗盘方位图1中。
步骤9:对罗盘方位图1中的颜色做线性插值变换,生成污染玫瑰图。
步骤10:添加浓度映射图4,将表示污染物浓度与颜色深度映射关系的浓度映射图4添加到污染玫瑰图中。
进一步的,由于污染物传感器检测到的监测数据上传至云数据库,实现了对监测数据的历史记录,在上述研究终端设置时间选项,使得研究人员直接通过研究终端选定时间段,从而生成相应时间段的污染玫瑰图。
进一步的,上述线性插值变换设置为双线性插值法,具体内容如下:
选定颜色添加位置(u0,v0),并选定颜色添加区域的4个邻点灰度值,4个邻点灰度值为相应时间点的灰度值,分别命名为(u,v)、(u+1,v)、(u,v+1)和(u+1,v+1)。
假定[S]表示不超过S的最大整数,设定:u=[u0]; v=[v0];α= u0-[u0]; β= v0-[v0] 。
S1:由f(u,v)、f(u+1,v)做第一次水平方向的插值计算,求得(u0,v)。
f(u0,v)= f(u,v)+α[f(u+1,v)- f(u,v)]
S2: 由f(u,v+1)、f(u+1,v+1)做第二次水平方向的插值计算,求得(u0,v+1)。
f(u0,v+1)= f(u,v+1)+α[f(u+1,v+1)- f(u,v+1)]
S3: 由(u0,v)和(u0,v+1)做第三次垂直方向的插值计算,求得(u0,v0)。
f(u0,v0)= f(u0,v)+β[f(u0,v+1)- f(u0,v)]
= f(u,v)(1-α)(1-β)+ f(u+1,v)α(1-β)
+ f(u,v+1)(1-α)β+ f(u+1,v+1)αβ
从而得到(u0,v0)的颜色。
进一步的,浓度映射图4优选为条形图,包括颜色由浅到深的颜色渐变条形图以及与颜色条形图中颜色深度对应的浓度标号(图中未示出)。
采用以上技术方案制作出的污染玫瑰图,不但能够通过风频廓线3表示风向,同时能够通过没一点的污染物浓度以及风向表示出风对污染物传播的影响,方便研究人员对污染物传播的研究和预算,对环境预警具有明显的帮助。
本具体实施例仅仅是对本发明的解释,其并不是对本发明的限制,本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改,但只要在本发明的权利要求范围内都受到专利法的保护。
Claims (6)
1.一种应用于污染分析的玫瑰图制作方法,其特征在于:包括以下步骤:
步骤1:选取预定地点的数据获取时间段,在所述时间段内选定若干测试时间点;
步骤2:获取每一时间点的风向数据、风速数据以及污染物浓度数据;
步骤3:以预定地点为中心,制作罗盘方位图(1);
步骤4:以罗盘方位图(1)的中心为圆心,制作风速线性差值图(2);
步骤5:依据每一时间点的风向数据在罗盘方位图(1)中生成风频廓线(3);
步骤6:依据所述时间点的风速数据和风向数据将每一时间点对应显示到所述罗盘方位图(1)中;
步骤7:将污染物浓度映射成颜色深度;
步骤8:将与每一时间点的污染物浓度映射的颜色叠加到所述罗盘方位图(1)中。
2.根据权利要求1所述的一种应用于污染分析的玫瑰图制作方法,其特征在于:还包括步骤9:对罗盘方位图(1)中的颜色做线性插值变换。
3.根据权利要求2所述的一种应用于污染分析的玫瑰图制作方法,其特征在于:所述线性插值变换设置为双线性插值法。
4.根据权利要求3所述的一种应用于污染分析的玫瑰图制作方法,其特征在于:所述罗盘方位图(1)为八罗盘方位图(1)。
5.根据权利要求3所述的一种应用于污染分析的玫瑰图制作方法,其特征在于:所述罗盘方位图(1)为十六罗盘方位图(1)。
6.根据权利要求3所述的一种应用于污染分析的玫瑰图制作方法,其特征在于:还包括有表示污染物浓度与颜色深度映射关系的浓度映射图(4)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810949653.2A CN109375285A (zh) | 2018-08-20 | 2018-08-20 | 一种应用于污染分析的玫瑰图制作方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810949653.2A CN109375285A (zh) | 2018-08-20 | 2018-08-20 | 一种应用于污染分析的玫瑰图制作方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109375285A true CN109375285A (zh) | 2019-02-22 |
Family
ID=65403733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810949653.2A Pending CN109375285A (zh) | 2018-08-20 | 2018-08-20 | 一种应用于污染分析的玫瑰图制作方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109375285A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109991683A (zh) * | 2019-04-12 | 2019-07-09 | 中国气象局沈阳大气环境研究所 | 风频数据处理方法及装置 |
CN110610279A (zh) * | 2019-09-27 | 2019-12-24 | 复旦大学 | 一种大气细颗粒物污染源识别方法及其应用 |
CN110687255A (zh) * | 2019-10-21 | 2020-01-14 | 软通动力信息技术有限公司 | 一种空气污染物的溯源方法、装置、设备和存储介质 |
CN110750257A (zh) * | 2019-09-28 | 2020-02-04 | 天津同阳科技发展有限公司 | 一种玫瑰图显示方法及*** |
CN110954658A (zh) * | 2019-12-02 | 2020-04-03 | 上海宝钢工业技术服务有限公司 | 工业园区挥发性有机物网格化溯源调查的方法 |
CN112051362A (zh) * | 2020-09-01 | 2020-12-08 | 杭州壹丈科技有限公司 | 风玫瑰预判污染物变化模型 |
CN113759441A (zh) * | 2021-09-08 | 2021-12-07 | 长春嘉诚信息技术股份有限公司 | 基于风玫瑰图及污染源监测的空气质量溯源方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08166438A (ja) * | 1994-12-15 | 1996-06-25 | Kaijo Corp | 風速測定用ドップラーレーダー |
CN102426399A (zh) * | 2011-09-01 | 2012-04-25 | 南京信息工程大学 | 高速公路横风预警方法 |
CN107656298A (zh) * | 2017-09-15 | 2018-02-02 | 深圳市卡普瑞环境科技有限公司 | 一种获取方位信息的方法及相关设备 |
CN108182491A (zh) * | 2017-12-27 | 2018-06-19 | 宇星科技发展(深圳)有限公司 | 大气细颗粒物(pm2.5)实时源解析定位方法 |
CN108416521A (zh) * | 2018-03-06 | 2018-08-17 | 深圳市卡普瑞环境科技有限公司 | 一种信息来源确定的方法以及信息来源确定装置 |
-
2018
- 2018-08-20 CN CN201810949653.2A patent/CN109375285A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08166438A (ja) * | 1994-12-15 | 1996-06-25 | Kaijo Corp | 風速測定用ドップラーレーダー |
CN102426399A (zh) * | 2011-09-01 | 2012-04-25 | 南京信息工程大学 | 高速公路横风预警方法 |
CN107656298A (zh) * | 2017-09-15 | 2018-02-02 | 深圳市卡普瑞环境科技有限公司 | 一种获取方位信息的方法及相关设备 |
CN108182491A (zh) * | 2017-12-27 | 2018-06-19 | 宇星科技发展(深圳)有限公司 | 大气细颗粒物(pm2.5)实时源解析定位方法 |
CN108416521A (zh) * | 2018-03-06 | 2018-08-17 | 深圳市卡普瑞环境科技有限公司 | 一种信息来源确定的方法以及信息来源确定装置 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109991683A (zh) * | 2019-04-12 | 2019-07-09 | 中国气象局沈阳大气环境研究所 | 风频数据处理方法及装置 |
CN110610279A (zh) * | 2019-09-27 | 2019-12-24 | 复旦大学 | 一种大气细颗粒物污染源识别方法及其应用 |
CN110750257A (zh) * | 2019-09-28 | 2020-02-04 | 天津同阳科技发展有限公司 | 一种玫瑰图显示方法及*** |
CN110687255A (zh) * | 2019-10-21 | 2020-01-14 | 软通动力信息技术有限公司 | 一种空气污染物的溯源方法、装置、设备和存储介质 |
CN110954658A (zh) * | 2019-12-02 | 2020-04-03 | 上海宝钢工业技术服务有限公司 | 工业园区挥发性有机物网格化溯源调查的方法 |
CN110954658B (zh) * | 2019-12-02 | 2023-09-29 | 上海宝钢工业技术服务有限公司 | 工业园区挥发性有机物网格化溯源调查的方法 |
CN112051362A (zh) * | 2020-09-01 | 2020-12-08 | 杭州壹丈科技有限公司 | 风玫瑰预判污染物变化模型 |
CN112051362B (zh) * | 2020-09-01 | 2022-06-24 | 杭州壹丈科技有限公司 | 风玫瑰预判污染物变化模型 |
CN113759441A (zh) * | 2021-09-08 | 2021-12-07 | 长春嘉诚信息技术股份有限公司 | 基于风玫瑰图及污染源监测的空气质量溯源方法 |
CN113759441B (zh) * | 2021-09-08 | 2022-04-22 | 长春嘉诚信息技术股份有限公司 | 基于风玫瑰图及污染源监测的空气质量溯源方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109375285A (zh) | 一种应用于污染分析的玫瑰图制作方法 | |
CN107917987B (zh) | 一种面向城市空气污染物溯源分析方法 | |
Zhong et al. | A city-scale estimation of rooftop solar photovoltaic potential based on deep learning | |
CN105550681B (zh) | 一种基于稳定映射点夜间灯光影像间相互辐射校正方法 | |
CN112989692A (zh) | 基于遥感数据的湖泊富营养化反演方法 | |
CN111879915B (zh) | 一种滨海湿地高分辨率的逐月土壤盐度监测方法及*** | |
CN110119556B (zh) | 一种区域水源涵养功能的时空演变分析方法 | |
Zhang et al. | Application of nonlinear land use regression models for ambient air pollutants and air quality index | |
CN106680833B (zh) | 基于Logistic大气消光廓线识别混合层厚度方法 | |
CN106886564B (zh) | 一种基于空间聚类订正nwp风能图谱的方法及装置 | |
Shao et al. | Impact of urban expansion on meteorological observation data and overestimation to regional air temperature in China | |
Zhou et al. | A multi-perspective study of atmospheric urban heat island effect in China based on national meteorological observations: Facts and uncertainties | |
Cheng et al. | Vegetation Cover Change and Relative Contributions of Associated Driving Factors in the Ecological Conservation and Development Zone of Beijing, China. | |
CN102288646A (zh) | 一种海水中泥沙浓度垂向分布的原位自动监测方法及*** | |
CN104198998B (zh) | 非均匀背景下基于聚类处理的恒虚警检测方法 | |
CN110132246B (zh) | 基于残差分区的水深遥感探测方法 | |
Liu et al. | Local temperature responses to actual land cover changes present significant latitudinal variability and asymmetry | |
CN113469097B (zh) | 一种基于ssd网络的水面漂浮物多相机实时检测方法 | |
CN111881566B (zh) | 基于实景仿真的滑坡位移检测方法及装置 | |
CN107038700B (zh) | 基于gssim的土壤湿度空间分布特征定量分析方法 | |
CN103336972A (zh) | 一种基于完备化局部三值模式的地基云图分类方法 | |
CN112166688B (zh) | 基于小卫星的沙漠与沙漠化土地监测方法 | |
CN117131441A (zh) | 夜间光污染监测方法、装置、计算机设备和存储介质 | |
Schlager et al. | Generation of high-resolution wind fields from the WegenerNet dense meteorological station network in southeastern Austria | |
CN116011881A (zh) | 一种基于像元尺度的黑土农田生产力评价方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Hu Jun Inventor after: Guan Zuguang Inventor after: Gao Fei Inventor after: Weng Yingjie Inventor before: Guan Zuguang Inventor before: Gao Fei Inventor before: Weng Yingjie |
|
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190222 |