CN102663524A - Land utilization class factor analysis method for urban and rural planning - Google Patents

Abstract

The invention discloses a method for quantitative zoning by utilizing land utilization class factors. The method comprises the following steps of: conducting investigation and measurement on multiple land single factors in a target area, and drawing a thematic map according to the result; carrying out vectorization on each land single factor by using a GIS (geographic information system), establishing a unified coordinated system, and drawing a picture for each land single factor; classifying each land single factor according to a preset classification standard; determining a weight value of each land single factor by using an analytic hierarchy process; superposing all the multiple land single factors by using the GIS; carrying out weighted average on each land single factor to obtain a land class factor analysis result; classify the land class factor analysis result by adopting a forced normal distribution method to finish a land class factor analysis chart, and thus carrying out quantitative zoning on the target area. According to the method for quantitative zoning by utilizing the land utilization class factors, provided by the embodiment of the invention, land utilization and planning can be more scientifically and objectively analyzed and researched..

Description

A kind of soil of urban and rural planning utilizes a type factor-analysis approach

Technical field

The present invention relates to environmental science and urban ecological planning field, particularly relate to and utilize the soil class factor to carry out the quantification partitioned method.

Background technology

Land resource is to carry unify material base, especially Human's production and the life of human settlements system of natural ecosystem to receive the very big restriction of land resource, thereby mode, method and the pattern of mankind's activity have been confirmed in the evaluation of land resource.The core of urban and rural planning is the soil utilization, thereby the land use class factorial analysis of town and country ecological planning is very important.

The land use class factor-analysis approach of town and country ecological planning is special a kind of space ecological analysis method of utilizing to the soil in the town and country ecological planning class factor-analysis approach, is directed against the spatial analysis of town and country ecological planning especially.Through analyze, research and estimate urban land and utilize various features key elements such as state and ecological effect as the ecological resources state; To have the future city ecological planning selects several leading factor to compose the weighted values stack in the object of influence; Quantitative evaluation; Thereby, also enriched simultaneously and estimated the technological means that the soil utilizes in the urban ecological planning for ecological The Study on Resources and ecological suitability evaluation provide foundation.

The soil utilizes the evaluation method representativeness to mainly contain following several kinds at present:

1, intensive utilization land valuation method: be one of main path of cities and towns sound development.Urban land intensive utilization notion is defined; From economic benefit target, social benefit target, ecological environment benefit target and four aspects of dynamic object, set up the assessment indicator system that the urban land utilizes intensity, according to the significance level of each index factor; Confirm the weight of evaluation index; On standardized basis, according to the index weight, with each desired value weighted sum to raw data; Obtain the soil and utilize intensive aggregative index (intensity index), confirm that quantitatively the urban land utilizes intensive level.

2, the soil utilizes dynamic remote monitoring area precision evaluation method: along with the application of remote sensing image in all kinds of planning in city; The research soil utilizes the evaluation method and the evaluation index of remote sensing monitoring figure spot area precision, is the survey of land and resources service of carrying out year by year.Research method: Monte Carlo method, the utilization of relative error theory and law of propagation of errors.Result of study: drawn the error profile rule, set up evaluation index, formed systematized variation monitoring spatial sampling and area precision evaluation method system.

3, folded diagram technology: in city planning and the ecological planning each is represented a natural resources element (such as soil and vegetation); Then many transparent drawings are fallen and analyzed inference together, on another drawing, drawn the pattern that new soil utilizes.All drawings are drawn according to the same engineer's scale.

Summary of the invention

The present invention one of is intended to solve the problems of the technologies described above at least to a certain extent or provides a kind of useful commerce to select at least.For this reason, one object of the present invention is to propose a kind of soil class factor of utilizing and carries out the quantification partitioned method.

The soil class factor of utilizing according to the embodiment of the invention is carried out the quantification partitioned method, may further comprise the steps: step S100, investigate measurement to a plurality of soils single-factor, and the result be depicted as thematic maps in the target area; Step S200 uses the GIS vector quantization with each said soil single-factor, sets up the unified coordinate system system, draws each said soil single-factor figure; Step S300 carries out classification with each said soil single-factor according to preset grade scale; Step S400 adopts analytical hierarchy process to confirm each monofactorial weighted value in said soil; Step S500 adopts GIS that said a plurality of soils single-factor is superposeed mutually; Step S600 obtains soil class results of factor analysis to each said soil single-factor weighted mean; With step S700, adopt and force normal distribution method that said soil class results of factor analysis is classified with completion soil class factorial analysis figure, thereby the quantification subregion is carried out in the target area.

