一种基于3S技术的区域蚊虫危害形势演示方法

摘要

本发明公开了一种基于3S技术的区域蚊虫危害形势演示方法，主要包括以下步骤：确定管理区域与蚊虫孳生、栖息、宿主动物相关的地物特征，利用灰度转换获取遥感影像图中每种地物特征的灰度区间；依据管理区域的遥感影像图，将环境类型分成若干个类型，在每个类型的环境中选择若干个点位,根据每种地物特征的灰度区间，利用机器判读的方式获得每个点位半径为250米缓冲区内的每种地物特征的大小；应用相似性分析方法，对所有点位的地物特征组成进行相似性分析；建立蚊虫数据矩阵；建立蚊虫优势种数据矩阵等。本发明的方法可以应用于不同地区的蚊虫危害形势演示。

权利要求

1.一种基于3S技术的区域蚊虫危害形势演示方法，其特征在于，包括以下步骤：第一步，确定管理区域与蚊虫孳生、栖息、宿主动物相关的地物特征，利用灰度转换获取遥感影像图中每种地物特征的灰度区间；

第二步，依据管理区域的遥感影像图，将环境类型分成若干个类型，在每个类型的环境中选择若干个点位,根据每种地物特征的灰度区间，利用机器判读的方式获得每个点位半径为250米缓冲区内的每种地物特征的大小；

第三步，应用相似性分析方法，对所有点位的地物特征组成进行相似性分析，在δ％相似性水平下，将所有点位分成相互之间没有交集的若干个组，计算每个组内所有点位每种地物特征的平均值，建立区间数据矩阵；

第四步，采用GB/T23797-2009中的某种或几种蚊虫监测方法调查所有点位的蚊虫群落数据，获得每个点位各蚊种的数据，将每个点位中优势度小于10％的蚊虫数据去除，建立蚊虫数据矩阵；

第五步，对蚊虫群落数据进行相似性分析，在ε％相似性水平下，将所有点位分成相互之间没有交集的若干个组，组数不大于第三步中地物特征的组数，计算获得每组蚊虫优势种排序，建立蚊虫优势种数据矩阵；

第六步，建立第三步与第五步中的组别对应关系；

第七步，根据第二步中的方法分析每个点位缓冲区内的地物构成，将每个点位的地物构成数据与第三步中的若干个组的地物构成数据进行相似性分析，与目标点位相似度最大的那个组所对应的蚊虫优势类群即为目标点位的优势蚊虫；

第八步，获得管理区域的遥感影响图，利用地理信息系统软件建立地物特征的分析点位、缓冲区和蚊虫像元；

第九步，将蚊虫优势类群相同的像元以相同的颜色显示，就能显示出管理区域中的蚊种分布示意图，进而揭示蚊虫潜在危害的分布形势。

说明书

一种基于3S技术的区域蚊虫危害形势演示方法

技术领域

[0001] 本发明涉及一种区域蚊虫危害形势判别方法，具体地说，涉及一种基于3S技术的区域蚊虫危害形势演示方法。

背景技术

[0002] 目前的蚊虫危害形势是通过具体的监测方法，如灯诱法、人诱法等来获得具体点位的蚊虫群落现状。这种方法，要消耗大量的人力、财力、物力才能了解一定区域的蚊虫现状，对于一个城市来讲，这种操作方法基本不可行。因此对于一定范围的区域来讲，无法从宏观层面弄清所辖范围的蚊虫危害形势，即蚊虫不同危害程度的分布状况，从而无法制定更加科学详细的蚊虫及蚊媒病防制预案。

[0003] 蚊虫的生活史包括卵、幼虫、蛹、成虫四个阶段，前三个阶段都生活在水中，适合的水体就是蚊虫的孳生地，孳生地的类型决定了孳生的蚊虫种类，孳生地的数量决定了蚊蚴的环境容量，即可以生存的蚊蚴最大数量。成虫生活在陆地上，靠刺吸植物汁液获得能量和营养生存，雌蚊靠吸血宿主动物血液来完成繁殖，因此陆地的地物特征类型也决定了蚊虫种类，地物特征数量决定了成蚊的环境容量。因此孳生地和陆地地物特征就决定了特定生境中的蚊虫优势种类和数量。3S技术可以通过高分辨率的遥感影像图来确定与蚊虫孳生、栖息、吸血相关的地物特征，因而可以通过3S技术来判别特定生境中的优势蚊种及潜在数量，从而可以判别蚊虫的危害形势。

发明内容

[0004] 本发明的目的在于克服上述技术存在的缺陷，提供一种基于3S技术的区域蚊虫危害形势演示方法。其具体技术方案为：

[0005] 一种基于3S技术的区域蚊虫危害形势演示方法，包括以下步骤：

[0006] 第一步，确定管理区域与蚊虫孳生、栖息、宿主动物相关的地物特征，利用灰度转换获取遥感影像图中每种地物特征的灰度区间；

