CN102254088B - Biohazard assessment visualization method based on Google map - Google Patents

Abstract

The invention relates to a biohazard assessment visualization method based on a Google map. The method comprises the following steps: 1) taking a Bernoulli curve function as a biohazard particle simulated aerosol diffusion model; 2) expressing the aerosol diffusion model as a polar coordinate form; 3) drawing an aerosol diffusion curve by using the GPolygon class of Google Maps API, and obtainingan actual hazard area corresponding to the aerosol diffusion curve by using getArea (); 4) dividing the actual hazard area into a plurality of areas, arbitrarily acquiring the longitude/latitude coordinates of a point in each area, then judging which county the area belongs to by a getLocations () function of Google Maps; 5) invoking the population density data of the county from a database of Google Maps, and multiplying the population density data by the area of the county so as to obtain the number of harmed people in the county; and 6) adding the numbers of harmed people in the divided areas so as to obtain a total number of harmed people. The method provided by the invention is high in accuracy and efficiency, and can be widely used in sudden biohazard assessment processes.

Description

A kind of biohazard assessment method for visualizing based on Google Maps

Technical field

The present invention relates to a kind of biohazard assessment method for visualizing, particularly about a kind of biohazard assessment method for visualizing based on Google Maps (Google Maps).

Background technology

" anthrax mail " incident after the U.S. " 911 " incident, indicate that biohazard becomes a reality once more, and, face new International Politics, economy and military environment, biohazard not only is confined to biological warfare, also comprises the harm that bio-terrorism, biological accident, New Development infectious disease, deadly infectious disease etc. bring for military operation, national security, social stability, economic development, national existence and health.

The biohazard incident early has precedent, and Japan " 731 " army attacks by the biological biology of implementing of media in World War II, and at that time the Chinese armed forces and the people have been caused huge damage; After this, the biohazard incident emerges in an endless stream, and the U.S. after the September 11 attacks " anthrax mail " incident indicates that biohazard develops into a new stage, and a series of ricin (WA) incident has taken place again thereafter, has caused great fear to the whole world.And the develop rapidly of modern molecular biology and related discipline, make enforcement attack of terrorism technology perfect, increased the possibility that extensive biohazard incident takes place greatly.Two prewar, biohazard is many from biological warfare, thereafter international community makes great efforts to advance the plan of control biological weapons always, just repeatedly held the international conference of anti-bio-terrorism as far back as scientists in 1998, also set up common people's biophylaxis tissue of forming by scientist, yet the research and development of current biohazard technology still are in to be prohibited and more than state, and is presenting to higher level, the gesture of high-stage development and diffusion more.High incidence and high mortality that modern biohazard has; Obtain easily; Producing cost is low; Can generate a panic; Diagnosis, treatment difficulty; Part can cause the human world to infect; Advantages such as the thrower is easy to run away obtain terrorist's favor.

After the biohazard incident takes place, take fast and effectively that emergency measure and countermeasure are very necessary, appropriate measures can remedy the loss that the biohazard incident is caused.Many countries have begun the Protective Research at burst biohazard incident, arms reaction as chemistry, burst biohazard incident is supported, biohazard incident medical science treatment training, rapid response force's tissue and communal facility biohazard are attacked the detection Study on Technology, and the analogue simulation of biohazard incident etc., biological warfare has also caused the great vigilance of various countries.

Summary of the invention

At the problems referred to above, the purpose of this invention is to provide a kind of biohazard incident assessment method for visualizing based on Google Maps (Google Maps), this method can obtain rapidly the to happen suddenly damaging range of biohazard incident and the number of being injured are to take the necessary condition of emergency measure and countermeasure fast and effectively.

For achieving the above object, the present invention takes following technical scheme: a kind of biohazard assessment method for visualizing based on Google Maps, and it may further comprise the steps: 1) adopt the model of Bernoulli Jacob's curvilinear function as the aerosol diffusion of simulation burst biohazard incident hazardous material:

[x ²+y ²] ²＝L ²[x ²-y ²](x＞0)(i＝0，1，·，n) (1)