The soil class factor of utilizing according to the embodiment of the invention is carried out the quantification partitioned method; Through GIS the indication information of soil each side is comprehensively superposeed; Utilize analytical hierarchy process to set up the quantification system, the analyzing of quantification, the science utilization and the planning in soil of analyzing and researching objectively more to the soil utilization class factor of target area; And then utilize class factor pair target area, soil to carry out quantitative ecological zoning, the auxiliary planning boundary of confirming the target area construction land.

In addition, according to the above embodiment of the present inventionly carry out the quantification partitioned method with the soil class factor and can also have following additional technical characterictic:

According to one embodiment of present invention, said a plurality of soils single-factor comprises land type and formation, fertility, the face of land, soil physical state.

According to one embodiment of present invention, said step S200 may further comprise the steps:

Step S201, with said each soil single-factor vector quantization, polar plot is continuously nonoverlapping unique territory, unique territory is meant that each key element of vector quantization interconnects but non-overlapping copies, need meet the following conditions

S1U?S2U...USn＝S

Si∩Sj＝0

S is the monofactorial scope in each soil, and S1, S2, Si, Sj, Sn etc. are meant a plurality of small area that each single-factor the inside, soil is divided, and n is the zonule sum; With

Step S202 sets up unified coordinate system system, to each the soil single-factor figure behind the vector quantization in GIS, define projection, conversion is projected as the Gauss projection coordinate system and draws each said soil single-factor figure.

According to one embodiment of present invention, said step S400 may further comprise the steps:

Step S401, said a plurality of soils single-factor is n, relatively n soil single-factor yz=(yz ₁, yz ₂... yz _n) to the influence of target zz (present status of land utilization), thus confirm their shared proportion in z, get two factor y at every turn _iAnd y _jUse a _IjExpression y _iWith y _jTo the ratio of the influence degree of z, measure a in 1～9 ratio scale _Ij, with numeral 1～9 and reciprocal, obtain weighted value by the difference of monofactorial importance value of each land use of quantification shown in the table 1 and influence property as scale,

Table 1 weight fiducial value

Judgment value	Comparison	Intensity
			1	Y i＝Y j	Equate
3	Y i＞∞Y j	Y iBe better than Y slightly j
			5	Y i＞＞Y j	Y iObviously be better than Y j
7	Y i＞＞＞Y j	Y iCompare Y jMuch better
			9	Y i＞Y j	Y iExtremely be better than Y j
1/3	Y i＜Y j	Y iSlightly inferior to Y j
			1/5	Y i＜＜Y j	Y iObviously inferior to Y j
1/7	Y i＜＜＜Y j	Y iCompare Y jMuch inferior
			1/9	Y i＜∞Y j	Y iDefinitely inferior to Y j

Annotate: 2,4,6,8 and 1/2,1/4,1/6,1/8 expression intensity is between corresponding adjacent rank;

Step S402, n judgment matrix A=(a that is compared (in pairs relatively) by one of soil single-factor formation relatively in twos _Ij) _{N * n}., obviously, judgment matrix has character:

$A=\left(\begin{array}{cccc}{a}_{11}& a_{12}& \·\·\·& a_{1n}\\ a_{21}& a_{22}& \·\·\·& a_{2n}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ a_{n1}& a_{n2}& \·\·\·& a_{\mathrm{nn}}\end{array}\right)$

a _ij＞0， a _ij＝1(i，j＝1，2，...n)

Step S403, weight coefficient calculates, the product Z of each row of judgment matrix _iThen press Formula calculates Z _iN th root W _i, to vectorial W _i=[W ₁, W ₂... W _n] ^TCarry out normalization and calculate, obtain proper vector, calculate maximum characteristic root $\mathrm{\λ}\mathrm{Max}={\mathrm{\Σ}}_{i=1}^n\frac{\left(\mathrm{KW}\right)i}{\mathrm{NWi}}$