[0007] 第二步，依据管理区域的遥感影像图，将环境类型分成若干个类型，在每个类型的环境中选择若干个点位,根据每种地物特征的灰度区间(包括色调区间和饱和度区间)，利用机器判读的方式获得每个点位半径为250米缓冲区内的每种地物特征的大小；

[0008] 第三步，应用相似性分析方法，对所有点位的地物特征组成进行相似性分析，在δ％相似性水平下，将所有点位分成相互之间没有交集的若干个组，计算每个组内所有点位每种地物特征的平均值，建立区间数据矩阵；

[0009] 第四步，采用GB/T23797-2009中的某种或几种蚊虫监测方法调查所有点位的蚊虫群落数据，获得每个点位各蚊种的数据，将每个点位中优势度小于10％的蚊虫数据去除，建立蚊虫数据矩阵；

[0010] 第五步，对蚊虫群落数据进行相似性分析，在ε％相似性水平下，将所有点位分成相互之间没有交集的若干个组，组数不大于第三步中地物特征的组数，计算获得每组蚊虫优势种排序，建立蚊虫优势种数据矩阵；

[0011] 第六步，建立第三步与第五步中的组别对应关系；

[0012] 第七步，根据第二步中的方法分析每个点位缓冲区内的地物构成，将每个点位的地物构成数据与第三步中的若干个组的地物构成数据进行相似性分析，与目标点位相似度最大的那个组所对应的蚊虫优势类群即为目标点位的优势蚊虫；

[0013] 第八步，获得管理区域的遥感影响图，利用地理信息系统软件建立地物特征的分析点位、缓冲区和蚊虫像元；

[0014] 第九步，将蚊虫优势类群相同的像元以相同的颜色显示，就能显示出管理区域中的蚊种分布示意图，进而揭示蚊虫潜在危害的分布形势。

[0015] 与现有技术相比，本发明的有益效果为：

[0016] 本发明所述方法可以从宏观层面一目了然的了解管理区域的蚊虫分布情况，从而掌握相应蚊虫传染病的传播风险区域，为划定准确的目标风险区域开展蚊虫和相关蚊媒病防制工作提供科学直观的依据。以前的方法是按照行政区域来划分，不能按照自然客观规律来划分风险区域，目标区域不准确。

具体实施方式

[0017] 为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解，下面结合具体实例进一步阐述本发明。