In the formula, the harm depth when L is aerosol diffusion:

L＝v×t×u×ω×q(i＝0，1，·，n) (2)

In the formula, v represents aerosol diffusion speed; T represents to spread the duration; U represents the topographic correction coefficient; ω represents the precipitation correction coefficient; Q represents the wind speed correction coefficient; 2) aerosol diffusion model (1) is expressed as polar form:

r ²＝L ²cos2θ

Wherein, L is the harm depth, and getting direct north is 0, the span of angle θ-90 °～90 °; 3) adopt the GPolygon class of Google Maps API to draw the aerosol diffusion curve: angle value θ from-90 ° of every variations of beginning once calculates the position of a point, equals 90 ° up to θ, obtains 181 points altogether, and wherein the coordinate of each point can be expressed as (r _i, θ _i), (i=1,2 ... 181), owing to the point on the Google Maps indicates with longitude and latitude, therefore, with polar coordinates (r _i, θ _i) be converted to latitude and longitude coordinates:

y＝arcsin(sin(oy)×cos(d)+cos(oy)×sin(d)×cos(h))

$x=\mathrm{ox}+\arctan \left(\frac{\sin \left(h\right)\×\sin \left(d\right)\×\cos \left(\mathrm{oy}\right)}{\cos \left(d\right)-\sin \left(y\right)}\right)$

In the formula, ox represents the radian value of aerosol diffusion starting point longitude, and oy represents the radian value of aerosol diffusion starting point latitude, and d represents to endanger the radian value of depth L, and h represents θ _iRadian value, (x, the y) (r that tries to achieve of expression _i, θ _i) corresponding latitude and longitude coordinates; After obtaining the latitude and longitude coordinates of 181 points, draw the aerosol diffusion curve by the GPolygon class of Google Maps API; 4) the function getArea () according to the closed curve reference area that utilizes the GPolygon class of Google Maps API to provide draws the area of the actual harm region of the aerosol diffusion curve correspondence of drawing in the step 3); 5) area with the actual harm region of the aerosol diffusion curve correspondence that obtains in the step 4) is divided into some, in each piece, obtain the latitude and longitude coordinates of a point arbitrarily, judge this piece zone by the getLocations () function of GoogleMaps and belong to which district; 6) call the density of population data in this district in the database of Google Maps, the area in itself and this district is multiplied each other, what obtain this district is subjected to the harm number; The number of the being injured addition in each the piece district that 7) will divide obtains total harm number that is subjected to, and takes emergency measure and countermeasure fast and effectively.

In the described step 1), the obtaining value method of each coefficient of formula (2) is as follows: 1. the aerosol diffusion speed v is the wind speed in the corresponding moment; 2. aerosol diffusion duration t, as shown in table 1:

10 half life period of table 1 biohazard material

Grade	Burst biohazard material title	10 half life period (hour)
			1	Bacillus anthracis	10
2	Yersinia pestis, bacterium tularense	5.0
			3	Brucella	2.8
4	Comma bacillus	2.5
			5	The hot Richettsia of Q	1.0
6	Chlamydia psittaci	0.5
			7	Virus, mycin class	0.3

3. the value of topographic correction coefficient u is as shown in table 2;

Table 2 topographic correction coefficient

Landform	Correction factor
		Plain, city	1
Thick forest, hills, mountain region	0.6
		Jungle, rainforest	0.2

4. precipitation correction coefficient ω: daily rainfall does not need to revise during less than 10 millimeters; Daily rainfall is just needing during greater than 10 millimeters, correction factor ω=0.2; 5. wind speed correction coefficient q: value is as shown in table 3;

Table 3 wind speed correction coefficient

In the described step 3):

D=L/6371;

, (lon, lat) the latitude, longitude coordinate of expression aerosol diffusion starting point.