Step S404 carries out consistency check to each the soil single-factor weighted value that obtains, and the check formula is CI=(λ max-n)/(n-1)

Step S405 obtains the monofactorial weight in each soil mean value of wi as a result,

$\mathrm{Wi}=\frac{{\mathrm{\Σ}}_{i=1}^m\mathrm{Di}}{m}$ And w1+w2+...+wn=1

Wherein wi is a land use single-factor weight, and m is for calculating the sum of weight, and di is the monofactorial weighted value of this land use.

According to one embodiment of present invention, the average weighted formula among the said step S600 is:

Lt＝f1*w1+f2*w2+...+fn*wn

Wherein Lt is meant that the soil utilizes the numerical value of type factor, and f1-fn is meant each soil single-factor grade separation value, and w1-wn is that each soil utilizes monofactorial weighted value.

Additional aspect of the present invention and advantage part in the following description provide, and part will become obviously from the following description, or recognize through practice of the present invention.

Description of drawings

Above-mentioned and/or additional aspect of the present invention and advantage are from obviously with easily understanding becoming the description of embodiment in conjunction with figs, wherein:

Fig. 1 carries out the process flow diagram of the method for quantitative ecological zoning according to the soil class factor of utilizing of the embodiment of the invention;

Fig. 2 carries out the synoptic diagram that in the method for quantitative ecological zoning a plurality of soils single-factor is superposeed mutually according to the soil class factor of utilizing of the embodiment of the invention; With

Fig. 3 is the soil class factorial analysis figure that adopts the method for utilizing the soil class factor to carry out quantitative ecological zoning according to the embodiment of the invention to obtain.

Embodiment

Describe embodiments of the invention below in detail, the example of said embodiment is shown in the drawings, and wherein identical from start to finish or similar label is represented identical or similar elements or the element with identical or similar functions.Be exemplary through the embodiment that is described with reference to the drawings below, be intended to be used to explain the present invention, and can not be interpreted as limitation of the present invention.

In description of the invention; It will be appreciated that; The orientation of indications such as term " " center ", " vertically ", " laterally ", " length ", " width ", " thickness ", " on ", D score, " preceding ", " back ", " left side ", " right side ", " vertically ", " level ", " top ", " end " " interior ", " clockwise ", " counterclockwise " or position relation are for based on orientation shown in the drawings or position relation; only be to describe with simplifying for the ease of describing the present invention; rather than the device or the element of indication or hint indication must have specific orientation, with specific azimuth configuration and operation, therefore can not be interpreted as limitation of the present invention.

In addition, term " first ", " second " only are used to describe purpose, and can not be interpreted as indication or hint relative importance or the implicit quantity that indicates indicated technical characterictic.Thus, one or more a plurality of these characteristics can be shown or impliedly comprised to the characteristic that is limited with " first ", " second " clearly.In description of the invention, except as otherwise noted, the implication of " a plurality of " is two or more, only if clear and definite qualification is arranged in addition.

In the present invention, only if clear and definite regulation and qualification are arranged in addition, broad understanding should be done in terms such as term " installation ", " linking to each other ", " connection ", " fixing ", for example, can be to be fixedly connected, and also can be to removably connect, or connect integratedly; Can be mechanical connection, also can be to be electrically connected; Can be directly to link to each other, also can link to each other indirectly through intermediary, can be the connection of two element internals.For those of ordinary skill in the art, can understand above-mentioned term concrete implication in the present invention as the case may be.

In the present invention; Only if clear and definite regulation and qualification are arranged in addition; First characteristic second characteristic it " on " or D score can comprise that first and second characteristics directly contact, can comprise that also first and second characteristics are not directly contacts but through the contact of the additional features between them.And, first characteristic second characteristic " on ", " top " and " above " comprise first characteristic directly over second characteristic and oblique upper, or only represent that the first characteristic level height is higher than second characteristic.First characteristic second characteristic " under ", " below " and " below " comprise first characteristic directly over second characteristic and oblique upper, or only represent that the first characteristic level height is less than second characteristic.