[0018] 第一步，确定管理区域与蚊虫孳生、栖息、宿主动物相关的地物特征，如建筑物、绿地、耕地、湿地、水体、圈养的牲畜，分别用x1、x2、x3、x4、x5、x6表示，利用灰度转换获取遥感影像图中每种地物特征的灰度区间，饱和度(ci，cj)和色调(hi，hj)，如表1(<0i

[0019] 表1

[0020]    x1 x2 x3 x4 x5
饱和度 (ci，cj) (85,125) (130,170) (45,85) (190,230) (200,240)
色调 (hi，hj) (140，180) (10，50) (20，50) (70，90) (70，110)

[0021] 第二步，依据管理区域的遥感影像图，将环境类型分成若干个类型，在每个类型的环境中选择若干个点位(s1、s2、s3、s4、s5、s6、s7、…)，如表。

[0022]点位序号 点位名称 纬度 经度
s1 简新新村 31°35.212′ 120°19.406′
s2 小娄巷9号 31°34.615′ 120°18.121′
s3 东林书院 31°34.754′ 120°18.248′
s4 映山河医院 31°34.840′ 120°17.948′
s5 五星家园 31°32.388′ 120°18.996′
s6 长欣小区 31°32.953′ 120°21.394′
s7 旺庄医院 31°33.148′ 120°21.324′
s8 七院 31°35.655′ 120°13.668′
s9 黄巷医院 31°36.228′ 120°16.613′
s10 湖父洑西村 31°12.232′ 119°45.466′
s11 大陈巷 31°37.125′ 120°34.616′
s12 费家塘21号 31°37.248′ 120°34.357′
s13 长广溪湿地 31°31.246′ 120°15.675′
s14 长广溪湿地 31°30.152′ 120°15.601′
s15 长广溪湿地 31°30.019′ 120°15.503′
s16 鼋头渚 31°31.702′ 120°12.962′
s17 马山湖边 31°29.206′ 120°07.836′
s18 桃园山庄 31°24.757′ 120°07.594′
s19 彭家塘 31°36.963′ 120°33.295′
s20 前张村 31°37.171′ 120°34.496′
s21 钱更上村 31°34.591′ 120°30.517′

[0023] 根据每种地物特征的灰度区间，利用机器判读的方式获得每个点位半径为250米缓冲区内的每种地物特征的大小。建立地物特征数据矩阵，如下表2。

[0024] 表2

[0025]  x1 x2 x3 x4 x5 x6
s1 146089.509 28052.296 0 0 349.7393 0
s2 176541.196 19633.751 0 0 75.0528 0
s3 171282.972 23046.689 0 0 1920.3397 0
s4 182135.545 11887.175 0 0 2227.2803 0
s5 140907.561 51640.235 0 0 3702.204 0
s6 138333.717 54079.943 0 0 3836.3399 0
s7 120257.584 73541.959 0 0 2450.4577 0
S8 153381.699 42868.301 0 0 0 0
S9 165930.562 23758.124 0 0 480.5601 0
S10 11105 94488 86905 0 0 50000
S11 51337.5429 12723.8 114064.0255 0 18124.632 20000
S12 36774.2276 18137.435 133990.7239 0 7347.6135 80000
S13 41893 72847 0 44566 1779 0
S14 46896.1283 118566.42 0 28039.2639 527.087 0
S15 50759.643 73498.179 0 19677.2518 2083.1021 0
S16 40161.1895 149269.38 0 0 6271.3061 0
S17 62807.2394 38638.454 0 0 8757.1927 0
S18 58204.4028 110725.82 25498.3498 551.0893 1270.3355 0
S19 47259.2342 33239.955 91621.7084 0 24129.103 0
S20 62010.0531 7770.6556 115519.9325 0 10949.359 0
S21 35499.7534 29592.835 128601.2737 0 2556.1378 0

[0026] 第三步，应用相似性分析方法，对所有点位的地物特征组成进行相似性分析，在δ％相似性水平下，将所有点位分成相互之间没有交集的若干个组(Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、…)，计算每个组内所有点位每种地物特征的平均值 建立数据矩阵，如表3所示。

[0027] 表3

[0028]

[0029] 第四步，采用GB/T23797-2009中的某种或几种蚊虫监测方法调查所有点位的蚊虫群落数据，获得每个点位各蚊种的数据，将每个点位中优势度小于5％的蚊虫数据去除，建立蚊虫数据矩阵，m1、m2、m3、m4、m5分别表示淡色库蚊、三带喙库蚊、白纹伊蚊、中华按蚊、骚扰阿蚊，如表4。

[0030] 表4

[0031]  m1 m2 m3 m4 m5
s1 1088 38 24 7 0
s2 128 6 41 23 0
s3 562 13 60 15 0
s4 441 44 8 8 0
s5 251 5 8 5 0
s6 180 7 6 5 0
s7 416 7 2 2 0
S8 317 9 13 4 0
S9 206 6 4 2 0
S10 27 46926 3 7514 4276
S11 6404 117110 12 6270 424
S12 2176 190868 0 9262 50
S13 2642 518 6 38 2
S14 1730 684 8 34 0
S15 1556 2618 4 182 16
S16 158 200 14 25 46
S17 245 421 10 72 4
S18 92 253 0 38 6
S19 652 104 0 206 182
S20 188 910 0 134 232
S21 256 743 0 97 201

[0032] 第五步，对蚊虫群落数据进行相似性分析，在ε％相似性水平下，将所有点位分成相互之间没有交集的若干个组(ⅰ、ⅱ、ⅲ、ⅳ、ⅴ、…)，组数不大于第三步中地物特征的组数。计算获得每组蚊虫优势种排序(M1>M2>M3>M4)，建立蚊虫优势种数据矩阵，如表5。

[0033] 表5

[0034]  ⅰ ⅱ ⅲ ⅳ ⅴ
包含点位 s1-s9 s10-s12 s13-s15 s16-s18 s19-s21
M1 淡色库蚊 三带喙库蚊 淡色库蚊 三带喙库蚊 三带喙库蚊
M2 白纹伊蚊 中华按蚊 三带喙库蚊 淡色库蚊 淡色库蚊
M3 - - - 中华按蚊 骚扰阿蚊
M4 - - - - 中华按蚊

[0035] 第六步，建立第三步与第五步中的组别对应关系，ⅰ、ⅱ、ⅲ、ⅳ、ⅴ分别与Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ相对应。

[0036] 第七步，获得管理区域的遥感影响图，利用地理信息系统软件建立地物特征的分析点位、缓冲区和蚊虫像元。在区域遥感影响图中，以 为边长建立最大数目的正方形网格图，每个正方形网格即为一个蚊虫像元。以每个像元的几何中心为圆心，以为半径建立缓冲区，每个像元的几何中心即为分析点位。

[0037] 第八步，根据第二步中的方法分析每个点位缓冲区内的地物构成，将每个点位的地物构成数据与第三步中的若干个组的地物构成数据进行相似性分析，与目标点位相似度最大的那个组所对应的蚊虫优势类群即为目标点位的优势蚊虫。

[0038] 第九步，将蚊虫优势类群相同的像元以相同的颜色显示，就能显示出管理区域中的蚊种分布示意图，进而揭示蚊虫潜在危害的分布形势。

[0039] 以上所述，仅为本发明最佳实施方式，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可显而易见地得到的技术方案的简单变化或等效替换均落入本发明的保护范围内。