The present invention is owing to take above technical scheme, it has the following advantages: 1, the present invention is compiling on the aerosol diffusion hazard evaluation model based, utilize the Bai Nuli curve to realize visual hazard evaluation in conjunction with Google Maps, can on map, assess hazard area rapidly, intuitively, have that accuracy height, computing velocity are fast, model evaluation characteristics intuitively in the animation mode.2, the present invention is aspect influenced number calculating, proposition is based on influenced people's number calculating method in geographic position, to in the past can only be according to the influenced number of national density of population data computation, bring up to and to calculate according to district density of population data accurately, further improved the degree of accuracy of assessment.The inventive method accuracy height, the efficient height can be widely used in the burst biohazard incident evaluation process.

Description of drawings

Fig. 1 is a schematic flow sheet of the present invention

Fig. 2 is the synoptic diagram of Bai Nuli simulation aerosol diffusion of the present invention

Fig. 3 is the function synoptic diagram of Bai Nuli curve of the present invention under polar coordinates

Fig. 4 is that the present invention utilizes Google Maps API to calculate the synoptic diagram of hazard area

Fig. 5 is that the present invention utilizes the piecemeal method to find the solution the synoptic diagram of influenced number

Embodiment

Below in conjunction with drawings and Examples the present invention is described in detail.

Google Maps API (Application Programming Interface, application programming interface) is a free services that provides by Google, allow the user to use Google Maps in the system of oneself, also provide a whole set of program development interface to be used for developing some application of oneself simultaneously, the present invention utilizes the drawing function of Google Maps to draw aerocolloidal range of scatter, utilizes the geographical information query function of Google Maps to calculate affected number in conjunction with demographic data.

As shown in Figure 1, the inventive method may further comprise the steps:

1, determines the mathematical model that hazard evaluation is required

The present invention utilizes the Bai Nuli curve to simulate the aerosol diffusion form of burst biohazard incident, makes people can describe and study the Polluted area of biohazard with the mode of mathematics.The Bai Nuli curve is at the figure in the difference moment (t1-t5), as shown in Figure 2.

1) adopt the aerosol diffusion model of Bernoulli Jacob's curvilinear function as simulation biohazard incident:

[x ²+y ²] ²＝L ²[x ²-y ²](x＞0)(i＝0，1，·，n) (1)

In the formula, the harm depth when L represents aerosol diffusion, computing formula is as follows:

L＝v×t×u×ω×q(i＝0，1，·，n) (2)

In the formula, v represents aerosol diffusion speed; T represents to spread the duration; U represents the topographic correction coefficient; ω represents the precipitation correction coefficient; Q represents the wind speed correction coefficient.

2) value of each parameter is analyzed as follows in the formula (2):

1. aerosol diffusion speed v: because gasoloid mainly relies on wind-force to propagate, therefore, the aerosol diffusion speed v is defaulted as the wind speed in this moment.

2. spread duration t:

After it has been generally acknowledged that 10 half life period, it is original 0.1% that the concentration value of biohazard particle has only, and can think to have lost appeal.Therefore, the duration with diffusion when calculating is decided to be 10 times half life period, and is the diffusion time of common biohazard particle, as shown in table 1.

10 half life period of the common biohazard particle of table 1

Grade	Biohazard particle title	10 half life period (hour)
			1	Bacillus anthracis	10
2	Yersinia pestis, bacterium tularense	5.0
			3	Brucella	2.8
4	Comma bacillus	2.5
			5	The hot Richettsia of Q	1.0
6	Chlamydia psittaci	0.5
			7	Virus, mycin class	0.3

3. the value of topographic correction coefficient u is as shown in table 2.

Table 2 topographic correction coefficient

Landform	Correction factor
		Plain, city	1
Thick forest, hills, mountain region	0.6
		Jungle, rainforest	0.2

4. precipitation correction coefficient ω: precipitation event has influence on the decline of biohazard particle.Consider the influence of rain factor, can revise with precipitation correction factor ω according to precipitation event.Daily rainfall does not need during less than 10 millimeters to revise; Daily rainfall is just needing during greater than 10 millimeters, correction factor ω=0.2.

(5) wind speed correction coefficient q:

Wind speed will influence the aerocolloidal range of scatter of biohazard particle, therefore need to increase the wind speed correction coefficient, overhead height H, atmospheric conditions (being divided into inversion and isothermal) are relevant during the aerosol diffusion of this coefficient and biohazard particle, and its value is as shown in table 3.

