CN109087011A - Assess method and device, storage medium and the processor of state of ecological environment - Google Patents
Assess method and device, storage medium and the processor of state of ecological environment Download PDFInfo
- Publication number
- CN109087011A CN109087011A CN201810905136.5A CN201810905136A CN109087011A CN 109087011 A CN109087011 A CN 109087011A CN 201810905136 A CN201810905136 A CN 201810905136A CN 109087011 A CN109087011 A CN 109087011A
- Authority
- CN
- China
- Prior art keywords
- target area
- mentioned
- area
- index
- ecological environment
- 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
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000003860 storage Methods 0.000 title claims abstract description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 100
- 230000006378 damage Effects 0.000 claims abstract description 49
- 239000003245 coal Substances 0.000 claims description 43
- 230000008859 change Effects 0.000 claims description 42
- 238000010606 normalization Methods 0.000 claims description 41
- 230000015556 catabolic process Effects 0.000 claims description 16
- 238000006731 degradation reaction Methods 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 13
- 238000007781 pre-processing Methods 0.000 claims description 7
- 239000002689 soil Substances 0.000 claims description 6
- 230000007613 environmental effect Effects 0.000 description 9
- 230000006870 function Effects 0.000 description 9
- 239000000047 product Substances 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06395—Quality analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of method and device, storage medium and processors for assessing state of ecological environment.Wherein, this method comprises: obtaining the weight parameter of the remotely-sensed data of target area and the evaluation index and above-mentioned evaluation index of initial environment data and above-mentioned target area;Above-mentioned remotely-sensed data and above-mentioned initial environment data are pre-processed, preprocessed data is obtained, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface area;According to above-mentioned evaluation index and above-mentioned weight parameter, above-mentioned preprocessed data is assessed, obtains the assessment result of the state of ecological environment for evaluating above-mentioned target area.The technical issues of it is narrow that the present invention solves the data surface that Ecology Environment Evaluation method in the prior art obtains, and evaluation index is single, leads to not accurate evaluation state of ecological environment.
Description
Technical field
The present invention relates to mining area ecological environment evaluation areas, in particular to a kind of side for assessing state of ecological environment
Method and device, storage medium and processor.
Background technique
Currently, still providing comprehensive, system evaluation method without being directed to state of ecological environment in the related technology, only focusing on
Big zone index numerical value, cannot examine state of ecological environment from part, have limitation, one-sidedness, to ecological environment shape
Condition evaluates the method system still without system, specification.
Specifically, in terms of above-mentioned limitation and one-sidedness are mainly reflected in following two: 1, the setting of existing evaluation index
It is assessed mainly for big Macro, but to the micro cells such as coal mine domain Ecology Environment Evaluation, the index
Lack operability support, some critical evaluation achievement datas can not obtain, as index, environmental restrictions index, the network of rivers are coerced in soil
Dnesity index etc.;2, existing evaluation index method lacks the assessment to ecological environment damage process and mode, the prior art
In do not consider different ecological environment examination elements, for the difference of coal mine pit mining and strip mining transformation mode, provide pair
Different index calculating methods is answered, the assessment report for causing assessment to obtain is inaccurate.
The data surface obtained for above-mentioned coal mine Ecology Environment Evaluation method in the prior art is narrow, and evaluation index list
One, the problem of leading to not accurate evaluation state of ecological environment, currently no effective solution has been proposed.
Summary of the invention
The embodiment of the invention provides it is a kind of assess state of ecological environment method and device, storage medium and processor,
It is narrow at least to solve the data surface that Ecology Environment Evaluation method in the prior art obtains, and evaluation index is single, leads to not
The technical issues of accurate evaluation state of ecological environment.
According to an aspect of an embodiment of the present invention, a kind of method for assessing state of ecological environment is provided, comprising: obtain
The power of the remotely-sensed data and initial environment data of target area and the evaluation index of above-mentioned target area and above-mentioned evaluation index
Weight parameter;Above-mentioned remotely-sensed data and above-mentioned initial environment data are pre-processed, preprocessed data is obtained, wherein is above-mentioned pre-
Processing data include at least: vegetative coverage area and destruction of surface area;It is right according to above-mentioned evaluation index and above-mentioned weight parameter
Above-mentioned preprocessed data is assessed, and the assessment result of the state of ecological environment for evaluating above-mentioned target area is obtained.
Further, above-mentioned assessment result includes: state of ecological environment, is joined according to above-mentioned evaluation index and above-mentioned weight
Number, assesses above-mentioned preprocessed data, obtains the assessment result of the state of ecological environment for evaluating above-mentioned target area,
Include: to be calculated by the way of weighted sum above-mentioned evaluation index and above-mentioned weight parameter, obtains above-mentioned target area
State of ecological environment index;It is assessed, is obtained above-mentioned according to the above-mentioned above-mentioned preprocessed data of state of ecological environment exponent pair
The state of ecological environment of target area, wherein above-mentioned state of ecological environment includes at least: the state of ecological environment of first level,
State of ecological environment, the state of ecological environment of third level, the state of ecological environment of fourth level, fifth level of second level
State of ecological environment;The amplitude of variation of above-mentioned state of ecological environment includes at least: the first amplitude of variation, the second amplitude of variation,
Third amplitude of variation, the 4th amplitude of variation;The rate of change of above-mentioned state of ecological environment includes at least: the first rate of change,
Two rate of changes, third rate of change, the 4th rate of change.
Further, above-mentioned evaluation index includes at least one following: eco-environmental quality index HQI, vegetation coverage
Index VCI, land degradation index LDI and destruction of surface index LSDI.
Further, above-mentioned eco-environmental quality index HQI:HQI=A is calculated by following formulabio×(Slw×Sl+Scw
×Sc+Ssw×Ss+Sgw×Sg+Sjw×Sj+Sww×Sw)/SCM;Wherein, AbioFor the eco-environmental quality index of above-mentioned target area
Normalization coefficient;SlFor the forest land area of above-mentioned target area;SlwFor the weighted value in above-mentioned forest land;ScFor above-mentioned target area
Grassland area;ScwFor the weighted value on above-mentioned meadow;SsFor the Wetland Area area of above-mentioned target area;SswIt is wet for above-mentioned waters
The weighted value on ground;SgFor the cultivated area of above-mentioned target area;SgwFor the weighted value in above-mentioned arable land;SjFor above-mentioned target area
Construction land area;SjwFor the weighted value of above-mentioned construction land;SwFor the unused land area of above-mentioned target area;SwwIt is above-mentioned
The weighted value of unused land;SCMFor the gross area of above-mentioned target area.
Further, above-mentioned vegetation coverage index VCI:VCI=A is calculated by following formulavfc×(Svfc1w×Svfc1+
Svfc2w×Svfc2+Svfc3w×Svfc3+Svfc4w×Svfc4+Svfc5w×Svfc5)/SCM;Wherein, AvfcFor the vegetation of above-mentioned target area
The normalization coefficient of coverage index;Svfc1For the high vegetation coverage area of above-mentioned target area;Svfc1wFor above-mentioned high vegetation
The weighted value of coverage;Svfc2For vegetation coverage area high in above-mentioned target area;Svfc2wFor high vegetative coverage among the above
The weighted value of degree;Svfc3For the middle vegetation coverage area of above-mentioned target area;Svfc3wFor the weight of vegetation coverage among the above
Value;Svfc4For the low vegetation coverage area of above-mentioned target area;Svfc4wFor the weighted value of above-mentioned low vegetation coverage;Svfc5For
The extremely low vegetation coverage area of above-mentioned target area;Svfc5wFor the weighted value of above-mentioned extremely low vegetation coverage;SCMFor above-mentioned mesh
Mark the gross area in region.
Further, above-mentioned land degradation index LDI:LDI=A is calculated by following formulaero×(Sero1w×Sero1+
Sero2w×Sero2+Sero3w×Sero3)/SCM;Wherein, AeroFor the normalization coefficient of the land degradation index of above-mentioned target area;
Sero1For the area of the slight eroded field of above-mentioned target area;Sero1wFor the weighted value of above-mentioned slight eroded field;Sero2For
The area of the moderate eroded field of above-mentioned target area;Sero2wFor the weighted value of above-mentioned moderate eroded field;Sero3For above-mentioned mesh
Mark the area of the severe eroded field in region;Sero3wFor the weighted value of above-mentioned severe eroded field;SCMFor above-mentioned target area
The gross area.
Further, above-mentioned destruction of surface index LSDI is included at least: the destruction of surface index LSDI1 of underground mine, is led to
Cross destruction of surface index LSDI1:LSDI1=(A × A that following formula calculates above-mentioned underground minesus×Ssus)/Scm+(B×Acmw
×Scmw)/Scm;Wherein, AsusFor the normalization coefficient of the surface collapse of above-mentioned target area;AcmwFor the row of above-mentioned target area
The normalization coefficient of cash field;SsusFor the area of the surface collapse of above-mentioned target area;ScmwFor the discharge refuse field of above-mentioned target area
Area;A is the weighted value of above-mentioned surface collapse;B is the weighted value of above-mentioned discharge refuse field;SCMFor total face of above-mentioned target area
Product.
Further, above-mentioned destruction of surface index LSDI is included at least: the table destructing index LSDI2 of open coal mine passes through
Following formula calculates the destruction of surface index LSDI2:LSDI2=(S of above-mentioned open coal mineldsw×Alds×Slds)/Scm+(Scmdw
×Acmd×Scmd)/Scm;Wherein, AldsFor the normalization coefficient of the excavating plant of above-mentioned target area;AcmdFor above-mentioned target area
Refuse dump normalization coefficient;SldsFor the area of the excavating plant of above-mentioned target area;SldswFor the power of above-mentioned excavating plant
Weight values;ScmdFor the area of the refuse dump of above-mentioned target area;ScmdwFor the weighted value of above-mentioned refuse dump;SCMFor above-mentioned target area
The gross area in domain.
According to another aspect of an embodiment of the present invention, a kind of device for assessing state of ecological environment is additionally provided, comprising: obtain
Modulus block, for obtain target area remotely-sensed data and initial environment data and the evaluation index of above-mentioned target area and
The weight parameter of above-mentioned evaluation index;Preprocessing module, it is pre- for being carried out to above-mentioned remotely-sensed data and above-mentioned initial environment data
Processing, obtains preprocessed data, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface area;It comments
Module is estimated, for assessing above-mentioned preprocessed data, obtaining for commenting according to above-mentioned evaluation index and above-mentioned weight parameter
The assessment result of the state of ecological environment of the above-mentioned target area of valence.
According to another aspect of an embodiment of the present invention, a kind of storage medium is additionally provided, above-mentioned storage medium includes storage
Program, wherein above procedure operation when control above-mentioned storage medium where equipment execute following steps: obtain target area
Remotely-sensed data and initial environment data and above-mentioned target area evaluation index and above-mentioned evaluation index weight parameter;
Above-mentioned remotely-sensed data and above-mentioned initial environment data are pre-processed, preprocessed data is obtained, wherein above-mentioned preprocessed data
It includes at least: vegetative coverage area and destruction of surface area;According to above-mentioned evaluation index and above-mentioned weight parameter, to above-mentioned pre- place
Reason data are assessed, and the assessment result of the state of ecological environment for evaluating above-mentioned target area is obtained.
According to another aspect of an embodiment of the present invention, a kind of processor is additionally provided, above-mentioned processor is used to run program,
Wherein, following steps are executed when above procedure is run: obtaining the remotely-sensed data and initial environment data and above-mentioned of target area
The weight parameter of the evaluation index of target area and above-mentioned evaluation index;To above-mentioned remotely-sensed data and above-mentioned initial environment data into
Row pretreatment, obtains preprocessed data, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface face
Product;According to above-mentioned evaluation index and above-mentioned weight parameter, above-mentioned preprocessed data is assessed, is obtained for evaluating above-mentioned mesh
Mark the assessment result of the state of ecological environment in region.
In embodiments of the present invention, pass through the remotely-sensed data and initial environment data of acquisition target area and above-mentioned mesh
Mark the evaluation index in region and the weight parameter of above-mentioned evaluation index;Above-mentioned remotely-sensed data and above-mentioned initial environment data are carried out
Pretreatment, obtains preprocessed data, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface area;
According to above-mentioned evaluation index and above-mentioned weight parameter, above-mentioned preprocessed data is assessed, is obtained for evaluating above-mentioned target
The assessment result of the state of ecological environment in region.
The state of ecological environment of accurate and comprehensive objective assessment target area is achieved the purpose that as a result, to realize
The technical effect of the accuracy of the assessment result of ecological environment is improved, and then solves Ecology Environment Evaluation in the prior art
The technical issues of data surface that method obtains is narrow, and evaluation index is single, leads to not accurate evaluation state of ecological environment.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, this hair
Bright illustrative embodiments and their description are used to explain the present invention, and are not constituted improper limitations of the present invention.In the accompanying drawings:
Fig. 1 is a kind of flow chart of method for assessing state of ecological environment according to an embodiment of the present invention;
Fig. 2 is a kind of flow chart of the method for optional assessment state of ecological environment according to an embodiment of the present invention;And
Fig. 3 is a kind of structural schematic diagram of device for assessing state of ecological environment according to an embodiment of the present invention.
Specific embodiment
In order to enable those skilled in the art to better understand the solution of the present invention, below in conjunction in the embodiment of the present invention
Attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only
The embodiment of a part of the invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people
The model that the present invention protects all should belong in member's every other embodiment obtained without making creative work
It encloses.
It should be noted that description and claims of this specification and term " first " in above-mentioned attached drawing, "
Two " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way
Data be interchangeable under appropriate circumstances, so as to the embodiment of the present invention described herein can in addition to illustrating herein or
Sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that cover
Cover it is non-exclusive include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to
Step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, product
Or other step or units that equipment is intrinsic.
Embodiment 1
According to embodiments of the present invention, a kind of embodiment of method for assessing state of ecological environment is provided, needs to illustrate
It is that step shown in the flowchart of the accompanying drawings can execute in a computer system such as a set of computer executable instructions,
Also, although logical order is shown in flow charts, and it in some cases, can be to be different from sequence execution herein
Shown or described step.
Fig. 1 is a kind of flow chart of method for assessing state of ecological environment according to an embodiment of the present invention, as shown in Figure 1,
This method comprises the following steps:
Step S102 obtains the remotely-sensed data of target area and the evaluation of initial environment data and above-mentioned target area
The weight parameter of index and above-mentioned evaluation index.
In above-mentioned steps S102, above-mentioned target area, which is seen, be can be, but not limited to as coal mine region;Above-mentioned remotely-sensed data can
Think that the remotely-sensed data in region locating for the coal mine obtained by remote sensing equipment, the remotely-sensed data are image data;Above-mentioned evaluation refers to
Mark is used for the state of ecological environment in evaluation goal region, can be, but not limited to are as follows: state of ecological environment evaluation index, wherein raw
State environmental aspect index (CMRI) is the evaluation index for reflecting coal field ecological environment integrated status.
As a kind of optional embodiment, above-mentioned evaluation index includes at least one following: eco-environmental quality index
HQI, vegetation coverage index VCI, land degradation index LDI and destruction of surface index LSDI.
Step S104 pre-processes above-mentioned remotely-sensed data and above-mentioned initial environment data, obtains preprocessed data,
In, above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface area.
In the embodiment of the present application, can be, but not limited to by remote sensing professional treatment software platform to above-mentioned remote sensing number
According to being handled, including but not limited to following processing mode: image radiation calibration, atmospheric correction, geometric accurate correction, visual fusion
Equal enhancings handle work;And existing difference between the remotely-sensed data and above-mentioned initial environment data obtained in real time by comparing,
Obtain above-mentioned preprocessed data.
Step S106 assesses above-mentioned preprocessed data, obtains according to above-mentioned evaluation index and above-mentioned weight parameter
For evaluating the assessment result of the state of ecological environment of above-mentioned target area.
In an alternative embodiment, by using the evaluation index and weight parameter got, to obtained pre- place
Reason data assessed, due in above-mentioned preprocessed data vegetative coverage area and destruction of surface area can be used for characterizing life
State environmental aspect, thus above-mentioned preprocessed data is assessed, namely state of ecological environment is assessed, and then can obtain
To assessment result, which can be used for evaluating the state of ecological environment in coal mine region.
In addition, after obtaining above-mentioned assessment result, the embodiment of the present application can also be exported above-mentioned assessment result to being
System platform or user terminal, in order to which staff obtains above-mentioned assessment result in time, and then can take corresponding reply
Measure.Also, the embodiment of the present application can be according to These parameters as a result, generating state of ecological environment remoteensing evaluation technological achievement point
Analysis report, and above-mentioned analysis report is sent to above-mentioned user.
In embodiments of the present invention, pass through the remotely-sensed data and initial environment data of acquisition target area and above-mentioned mesh
Mark the evaluation index in region and the weight parameter of above-mentioned evaluation index;Above-mentioned remotely-sensed data and above-mentioned initial environment data are carried out
Pretreatment, obtains preprocessed data, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface area;
According to above-mentioned evaluation index and above-mentioned weight parameter, above-mentioned preprocessed data is assessed, is obtained for evaluating above-mentioned target
The assessment result of the state of ecological environment in region.
The state of ecological environment of accurate and comprehensive objective assessment target area is achieved the purpose that as a result, to realize
The technical effect of the accuracy of the assessment result of ecological environment is improved, and then solves Ecology Environment Evaluation in the prior art
The technical issues of data surface that method obtains is narrow, and evaluation index is single, leads to not accurate evaluation state of ecological environment.
In an alternative embodiment, the method for assessment state of ecological environment provided by the present application can be, but not limited to answer
In assessment for the state of ecological environment of coal mine, the essential core of the method for the assessment state of ecological environment due to restricting coal mine
Be the different ecologicals element such as the land use for how determining coal mine, vegetative coverage, the soil erosion and landscape ecological the situation factor,
Evaluation selecting index and weight are determined and index calculating method.
Under the requirement of above-mentioned core point, the application is adopted under the support of geographic information system technology using remotely-sensed data
The mode combined with man computer interactive interpretation with computer automated information retrieval, indoor comprehensive research are mutually tied with on-site inspection verifying
The operating mode of conjunction realizes the assessment of state of ecological environment.
In an alternative embodiment, Fig. 2 is a kind of optional assessment ecological environment shape according to an embodiment of the present invention
The flow chart of the method for condition, as shown in Fig. 2, above-mentioned assessment result includes: state of ecological environment, according to above-mentioned evaluation index and upper
Weight parameter is stated, above-mentioned preprocessed data is assessed, obtains the state of ecological environment for evaluating above-mentioned target area
Assessment result, comprising:
Step S202 is calculated above-mentioned evaluation index and above-mentioned weight parameter by the way of weighted sum, obtained
The state of ecological environment index of above-mentioned target area.
In the embodiment of the present application, the calculation method of above-mentioned state of ecological environment index can be, but not limited to ask using weighting
And method, the corresponding weight parameter of indices is set according to expert evaluation, when being calculated, is obtained corresponding with evaluation index
Weight parameter, wherein the weighted value of every evaluation number can be, but not limited to as shown in table 1:
Table 1
In the embodiment of the present application, above-mentioned state of ecological environment formula of index can be, but not limited to as follows:
CMEI=0.35 × HQI+0.25 × VCI+0.25 × (100-LDI)+0.15 × (100-LSDI);
Wherein, CMEI is the state of ecological environment index of above-mentioned target area;HQI is the ecological environment of above-mentioned target area
Performance figure;VCI is the vegetation coverage index of above-mentioned target area;LDI is the land degradation index of above-mentioned target area;
LSDI is the destruction of surface index of above-mentioned target area.
Step S204 is assessed according to the above-mentioned above-mentioned preprocessed data of state of ecological environment exponent pair, obtains above-mentioned mesh
Mark the state of ecological environment in region.
In an alternative embodiment, above-mentioned state of ecological environment includes at least: the state of ecological environment of first level,
State of ecological environment, the state of ecological environment of third level, the state of ecological environment of fourth level, fifth level of second level
State of ecological environment.
According to coal mine state of ecological environment index, coal mine ecological environment is divided into five grades, may include: first level
State of ecological environment (grade be excellent), the state of ecological environment of second level (grade is good), third level ecological environment
State of ecological environment (etc. of situation (grade is general), the state of ecological environment of fourth level (grade is poor), fifth level
Grade is poor).Wherein, the classification of coal mine state of ecological environment, evaluation number range and descriptive grade etc. can be, but not limited to such as table 2
It is shown:
Table 2
In an alternative embodiment, the amplitude of variation of above-mentioned coal mine state of ecological environment is included at least: first becomes
Change amplitude, the second amplitude of variation, third amplitude of variation, the 4th amplitude of variation;The rate of change of above-mentioned coal mine state of ecological environment
It includes at least: the first rate of change, the second rate of change, third rate of change, the 4th rate of change.
In the above embodiments of the present application, the amplitude of variation of above-mentioned coal mine state of ecological environment can be △ CMEI, wherein
The amplitude of variation △ CMEI of coal mine state of ecological environment is the difference of current value and a reference value, can be incited somebody to action in the embodiment of the present application
△ CMEI points are four grades, comprising: the first amplitude of variation (no significant change), the second amplitude of variation (being slightly changed), third
Amplitude of variation (significant change), the 4th amplitude of variation (significant changes), the classification of coal mine state of ecological environment amplitude of variation, variation width
Degree and variation characteristic can be as shown in Figure 3:
Table 3
In the above embodiments of the present application, the rate of change (V of above-mentioned coal mine state of ecological environment△CMEI) it is status benchmark
The ratio of difference and a reference value, by V△CMEIIt is divided into level Four, i.e. the first rate of change (rate is without significant change), the second variation speed
Rate (rate is slightly changed), third rate of change (rate significant change), the 4th rate of change (rate significant changes), above-mentioned coal
The rate of change of mine state of ecological environment is classified, rate of change and variation characteristic are as shown in table 4:
Table 4
In an alternative embodiment, above-mentioned eco-environmental quality index HQI is calculated by following formula:
HQI=Abio×(Slw×Sl+Scw×Sc+Ssw×Ss+Sgw×Sg+Sjw×Sj+Sww×Sw)/SCM;
Wherein, AbioFor the normalization coefficient of the eco-environmental quality index of above-mentioned target area;SlFor above-mentioned target area
Forest land area;SlwFor the weighted value in above-mentioned forest land;ScFor the grassland area of above-mentioned target area;ScwFor the weight on above-mentioned meadow
Value;SsFor the Wetland Area area of above-mentioned target area;SswFor the weighted value of above-mentioned Wetland Area;SgFor above-mentioned target area
Cultivated area;SgwFor the weighted value in above-mentioned arable land;SjFor the construction land area of above-mentioned target area;SjwFor above-mentioned construction land
Weighted value;SwFor the unused land area of above-mentioned target area;SwwFor the weighted value of above-mentioned unused land;SCMFor above-mentioned target
The gross area in region.
In the above-described embodiments, above-mentioned Sl、Sc、Ss、Sg、Sw、SCMUnit all can be hm2。
In the embodiment of the present application, it can be, but not limited to that above-mentioned eco-environmental quality is calculated using weighted sum method
Index HQI, the type of ecological environment and corresponding weighted value are as shown in table 5:
Table 5
Type | Forest land | Meadow | Wetland Area | Arable land | Construction land | Unused land |
Weighted value | 0.35 | 0.30 | 0.20 | 0.10 | 0.04 | 0.01 |
As a kind of optional embodiment, above-mentioned SlwIt is 0.35, above-mentioned ScwIt is 0.30, above-mentioned SswIt is 0.20, above-mentioned SgwFor
0.10, above-mentioned SjwIt is 0.04, above-mentioned SwwIt is 0.01.
In an alternative embodiment, above-mentioned vegetation coverage index VCI is calculated by following formula:
VCI=Avfc×(Svfc1w×Svfc1+Svfc2w×Svfc2+Svfc3w×Svfc3+Svfc4w×Svfc4+Svfc5w×Svfc5)/
SCM;
Wherein, AvfcFor the normalization coefficient of the vegetation coverage index of above-mentioned target area;Svfc1For above-mentioned target area
High vegetation coverage area;Svfc1wFor the weighted value of above-mentioned high vegetation coverage;Svfc2For the middle Gao Zhi of above-mentioned target area
Coating cover degree area;Svfc2wFor the weighted value of high vegetation coverage among the above;Svfc3For the middle vegetative coverage of above-mentioned target area
Spend area;Svfc3wFor the weighted value of vegetation coverage among the above;Svfc4For the low vegetation coverage area of above-mentioned target area;
Svfc4wFor the weighted value of above-mentioned low vegetation coverage;Svfc5For the extremely low vegetation coverage area of above-mentioned target area;Svfc5wFor
The weighted value of above-mentioned extremely low vegetation coverage;SCMFor the gross area of above-mentioned target area.
In the above-described embodiments, above-mentioned Svfc1、Svfc2、Svfc3、Svfc4、Svfc5、SCMUnit all can be hm2.In this Shen
It please can be, but not limited to that above-mentioned vegetation coverage index VCI, vegetative coverage is calculated using weighted sum method in embodiment
The type of degree and corresponding weighted value are as shown in table 6:
Table 6
As a kind of optional embodiment, above-mentioned Svfc1wIt is 0.30, above-mentioned Svfc2wIt is 0.25, above-mentioned Svfc3wIt is 0.20, on
State Svfc4wIt is 0.15, above-mentioned Svfc1wIt is 0.10.
In an alternative embodiment, above-mentioned land degradation index LDI is calculated by following formula:
LDI=Aero×(Sero1w×Sero1+Sero2w×Sero2+Sero3w×Sero3)/SCM;
Wherein, AeroFor the normalization coefficient of the land degradation index of above-mentioned target area;Sero1For above-mentioned target area
The area of slight eroded field;Sero1wFor the weighted value of above-mentioned slight eroded field;Sero2Moderate for above-mentioned target area is invaded
Lose the area in soil;Sero2wFor the weighted value of above-mentioned moderate eroded field;Sero3For the severe eroded field of above-mentioned target area
Area;Sero3wFor the weighted value of above-mentioned severe eroded field;SCMFor the gross area of above-mentioned target area.
In the above-described embodiments, above-mentioned Sero1、Sero2、Sero3、SCMUnit all can be hm2.In the embodiment of the present application
In, it can be, but not limited to that above-mentioned land degradation index LDI, the type of land deterioration and institute is calculated using weighted sum method
Corresponding weighted value is as shown in table 7:
Table 7
As a kind of optional embodiment, above-mentioned Sero1wIt is 0.05, above-mentioned Sero2wIt is 0.25, above-mentioned Sero3wIt is 0.70.
In an alternative embodiment, above-mentioned destruction of surface index LSDI is included at least: the destruction of surface of underground mine
Index LSDI1 calculates the destruction of surface index LSDI1 of above-mentioned underground mine by following formula:
LSDI1=(A × Asus×Ssus)/Scm+(B×Acmw×Scmw)/Scm;
Wherein, AsusFor the normalization coefficient of the surface collapse of above-mentioned target area;AcmwFor the discharge refuse of above-mentioned target area
The normalization coefficient of field;SsusFor the area of the surface collapse of above-mentioned target area;ScmwFor the discharge refuse field of above-mentioned target area
Area;A is the weighted value of above-mentioned surface collapse;B is the weighted value of above-mentioned discharge refuse field;SCMFor the gross area of above-mentioned target area.
In the above-described embodiments, above-mentioned Ssus、Scmw、SCMUnit all can be hm2.It in the embodiment of the present application, can be with
But it is not limited to that the destruction of surface index LSDI1 of above-mentioned underground mine, the earth's surface of underground mine is calculated using weighted sum method
The type of destruction and corresponding weighted value are as shown in table 8:
Table 8
Destruction of surface type | Surface collapse/excavating plant | Discharge refuse field/refuse dump |
Weighted value | 0.50 | 0.50 |
As a kind of optional embodiment, the weighted value A of above-mentioned surface collapse is 0.50, the weighted value B of above-mentioned discharge refuse field
It is 0.50.
In an alternative embodiment, above-mentioned destruction of surface index LSDI is included at least: the table destruction of open coal mine refers to
Number LSDI2, the destruction of surface index LSDI2 of above-mentioned open coal mine is calculated by following formula:
LSDI2=(Sldsw×Alds×Slds)/Scm+(Scmdw×Acmd×Scmd)/Scm;
Wherein, AldsFor the normalization coefficient of the excavating plant of above-mentioned target area;AcmdFor the casting of above-mentioned target area
The normalization coefficient of field;SldsFor the area of the excavating plant of above-mentioned target area;SldswFor the weighted value of above-mentioned excavating plant;
ScmdFor the area of the refuse dump of above-mentioned target area;ScmdwFor the weighted value of above-mentioned refuse dump;SCMFor the total of above-mentioned target area
Area.
In the above-described embodiments, above-mentioned Slds、Scmd、SCMUnit all can be hm2.It in the embodiment of the present application, can be with
But it is not limited to that the destruction of surface index LSDI2 of above-mentioned open coal mine, the earth's surface of open coal mine is calculated using weighted sum method
Shown in the type of destruction and corresponding weighted value table 8 as above:
As a kind of optional embodiment, the weighted value S of above-mentioned surface collapseldswIt is 0.50, the weight of above-mentioned discharge refuse field
Value ScmdwIt is 0.50.
It should be noted that before land reclamation or below standard after land reclamation being included in brokenly when destruction of surface index calculates
Bad area, it is up to standard after reclaiming to be not counted in damage area;What opencast coal mine dump was arranged in realizing, when excavating plant areal calculation
It does not include internal dumping area.
The embodiment of the present application can provide comprehensive, system evaluation side for the state of ecological environment of the zonules such as coal mine
Law system provides technical support for comprehensive, objective appraisal state of ecological environment;Also, the embodiment of the present application is from the reality of coal mine
The principal element that border production practice process and environment influence considers, has combed 4 sport, 16 points of indexs and weight, and earth's surface is broken
Bad index index is subdivided into underground mine and open coal mine, be it is perfect to the supplement of current evaluation method system so that commenting
Valence result is more objective comprehensively.
Also, the embodiment of the present application considers from the Spatial dimensionality angle of coal mine ecological environment dynamic change, by coal mine ecology
Environmental aspect classification is included in appraisement system, coal mine state of ecological environment amplitude of variation and rate index, embodies coal mine ecology
Variation and lasting comparison of the environmental aspect in Spatial dimensionality;In addition, the embodiment of the present application is applicable not only to coal mine ecological environment shape
Condition remoteensing evaluation is also applied for non-coal mine state of ecological environment remoteensing evaluation, has general applicability and generalization.
Embodiment 2
According to embodiments of the present invention, it additionally provides a kind of for implementing the device of the method for above-mentioned assessment state of ecological environment
Embodiment, Fig. 3 is a kind of structural schematic diagram of device for assessing state of ecological environment according to an embodiment of the present invention, such as Fig. 3 institute
Show, the device of above-mentioned assessment state of ecological environment, comprising: obtain module 30, preprocessing module 32 and evaluation module 34, in which:
Module 30 is obtained, for obtaining remotely-sensed data and initial environment data and the above-mentioned target area of target area
Evaluation index and above-mentioned evaluation index weight parameter;Preprocessing module 32, for above-mentioned remotely-sensed data and above-mentioned initial
Environmental data is pre-processed, and preprocessed data is obtained, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and
Destruction of surface area;Evaluation module 34, for according to above-mentioned evaluation index and above-mentioned weight parameter, to above-mentioned preprocessed data into
Row assessment, obtains the assessment result of the state of ecological environment for evaluating above-mentioned target area.
It should be noted that above-mentioned modules can be realized by software or hardware, for example, for the latter,
Can be accomplished by the following way: above-mentioned modules can be located in same processor;Alternatively, above-mentioned modules are with any
Combined mode is located in different processors.
Herein it should be noted that above-mentioned acquisition module 30, preprocessing module 32 and evaluation module 34 correspond to embodiment 1
In step S102 to step S106, the example and application scenarios that above-mentioned module and corresponding step are realized be identical but unlimited
In 1 disclosure of that of above-described embodiment.It should be noted that above-mentioned module may operate in calculating as a part of device
In machine terminal.
It should be noted that the optional or preferred embodiment of the present embodiment may refer to the associated description in embodiment 1,
Details are not described herein again.
The device of above-mentioned assessment state of ecological environment can also include processor and memory, above-mentioned assessment ecological environment
Device of situation etc. stores in memory as program unit, executes above procedure stored in memory by processor
Unit realizes corresponding function.
Include kernel in processor, is gone in memory to transfer corresponding program unit by kernel, above-mentioned kernel can be set
One or more.Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM)
And/or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM), memory includes at least one
Storage chip.
According to the embodiment of the present application, a kind of storage medium embodiment is additionally provided.Optionally, in the present embodiment, above-mentioned
Storage medium includes the program of storage, wherein equipment where controlling above-mentioned storage medium in above procedure operation executes above-mentioned
The method of any one assessment state of ecological environment.
Optionally, in the present embodiment, above-mentioned storage medium can be located in computer network in computer terminal group
In any one terminal, or in any one mobile terminal in mobile terminal group, above-mentioned storage medium packet
Include the program of storage.
Optionally, when program is run, equipment where control storage medium executes following functions: obtaining the distant of target area
Feel the weight parameter of the evaluation index and above-mentioned evaluation index of data and initial environment data and above-mentioned target area;To upper
It states remotely-sensed data and above-mentioned initial environment data is pre-processed, obtain preprocessed data, wherein above-mentioned preprocessed data is at least
It include: vegetative coverage area and destruction of surface area;According to above-mentioned evaluation index and above-mentioned weight parameter, to above-mentioned pretreatment number
According to being assessed, the assessment result of the state of ecological environment for evaluating above-mentioned target area is obtained.
Optionally, when program is run, equipment where control storage medium executes following functions: to above-mentioned evaluation index and
Above-mentioned weight parameter is calculated by the way of weighted sum, obtains the state of ecological environment index of above-mentioned target area;According to
It is assessed according to the above-mentioned above-mentioned preprocessed data of state of ecological environment exponent pair, obtains the ecological environment shape of above-mentioned target area
Condition, wherein above-mentioned state of ecological environment includes at least: the ecological environment shape of the state of ecological environment of first level, second level
Condition, the state of ecological environment of third level, the state of ecological environment of fourth level, the state of ecological environment of fifth level;It is above-mentioned
The amplitude of variation of state of ecological environment includes at least: the first amplitude of variation, the second amplitude of variation, third amplitude of variation, the 4th become
Change amplitude;The rate of change of above-mentioned state of ecological environment includes at least: the first rate of change, the second rate of change, third variation
Rate, the 4th rate of change.
Optionally, when program is run, equipment where control storage medium executes following functions: being calculated by following formula
Above-mentioned eco-environmental quality index HQI:HQI=Abio×(Slw×Sl+Scw×Sc+Ssw×Ss+Sgw×Sg+Sjw×Sj+Sww×
Sw)/SCM;Wherein, AbioFor the normalization coefficient of the eco-environmental quality index of above-mentioned target area;SlFor above-mentioned target area
Forest land area;SlwFor the weighted value in above-mentioned forest land;ScFor the grassland area of above-mentioned target area;ScwFor the weight on above-mentioned meadow
Value;SsFor the Wetland Area area of above-mentioned target area;SswFor the weighted value of above-mentioned Wetland Area;SgFor above-mentioned target area
Cultivated area;SgwFor the weighted value in above-mentioned arable land;SjFor the construction land area of above-mentioned target area;SjwFor above-mentioned construction land
Weighted value;SwFor the unused land area of above-mentioned target area;SwwFor the weighted value of above-mentioned unused land;SCMFor above-mentioned target
The gross area in region.
Optionally, when program is run, equipment where control storage medium executes following functions: being calculated by following formula
Above-mentioned vegetation coverage index VCI:VCI=Avfc×(Svfc1w×Svfc1+Svfc2w×Svfc2+Svfc3w×Svfc3+Svfc4w×Svfc4+
Svfc5w×Svfc5)/SCM;Wherein, AvfcFor the normalization coefficient of the vegetation coverage index of above-mentioned target area;Svfc1It is above-mentioned
The high vegetation coverage area of target area;Svfc1wFor the weighted value of above-mentioned high vegetation coverage;Svfc2For above-mentioned target area
In high vegetation coverage area;Svfc2wFor the weighted value of high vegetation coverage among the above;Svfc3For in above-mentioned target area
Vegetation coverage area;Svfc3wFor the weighted value of vegetation coverage among the above;Svfc4For the low vegetative coverage of above-mentioned target area
Spend area;Svfc4wFor the weighted value of above-mentioned low vegetation coverage;Svfc5For the extremely low vegetation coverage area of above-mentioned target area;
Svfc5wFor the weighted value of above-mentioned extremely low vegetation coverage;SCMFor the gross area of above-mentioned target area.
Optionally, when program is run, equipment where control storage medium executes following functions: being calculated by following formula
Above-mentioned land degradation index LDI:LDI=Aero×(Sero1w×Sero1+Sero2w×Sero2+Sero3w×Sero3)/SCM;Wherein, Aero
For the normalization coefficient of the land degradation index of above-mentioned target area;Sero1For the face of the slight eroded field of above-mentioned target area
Product;Sero1wFor the weighted value of above-mentioned slight eroded field;Sero2For the area of the moderate eroded field of above-mentioned target area;Sero2w
For the weighted value of above-mentioned moderate eroded field;Sero3For the area of the severe eroded field of above-mentioned target area;Sero3wIt is above-mentioned
The weighted value of severe eroded field;SCMFor the gross area of above-mentioned target area.
Optionally, when program is run, equipment where control storage medium executes following functions: being calculated by following formula
Destruction of surface index LSDI1:LSDI1=(A × A of above-mentioned underground minesus×Ssus)/Scm+(B×Acmw×Scmw)/Scm;Its
In, AsusFor the normalization coefficient of the surface collapse of above-mentioned target area;AcmwFor the normalization of the discharge refuse field of above-mentioned target area
Coefficient;SsusFor the area of the surface collapse of above-mentioned target area;ScmwFor the area of the discharge refuse field of above-mentioned target area;A is upper
State the weighted value of surface collapse;B is the weighted value of above-mentioned discharge refuse field;SCMFor the gross area of above-mentioned target area.
Optionally, when program is run, equipment where control storage medium executes following functions: being calculated by following formula
The destruction of surface index LSDI2:LSDI2=(S of above-mentioned open coal mineldsw×Alds×Slds)/Scm+(Scmdw×Acmd×Scmd)/
Scm;Wherein, AldsFor the normalization coefficient of the excavating plant of above-mentioned target area;AcmdReturn for the refuse dump of above-mentioned target area
One changes coefficient;SldsFor the area of the excavating plant of above-mentioned target area;SldswFor the weighted value of above-mentioned excavating plant;ScmdIt is upper
State the area of the refuse dump of target area;ScmdwFor the weighted value of above-mentioned refuse dump;SCMFor the gross area of above-mentioned target area.
According to the embodiment of the present application, a kind of processor embodiment is additionally provided.Optionally, in the present embodiment, above-mentioned place
Reason device is for running program, wherein above procedure executes any one of the above assessment state of ecological environment method when running.
The embodiment of the present application provides a kind of equipment, equipment include processor, memory and storage on a memory and can
The program run on a processor, processor execute program when perform the steps of obtain target area remotely-sensed data and just
The weight parameter of the evaluation index and above-mentioned evaluation index of beginning environmental data and above-mentioned target area;To above-mentioned remotely-sensed data
It is pre-processed with above-mentioned initial environment data, obtains preprocessed data, wherein above-mentioned preprocessed data includes at least: vegetation
Area coverage and destruction of surface area;According to above-mentioned evaluation index and above-mentioned weight parameter, above-mentioned preprocessed data is commented
Estimate, obtains the assessment result of the state of ecological environment for evaluating above-mentioned target area.
Optionally, when above-mentioned processor executes program, above-mentioned evaluation index and above-mentioned weight parameter can also be used and is added
The mode of power summation is calculated, and the state of ecological environment index of above-mentioned target area is obtained;According to above-mentioned state of ecological environment
The above-mentioned preprocessed data of exponent pair is assessed, and the state of ecological environment of above-mentioned target area is obtained, wherein above-mentioned ecological environment
Situation includes at least: the ecological environment of the state of ecological environment of first level, the state of ecological environment of second level, third level
Situation, the state of ecological environment of fourth level, the state of ecological environment of fifth level;The amplitude of variation of above-mentioned state of ecological environment
It includes at least: the first amplitude of variation, the second amplitude of variation, third amplitude of variation, the 4th amplitude of variation;Above-mentioned state of ecological environment
Rate of change include at least: the first rate of change, the second rate of change, third rate of change, the 4th rate of change.
Optionally, when above-mentioned processor executes program, above-mentioned eco-environmental quality can also be calculated by following formula and referred to
Number HQI:HQI=Abio×(Slw×Sl+Scw×Sc+Ssw×Ss+Sgw×Sg+Sjw×Sj+Sww×Sw)/SCM;Wherein, AbioIt is above-mentioned
The normalization coefficient of the eco-environmental quality index of target area;SlFor the forest land area of above-mentioned target area;SlwFor above-mentioned woods
The weighted value on ground;ScFor the grassland area of above-mentioned target area;ScwFor the weighted value on above-mentioned meadow;SsFor above-mentioned target area
Wetland Area area;SswFor the weighted value of above-mentioned Wetland Area;SgFor the cultivated area of above-mentioned target area;SgwFor above-mentioned arable land
Weighted value;SjFor the construction land area of above-mentioned target area;SjwFor the weighted value of above-mentioned construction land;SwFor above-mentioned target
The unused land area in region;SwwFor the weighted value of above-mentioned unused land;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned processor executes program, above-mentioned vegetation coverage index can also be calculated by following formula
VCI:VCI=Avfc×(Svfc1w×Svfc1+Svfc2w×Svfc2+Svfc3w×Svfc3+Svfc4w×Svfc4+Svfc5w×Svfc5)/SCM;Its
In, AvfcFor the normalization coefficient of the vegetation coverage index of above-mentioned target area;Svfc1It is covered for the high vegetation of above-mentioned target area
Cover degree area;Svfc1wFor the weighted value of above-mentioned high vegetation coverage;Svfc2For vegetation coverage face high in above-mentioned target area
Product;Svfc2wFor the weighted value of high vegetation coverage among the above;Svfc3For the middle vegetation coverage area of above-mentioned target area;Svfc3w
For the weighted value of vegetation coverage among the above;Svfc4For the low vegetation coverage area of above-mentioned target area;Svfc4wIt is above-mentioned low
The weighted value of vegetation coverage;Svfc5For the extremely low vegetation coverage area of above-mentioned target area;Svfc5wFor above-mentioned extremely low vegetation
The weighted value of coverage;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned processor executes program, above-mentioned land degradation index can also be calculated by following formula
LDI:LDI=Aero×(Sero1w×Sero1+Sero2w×Sero2+Sero3w×Sero3)/SCM;Wherein, AeroFor above-mentioned target area
The normalization coefficient of land degradation index;Sero1For the area of the slight eroded field of above-mentioned target area;Sero1wIt is above-mentioned light
Spend the weighted value of eroded field;Sero2For the area of the moderate eroded field of above-mentioned target area;Sero2wFor the erosion of above-mentioned moderate
The weighted value in soil;Sero3For the area of the severe eroded field of above-mentioned target area;Sero3wFor above-mentioned severe eroded field
Weighted value;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned processor executes program, the earth's surface of above-mentioned underground mine can also be calculated by following formula
Destructing index LSDI1:LSDI1=(A × Asus×Ssus)/Scm+(B×Acmw×Scmw)/Scm;Wherein, AsusFor above-mentioned target area
The normalization coefficient of the surface collapse in domain;AcmwFor the normalization coefficient of the discharge refuse field of above-mentioned target area;SsusFor above-mentioned target
The area of the surface collapse in region;ScmwFor the area of the discharge refuse field of above-mentioned target area;A is the weighted value of above-mentioned surface collapse;
B is the weighted value of above-mentioned discharge refuse field;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned processor executes program, the earth's surface of above-mentioned open coal mine can also be calculated by following formula
Destructing index LSDI2:LSDI2=(Sldsw×Alds×Slds)/Scm+(Scmdw×Acmd×Scmd)/Scm;Wherein, AldsFor above-mentioned mesh
Mark the normalization coefficient of the excavating plant in region;AcmdFor the normalization coefficient of the refuse dump of above-mentioned target area;SldsIt is above-mentioned
The area of the excavating plant of target area;SldswFor the weighted value of above-mentioned excavating plant;ScmdFor the refuse dump of above-mentioned target area
Area;ScmdwFor the weighted value of above-mentioned refuse dump;SCMFor the gross area of above-mentioned target area.
Present invention also provides a kind of computer program products, when executing on data processing equipment, are adapted for carrying out just
The program of beginningization there are as below methods step: the remotely-sensed data and initial environment data and above-mentioned target area of target area are obtained
The weight parameter of the evaluation index in domain and above-mentioned evaluation index;Above-mentioned remotely-sensed data and above-mentioned initial environment data are located in advance
Reason, obtains preprocessed data, wherein above-mentioned preprocessed data includes at least: vegetative coverage area and destruction of surface area;Foundation
Above-mentioned evaluation index and above-mentioned weight parameter, assess above-mentioned preprocessed data, obtain for evaluating above-mentioned target area
State of ecological environment assessment result.
Optionally, when above-mentioned computer program product executes program, above-mentioned evaluation index and above-mentioned weight can also be joined
Number is calculated by the way of weighted sum, obtains the state of ecological environment index of above-mentioned target area;According to above-mentioned ecology
The above-mentioned preprocessed data of environmental aspect exponent pair is assessed, and obtains the state of ecological environment of above-mentioned target area, wherein above-mentioned
State of ecological environment includes at least: the state of ecological environment of first level, the state of ecological environment of second level, third level
State of ecological environment, the state of ecological environment of fourth level, the state of ecological environment of fifth level;Above-mentioned state of ecological environment
Amplitude of variation includes at least: the first amplitude of variation, the second amplitude of variation, third amplitude of variation, the 4th amplitude of variation;Above-mentioned ecology
The rate of change of environmental aspect includes at least: the first rate of change, the second rate of change, third rate of change, the 4th variation speed
Rate.
Optionally, when above-mentioned computer program product executes program, above-mentioned ecological ring can also be calculated by following formula
Border performance figure HQI:HQI=Abio×(Slw×Sl+Scw×Sc+Ssw×Ss+Sgw×Sg+Sjw×Sj+Sww×Sw)/SCM;Wherein,
AbioFor the normalization coefficient of the eco-environmental quality index of above-mentioned target area;SlFor the forest land area of above-mentioned target area;Slw
For the weighted value in above-mentioned forest land;ScFor the grassland area of above-mentioned target area;ScwFor the weighted value on above-mentioned meadow;SsFor above-mentioned mesh
Mark the Wetland Area area in region;SswFor the weighted value of above-mentioned Wetland Area;SgFor the cultivated area of above-mentioned target area;SgwFor
The weighted value in above-mentioned arable land;SjFor the construction land area of above-mentioned target area;SjwFor the weighted value of above-mentioned construction land;SwFor
The unused land area of above-mentioned target area;SwwFor the weighted value of above-mentioned unused land;SCMFor total face of above-mentioned target area
Product.
Optionally, when above-mentioned computer program product executes program, above-mentioned vegetation can also be calculated by following formula and covered
Cover degree index VCI:VCI=Avfc×(Svfc1w×Svfc1+Svfc2w×Svfc2+Svfc3w×Svfc3+Svfc4w×Svfc4+Svfc5w×
Svfc5)/SCM;Wherein, AvfcFor the normalization coefficient of the vegetation coverage index of above-mentioned target area;Svfc1For above-mentioned target area
The high vegetation coverage area in domain;Svfc1wFor the weighted value of above-mentioned high vegetation coverage;Svfc2For the middle height of above-mentioned target area
Vegetation coverage area;Svfc2wFor the weighted value of high vegetation coverage among the above;Svfc3It is covered for the middle vegetation of above-mentioned target area
Cover degree area;Svfc3wFor the weighted value of vegetation coverage among the above;Svfc4For the low vegetation coverage area of above-mentioned target area;
Svfc4wFor the weighted value of above-mentioned low vegetation coverage;Svfc5For the extremely low vegetation coverage area of above-mentioned target area;Svfc5wFor
The weighted value of above-mentioned extremely low vegetation coverage;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned computer program product executes program, above-mentioned soil can also be calculated by following formula and moved back
Change index LDI:LDI=Aero×(Sero1w×Sero1+Sero2w×Sero2+Sero3w×Sero3)/SCM;Wherein, AeroFor above-mentioned target
The normalization coefficient of the land degradation index in region;Sero1For the area of the slight eroded field of above-mentioned target area;Sero1wFor
The weighted value of above-mentioned slight eroded field;Sero2For the area of the moderate eroded field of above-mentioned target area;Sero2wFor among the above
Spend the weighted value of eroded field;Sero3For the area of the severe eroded field of above-mentioned target area;Sero3wFor the erosion of above-mentioned severe
The weighted value in soil;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned computer program product executes program, above-mentioned well work coal can also be calculated by following formula
Destruction of surface index LSDI1:LSDI1=(A × A of minesus×Ssus)/Scm+(B×Acmw×Scmw)/Scm;Wherein, AsusIt is upper
State the normalization coefficient of the surface collapse of target area;AcmwFor the normalization coefficient of the discharge refuse field of above-mentioned target area;SsusFor
The area of the surface collapse of above-mentioned target area;ScmwFor the area of the discharge refuse field of above-mentioned target area;A is above-mentioned surface collapse
Weighted value;B is the weighted value of above-mentioned discharge refuse field;SCMFor the gross area of above-mentioned target area.
Optionally, when above-mentioned computer program product executes program, above-mentioned outdoor coal can also be calculated by following formula
The destruction of surface index LSDI2:LSDI2=(S of mineldsw×Alds×Slds)/Scm+(Scmdw×Acmd×Scmd)/Scm;Wherein, Alds
For the normalization coefficient of the excavating plant of above-mentioned target area;AcmdFor the normalization coefficient of the refuse dump of above-mentioned target area;
SldsFor the area of the excavating plant of above-mentioned target area;SldswFor the weighted value of above-mentioned excavating plant;ScmdFor above-mentioned target area
The area of the refuse dump in domain;ScmdwFor the weighted value of above-mentioned refuse dump;SCMFor the gross area of above-mentioned target area.
The serial number of the above embodiments of the invention is only for description, does not represent the advantages or disadvantages of the embodiments.
In the above embodiment of the invention, it all emphasizes particularly on different fields to the description of each embodiment, does not have in some embodiment
The part of detailed description, reference can be made to the related descriptions of other embodiments.
In several embodiments provided herein, it should be understood that disclosed technology contents can pass through others
Mode is realized.Wherein, the apparatus embodiments described above are merely exemplary, such as the division of the unit, Ke Yiwei
A kind of logical function partition, there may be another division manner in actual implementation, for example, multiple units or components can combine or
Person is desirably integrated into another system, or some features can be ignored or not executed.Another point, shown or discussed is mutual
Between coupling, direct-coupling or communication connection can be through some interfaces, the INDIRECT COUPLING or communication link of unit or module
It connects, can be electrical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
On unit.It can some or all of the units may be selected to achieve the purpose of the solution of this embodiment according to the actual needs.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit
It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list
Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product
When, it can store in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially
The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words
It embodies, which is stored in a storage medium, including some instructions are used so that a computer
Equipment (can for personal computer, server or network equipment etc.) execute each embodiment the method for the present invention whole or
Part steps.And storage medium above-mentioned includes: that USB flash disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), mobile hard disk, magnetic or disk etc. be various to can store program code
Medium.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (11)
1. a kind of method for assessing state of ecological environment characterized by comprising
The remotely-sensed data and initial environment data and the evaluation index and the evaluation of the target area of acquisition target area
The weight parameter of index;
The remotely-sensed data and the initial environment data are pre-processed, preprocessed data is obtained, wherein the pretreatment
Data include at least: vegetative coverage area and destruction of surface area;
According to the evaluation index and the weight parameter, the preprocessed data is assessed, is obtained described for evaluating
The assessment result of the state of ecological environment of target area.
2. the method according to claim 1, wherein the assessment result includes: state of ecological environment, according to institute
Evaluation index and the weight parameter are stated, the preprocessed data is assessed, is obtained for evaluating the target area
The assessment result of state of ecological environment, comprising:
The evaluation index and the weight parameter are calculated by the way of weighted sum, obtain the target area
State of ecological environment index;
It is assessed according to preprocessed data described in the state of ecological environment exponent pair, obtains the ecological ring of the target area
Border situation, wherein the state of ecological environment includes at least: the ecological ring of the state of ecological environment of first level, second level
Border situation, the state of ecological environment of third level, the state of ecological environment of fourth level, the state of ecological environment of fifth level;
The amplitude of variation of the state of ecological environment includes at least: the first amplitude of variation, the second amplitude of variation, third change width
Degree, the 4th amplitude of variation;The rate of change of the state of ecological environment includes at least: the first rate of change, the second rate of change,
Third rate of change, the 4th rate of change.
3. the method according to claim 1, wherein the evaluation index includes at least one following: ecological ring
Border performance figure HQI, vegetation coverage index VCI, land degradation index LDI and destruction of surface index LSDI.
4. according to the method described in claim 3, it is characterized in that, calculating the eco-environmental quality index by following formula
HQI:
HQI=Abio×(Slw×Sl+Scw×Sc+Ssw×Ss+Sgw×Sg+Sjw×Sj+Sww×Sw)/SCM;
Wherein, AbioFor the normalization coefficient of the eco-environmental quality index of the target area;SlFor the woods of the target area
Ground area;SlwFor the weighted value in the forest land;ScFor the grassland area of the target area;ScwFor the weighted value on the meadow;
SsFor the Wetland Area area of the target area;SswFor the weighted value of the Wetland Area;SgFor ploughing for the target area
Ground area;SgwFor the weighted value in the arable land;SjFor the construction land area of the target area;SjwFor the construction land
Weighted value;SwFor the unused land area of the target area;SwwFor the weighted value of the unused land;SCMFor the target area
The gross area in domain.
5. according to the method described in claim 3, it is characterized in that, calculating the vegetation coverage index by following formula
VCI:
VCI=Avfc×(Svfc1w×Svfc1+Svfc2w×Svfc2+Svfc3w×Svfc3+Svfc4w×Svfc4+Svfc5w×Svfc5)/SCM;
Wherein, AvfcFor the normalization coefficient of the vegetation coverage index of the target area;Svfc1For the height of the target area
Vegetation coverage area;Svfc1wFor the weighted value of the high vegetation coverage;Svfc2It is covered for high vegetation in the target area
Cover degree area;Svfc2wFor the weighted value of the middle high vegetation coverage;Svfc3For the middle vegetation coverage face of the target area
Product;Svfc3wFor the weighted value of the middle vegetation coverage;Svfc4For the low vegetation coverage area of the target area;Svfc4wFor
The weighted value of the low vegetation coverage;Svfc5For the extremely low vegetation coverage area of the target area;Svfc5wFor the pole
The weighted value of low vegetation coverage;SCMFor the gross area of the target area.
6. according to the method described in claim 3, it is characterized in that, calculating the land degradation index LDI by following formula:
LDI=Aero×(Sero1w×Sero1+Sero2w×Sero2+Sero3w×Sero3)/SCM;
Wherein, AeroFor the normalization coefficient of the land degradation index of the target area;Sero1For the slight of the target area
The area of eroded field;Sero1wFor the weighted value of the slight eroded field;Sero2For the moderate eroded soil of the target area
The area on ground;Sero2wFor the weighted value of the moderate eroded field;Sero3For the face of the severe eroded field of the target area
Product;Sero3wFor the weighted value of the severe eroded field;SCMFor the gross area of the target area.
7. according to the method described in claim 3, it is characterized in that, the destruction of surface index LSDI is included at least: well work coal
The destruction of surface index LSDI1 of mine calculates the destruction of surface index LSDI1 of the underground mine by following formula:
LSDI1=(A × Asus×Ssus)/Scm+(B×Acmw×Scmw)/Scm;
Wherein, AsusFor the normalization coefficient of the surface collapse of the target area;AcmwFor the discharge refuse field of the target area
Normalization coefficient;SsusFor the area of the surface collapse of the target area;ScmwFor the area of the discharge refuse field of the target area;
A is the weighted value of the surface collapse;B is the weighted value of the discharge refuse field;SCMFor the gross area of the target area.
8. the method according to the description of claim 7 is characterized in that the destruction of surface index LSDI is included at least: outdoor coal
The table destructing index LSDI2 of mine calculates the destruction of surface index LSDI2 of the open coal mine by following formula:
LSDI2=(Sldsw×Alds×Slds)/Scm+(Scmdw×Acmd×Scmd)/Scm;
Wherein, AldsFor the normalization coefficient of the excavating plant of the target area;AcmdFor the refuse dump of the target area
Normalization coefficient;SldsFor the area of the excavating plant of the target area;SldswFor the weighted value of the excavating plant;ScmdFor
The area of the refuse dump of the target area;ScmdwFor the weighted value of the refuse dump;SCMFor the gross area of the target area.
9. a kind of device for assessing state of ecological environment characterized by comprising
Module is obtained, for obtaining the remotely-sensed data of target area and the evaluation of initial environment data and the target area
The weight parameter of index and the evaluation index;
Preprocessing module obtains preprocessed data for pre-processing to the remotely-sensed data and the initial environment data,
Wherein, the preprocessed data includes at least: vegetative coverage area and destruction of surface area;
Evaluation module, for assessing the preprocessed data, obtaining according to the evaluation index and the weight parameter
For evaluating the assessment result of the state of ecological environment of the target area.
10. a kind of storage medium, which is characterized in that the storage medium includes the program of storage, wherein run in described program
When control the storage medium where equipment execute following steps: obtain the remotely-sensed data and initial environment data of target area,
And the target area evaluation index and the evaluation index weight parameter;To the remotely-sensed data and the initial ring
Border data are pre-processed, and preprocessed data is obtained, wherein the preprocessed data includes at least: vegetative coverage area and ground
Table damage area;According to the evaluation index and the weight parameter, the preprocessed data is assessed, is obtained for commenting
The assessment result of the state of ecological environment of target area described in valence.
11. a kind of processor, which is characterized in that the processor is for running program, wherein executed when described program is run with
Lower step: the remotely-sensed data of target area and the evaluation index of initial environment data and the target area and described are obtained
The weight parameter of evaluation index;The remotely-sensed data and the initial environment data are pre-processed, preprocessed data is obtained,
Wherein, the preprocessed data includes at least: vegetative coverage area and destruction of surface area;According to the evaluation index and described
Weight parameter assesses the preprocessed data, obtains commenting for the state of ecological environment for evaluating the target area
Estimate result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810905136.5A CN109087011A (en) | 2018-08-09 | 2018-08-09 | Assess method and device, storage medium and the processor of state of ecological environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810905136.5A CN109087011A (en) | 2018-08-09 | 2018-08-09 | Assess method and device, storage medium and the processor of state of ecological environment |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109087011A true CN109087011A (en) | 2018-12-25 |
Family
ID=64834566
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810905136.5A Pending CN109087011A (en) | 2018-08-09 | 2018-08-09 | Assess method and device, storage medium and the processor of state of ecological environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109087011A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110553980A (en) * | 2019-07-03 | 2019-12-10 | 中南大学 | multi-index monitoring and evaluating method for rocky desertification control effect based on multi-source remote sensing data |
CN110659838A (en) * | 2019-09-27 | 2020-01-07 | 佛山科学技术学院 | Ecological environment quality assessment method based on remote sensing and GIS |
CN110837875A (en) * | 2019-11-18 | 2020-02-25 | 国家基础地理信息中心 | Method and device for judging quality abnormity of earth surface coverage data |
CN111199355A (en) * | 2019-12-31 | 2020-05-26 | 国家卫星气象中心(国家空间天气监测预警中心) | Evaluation method and device for ecological environment and electronic system |
CN112100579A (en) * | 2020-09-17 | 2020-12-18 | 国家卫星气象中心(国家空间天气监测预警中心) | Remote sensing ecological environment quality evaluation method and device and electronic equipment |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049655A (en) * | 2012-12-18 | 2013-04-17 | 中国矿业大学(北京) | Mining area land ecological damage information acquisition method based on multiple data sources |
CN104103016A (en) * | 2014-05-13 | 2014-10-15 | 杭州师范大学 | Comprehensive evaluation method for ecosystem health of wetland on the basis of remote sensing technology |
CN108229821A (en) * | 2018-01-02 | 2018-06-29 | 中国神华能源股份有限公司 | Appraisal procedure, device, storage medium and the system of mining area ecological environment |
-
2018
- 2018-08-09 CN CN201810905136.5A patent/CN109087011A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103049655A (en) * | 2012-12-18 | 2013-04-17 | 中国矿业大学(北京) | Mining area land ecological damage information acquisition method based on multiple data sources |
CN104103016A (en) * | 2014-05-13 | 2014-10-15 | 杭州师范大学 | Comprehensive evaluation method for ecosystem health of wetland on the basis of remote sensing technology |
CN108229821A (en) * | 2018-01-02 | 2018-06-29 | 中国神华能源股份有限公司 | Appraisal procedure, device, storage medium and the system of mining area ecological environment |
Non-Patent Citations (1)
Title |
---|
徐嘉兴等: "矿区土地生态质量评价及动态变化", 《煤炭学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110553980A (en) * | 2019-07-03 | 2019-12-10 | 中南大学 | multi-index monitoring and evaluating method for rocky desertification control effect based on multi-source remote sensing data |
CN110553980B (en) * | 2019-07-03 | 2022-01-14 | 中南大学 | Multi-index monitoring and evaluating method for rocky desertification control effect based on multi-source remote sensing data |
CN110659838A (en) * | 2019-09-27 | 2020-01-07 | 佛山科学技术学院 | Ecological environment quality assessment method based on remote sensing and GIS |
CN110837875A (en) * | 2019-11-18 | 2020-02-25 | 国家基础地理信息中心 | Method and device for judging quality abnormity of earth surface coverage data |
CN111199355A (en) * | 2019-12-31 | 2020-05-26 | 国家卫星气象中心(国家空间天气监测预警中心) | Evaluation method and device for ecological environment and electronic system |
CN111199355B (en) * | 2019-12-31 | 2023-09-05 | 国家卫星气象中心(国家空间天气监测预警中心) | Ecological environment evaluation method, device and electronic system |
CN112100579A (en) * | 2020-09-17 | 2020-12-18 | 国家卫星气象中心(国家空间天气监测预警中心) | Remote sensing ecological environment quality evaluation method and device and electronic equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109087011A (en) | Assess method and device, storage medium and the processor of state of ecological environment | |
Nussbaum et al. | Evaluation of digital soil mapping approaches with large sets of environmental covariates | |
Morandi et al. | Delimitation of ecological corridors between conservation units in the Brazilian Cerrado using a GIS and AHP approach | |
Di Febbraro et al. | Expert-based and correlative models to map habitat quality: Which gives better support to conservation planning? | |
Mäkelä et al. | Estimation of forest stand volumes by Landsat TM imagery and stand-level field-inventory data | |
Dettmers et al. | A GIS modeling method applied to predicting forest songbird habitat | |
Saalfeld et al. | Distribution and abundance of the feral camel (Camelus dromedarius) in Australia | |
Riffell et al. | Effects of the Conservation Reserve Program on northern bobwhite and grassland birds | |
Ghadimi et al. | Possibility of cooperative management in groundwater resources using an evolutionary hydro-economic simulation-optimization model | |
Ungaro et al. | Assessing and mapping topsoil organic carbon stock at regional scale: a scorpan kriging approach conditional on soil map delineations and land use | |
Robinson et al. | Comparison of alternative strategies for invasive species distribution modeling | |
Bicknell et al. | Designing protected area networks that translate international conservation commitments into national action | |
Schulp et al. | A quantitative assessment of policy options for no net loss of biodiversity and ecosystem services in the European Union | |
Bertrand et al. | Sensitivity analysis of environmental changes associated with riverscape evolutions following sediment reintroduction: geomatic approach on the Drôme River network, France | |
Beverly et al. | A simple metric of landscape fire exposure | |
Alam et al. | Understanding human–water feedbacks of interventions in agricultural systems with agent based models: a review | |
Stoebner et al. | Geographic determinants of rural land covers and the agricultural margin in the Central United States | |
Yang et al. | Reforestation improves vegetation coverage and biomass, but not spatial structure, on semi-arid mine dumps | |
Rathore et al. | Efficacy of species distribution models (SDMs) for ecological realms to ascertain biological conservation and practices | |
Kornse et al. | Assessing restoration potential for beaver (Castor canadensis) in the semiarid foothills of the Southern Rockies, USA | |
Groeneveld | Species-specific spatial characteristics in reserve site selection | |
Brinley Buckley et al. | Assessing the hydroregime of an archetypal riverine wet meadow in the central Great Plains using time‐lapse imagery | |
Mango | Modeling the effect of land use and climate change scenarios on the water flux of the Upper Mara River flow, Kenya | |
Smith et al. | Cost-effective targeting for reducing soil erosion in a large agricultural watershed | |
Mendoza et al. | Integrating multi-criteria analysis and GIS for land condition assessment: Part 2—Allocation of military training areas |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181225 |