Describe in detail according to the soil class factor of utilizing of the embodiment of the invention with reference to accompanying drawing below and carry out the quantification partitioned method.

To shown in Figure 3, carry out the quantification partitioned method like Fig. 1, may further comprise the steps according to the soil class factor of utilizing of the embodiment of the invention:

Step S100 investigates measurement to a plurality of soils single-factor, and the result is depicted as thematic maps in the target area.Specifically, the soil of collecting the target area utilizes information, investigates, measures, analyzes, and is depicted as thematic maps.The information content should contain: land type and formation (comprising forest, grassland, farmland, water body, vegetation, urban land etc.), fertility (soil types, fertility etc.), the face of land, the soil physical state (gradient, aspect, elevation etc.) of utilizing the GIS technical Analysis to draw.

Step S200 uses the GIS vector quantization with each said soil single-factor, sets up the unified coordinate system system, draws each said soil single-factor figure.

For example, select the main soil single-factor in three types of land type and formation, fertility, the face of land, the soil physical states, utilize single-factor to use the GIS vector quantization in each soil, set up unified coordinate system, drafting soil single-factor figure.

Particularly, with each soil single-factor vector quantization, it is continuously nonoverlapping unique territory that polar plot requires.Unique territory is meant that each key element of vector quantization interconnects but non-overlapping copies.This be since in city planning every land used only corresponding to a kind of soil utilize character be determine.Need meet the following conditions

S1U?S2U...USn＝S

Si∩Sj＝0

Wherein, S is certain monofactorial scope, and S1, S2, Si, Sj, Sn etc. are meant a plurality of small area that certain single-factor the inside, soil is divided, and n is the zonule sum.

Then, set up unified coordinate system system, utilize single-factor figure to be unified in to the soil behind the vector quantization and define projection among the GIS, conversion is projected as Gauss projection coordinate system commonly used.

Step S300 carries out classification with each said soil single-factor according to preset grade scale.

To each soil single-factor of choosing according to city planning with build susceptibility and importance divided rank.Grade is low more, and susceptibility is strong more, and uncomfortable more Yicheng City is built; Otherwise grade is high more, and is suitable constructive strong more, meets the requirement of urban construction more.

Step S400 adopts analytical hierarchy process to confirm each monofactorial weighted value in said soil.

Particularly, a plurality of soils single-factor is n, relatively n soil single-factor yz=(yz ₁, yz ₂... yz _n) to the influence of target zz (present status of land utilization), thus confirm their shared proportion in z, get two factor y at every turn _iAnd y _jUse a _IjExpression y _iWith y _jTo the ratio of the influence degree of z, measure a in 1～9 ratio scale _IjWith numeral 1～9 and reciprocal, obtain weighted value by the difference of monofactorial importance value in the different soils of quantification shown in the table 1 and influence property as scale.

Table 1

Judgment value

Comparison

Intensity

1	Y i＝Y j	Equate
			3	Y i＞∞Y j	Y iBe better than Y slightly j
5	Y i＞＞Y j	Y iObviously be better than Y j
			7	Y i＞＞＞Y j	Y iCompare Y jMuch better
9	Y i＞Y j	Y iExtremely be better than Y j
			1/3	Y i＜Y j	Y iSlightly inferior to Y j
1/5	Y i＜＜Y j	Y iObviously inferior to Y j
			1/7	Y i＜＜＜Y j	Y iCompare Y jMuch inferior
1/9	Y i＜∞Y j	Y iDefinitely inferior to Y j

Annotate: 2,4,6,8 and 1/2,1/4,1/6,1/8 expression intensity is between corresponding adjacent rank.

N is constituted one by element relatively and compares the judgment matrix A=(a of (in pairs relatively) in twos _Ij) _{N * n}. obviously, judgment matrix has character:

$A=\left(\begin{array}{cccc}{a}_{11}& a_{12}& \·\·\·& a_{1n}\\ a_{21}& a_{22}& \·\·\·& a_{2n}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ a_{n1}& a_{n2}& \·\·\·& a_{\mathrm{nn}}\end{array}\right)$

a _ij＞0，

a _ij＝1(i，j＝1，2，...n)

Weight coefficient calculates.The product Z of each row of judgment matrix _iThen press

