CN107909062A - A kind of Algorithm for Discretization of Remote Sensing Image Attribute and system based on comentropy - Google Patents

A kind of Algorithm for Discretization of Remote Sensing Image Attribute and system based on comentropy Download PDF

Info

Publication number
CN107909062A
CN107909062A CN201711311000.3A CN201711311000A CN107909062A CN 107909062 A CN107909062 A CN 107909062A CN 201711311000 A CN201711311000 A CN 201711311000A CN 107909062 A CN107909062 A CN 107909062A
Authority
CN
China
Prior art keywords
discretization
comentropy
section
image information
attribute
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.)
Granted
Application number
CN201711311000.3A
Other languages
Chinese (zh)
Other versions
CN107909062B (en
Inventor
黄梦醒
陈琼
张雨
冯文龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hainan University
Original Assignee
Hainan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hainan University filed Critical Hainan University
Priority to CN201711311000.3A priority Critical patent/CN107909062B/en
Publication of CN107909062A publication Critical patent/CN107909062A/en
Application granted granted Critical
Publication of CN107909062B publication Critical patent/CN107909062B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/243Classification techniques relating to the number of classes
    • G06F18/2433Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection

Abstract

The present invention relates to a kind of Algorithm for Discretization of Remote Sensing Image Attribute and system based on comentropy, the method comprising the steps of:Establish the image information decision table based on pixel value and generic;Section based on comentropy separates model computation interval comentropy, and interal separation is carried out based on comentropy;Discretization is carried out to the image information decision table, and calculates the Indiscernible relation diversity factor of the forward and backward image information table of discretization respectively;The result of discretization is assessed, selects optimal discretization scheme.The method of the present invention and system can reduce discretization section number, avoid the redundancy of block information, alleviate discretization data distortion problem.

Description

A kind of Algorithm for Discretization of Remote Sensing Image Attribute and system based on comentropy
Technical field
The present invention relates to characteristics of remote sensing image extractive technique field, more particularly to a kind of remote sensing image based on comentropy is special Levy discretization method and system.
Background technology
Environmental monitoring is the important component of marine informatization construction, is the environmental protection and management of ocean, and resource is opened Send out and utilize and provide strong scientific basis.Seasat effectively compensate for the deficiency of Conventional marine observation method, based on more Kind of the remote sensor continuously observation to ocean, makes the mankind greatly deepen the understanding to ocean, in preventing and reducing natural disasters for Oceanic disasters, Development of resources, ocean right-safeguarding, the numerous areas such as marine ecology and environmental protection play an important role.Remote sensing that will be to collecting Image data is analyzed and processed it is necessary to first carry out feature extraction to it, these image features are often continuous data.However, work as Preceding most knowledge extraction (Knowledge Extraction) and data mining (Data Mining) algorithm can only handle from Type data are dissipated, although some algorithms can also handle continuous data, poor-performing.Therefore, it is necessary to be converted into continuous feature Discrete features reduce time and the space expense of algorithm, improve cluster energy of the system to sample to adapt to these intellectualized algorithms Power, strengthening system anti-noise ability, and the study precision of algorithm is improved, expand application range.
A kind of data reducti techniques of the discretization (Discretization) as data preprocessing phase, were subject in recent years Extensive concern and research, and achieve plentiful and substantial achievement in research.The essence of discretization is exactly to determine that selection is more briefly Few cut-point and definite cut-point position.Whether the method for discretization is very much, can be it comprising classification information according to data Be divided into and have supervision discretization and unsupervised discretization, there is supervision discretization to need to consider classification information, and unsupervised discretization Then need not.
At this stage, the unsupervised method of comparative maturity is also fewer, in the case of no category information, to be got well Discretization results are relatively difficult, and the result of discretization is also difficult to weigh.And in having supervision discretization algorithm, single argument has Supervise discretization algorithm and do not have the relation of interdependence considered between attribute, simply from unsupervised discretization mode development to orphan On the spot consider some attribute to the unidirectional dependence or attribute of classification and the relation of interdependence of classification.When it is such from When dispersion algorithm is applied to assorting process, the missing of substantial amounts of redundancy breakpoint and necessary breakpoint will cause overall process between attribute Failure.And existing multivariable has supervision discretization algorithm (Multi Variable Supervised Discretization) Although compensate for above-mentioned single argument to a certain extent has the deficiency of supervision discretization algorithm, also fail to reach and rejecting The Indiscernible relation of decision table is kept while redundancy breakpoint.It is when selecting breakpoint, not in view of between attribute and belonging to Property internal breakpoints alternative, this obviously be not suitable for processing with multiple features remote sensing image data.
The content of the invention
It is an object of the invention to improve the deficiency in the presence of the prior art, there is provided a kind of remote sensing shadow based on comentropy As feature discretization method and system.
In order to realize foregoing invention purpose, an embodiment of the present invention provides following technical scheme:
A kind of Algorithm for Discretization of Remote Sensing Image Attribute based on comentropy, comprises the following steps:
List and pixel value and generic of the sample in each wave band are selected in remote sensing image, and establish and be based on pixel value With the image information decision table of generic;
Section based on comentropy separates model computation interval comentropy, and interal separation is carried out based on comentropy;
To the image information decision table carry out discretization, and respectively calculate the forward and backward image information table of discretization can not Resolution relation diversity factor;
The result of discretization is assessed according to the Indiscernible relation diversity factor, selects optimal discretization side Case.
On the other hand, the embodiment of the present invention additionally provides a kind of characteristics of remote sensing image discretized system based on comentropy, Including:
Image information decision table establishes unit, and pixel value of the sample in each wave band is selected in remote sensing image for listing And generic, and establish the image information decision table based on pixel value and generic;
Comentropy computing unit, separates model computation interval comentropy for the section based on comentropy, and is based on information Entropy carries out interal separation;
Discretization unit, for carrying out discretization to the image information decision table, and calculates the forward and backward shadow of discretization respectively As the Indiscernible relation diversity factor of information table;
Assessment unit, for being assessed according to the Indiscernible relation diversity factor the result of discretization, selection is most Excellent discretization scheme.
Compared with prior art, beneficial effects of the present invention:
Firstly, since the present invention is the interal separation modelling technique based on comentropy, the division of discrete segment is built upon Carried out on the basis of the information gain of attribute independent, and existing other methods be substantially foundation in the generic attribute degree of correlation and Carried out on the basis of similarity between intervals, by experimental data below it can be confirmed that the method for the present invention can reduce discretization Section number;Secondly, when being selected to Candidate point, calculated by section entropy, avoid the redundancy of block information Property, improve the calculating accuracy rate that subsequent classification prediction is handled;Finally, due to it is the calculating based on Indiscernible relation, therefore The data distortion problem after discretization is alleviated to a certain extent, and this method is concentrated use in more long, acquisition in data Discretization scheme it is more excellent.
Brief description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, therefore be not construed as pair The restriction of scope, for those of ordinary skill in the art, without creative efforts, can also be according to this A little attached drawings obtain other relevant attached drawings.
Fig. 1 is the flow chart of the Algorithm for Discretization of Remote Sensing Image Attribute of the invention based on comentropy.
Fig. 2 divides schematic diagram for approximate set.
Fig. 3 is the comprehensive schematic diagram of approximate set.
Fig. 4 is the performance comparison figure of several algorithms.
Fig. 5 is the influence schematic diagram of Indiscernible relation diversity factor.
Fig. 6 analyzes schematic diagram for data consistency.
Fig. 7 is the composition frame chart of the characteristics of remote sensing image discretized system of the invention based on comentropy.
Embodiment
Below in conjunction with attached drawing in the embodiment of the present invention, the technical solution in the embodiment of the present invention is carried out clear, complete Ground describes, it is clear that described embodiment is only part of the embodiment of the present invention, instead of all the embodiments.Usually exist The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and designed with a variety of configurations herein.Cause This, the detailed description of the embodiment of the present invention to providing in the accompanying drawings is not intended to limit claimed invention below Scope, but it is merely representative of the selected embodiment of the present invention.Based on the embodiment of the present invention, those skilled in the art are not doing Go out all other embodiments obtained on the premise of creative work, belong to the scope of protection of the invention.
Referring to Fig. 1, the Algorithm for Discretization of Remote Sensing Image Attribute based on comentropy provided in the present embodiment, including with Lower step:
S101, by carrying out radiation calibration, atmospheric correction, sample selection processing to remote sensing image, is listed in remote sensing image Pixel value and generic of the sample in each wave band are selected, and establishes the image information based on pixel value and generic and determines Plan table.
In the present embodiment, specifically, this step can be realized in the following ways:
Training dataset of the significant region of several block features as discretization is selected from remote sensing image, is denoted as TDS;Together When, wave band attribute is isolated in TDS, a newly-increased row are used as category attribute, and respectively every piece of region is labeled.It is right herein Remote sensing images carry out information extraction, isolate the method for the various information such as geographic coordinate information, wave band attribute can use it is existing Technology realizes that this discrete method does not improve this, therefore to simplify description, this is not stated carefully herein.
Multiple samples are chosen by image sampling, and obtain the pixel value of each selected sample, and at each Corresponding Land cover types are marked in sample, all samples form sample set.For example, according to area marking in region Sample assignment in category attribute row;Then, the wave band attribute of each sample is selected in the much information obtained after TDS separation With category attribute (other information give up should not), merge wave band attribute and category attribute, to build decision table matrix D S;Most Afterwards, the value of category attribute column is ranked up sample, sample set of the generation with pixel value and land type label.
According to the sample set with pixel value and land type label, establish using wave band as conditional attribute, with type of ground objects For the decision table matrix of decision attribute.Matrix pattern is as follows:
Each of which row represents 1 sample item, sample set U={ u1,u2,...,un, i.e., sample is concentrated with n sample, with Wave band is the attribute column C={ c of conditional attribute1,c2,...,cm, represent pixel value of the sample in m wave band, last row is Attribute column D using type of ground objects as decision attribute, identifies the classification information of sample.Each sample item is by sample sequence number, wave band Attribute and category attribute composition.The value range of wave band attribute is 0≤cij≤ 1, wherein cijIt is i-th of sample in j-th of wave band Pixel value;Category attribute value is represented with nature numerical value, and value range is determined according to the classification number of definition, such as:Assuming that definition Classification number be 5, then value range is D={ 1,2,3,4,5 }.
S102, the section based on comentropy separate model computation interval comentropy, and interal separation is carried out based on comentropy.
In the present embodiment, specifically, this step can be accomplished by the following way:
A) according to the status information in current all sections, (status information refers to maximum, minimum value in the entropy table of section And entropy) establish or improve section entropy table;
In order to be quickly found the subinterval of entropy maximum, it is necessary to establish a table to record the entropy in current all sections Value, i.e. section entropy table (IET).For IET altogether comprising 3 row, the 1st row preserve the minimum value in section, and the 2nd row preserve the maximum in section, 3rd row are corresponding entropy obtained by calculation, as shown in table 1.
Table 1 section entropy table (IET) structure
Every a line in table corresponds to a subinterval, originally, a line is comprised only in IET, i.e., whole connection attribute section, with The operation of algorithm, starts to divide (b below) step), by increasing new a line and renewal when predivision section position The minimum and maximum value in two sections, realizes preservations of the IET to all sections.
B) current all corresponding comentropies in section are calculated, judge whether the comentropy in each section is more than given threshold, If it is the section is continued to split, until the entropy in the section is less than given threshold, updates section entropy table.
Sample size according to contained by section, calculates the comentropy in section in above-mentioned IET tables successively, and the calculating of comentropy is public Formula is as follows:
Wherein, S is object set;K is classification number;CiExpression classification in object set S is the example number of i;P (Ci, S) and represent the probability that classification i occurs in object set S, i.e. and the example number of class i accounts for the percentage of whole object set S Than.
Try to achieve in IET tables after the comentropy in all sections, arranged from small to large according to entropy.Every time since last Search entropy it is maximum can by stages, can by stages refer to including at least two breakpoints (upper bound in section and lower bound are unequal) And entropy is more than the section of given threshold value.
A), b) process can be merged understanding, that is, that completes section entropy table establishes process, if the entropy of current interval More than given threshold, then current interval is split, obtain two sections, and update the maximum in two sections, minimum Value, threshold value, then the entropy in each section is calculated, and progress interal separation is determined whether according to entropy, so constantly divided Cut, until the entropy in all sections is respectively less than or equal to given threshold.
C) the segmentation section of each wave band is represented with Candidate point successively, and preserved with a matrix type.
Section can be represented using the lower bound in section as break value, the Candidate point of all wave bands is protected with a matrix type Deposit, as shown in table 2.
The Candidate point matrix of 2 wave band of table
Wherein, BandjRepresent j-th of wave band, breakpointijRepresent i-th of Candidate point value in j-th of wave band, totalkRepresent the Candidate point number that k-th of wave band includes.
S103, discretization is carried out to the raw video information decision table obtained in step S101, and calculates discretization respectively The Indiscernible relation diversity factor of forward and backward image information table.
Specifically, in this embodiment, this step can be accomplished by the following way:
1) discretization is carried out to raw video information decision table using the Candidate point of each wave band successively.
For example, for wave band j, if the pixel value bandvalue of i-th of sampleiK-th of section of this wave band is fallen into, That is, Lk≤bandvaluei≤Uk, then, and bandvaluei=Lk=breakpointkj.Wherein, breakpointkjIt is in wave band j K-th of Candidate point being sorted in ascending order, LkAnd UkIt is lower bound and the upper bound in the segmentation section that this Candidate point represents respectively.Will All samples, by aforesaid operations, complete the sliding-model control to raw information decision table in each wave band.
2) lower and upper approximations on decision attribute of the forward and backward image information decision table of discretization are solved respectively.
Given decision table S=(U, R, V, f), wherein, U is limited object set, i.e. domain, and R is attribute set, bag Collection containing conditional attribute C and decision kind set D.Binary crelation can not be differentiated:For each attribute setDefinition can not Binary crelation IND (A) is differentiated, i.e.,:
Equivalence class in domain U on attribute set A:
To each subsetAnd the equivalence class in domain U on attribute set A, the lower and upper approximations of X respectively by It is defined as follows:
Here A=C, X ∈ U are taken | IND (d), d are the decision attribute of decision kind set D, before being calculated discretization respectively, The lower aprons collection C of each decision attribute values is corresponded in the decision attribute d of image information decision table afterwards-(dX), C-(dX) ' and it is upper Approximate set C-(dX), C-(dX) ', C-(dX) be discretization before lower aprons collection, C-(dX) ' be discretization after lower aprons collection, C- (dX) be discretization before upper approximate set, C-(dX) ' be discretization after upper approximate set.
3) the Indiscernible relation diversity factor of the forward and backward image information decision table of discretization is calculated.
, can be in the hope of the difference number N of the forward and backward lower and upper approximations of discretization according to above-mentioned definitiond=Nl+Nu, wherein, Nl For the difference number of lower aprons collection, NuFor the difference number of upper approximate set, Nd, NlAnd Nu0 is initialized as respectively.Discretization is forward and backward The difference number of lower and upper approximations is the Indiscernible relation diversity factor of the forward and backward image information decision table of discretization.NdCalculating Step is as follows:
Step 1:Element d is chosen from the decision attribute d of raw video information decision tablei(i=1,2 ..., n), its In, n is the number of different values of the decision attribute d on domain U, i.e. classification number;
Step 2:Calculate diLower aprons collection C-(di) and upper approximate set C-(di), if also having element not in decision attribute d Calculated, then return to Step 1 and continue to execute, otherwise, into next step Step 3;
Step 3:With the IET ultimately generated to raw video information decision table discretization, new image information decision-making is obtained Table SE=(U, R, VE,fE);
Step 4:Element d is chosen from the decision attribute d of new image information decision tablei(i=1,2 ..., n), calculate diLower aprons collection C-(di) ' and upper approximate set C- (di)';
Step 5:The forward and backward d of discretization is judged respectivelyiLower and upper approximations it is whether equal, if C-(di)'≠C-(di), Then Nl=Nl+ 1, if C- (di)'≠C-(di), then Nu=Nu+1;If also having element not calculated in decision attribute d, return Step 4 is returned to continue to execute, otherwise, Nd=Nl+Nu, terminate.
S104, assesses the result of discretization according to the Indiscernible relation diversity factor, selects optimal discrete Change scheme.
Wherein, this step assessment is carried out to the result of discretization can be in the following way:
The deliberated index of the result of discretization is essentially from nicety of grading, the uniformity of section number and sample, and classifies Precision and the uniformity of sample depend on the size that Indiscernible relation changes degree again, section number then by the threshold value that sets and Iterations controls.Based on the Indiscernible relation diversity factor being calculated in step 103, different threshold values and iterations are set The discretization results fixed score, and formula is as follows:
Wherein, NdFor the Indiscernible relation diversity factor of the forward and backward image information decision table of discretization, | d | it is the class of sample Other number, each classification correspond to a pair of of lower and upper approximations;It is for the Candidate point sum of all wave bands, i.e., discrete The section number obtained after change,For the breakpoint sum of raw video information decision table.
Select optimal discretization scheme to recommend targeted customer described in this step, can by the different threshold values of setting and Discretization assessment result under iterations is ranked up, and selects the highest discretization scheme of score value to recommend targeted customer.
The step of Algorithm for Discretization of Remote Sensing Image Attribute of the present invention based on comentropy is described above, discretization of the present invention Method has the redundancy that can be reduced discretization section number, avoid block information, alleviates discretization data distortion problem.Under Face compared with the method for the present invention, will be verified with several discretization methods of current mainstream in terms of theoretical and test data two The superiority of the method for the present invention.
1) generic information calculates
There is the discretization algorithm of supervision due to the calculating based on generic information, have and make full use of class label and target The superiority of attribute information, can be more easily found suitable breakpoint location than unsupervised discretization algorithm.Have supervision from Dispersion method has very much, wherein more typical generic information calculation formula mainly has following three kinds.
The a generic degrees of correlation
Wherein, CAIR is class-attribute relevant redundancy degree, and I (C, D | F) represents class-attribute of generic attribute C and discrete division D Shared mutual information, H (C, D | F) represent the combination entropy between generic attribute C and discrete division D.
B similarity between intervals
Wherein, χ2Represent the departure degree between observed value and theoretical value, AijRepresent the sample for belonging to class j in the i-th section This number, EijRepresent the expecterd frequency to the class j in the i-th section.
C information gain degrees
Wherein, Gain (S) is the information gain amount before and after interval division, and Ent (S) represents the information content of section S, S1And S2 Represent respectively by two sections of breakpoint division generation;| S | represent the example number that section S is included.
The traditional generic information computational methods of above three have more complete theoretical foundation.Calculation based on the generic degree of correlation Method mainly has CADD, CAIM and CACC etc..CADD algorithms need user to specify interval number in initialization, and CAIM need not be pre- First set interval number, and make certain improvements in the calculating of degree of correlation standard.Although CAIM compensate for some of CADD Deficiency, but since warp is frequently with interval number approximate target attribute value number, and for a section, the discrete differentiation provided Formula only considers that the section includes the generic attribute of most samples, have ignored the distribution situation of other generic attributes, eventually produces excessive Section and cause over training, cause overfitting.CACC is similar with CAIM, and simply discretization evaluation criterion is repaiied Change, determine whether section needs to split by the coefficient of association cacc of computation interval, and replace in CAIR standards using log (n) N, so as to compensate for the deficiency of CAIM.Although CACC can accelerate the process of discretization, overfitting is prevented, to single category Performance obtains the preferable breakpoint of quality, but since CACC accepted standards are to maximize the generic degree of correlation, only to single attribute Realize that Interval Discrete result is optimal, lack the description to data Global Information, also do not account for the uniformity of data before and after the processing, It will necessarily cause the loss of initial data important information amount.
ChiMerge, Chi2 and Extended Chi2 scheduling algorithms employ the generic information based on similarity between intervals and calculate Method carries out differentiation merging to adjacent interval.These methods apply the Pearson came statistic in statistics to differentiate current breakpoint Whether should be removed, i.e., adjacent with the breakpoint two interval whether the merging.Can by the Pearson came theorem in mathematical statistics Know, statistic χ2Progressive distribution be the free degree be k-1 χ2Distribution, as given level of signifiance α, it may be determined that corresponding critical Value χα.ChiMerge is using fixed level of signifiance α, and Chi2 is then constantly to decline level of signifiance α, so that critical value χαConstantly increase, passes through mergingMaximum Candidate point, improves the utilization rate of computational efficiency and information.But It is that Chi2 uses total classification number that fixed free degree v=k-1, wherein k are systems.In fact, the selection of the free degree should The classification number that the foundation section is associated, total classification number without that should be system, therefore, ExtendedChi2 algorithms are in Chi2 On improved, useInstead of in Chi2 algorithmsCalculated.However, Extended Chi2 are calculated Used in methodAlthough the discrimination standard as breakpoint importance is than usingIt is rational more, but be a lack of Corresponding theoretical foundation, at the same it is computationally also inaccurate.
The most common discretization method calculated using information gain generic information, is based on most short description length principle (MDLP) method.It is after connection attribute according to value size sequence, and the border between different target class is set to candidate's division points, Then being found out in candidate's division points makes the division border T of information gain degree maximum as two points of divergent boundaries, successively iteration time Definite potential division points (Potential cut point) are gone through until meeting minimum description length principle.The algorithm is due to using The measurement standard of information gain, for more above-mentioned algorithm, can largely ensure the uniformity of section internal specimen generic, The situation that suitable Target Attribute values are evenly distributed.But in remote sensing image processing, the sample taken often be subject to noise and The influence of other impurity, Target Attribute values be distributed after sequence it is still more scattered, to each attribute need carry out N-1 times Information gain degree calculates, and wherein N is the number of samples with different attribute value, therefore, when data scale is very big, can be caused non- Often big time overhead.
2) approximately equivalent model up and down
Information system is four-tuple (U, Q, V, f).Wherein, U is object set, also referred to as domain;Q is attribute set, bag Collection containing conditional attribute C and decision kind set D;V is the codomain of attribute;F is mapping function, represents that the attribute of object is mapped to attribute Some value of codomain.Knowledge is referred to as according to the ability that the characteristic difference of things is classified.It can not be differentiated each other in domain Object composition set be form knowledge particle (granule).Knowledge has granularity, and the smaller explanation of granularity can accurate table The concept reached is more, and granularity uses Indiscernible relation, i.e. equivalence class is described, and is the least unit for representing knowledge.
In assorting process, the individual subject being not much different is attributed to same class, and the relation formed between them is exactly not Distinguishable relation, then be defined as below:
If IND (P)={ (x, y) ∈ U × U:F (x, a)=f (y, a), a ∈ P }, be on any one attribute set P can not Resolution relation.Judge whether an object a belongs to set X according to the Indiscernible relation cluster in knowledge base, three kinds of situations can be divided: A necessarily belongs to set X;A is necessarily not belonging to set X;A, which may belong to, may also be not belonging to set X.So, set X is on can not Differentiate the lower aprons (Lower approximation) of relations IIt is exactly according to I by those Judgement belongs to the maximum set that the object of X is formed certainly, also referred to as the positive domain (positive region) of X, is denoted as POS (X);Upper approximate (Upper approximations) of the set X on Indiscernible relation IIt is by the union of all equivalence class I (X) for intersecting non-NULL with X, is that those may belong to In the minimal set that the object of X forms.If approximation set equality up and down, X is a precise set, is otherwise one coarse Set, wherein lower aprons are known as the positive region of the concept, up and down it is approximate it is poor be referred to as border, it is upper it is approximate beyond region be known as Negative region (Negative region), is denoted as NEG (x), and the model thus established is known as approximately equivalent model up and down.
It is assumed that information system S=(U, Q, V, f), there is two attribute Q1And Q2, attribute Q1There are 5 values, attribute Q2There are 6 values. Given approximation set X, uses attribute Q1And Q2Divided, situation is as shown in Figure 2.
Such as Fig. 2, attribute Q1Domain U 5 row, attribute Q have been divided into2Domain U 6 rows, Q have been divided into1×Q2Constitute one A two dimensional character, is divided into 30 lattices by domain U, that is, generates 30 knowledge particles altogether.Circle represents approximate set X, The region of lattice composition disjoint with circle is the negative region on X, its expression formula is:
NEG (X)=U-I-(X) (11)
Lower aprons I-(X) and upper approximation I-(X) extracted respectively from the knowledge particles of division, as shown in Figure 3.Can To see, the region with five-pointed star mark is the positive region on X, and the region with No. # mark is the frontier district of approximate set X Domain, its expression formula are:
BND (X)=I-(X)-I-(X) (12)
The concept such as lower aprons, upper approximate and borderline region carves the approximation properties for having changed an ambiguous set in border.Coarse journey Degree is calculated as follows:
In formula | I- (X) | and | I- (X) | the radix or gesture of lower aprons collection and upper approximate set are represented respectively, that is, are wrapped in gathering The element number contained.Obvious 0≤α1(X)≤1, if α1(X)=1, then set X is claimed relative to Indiscernible relation I to be clear , if α1(X) < 1, then title X is coarse relative to I.α1(X) it is regarded as approaching set X's under Indiscernible relation I Precision.
Ambiguous set can not judge whether some elements belong to this without clearly border according to existing knowledge in domain Set.In rough set, uncertain concept is for element is under the jurisdiction of the degree of set.By using pass can not be differentiated System, it is as follows to the coarse membership function of set X, formula under equivalence relation I can to define element x:
Indiscernible relation changes the importance for the conditional attribute set C that degree is actually reflected under decision attribute D Change, i.e. decision attribute set D is formulated as follows the degree of dependence of conditional attribute set C:
Wherein, POSc(D) it is positive regions of the attribute set C in U/IND (D).As can be seen from the above equation, degree of dependence is just It is ratio of the lower aprons set element with whole domain radix of the corresponding equivalence class of all categories value in decision attribute D. Since each approximation set X correspond to a lower aprons collection I of given equivalence relation I-(X) and one is gone up approximate set I- (X), the change of degree of dependence means that approximate set is also changed, and can replace asking by calculating the variable number of approximate set The ratio of positive domain and domain radix carrys out the intensity of variation of reaction information system.In fact, the change of information system is exactly by knowledge Caused by the equivalence relation in storehouse changes, i.e., directly portray the I of set X-(X) changed with I- (X).And γC(D) Simply reflect the ratio of element number, identical γC(D) different I may be correspond to-(X) and I-(X) change number, because This, directly uses I-(X) and I-(X) change number reflects that the intensity of variation of information system is relatively reasonable.According to above-mentioned analysis, Construction solves the function f of Indiscernible relation diversity factor:
I. the input of function is two decision tables before and after discretization, is represented respectively with DT and NDT;
Ii. all approximate sets, and being compared up and down of DT and NDT on decision attribute are obtained respectively;
Iii. the difference number of lower aprons collection and the difference number of upper approximate set that the output of function is DT and NDT.
Diversity factor function f reflects the intensity of variation of knowledge system, and lower aprons set and the upper approximate number gathered are by sample The value number of this class number, i.e. decision attribute determines.
The principle of the present invention is that the change for dividing front and rear information gain according to section decides whether to more than threshold value Section continues to divide, and entropy, the linear module as information uncertainty can reflect the steady of section generic well Determine degree.The knowledge base of information system is made of internal equivalence relation at the same time, so, closed using the equivalence of decision table System's change is to reflect the otherness between information system, and utilizing the relation diversity factor calculated, simultaneously combining information gain-splitting is calculated Method, the final discretization scheme for obtaining image feature.The method of the present invention introduces in the section fission process larger to entropy The equivalence model of rough set is compared calculating to the decision table before and after discretization, by the difference number of obtained approximate set up and down Mesh knows the change degree of Indiscernible relation, so as to obtaining optimal discretization interval number.
When being verified by experiment to the method for the present invention, it is South Sea inshore region to test the test data set used A width GF-2 satellite images, which includes 4 wave bands, the atural object on image is divided into settlement place in experiment, forest land is naked Ground, arable land and five major class of water body.Local sampling is carried out to imagery zone, selects the significant imagery zone of several block features as discrete The training dataset of change, the sample set with pixel value and class label finally obtained include 8379 samples, 4 wave bands altogether Initial breakpoint number be 6328,6596,6720 and 6611 respectively, sum be 26255.
The visualization of this experiment, programs, emulation, and test and numerical computations processing are in MATLAB (R2016a versions) environment Lower realization;The radiation calibration of image, air calibration, and comparison of classification is carried out in ENVI5.3 rings to the result before and after discretization Completed under border.First by carrying out radiation calibration to remote sensing image, atmospheric correction, sample selection processing, lists in image and is chosen Sample and establishes the information model of image, i.e. image information decision table in the pixel value and generic of each wave band.Image information Decision table is altogether comprising 5 row, and the 1st is classified as the numbering of sample, and the 2nd to the 5th is classified as wave band attribute, i.e. conditional attribute, and the 6th is classified as classification Attribute, i.e. decision attribute.This input data of image information decision table as discretization algorithm performs, every time with one-dimension array Form is exported as a result, last, and all results are preserved in the form of two-dimensional matrix, matrix totally 5 row, the 1st to 4 row respectively from The breakpoint number of each wave band after dispersion, the 5th is classified as Indiscernible relation difference number.
In order to verify the validity for proposing algorithm, experiment is accurate by section number, Indiscernible relation diversity factor, data Rate and run time this several respect analyze implementing result.And there is supervision with the former algorithm based on entropy and current main-stream Discretization algorithm 1R, ChiMerge are compared.
Figure 4, it is seen that the run time of four kinds of algorithms is similar.Since classification number is set as 5 major classes, so shadow As information decision table on the lower and upper approximations of category attribute sum be 10.Although former algorithm and side of the present invention based on entropy Method obtains almost consistent breakpoint number, but since Indiscernible relation difference number is 8, the number obtained than the method for the present invention 6 is high, and therefore, the accuracy 78.5% classified is similar 10 percentage points low compared with the method for the present invention 90.8%. The Indiscernible relation difference number of ChiMerge algorithms is consistent with the method for the present invention, but is obtained more than the method for the present invention Breakpoint number, and the accuracy classified also is slightly less than the method for the present invention.1R algorithms are relatively other due to the process in splitting section Algorithm is simple, and the breakpoint number of minimum is achieved in 4 wave bands, but the difference number of Indiscernible relation is maximum, reaches 10, Mean that the distortion level of decision table maximizes, it is clear that the accuracy rate of classification has also reached minimum 53.6% in four algorithms. Performance Evaluation is carried out to four kinds of algorithms using DME formula, the assessed value for trying to achieve four kinds of algorithms is respectively:0.2896,0.1441, 0.1828 and 0.Where it can be seen that the overall performance of the method for the present invention is more preferably.
By the analysis to Fig. 4, while accuracy rate generation of the Indiscernible relation diversity factor to classification of drawing a conclusion is important Influence, i.e. classification accuracy under different relation diversity factoies is different, and is reduced with the increase of difference number.
From figure 5 it can be seen that by varying the discretization flow iterations of the method for the present invention threshold different with setting Value, can obtain different relation diversity factoies, the value in [0,2,4,6,8,10] scope.Number of breakpoints is with difference number Reduce and increase, when heteromerism mesh of being on duty is 0, the number of breakpoints and original breakpoint sum that discretization obtains are basically identical, occur Over-fitting.Calculated with crossing, obtain the DME values under this 6 diversity factoies, be respectively 0.0251,0.0631,0.1253, 0.2896,0.1629 and 0.As can be seen that when diversity factor is 6, the DME values of maximum are obtained.Draw a conclusion, although diversity factor reduces The compatibility of system after discretization can be kept, but also brings breakpoint enormous amount, over-fitting occurs, obtained DME values It is relatively low.Therefore, the method for the present invention can therefrom obtain optimum discretization scheme under balance number of breakpoints and relation diversity factor.
Mainly the uniformity of data after each wave band discretization is analyzed in figure 6.Here uniformity refers to The accuracy rate of sample, i.e., after discretization, the sample in same discrete segment has identical category attribute.As can be seen that As difference number constantly increases, the accuracy rate of sample gradually reduces in each wave band.Wave band 1 difference number in [0,2,4, 6,8,10] during value, the accuracy rate of sample is respectively [0.985 0.965 0.913 0.895 0.767 0.665];Wave band 2 Sample accuracy rate is [0.963 0.957 0.935 0.925 0.793 0.632];The sample accuracy rate of wave band 3 is [0.956 0.924 0.906 0.894 0.802 0.763];The sample accuracy rate of wave band 4 is [0.972 0.943 0.883 0.856 0.757 0.602].This illustrates that the change of decision table Indiscernible relation can produce large effect to data consistency.And Experiment simulation and comprehensive analysis more than demonstrate the validity of the method for the present invention, and more other algorithms can obtain less area Between the classification accuracy of number and higher, while also alleviate the data distortion problem of system to a certain extent, improve discrete Change scheme recommends efficiency and the scalability of algorithm.
Referring to Fig. 7, the present invention gives the characteristics of remote sensing image discretized system based on comentropy at the same time, in the system Including:
Image information decision table establishes unit, and pixel value of the sample in each wave band is selected in remote sensing image for listing And generic, and establish the image information decision table based on pixel value and generic;
Comentropy computing unit, separates model computation interval comentropy for the section based on comentropy, and is based on information Entropy carries out interal separation;
Discretization unit, for carrying out discretization to the image information decision table, and calculates the forward and backward shadow of discretization respectively As the Indiscernible relation diversity factor of information table;
Assessment unit, for being assessed according to the Indiscernible relation diversity factor the result of discretization, selection is most Excellent discretization scheme.
Wherein, in the present embodiment, image information decision table is established unit and is specifically included:Sample choose module, for from Multiple samples are chosen in remote sensing image;Sample set establishes module, for obtaining each selected sample by image sampling Pixel value, and corresponding Land cover types are marked in each sample, all samples form sample set;Decision table square Battle array establishes module, for according to the sample set with pixel value and land type label, establishing using wave band as conditional attribute, with ground Species type is the decision table matrix of decision attribute.
Comentropy computing unit specifically includes:Section entropy table establishes module, for the state letter according to current all sections Breath establishes section entropy table;Comentropy computing module, for calculating current all corresponding comentropies in section, is obtained based on comentropy The position in next section to be split;Preserving module, for successively representing in the segmentation section of each wave band simultaneously with Candidate point Preserve with a matrix type.
Discretization unit specifically includes:Descretization module, for utilizing the Candidate point of each wave band successively to the shadow As information decision table carries out discretization;Approximate set asks for module, for image information table before solving discretization respectively on decision-making Upper approximate set, the lower aprons collection of attribute, and upper approximate set, lower aprons collection of the image information table on decision attribute after discretization; Diversity factor computing module, for calculating the Indiscernible relation diversity factor of the forward and backward image information table of discretization respectively.
Only the characteristics of remote sensing image discretized system based on comentropy is briefly described herein, is described in more detail It may refer to the associated description in the above-mentioned Algorithm for Discretization of Remote Sensing Image Attribute based on comentropy.
The above description is merely a specific embodiment, but protection scope of the present invention is not limited thereto, any Those familiar with the art the invention discloses technical scope in, change or replacement can be readily occurred in, should all be contained Cover within protection scope of the present invention.