Table 3 wind speed correction coefficient

2, the calculating of influence area

Google Maps API provide GPolygon function, and it can be formed according to number of vertex and build polygon.The present invention adopts the method curve plotting of segmented fitting.Graphing adopts the polar coordinates marker method of Berkeley curve to draw for convenience, the figure under the polar coordinates, as shown in Figure 3.

Equation is as follows:

r ²＝2a ²cos2θ

Getting direct north is 0, angle span-90 °～90 °, and the Berkeley curve can be expressed as so:

r ²＝L ²cos2θ

Wherein, L is the harm depth, and the θ span is-90 °～90 °.For matched curve, angle value from-90 ° of every variations of beginning once calculates the position of a point, equals 90 ° up to θ, obtains 181 points so altogether.Wherein the coordinate of each point can be expressed as (r _i, θ _i), (i=1,2 ... 181) because the point on the Google Maps indicates with longitude and latitude, therefore need carry out the conversion of polar coordinates to latitude and longitude coordinates.The latitude and longitude coordinates of aerosol diffusion starting point (lon, lat) known, conversion formula is as follows:

$\mathrm{ox}=\frac{\mathrm{lon}\×\mathrm{\π}}{180}$

$\mathrm{oy}=\frac{\mathrm{lat}\×\mathrm{\π}}{180}$

$d=\frac{L}{6371}$

$h=\frac{{\mathrm{\θ}}_i\×\mathrm{\π}}{180}$

y＝arcsin(sin(oy)×cos(d)+cos(oy)×sin(d)×cos(h))

$x=\mathrm{ox}+\arctan \left(\frac{\sin \left(h\right)\×\sin \left(d\right)\×\cos \left(\mathrm{oy}\right)}{\cos \left(d\right)-\sin \left(y\right)}\right)$

In the formula, ox represents the radian value of aerosol diffusion starting point longitude, and oy represents the radian value of aerosol diffusion starting point latitude, and d represents to endanger the radian value of depth L, and h represents θ _iRadian value, (the r that tries to achieve of (x y) expression _i, θ _i) corresponding latitude and longitude coordinates.

After obtaining the latitude and longitude coordinates of 181 points, by the GPolygon class of the Google Maps API aerosol diffusion curve (as shown in Figure 4) that just can draw.

GPolygon also provides the method according to the closed curve reference area: getArea () roughly the same the time, by calling the actual affected area that this method can easily draw aerosol diffusion curve correspondence.This method specifies as shown in table 4.

The explanation of table 4getArea () method

3, influenced number determines

In order to obtain affected effective strength, general method is directly to multiply each other with the hazard area and the national density of population, but because the density of population of different provinces and cities is greatly different, the degree of accuracy as a result that this mode is tried to achieve is not high, the present invention adopts the way of block sampling, multiply each other with actual district density of population data and corresponding region area, thereby obtain comparatively accurate data, specific embodiment is as follows:

1) hazard area is divided into 8, as shown in Figure 5, in each piece, obtains the latitude and longitude coordinates of a point arbitrarily,, can judge this piece zone and belong to which district by the getLocations () method (the method explanation is as shown in table 5) of Google Maps.

2) call the density of population data in this district in the database of Google Maps, itself and this district area are multiplied each other, what obtain this district is subjected to the harm number.

That 3) calculates whole 8 district areas is subjected to the harm number, draws total harm number that is subjected to, and takes emergency measure and countermeasure fast and effectively then.

The explanation of table 5getLocations () method

In the foregoing description, the various functions that Google Maps provides, and the getLocations () that provides of Google Maps API, the using method of getLocations functions such as (), belong to prior art, specifically use, method of operating no longer describes in detail at this.

The various embodiments described above only are used to illustrate the present invention, and every equivalents of carrying out on the basis of technical solution of the present invention and improvement all should not got rid of outside protection scope of the present invention.

Claims (2)

1. method for visualizing is assessed in the biohazard based on Google Maps, and it may further comprise the steps:

1) adopt the model of Bernoulli Jacob's curvilinear function as the aerosol diffusion of simulation burst biohazard incident hazardous material:

[x ²+ y ²] ²=L ²[x ²One y ²] (x＞0) (i=0,1 ..., n) (1)

In the formula, the harm depth when L is aerosol diffusion:

L＝v×t×u×ω×q(i＝0,1,…,n) （2）

In the formula, v represents aerosol diffusion speed; T represents to spread the duration; U represents the topographic correction coefficient; ω represents the precipitation correction coefficient; Q represents the wind speed correction coefficient;

2) aerosol diffusion model (1) is expressed as polar form:

r ²=L ²cos2θ

Wherein, L is the harm depth, and getting direct north is 0, the span of angle θ-90 °～90 °;

3) adopt the GPolygon class of Google Maps API to draw the aerosol diffusion curve:

Angle value θ from-90 ° of every variations of beginning once calculates the position of a point, equals 90 ° up to θ, obtains 181 points altogether, and wherein the polar coordinates of each point can be expressed as (r _i, θ _i), (i=1,2 ... 181), owing to the point on the Google Maps indicates with longitude and latitude, therefore, with polar coordinates (r _i, θ _i) be converted to latitude and longitude coordinates:

y＝arcsin(sin(oy)×cos(d)+cos(oy)×sin(d)×cos(h))

$x=\mathrm{ox}+\arctan \left(\frac{\sin \left(h\right)\×\sin \left(d\right)\×\cos \left(\mathrm{oy}\right)}{\cos \left(d\right)-\sin \left(y\right)}\right)$

In the formula, ox represents the radian value of aerosol diffusion starting point longitude,

Oy represents the radian value of aerosol diffusion starting point latitude,

D represents to endanger the radian value of depth L, d=L/6371; H represents θ _iRadian value,

(x, the y) (r that tries to achieve of expression _i, θ _i) corresponding latitude and longitude coordinates; Wherein, (lon, lat) the latitude, longitude coordinate of expression aerosol diffusion starting point;

After obtaining the latitude and longitude coordinates of 181 points, draw the aerosol diffusion curve by the GPolygon class of Google Maps API;

4) the function getArea () according to the closed curve reference area that utilizes the GPolygon class of Google Maps API to provide draws the area of the actual harm region of the aerosol diffusion curve correspondence of drawing in the step 3);

5) area with the actual harm region of the aerosol diffusion curve correspondence that obtains in the step 4) is divided into some, in each piece, obtain the latitude and longitude coordinates of a point arbitrarily, judge this piece zone by the getLocations () function of Google Maps and belong to which district;

6) call the density of population data in this district in the database of Google Maps, the area in itself and this district is multiplied each other, what obtain this district is subjected to the harm number;

The number that the endangered addition in each the piece district that 7) will divide obtains total harm number that is subjected to, and takes emergency measure and countermeasure fast and effectively.

2. a kind of biohazard assessment method for visualizing as claimed in claim 1 based on Google Maps, it is characterized in that: in the described step 1), the obtaining value method of each coefficient of formula (2) is as follows:

1. the aerosol diffusion speed v is the wind speed in the corresponding moment;

2. aerosol diffusion duration t, as shown in table 1:

10 half life period of table 1 biohazard material

Grade Burst biohazard material title 10 half life period (hour) 1 Bacillus anthracis 10 2 Yersinia pestis, bacterium tularense 5.0 3 Brucella 2.8 4 Comma bacillus 2.5 5 The hot Richettsia of Q 1.0 6 Chlamydia psittaci 0.5 7 Virus, mycin class 0.3

3. the value of topographic correction coefficient u is as shown in table 2;

Table 2 topographic correction coefficient

Landform Correction coefficient Plain, city 1 Thick forest, hills, mountain region 0.6 Jungle, rainforest 0.2

4. precipitation correction coefficient ω: daily rainfall does not need to proofread and correct during less than 10 millimeters; Daily rainfall needs during greater than 10 millimeters to proofread and correct correction coefficient ω=0.2;

5. wind speed correction coefficient q: value is as shown in table 3;

Table 3 wind speed correction coefficient

Images