Formula calculates Z _iN th root W _iTo vectorial W _i=[W ₁, W ₂... W _n] ^TCarry out normalization and calculate, obtain proper vector.

Calculate maximum characteristic root λ max: $\mathrm{\λ}\mathrm{Max}={\mathrm{\Σ}}_{i=1}^n\frac{\left(\mathrm{KW}\right)i}{\mathrm{NWi}}$

Utilize the single-factor weighted value to carry out consistency check to the soil that obtains.

In order to check the consistency problem of judgment matrix; The consistance that needs to calculate coincident indicator CI judgment matrix also has randomness; This consistance at random can represent that the value of RI is relevant with the dimension size of matrix with mean random coincident indicator RI, relative uniformity index CR=CI/RI.It is generally acknowledged that when CR≤0.1, judgment matrix meets the crash consistency condition basically; When CR＞0.1, given judgment matrix does not meet the crash consistency condition, need adjust and revise.

CI＝(λmax-n)/(n-1)

Popular, if the soil utilizes 1 pair of single-factor of single-factor, 2 influence value a＞1 (being that single-factor 1 is bigger than single-factor 2), 2 pairs of single-factors of single-factor, 3 influence value b＞1, then in logic, 1 pair of single-factor 3 influence value of single-factor are c=ab＞1 certainly.But it is impossible when being configured to the comparison matrix, requiring to satisfy c=ab.Therefore move back and require paired comparator matrix that certain consistance is arranged.If c＜1, then the result is inconsistent, must recomputate weighted value.

Obtain each soil and utilize monofactorial weight w as a result _iMean value.

$\mathrm{Wi}=\frac{{\mathrm{\Σ}}_{i=1}^m\mathrm{Di}}{m}$ And w1+w2+...+wn=1

Wherein wi is a certain single-factor weight, and m is for calculating the sum of weight, and di is that this soil utilizes monofactorial some weighted values.

Step S500 adopts GIS that said a plurality of soils single-factor is superposeed mutually.As shown in Figure 2, utilize the GIS geographical information technology to merge stack each other each soil single-factor, stacked system is represented two different single-factors respectively for the key element of input, and the figure layer of output has kept two monofactorial all characteristics.

Step S600 obtains soil class results of factor analysis to each said soil single-factor weighted mean.

Lt＝f1*w1+f2*w2+...+fn*wn

Wherein Lt is meant that the soil utilizes the numerical value of type factor, and f1-fn is meant that each soil utilizes single-factor grade separation value, and w1-wn is that each soil utilizes monofactorial weighted value.

Step S700 adopts and forces normal distribution method that said soil class results of factor analysis is classified with completion soil class factorial analysis figure, thereby the quantification subregion is carried out in the target area.

Forcing normal distribution is according to the normal distribution principle; The regularity of distribution of " broad in the middle, two is little " that promptly is commonly called as; Confirm opinion rating and each grade shared number percent in sum in advance, confirm the protection border of the plant class factor then according to the analysis result of the plant class factor.

The soil class factor of utilizing according to the embodiment of the invention is carried out the quantification partitioned method; Through GIS the indication information of soil each side is comprehensively superposeed; Utilize analytical hierarchy process to set up the quantification system, the analyzing of quantification, the science utilization and the planning in soil of analyzing and researching objectively more to the soil utilization class factor of target area; And then utilize class factor pair target area, soil to carry out quantitative ecological zoning, the auxiliary planning boundary of confirming the target area construction land.Through using the present invention; Overcome soil in the ecological planning of town and country utilize present status of land utilization evaluation that type factor exists in estimating comprehensively with the direct defective of guiding plan construction; The soil that has proposed quantification utilizes a type factor evaluation method, has enriched the city planning method.

Be that example is sketched according to the soil class factor of utilizing of the embodiment of the invention and carried out the quantification partitioned method below with the Huairou District, Beijing City.

The target principle that land use is analyzed is to confirm that scope, intensity and the reasonable space thereof of mankind's activity distribute.Land resource map, land type map, low yield and bad land type map, land use degree and present landuse map that explanation land resource quality is arranged in its basic factor are referring to table 2.Select the utilization of comprehensive soil, land use, soil erosion property and land use degree dominant factor for participating in estimating.