Claims (10)

1. a kind of Algorithm for Discretization of Remote Sensing Image Attribute based on comentropy, it is characterised in that comprise the following steps:
List and pixel value and generic of the sample in each wave band are selected in remote sensing image, and establish and be based on pixel value and institute Belong to the image information decision table of classification;
Section based on comentropy separates model computation interval comentropy, and interal separation is carried out based on comentropy;
Discretization is carried out to the image information decision table, and calculates can not differentiating for the forward and backward image information table of discretization respectively Relation diversity factor;
The result of discretization is assessed according to the Indiscernible relation diversity factor, selects optimal discretization scheme.
2. according to the method described in claim 1, it is characterized in that, described list is selected sample in each ripple in remote sensing image The pixel value and generic of section, and the image information decision table based on pixel value and generic is established, including step:
Multiple samples are chosen from remote sensing image;
The pixel value of each selected sample is obtained by image sampling, and corresponding soil is marked in each sample Cover type, all samples form sample set;
According to the sample set with pixel value and land type label, establish using wave band as conditional attribute, using type of ground objects as certainly The decision table matrix of plan attribute.
3. according to the method described in claim 1, it is characterized in that, the section based on comentropy separates model computation interval Comentropy, interal separation, including step are carried out based on comentropy:
Section entropy table is established according to the status information in current all sections;
Current all corresponding comentropies in section are calculated, judge whether the comentropy in each section is more than given threshold, if Then the section is continued to split, until the entropy in the section is less than given threshold, update section entropy table;
The segmentation section of each wave band is represented with Candidate point successively and is preserved with a matrix type.
4. according to the method described in claim 3, it is characterized in that, it is described to the image information decision table carry out discretization, And the Indiscernible relation diversity factor of the forward and backward image information table of discretization is calculated respectively, including step:
Discretization is carried out to the image information decision table using the Candidate point of each wave band successively;
Upper approximate set, lower aprons collection of the image information table on decision attribute before solution discretization respectively, and image after discretization Upper approximate set, lower aprons collection of the information table on decision attribute;
The Indiscernible relation diversity factor of the forward and backward image information table of discretization is calculated respectively.
5. according to the method described in claim 1, it is characterized in that, it is described according to the Indiscernible relation diversity factor to discrete The result of change is assessed, and selects optimal discretization scheme, including step:
Scored using equation below the discrete results under different set threshold value and iterations:Wherein, NdFor the forward and backward image information table of discretization Indiscernible relation diversity factor, | d | be the classification number of sample, each classification correspond to a pair of of lower and upper approximations;It is total for the Candidate point of all wave bands,Breakpoint for raw video information table is total Number, m is wave band number, and i, j, m is integer;
Choose the highest discretization scheme of score value.
A kind of 6. characteristics of remote sensing image discretized system based on comentropy, it is characterised in that including:
Image information decision table establishes unit, and pixel value and institute of the sample in each wave band are selected in remote sensing image for listing Belong to classification, and establish the image information decision table based on pixel value and generic;
Comentropy computing unit, for based on comentropy section separate model computation interval comentropy, and based on comentropy into Row interal separation;
Discretization unit, for carrying out discretization to the image information decision table, and calculates the forward and backward image letter of discretization respectively Cease the Indiscernible relation diversity factor of table;
Assessment unit, for being assessed according to the Indiscernible relation diversity factor the result of discretization, selects optimal Discretization scheme.
7. system according to claim 6, it is characterised in that image information decision table, which establishes unit, to be included:
Sample chooses module, for choosing multiple samples from remote sensing image;
Sample set establishes module, for obtaining the pixel value of each selected sample by image sampling, and at each Corresponding Land cover types are marked in sample, all samples form sample set;
Decision table matrix establishes module, for according to the sample set with pixel value and land type label, establish using wave band as Conditional attribute, the decision table matrix using type of ground objects as decision attribute.
8. system according to claim 6, it is characterised in that comentropy computing unit includes:
Section entropy table establishes module, for establishing section entropy table according to the status information in current all sections;
Comentropy computing module, for calculating current all corresponding comentropies in section, judge each section comentropy whether More than given threshold, if it is the section is continued to split, until the entropy in the section is less than given threshold, update section Entropy table;
Preserving module, for being represented the segmentation section of each wave band with Candidate point and being preserved with a matrix type successively.
9. system according to claim 8, it is characterised in that discretization unit includes:
Descretization module, discretization is carried out for the Candidate point successively using each wave band to the image information decision table;
Approximate set asks for module, for solving the upper approximate set, lower near of image information table before discretization on decision attribute respectively Like upper approximate set of the image information table after collection, and discretization on decision attribute, lower aprons collection;
Diversity factor computing module, for calculating the Indiscernible relation diversity factor of the forward and backward image information table of discretization respectively.
10. system according to claim 6, it is characterised in that assessment unit is specifically used for:Using equation below to difference Discrete results under given threshold and iterations score, and choose the highest discretization scheme of score value:
Wherein, NdFor the forward and backward image letter of discretization The Indiscernible relation diversity factor of table is ceased, | d | it is the classification number of sample, each classification correspond to a pair of of lower and upper approximations;It is total for the Candidate point of all wave bands,Breakpoint for raw video information table is total Number, m is wave band number, and i, j, m is integer.
CN201711311000.3A 2017-12-11 2017-12-11 Information entropy-based remote sensing image feature discretization method and system Active CN107909062B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711311000.3A CN107909062B (en) 2017-12-11 2017-12-11 Information entropy-based remote sensing image feature discretization method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711311000.3A CN107909062B (en) 2017-12-11 2017-12-11 Information entropy-based remote sensing image feature discretization method and system

Publications (2)

Publication Number Publication Date
CN107909062A true CN107909062A (en) 2018-04-13
CN107909062B CN107909062B (en) 2021-05-07

Family

ID=61865137

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711311000.3A Active CN107909062B (en) 2017-12-11 2017-12-11 Information entropy-based remote sensing image feature discretization method and system

Country Status (1)

Country Link
CN (1) CN107909062B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109076761A (en) * 2018-08-08 2018-12-25 安徽省现代农业装备产业技术研究院有限公司 One kind being based on remote sensing images variable fertilization method
CN112580483A (en) * 2020-12-15 2021-03-30 海南大学 Remote sensing image characteristic discretization method based on rough fuzzy model
CN113743464A (en) * 2021-08-02 2021-12-03 昆明理工大学 Continuous characteristic discretization loss information compensation method and application thereof
CN117075884A (en) * 2023-10-13 2023-11-17 南京飓风引擎信息技术有限公司 Digital processing system and method based on visual script

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1652146A1 (en) * 2003-07-21 2006-05-03 Tubitak-Bilten ( Turkiye Bilimsel Ve Teknik Arastirma Kurumu-Bilgi Teknolojileri Ve Elektronik Arastirma Enstitusu) Implementation of the jpeg2000 compression algorithm in hardware
CN101702172A (en) * 2009-11-10 2010-05-05 大连理工大学 Data discretization method based on category-attribute relation dependency
CN102879677A (en) * 2012-09-24 2013-01-16 西北工业大学 Intelligent fault diagnosis method based on rough Bayesian network classifier
CN102930275A (en) * 2012-10-29 2013-02-13 福州大学 Remote sensing image feature selection method based on Cramer's V index
CN104573050A (en) * 2015-01-20 2015-04-29 安徽科力信息产业有限责任公司 Continuous attribute discretization method based on Canopy clustering and BIRCH hierarchical clustering
CN105139021A (en) * 2015-07-08 2015-12-09 Tcl集团股份有限公司 Method and system for realizing television user rapid classification based on rough set theory
CN106971204A (en) * 2017-03-31 2017-07-21 上海海洋大学 A kind of dimension reduction method of the new high-spectrum remote sensing data based on rough set
TWI599896B (en) * 2016-10-21 2017-09-21 嶺東科技大學 Multiple decision attribute selection and data discretization classification method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1652146A1 (en) * 2003-07-21 2006-05-03 Tubitak-Bilten ( Turkiye Bilimsel Ve Teknik Arastirma Kurumu-Bilgi Teknolojileri Ve Elektronik Arastirma Enstitusu) Implementation of the jpeg2000 compression algorithm in hardware
CN101702172A (en) * 2009-11-10 2010-05-05 大连理工大学 Data discretization method based on category-attribute relation dependency
CN102879677A (en) * 2012-09-24 2013-01-16 西北工业大学 Intelligent fault diagnosis method based on rough Bayesian network classifier
CN102930275A (en) * 2012-10-29 2013-02-13 福州大学 Remote sensing image feature selection method based on Cramer's V index
CN104573050A (en) * 2015-01-20 2015-04-29 安徽科力信息产业有限责任公司 Continuous attribute discretization method based on Canopy clustering and BIRCH hierarchical clustering
CN105139021A (en) * 2015-07-08 2015-12-09 Tcl集团股份有限公司 Method and system for realizing television user rapid classification based on rough set theory
TWI599896B (en) * 2016-10-21 2017-09-21 嶺東科技大學 Multiple decision attribute selection and data discretization classification method
CN106971204A (en) * 2017-03-31 2017-07-21 上海海洋大学 A kind of dimension reduction method of the new high-spectrum remote sensing data based on rough set

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FENGJIANG ET AL: "A novel approach for discretization of continuous attributes in rough set theory", 《KNOWLEDGE-BASED SYSTEMS》 *
STEPHEN D. BAY ET AL: "Multivariate discretization of continuous variables for set mining", 《KDD "00: PROCEEDINGS OF THE SIXTH ACM SIGKDD INTERNATIONAL CONFERENCE ON KNOWLEDGE DISCOVERY AND DATA MINING》 *
梁静等: "基于信息熵的遥感影像特征离散化方法", 《地理空间信息》 *
王举范等: "基于信息熵的粗糙集连续属性多变量离散化算法", 《青岛科技大学学报(自然科学版)》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109076761A (en) * 2018-08-08 2018-12-25 安徽省现代农业装备产业技术研究院有限公司 One kind being based on remote sensing images variable fertilization method
CN109076761B (en) * 2018-08-08 2021-01-08 安徽省现代农业装备产业技术研究院有限公司 Variable fertilization method based on remote sensing image
CN112580483A (en) * 2020-12-15 2021-03-30 海南大学 Remote sensing image characteristic discretization method based on rough fuzzy model
WO2022127075A1 (en) * 2020-12-15 2022-06-23 海南大学 Feature discretization method for remote sensing image on the basis of rough fuzzy model
CN112580483B (en) * 2020-12-15 2023-12-01 海南大学 Remote sensing image feature discretization method based on rough fuzzy model
CN113743464A (en) * 2021-08-02 2021-12-03 昆明理工大学 Continuous characteristic discretization loss information compensation method and application thereof
CN113743464B (en) * 2021-08-02 2023-09-05 昆明理工大学 Continuous characteristic discretization loss information compensation method and application thereof
CN117075884A (en) * 2023-10-13 2023-11-17 南京飓风引擎信息技术有限公司 Digital processing system and method based on visual script
CN117075884B (en) * 2023-10-13 2023-12-15 南京飓风引擎信息技术有限公司 Digital processing system and method based on visual script