Table 2 soil class single-factor complete list

Set up each soil and utilized monofactorial evaluation criterion.As:

Be suitable for water system, the water surface of nature, vegetable plot, forest land etc. is evaluated as 1, and irrigation canals and ditches nursery etc. is 2, and cities and towns industrial and mineral, bare area that is suitable for the human settlements system etc. is 5, and wasteland, nonirrigated farmland etc. is 4, and all the other are 3.

The evaluation criterion of low yield and bad land use degree is that semifixed dust storm soil, meadow bog soil, waterloggogenic paddy soil etc. are 1, and the mountain area is corroded ground and spent moderate to be respectively 2-4, flat sandy soil ground etc. be 3 by slight, and moist soil is 5.

The evaluation criterion of soil erosion property is to corrode slight corrosion area according to soil erosion degree and erosion grade by moderate to be respectively grade 1-5.

For land resource, should agricultural two grades also wait land suitability wide, yield-power is high, and land quality is better, is respectively 1; Second-class ground suitable for afforestation land suitability and yield-power are higher for being 2, and be should farming third-class and should to herd third-class be 3; Betterment of land difficulty be 4, very poor should not agro-farming be 5.

Adopt analytical hierarchy process to carry out weight calculation, the result sees table 3.Calculate through carrying out the spatial analysis stack in the GIS-Geographic Information System, ask for evaluation result.Formulate the suitable classification of evaluation set, the construction soil utilization class factor reflects that each regional soil utilizes the quality and power of ecological basal status.

Table 3 soil class factor expert weight is estimated consolidated statement

From the present status of land utilization evaluation result, see that from space layout the zone of suitable human settlements concentrates on local wide river valley in Plain and mountain area and cheuch area clearly, the especially suitable human settlements in wherein existing settlement are for preferred human settlements ground, like the a-quadrant among Fig. 3.The human settlements district in Plain zone, topology layout is with the gathering of existing human settlements ground, along Caobai River in poplar Song and north house dispersed placement.The zone that is suitable for nature is large-area massif and water system in the mountain area, and is stretched over the zone, Plain by the mountain area with passing through existing human settlements, like the zone of the B among Fig. 3.

In the description of this instructions, the description of reference term " embodiment ", " some embodiment ", " example ", " concrete example " or " some examples " etc. means the concrete characteristic, structure, material or the characteristics that combine this embodiment or example to describe and is contained at least one embodiment of the present invention or the example.In this manual, the schematic statement to above-mentioned term not necessarily refers to identical embodiment or example.And concrete characteristic, structure, material or the characteristics of description can combine with suitable manner in any one or more embodiment or example.

Although illustrated and described embodiments of the invention above; It is understandable that; The foregoing description is exemplary; Can not be interpreted as limitation of the present invention, those of ordinary skill in the art can change the foregoing description under the situation that does not break away from principle of the present invention and aim within the scope of the invention, modification, replacement and modification.

Claims (5)

1. one kind is utilized the soil class factor to carry out the quantification partitioned method, it is characterized in that, may further comprise the steps:

Step S100 investigates measurement to a plurality of soils single-factor, and the result is depicted as thematic maps in the target area;

Step S200 uses the GIS vector quantization with each said soil single-factor, sets up the unified coordinate system system, draws each said soil single-factor figure;

Step S300 carries out classification with each said soil single-factor according to preset grade scale;

Step S400 adopts analytical hierarchy process to confirm each monofactorial weighted value in said soil;

Step S500 adopts GIS that said a plurality of soils single-factor is superposeed mutually;

Step S600 obtains soil class results of factor analysis to each said soil single-factor weighted mean; With

Step S700 adopts and forces normal distribution method that said soil class results of factor analysis is classified with completion soil class factorial analysis figure, thereby the quantification subregion is carried out in the target area.

2. the soil class factor of utilizing according to claim 1 is carried out the quantification partitioned method, it is characterized in that, said a plurality of soils single-factor comprises land type and formation, fertility, the face of land, soil physical state.

3. the soil class factor of utilizing according to claim 1 is carried out the quantification partitioned method, it is characterized in that said step S200 may further comprise the steps:

Step S201, with said each soil single-factor vector quantization, polar plot is continuously nonoverlapping unique territory, unique territory is meant that each key element of vector quantization interconnects but non-overlapping copies, need meet the following conditions

S1U?S2U...USn＝S

Si∩Sj＝0

S is the monofactorial scope in each soil, and S1, S2, Si, Sj, Sn etc. are meant a plurality of small area that each single-factor the inside, soil is divided, and n is the zonule sum; With

Step S202 sets up unified coordinate system system, to each the soil single-factor figure behind the vector quantization in GIS, define projection, conversion is projected as the Gauss projection coordinate system and draws each said soil single-factor figure.

4. the soil class factor of utilizing according to claim 1 is carried out the quantification partitioned method, it is characterized in that said step S400 may further comprise the steps:

Step S401, said a plurality of soils single-factor is n, relatively n soil single-factor yz=(yz ₁, yz ₂... yz _n) to the influence of target zz (present status of land utilization), thus confirm their shared proportion in z, get two factor y at every turn _iAnd y _jUse a _IjExpression y _iWith y _jTo the ratio of the influence degree of z, measure a in 1～9 ratio scale _Ij, with numeral 1～9 and reciprocal, obtain weighted value by the difference of monofactorial importance value of each land use of quantification shown in the table 1 and influence property as scale,

Table 1 weight fiducial value

Judgment value Comparison Intensity 1 Y _i＝Y _j Equate 3 Y _i＞∞Y _j Y _iBe better than Y slightly _j 5 Y _i＞＞Y _j Y _iObviously be better than Y _j 7 Y _i＞＞＞Y _j Y _iCompare Y _jMuch better 9 Y _i＞Y _j Y _iExtremely be better than Y _j 1/3 Y _i＜Y _j Y _iSlightly inferior to Y _j 1/5 Y _i＜＜Y _j Y _iObviously inferior to Y _j 1/7 Y _i＜＜＜Y _j Y _iCompare Y _jMuch inferior 1/9 Y _i＜∞Y _j Y _iDefinitely inferior to Y _j

Annotate: 2,4,6,8 and 1/2,1/4,1/6,1/8 expression intensity is between corresponding adjacent rank;

Step S402, n judgment matrix A=(a that is compared (in pairs relatively) by one of soil single-factor formation relatively in twos _Ij) _{N * n}., obviously, judgment matrix has character:

$A=\left(\begin{array}{cccc}{a}_{11}& a_{12}& \·\·\·& a_{1n}\\ a_{21}& a_{22}& \·\·\·& a_{2n}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ a_{n1}& a_{n2}& \·\·\·& a_{\mathrm{nn}}\end{array}\right)$

a _ij＞0，

a _ij＝1(i，j＝1，2，...n)

Step S403, weight coefficient calculates, the product Z of each row of judgment matrix _iThen press

Formula calculates Z _iN th root W _i, to vectorial W _i=[W ₁, W ₂... W _n] ^TCarry out normalization and calculate, obtain proper vector, calculate maximum characteristic root $\mathrm{\λ}\mathrm{Max}={\mathrm{\Σ}}_{i=1}^n\frac{\left(\mathrm{KW}\right)i}{\mathrm{NWi}}$

Step S404 carries out consistency check to each the soil single-factor weighted value that obtains, and the check formula is CI=(λ max-n)/(n-1)

Step S405 obtains the monofactorial weight in each soil mean value of wi as a result,

$\mathrm{Wi}=\frac{{\mathrm{\Σ}}_{i=1}^m\mathrm{Di}}{m}$ And w1+w2+...+wn=1

Wherein wi is a land use single-factor weight, and m is for calculating the sum of weight, and di is the monofactorial weighted value of this land use.

5. the soil class factor of utilizing according to claim 1 is carried out the quantification partitioned method, it is characterized in that the average weighted formula among the said step S600 is:

Lt＝f1*w1+f2*w2+...+fn*wn

Wherein Lt is meant that the soil utilizes the numerical value of type factor, and f1-fn is meant each soil single-factor grade separation value, and w1-wn is that each soil utilizes monofactorial weighted value.

Images