Also Published As

Publication number Publication date
CN107909062B (en) 2021-05-07

Similar Documents

Publication Publication Date Title
WO2019149147A1 (en) Method for dividing ecological and geological environment types based on coal resource development
CN107909062A (en) A kind of Algorithm for Discretization of Remote Sensing Image Attribute and system based on comentropy
CN107067405B (en) Remote sensing image segmentation method based on scale optimization
CN102930275B (en) Based on the characteristics of remote sensing image system of selection of Cramer ' s V index
CN103336801A (en) Multi-feature locality sensitive hashing (LSH) indexing combination-based remote sensing image retrieval method
CN108764527B (en) Screening method for soil organic carbon library time-space dynamic prediction optimal environment variables
CN111611960B (en) Large-area ground surface coverage classification method based on multilayer perceptive neural network
CN112087316B (en) Network anomaly root cause positioning method based on anomaly data analysis
CN107392887A (en) A kind of heterogeneous method for detecting change of remote sensing image based on the conversion of homogeneity pixel
Lunga et al. Resflow: A remote sensing imagery data-flow for improved model generalization
CN116129262A (en) Cultivated land suitability evaluation method and system for suitable mechanized transformation
Cheng et al. Multi-view graph convolutional network with spectral component decompose for remote sensing images classification
Zhang et al. Geochemically constrained prospectivity mapping aided by unsupervised cluster analysis
Nong et al. Airborne LiDAR point cloud classification using PointNet++ network with full neighborhood features
Wu et al. Geo-parcel-based geographical thematic mapping using C5. 0 decision tree: A case study of evaluating sugarcane planting suitability
CN117093890A (en) Comprehensive evaluation method for ecological environment of energy resource development area
Nguyen Land cover change detection in northwestern Vietnam using Landsat images and Google Earth Engine
Du et al. An ensemble learning approach for land use/land cover classification of arid regions for climate simulation: A case study of Xinjiang, Northwest China
CN115719453B (en) Deep learning-based remote sensing extraction method for rice planting structure
Roudier et al. Soilscapes of New Zealand: Pedologic diversity as organised along environmental gradients
Maurya et al. Estimation of major agricultural crop with effective yield prediction using data mining
Ji et al. Multicascaded feature fusion-based deep learning network for local climate zone classification based on the So2Sat LCZ42 benchmark dataset
Jiang et al. Remote sensing recognition and classification of forest vegetation based on image feature depth learning
Zhang et al. A hierarchical growth method for extracting 3D phenotypic trait of apple tree branch in edge computing
Li et al. Feature-guided dynamic graph convolutional network for wetland hyperspectral image classification

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
GR01 Patent grant
GR01 Patent grant