Summary of the invention
In order to solve at least one above-mentioned technical problem, the invention proposes a kind of ground utilizations based on multi-source data
Quantitative estimation method, system and storage medium.
To achieve the goals above, first aspect present invention proposes a kind of ground utilization quantization based on multi-source data
Appraisal procedure, comprising:
Road network information is obtained, is refined area unit by region division to be assessed according to the road network information;
Government affairs geography information is obtained, geocoding is generated according to the government affairs geography information, and the geocoding is turned
It is changed to coordinate information;
Arest neighbors spatial match is carried out according to the refined area unit and coordinate information, the coordinate letter after being optimized
Breath;
Efficiency index calculating is carried out according to the corresponding plot region of coordinate information after the optimization, obtains calculated result;
The calculated result is shown by spatialization form.
In the present solution, it is described the geocoding is converted to coordinate information after, further includes:
According to the government affairs geography information, multiple operation information of target object are obtained;
The multiple operation information is subjected to geocoding, obtains multiple latitude and longitude informations and corresponding accuracy class;
The high latitude and longitude information of choice accuracy grade is as final target object location information.
In the present solution, described carry out arest neighbors spatial match according to the refined area unit and coordinate information, specifically
Are as follows:
The coordinate information for obtaining target object, judges whether the coordinate information is in preset plot region;
If being not in preset plot region, the coordinate in the plot region nearest apart from the coordinate information is obtained
Information;
Using the coordinate information in this plot region as the coordinate information after optimization.
In the present solution, in the progress arest neighbors spatial match, after the coordinate information after being optimized, further includes:
Plot area information is obtained, and the plot region will be expanded into level-one according to default rule, after obtaining expansion
Plot area information;
The index value total amount in the plot region after calculating the expansion;
Index value is distributed to multiple plot regions in the plot region to expansion according to default rule.
In the present solution, the corresponding plot region of coordinate information according to after the optimization carries out efficiency index calculating,
Calculated result is obtained, specifically:
The plot region is divided into first order region and second level region;
The ratio that the commercial power and industry for calculating first order region are paid taxes, obtains input and output index value;
The equal power consumption index of settlement place for calculating second level region is worth by the ratio of residential electricity consumption total amount and block area
To unit electricity consumption density.
In the present solution, further include:
Obtain first order region and the density of population in second level region, input and output index, unit electricity consumption density index;
By the density of population in the first order region and second level region, input and output index, unit electricity consumption density index
It is normalized;
The unit electricity consumption density criterion value in first order region is calculated multiplied by density of registered inhabitants standard value, then takes radical sign 2
Value;
The input and output index value in second level region is calculated multiplied by working population's density criterion value, then takes the value of radical sign 2.
Second aspect of the present invention also provides a kind of ground utilization quantitative evaluation system based on multi-source data, described to be based on
The ground utilization quantitative evaluation system of multi-source data includes: memory and processor, includes a kind of in the memory based on more
The ground utilization quantitative estimation method program of source data, the ground utilization quantitative estimation method program based on multi-source data
Following steps are realized when being executed by the processor:
Road network information is obtained, is refined area unit by region division to be assessed according to the road network information;
Government affairs geography information is obtained, geocoding is generated according to the government affairs geography information, and the geocoding is turned
It is changed to coordinate information;
Arest neighbors spatial match is carried out according to the refined area unit and coordinate information, the coordinate letter after being optimized
Breath;
Efficiency index calculating is carried out according to the corresponding plot region of coordinate information after the optimization, obtains calculated result;
The calculated result is shown by spatialization form.
In the present solution, it is described the geocoding is converted to coordinate information after, further includes:
According to the government affairs geography information, multiple operation information of target object are obtained;
The multiple operation information is subjected to geocoding, obtains multiple latitude and longitude informations and corresponding accuracy class;
The high latitude and longitude information of choice accuracy grade is as final target object location information.
In the present solution, described carry out arest neighbors spatial match according to the refined area unit and coordinate information, specifically
Are as follows:
The coordinate information for obtaining target object, judges whether the coordinate information is in preset plot region;
If being not in preset plot region, the coordinate in the plot region nearest apart from the coordinate information is obtained
Information;
Using the coordinate information in this plot region as the coordinate information after optimization.
In the present solution, in the progress arest neighbors spatial match, after the coordinate information after being optimized, further includes:
Plot area information is obtained, and the plot region will be expanded into level-one according to default rule, after obtaining expansion
Plot area information;
The index value total amount in the plot region after calculating the expansion;
Index value is distributed to multiple plot regions in the plot region to expansion according to default rule.
In the present solution, the corresponding plot region of coordinate information according to after the optimization carries out efficiency index calculating,
Calculated result is obtained, specifically:
The plot region is divided into first order region and second level region;
The ratio that the commercial power and industry for calculating first order region are paid taxes, obtains input and output index value;
The equal power consumption index of settlement place for calculating second level region is worth by the ratio of residential electricity consumption total amount and block area
To unit electricity consumption density.
In the present solution, further include:
Obtain first order region and the density of population in second level region, input and output index, unit electricity consumption density index;
By the density of population in the first order region and second level region, input and output index, unit electricity consumption density index
It is normalized;
The unit electricity consumption density criterion value in first order region is calculated multiplied by density of registered inhabitants standard value, then takes radical sign 2
Value;
The input and output index value in second level region is calculated multiplied by working population's density criterion value, then takes the value of radical sign 2.
Third aspect present invention also provides a kind of computer readable storage medium, in the computer readable storage medium
Including a kind of ground utilization quantitative estimation method program based on multi-source data, the ground utilization amount based on multi-source data
When change appraisal procedure program is executed by processor, a kind of land used effect based on multi-source data as described in any one of the above embodiments is realized
The step of rate quantitative estimation method.
Additional aspect and advantage of the invention will provide in following description section, partially will from the following description
Become obvious, or practice through the invention is recognized.
Specific embodiment
To better understand the objects, features and advantages of the present invention, with reference to the accompanying drawing and specific real
Applying mode, the present invention is further described in detail.It should be noted that in the absence of conflict, the reality of the application
The feature applied in example and embodiment can be combined with each other.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, still, the present invention may be used also
To be implemented using other than the one described here other modes, therefore, protection scope of the present invention is not by described below
Specific embodiment limitation.
Fig. 1 shows the flow chart of the ground utilization quantitative estimation method the present invention is based on multi-source data.
As shown in Figure 1, to achieve the goals above, first aspect present invention proposes a kind of use based on multi-source data
Ground efficiency quantitative estimation method, comprising:
S102 obtains road network information, is refined area unit by region division to be assessed according to the road network information.
S104 obtains government affairs geography information, generates geocoding according to the government affairs geography information, and by the geography
Code conversion is coordinate information.
S106 carries out arest neighbors spatial match according to the refined area unit and coordinate information, the seat after being optimized
Mark information.
S108 carries out efficiency index calculating according to the corresponding plot region of coordinate information after the optimization, is calculated
As a result.
S110 shows the calculated result by spatialization form.
It is obtained it should be noted that obtaining road network information and can be by the data of some navigation softwares or backstage of navigating
Road network information is taken, the territory status road network information obtained is also possible to.For example, obtaining road network by Amap navigation software
Information.Then the village area Nei Jiu will be studied by fine navigation road network and old factory site block carries out piecemeal, that is evaluated is fine
Unit.
It further include the process manually checked specifically, being in refined area unit in region division to be assessed, tool
Body includes:
The first refined area unit is generated according to road network information;
Receive the parameter of the first refined area unit of change;
According to the first refined area unit described in the parameter modification.
It should be noted that can manually be proofreaded after having obtained the first refined area unit, it is fine by first
The parameter of territory element changes, and the first new refined area unit is obtained, as final refined area unit.
It should be noted that obtaining government affairs geography information, geocoding is generated according to the government affairs geography information, and will
The geocoding is converted to coordinate information, can be supplemented by government affairs geodata, by navigation software, for example, high moral
Figure realizes geocoding API and carries out coordinate transformation to literal address.Government affairs geography information therein can be government affairs backstage
Acquisition can also be by technical staff's typing.
It should be noted that arest neighbors spatial match is carried out according to the refined area unit and coordinate information in step,
In coordinate information after being optimized.The problem of due to Position disturbance or geocoding deviation, will lead to the presence of partial coordinates point
Certain random offset error, small part point can deviate evaluation unit.So being entangled by the space that this step carries out deviation point
Just.
It should be noted that after obtaining the coordinate information after optimization, then according to the coordinate information after the optimization
Corresponding plot region carries out efficiency index calculating, obtains calculated result.Wherein the plot region is object to be measured region,
It can be a street, a village, small towns etc. can specifically be determined to be measured with the actual needs of those skilled in the art
Target area.After optimizing coordinate information, efficiency index calculating can be carried out for these regions, the efficiency refers to
Mark can not be directed to the index and specifically be limited for indexs, the present invention such as electricity consumption, financial data, resident's density, but appoint
The method what carries out quantitative statistics according to the technical solution of the present invention is fallen in the scope of the present invention.Calculating
After stating result, the calculated result can be shown by way of spatialization.Wherein, the spatialization form can be with
For 2D map, 3D map, data form display form.Those skilled in the art can spatialization determine according to actual needs
The form of expression.
According to embodiments of the present invention, after the geocoding to be converted to coordinate information described, further includes:
According to the government affairs geography information, multiple operation information of target object are obtained;
The multiple operation information is subjected to geocoding, obtains multiple latitude and longitude informations and corresponding accuracy class;
The high latitude and longitude information of choice accuracy grade is as final target object location information.
It should be noted that obtaining multiple operation information of government affairs information, is then determined and corresponded to by these operation information
The address of record, obtains geocoding.Wherein each information has corresponded to latitude and longitude information, that is, coordinate information, and corresponding
Accuracy class.Accuracy class therein is the specific location class information for indicating geographical location, and higher grade, indicates address
More more detailed, more specific.Finally, the latitude and longitude information for selecting choice accuracy grade high is as final target object location information.
Namely select location information of the more detailed coordinate in address as target object.
For example, each record often has two address informations in government affairs geodata, paying taxes in data such as industry includes
The two attributes of taxpayer's title and production and operation address, same record all include the description information in geographical location, together
When two addresses are subjected to geocoding respectively, the latitude and longitude coordinates and positioning accuracy of two interpretations are obtained, by positioning accuracy
It is high as final positioning result.
According to embodiments of the present invention, described that arest neighbors space is carried out according to the refined area unit and coordinate information
Match, specifically:
The coordinate information for obtaining target object, judges whether the coordinate information is in preset plot region;
If being not in preset plot region, the coordinate in the plot region nearest apart from the coordinate information is obtained
Information;
Using the coordinate information in this plot region as the coordinate information after optimization.
It should be noted that geocoding result is sometimes survey office, state coordinate (GCJ02), there are encrypted location interference,
Cause partial coordinates point there are certain random offset error, small part point can deviate evaluation unit, so needing to carry out nearest
Adjacent spatial match.Firstly, obtaining the coordinate information of target object, judge whether the coordinate information is in preset plot area
In domain.Preset plot region can be the regions such as non-rice habitats region, lake, massif, for example, if coordinate information is in certain master
On arterial highway, then show to be not in predeterminable area.If being not in preset plot region, obtains and believe apart from the coordinate
Cease the coordinate information in nearest plot region.Wherein the plot region needs to be the plot for meeting predeterminable area.It obtains
After the region of plot, using the coordinate information in this plot region as the coordinate information after optimization.Since plot region is one
Regional scope, so there are multiple coordinate points, it is preferred, therefore, that choosing the coordinate of the central point in the plot region as excellent
Coordinate information after change.Certainly, those skilled in the art can set the coordinate information position in plot region according to actual needs.
According to embodiments of the present invention, after the progress arest neighbors spatial match, the coordinate information after being optimized,
Further include:
Plot area information is obtained, and the plot region will be expanded into level-one according to default rule, after obtaining expansion
Plot area information;
The index value total amount in the plot region after calculating the expansion;
Index value is distributed to multiple plot regions in the plot region to expansion according to default rule.
It should be noted that since the index after the geocoding of part has coverage rate deficiency, it is of the invention thus
This is solved the problems, such as using the distribution method of plot group, expands level-one upwards especially by the basis of Land unit.It is first
First, plot area information is obtained, and the plot region will be expanded into level-one according to default rule, the ground after being expanded
Block area information.For example, if obtain plot area information be street, can expand level-one to commercial circle range, it is therein
Commercial circle range may include multiple plot regions, that is to say, that may include multiple streets.Then, after calculating the expansion
Plot region index value total amount.Index value distributed according still further to default rule more in the plot region to expansion
A plot region.Wherein, the plot region after expansion contains multiple plot regions.In addition, default rule can be technology
Personnel are preset, for example, index value total amount can be averagely divided in multiple plot regions, according to above-described embodiment,
Then the index value total amount of commercial circle can be averagely divided into each street in commercial circle.
According to embodiments of the present invention, the corresponding plot region of coordinate information according to after the optimization carries out efficiency
Index calculates, and obtains calculated result, specifically:
The plot region is divided into first order region and second level region;
The ratio that the commercial power and industry for calculating first order region are paid taxes, obtains input and output index value;
The equal power consumption index of settlement place for calculating second level region is worth by the ratio of residential electricity consumption total amount and block area
To unit electricity consumption density.
It should be noted that the first order region is the region Jiu Cun, second area is old plant area domain.Then, the is calculated
The ratio that the commercial power of level-1 area and industry are paid taxes, obtains input and output index value;Calculate the settlement place in second level region
Equal power consumption index obtains unit electricity consumption density by the ratio of residential electricity consumption total amount and block area.Those skilled in the art can
To calculate different index values according to actual needs, but the scheme that any quantitative estimation method through the invention is calculated is all
It will fall into the scope of the present invention.
According to embodiments of the present invention, further includes:
Obtain first order region and the density of population in second level region, input and output index, unit electricity consumption density index;
By the density of population in the first order region and second level region, input and output index, unit electricity consumption density index
It is normalized;
The unit electricity consumption density criterion value in first order region is calculated multiplied by density of registered inhabitants standard value, then takes radical sign 2
Value;
The input and output index value in second level region is calculated multiplied by working population's density criterion value, then takes the value of radical sign 2.
It should be noted that specific formula for calculation is as follows when carrying out first order region and the calculating of second level region:
Old station-service:
The ground utilization of old factory needs to take into account input and output and use intensity, and the inefficient of any index can all use ground
Efficiency decline.
Old village's land used:
The ground utilization in old village needs to take into account vacancy rate and use intensity, and the inefficient of any index can all use ground effect
Rate decline.
It should be noted that influencing to eliminate the dimension of different indexs and extreme value, the village Jiang Jiu and the population of old factory are close
Degree, input and output index, unit electricity consumption density index carry out segmentation normalized, and formula is as follows:
IQR=(Q3-Q1)
xlow_outlier=Q3-1.5×IQR
xup_outlier=Q3+1.5×IQR
Xnorm--- criterion value
xlow_outlier--- it is the interior limit lower limit of index
xup_outlier--- it is the interior limit upper limit of index
IQR --- quartile spacing
Q3--- index is compared with large quantile
Q1--- the smaller quantile of index
Fig. 2 shows the present invention is based on the block diagrams of the ground utilization quantitative evaluation system of multi-source data.
As shown in Fig. 2, second aspect of the present invention also provides a kind of ground utilization quantitative evaluation system based on multi-source data
System, the ground utilization quantitative evaluation system based on multi-source data include: memory 21, processor 22, in the memory
Including a kind of ground utilization quantitative estimation method program based on multi-source data, the ground utilization amount based on multi-source data
Change when appraisal procedure program is executed by the processor and realize following steps:
Road network information is obtained, is refined area unit by region division to be assessed according to the road network information;
Government affairs geography information is obtained, geocoding is generated according to the government affairs geography information, and the geocoding is turned
It is changed to coordinate information;
Arest neighbors spatial match is carried out according to the refined area unit and coordinate information, the coordinate letter after being optimized
Breath;
Efficiency index calculating is carried out according to the corresponding plot region of coordinate information after the optimization, obtains calculated result;
The calculated result is shown by spatialization form.
It is obtained it should be noted that obtaining road network information and can be by the data of some navigation softwares or backstage of navigating
Road network information is taken, the territory status road network information obtained is also possible to.For example, obtaining road network by Amap navigation software
Information.Then the village area Nei Jiu will be studied by fine navigation road network and old factory site block carries out piecemeal, that is evaluated is fine
Unit.
It further include the process manually checked specifically, being in refined area unit in region division to be assessed, tool
Body includes:
The first refined area unit is generated according to road network information;
Receive the parameter of the first refined area unit of change;
According to the first refined area unit described in the parameter modification.
It should be noted that can manually be proofreaded after having obtained the first refined area unit, it is fine by first
The parameter of territory element changes, and the first new refined area unit is obtained, as final refined area unit.
It should be noted that obtaining government affairs geography information, geocoding is generated according to the government affairs geography information, and will
The geocoding is converted to coordinate information, can be supplemented by government affairs geodata, by navigation software, for example, high moral
Figure realizes geocoding API and carries out coordinate transformation to literal address.Government affairs geography information therein can be government affairs backstage
Acquisition can also be by technical staff's typing.
It should be noted that arest neighbors spatial match is carried out according to the refined area unit and coordinate information in step,
In coordinate information after being optimized.The problem of due to Position disturbance or geocoding deviation, will lead to the presence of partial coordinates point
Certain random offset error, small part point can deviate evaluation unit.So being entangled by the space that this step carries out deviation point
Just.
It should be noted that after obtaining the coordinate information after optimization, then according to the coordinate information after the optimization
Corresponding plot region carries out efficiency index calculating, obtains calculated result.Wherein the plot region is object to be measured region,
It can be a street, a village, small towns etc. can specifically be determined to be measured with the actual needs of those skilled in the art
Target area.After optimizing coordinate information, efficiency index calculating can be carried out for these regions, the efficiency refers to
Mark can not be directed to the index and specifically be limited for indexs, the present invention such as electricity consumption, financial data, resident's density, but appoint
The method what carries out quantitative statistics according to the technical solution of the present invention is fallen in the scope of the present invention.Calculating
After stating result, the calculated result can be shown by way of spatialization.Wherein, the spatialization form can be with
For 2D map, 3D map, data form display form.Those skilled in the art can spatialization determine according to actual needs
The form of expression.
According to embodiments of the present invention, after the geocoding to be converted to coordinate information described, further includes:
According to the government affairs geography information, multiple operation information of target object are obtained;
The multiple operation information is subjected to geocoding, obtains multiple latitude and longitude informations and corresponding accuracy class;
The high latitude and longitude information of choice accuracy grade is as final target object location information.
It should be noted that obtaining multiple operation information of government affairs information, is then determined and corresponded to by these operation information
The address of record, obtains geocoding.Wherein each information has corresponded to latitude and longitude information, that is, coordinate information, and corresponding
Accuracy class.Accuracy class therein is the specific location class information for indicating geographical location, and higher grade, indicates address
More more detailed, more specific.Finally, the latitude and longitude information for selecting choice accuracy grade high is as final target object location information.
Namely select location information of the more detailed coordinate in address as target object.
For example, each record often has two address informations in government affairs geodata, paying taxes in data such as industry includes
The two attributes of taxpayer's title and production and operation address, same record all include the description information in geographical location, together
When two addresses are subjected to geocoding respectively, the latitude and longitude coordinates and positioning accuracy of two interpretations are obtained, by positioning accuracy
It is high as final positioning result.
According to embodiments of the present invention, described that arest neighbors space is carried out according to the refined area unit and coordinate information
Match, specifically:
The coordinate information for obtaining target object, judges whether the coordinate information is in preset plot region;
If being not in preset plot region, the coordinate in the plot region nearest apart from the coordinate information is obtained
Information;
Using the coordinate information in this plot region as the coordinate information after optimization.
It should be noted that geocoding result is sometimes survey office, state coordinate (GCJ02), there are encrypted location interference,
Cause partial coordinates point there are certain random offset error, small part point can deviate evaluation unit, so needing to carry out nearest
Adjacent spatial match.Firstly, obtaining the coordinate information of target object, judge whether the coordinate information is in preset plot area
In domain.Preset plot region can be the regions such as non-rice habitats region, lake, massif, for example, if coordinate information is in certain master
On arterial highway, then show to be not in predeterminable area.If being not in preset plot region, obtains and believe apart from the coordinate
Cease the coordinate information in nearest plot region.Wherein the plot region needs to be the plot for meeting predeterminable area.It obtains
After the region of plot, using the coordinate information in this plot region as the coordinate information after optimization.Since plot region is one
Regional scope, so there are multiple coordinate points, it is preferred, therefore, that choosing the coordinate of the central point in the plot region as excellent
Coordinate information after change.Certainly, those skilled in the art can set the coordinate information position in plot region according to actual needs.
According to embodiments of the present invention, after the progress arest neighbors spatial match, the coordinate information after being optimized,
Further include:
Plot area information is obtained, and the plot region will be expanded into level-one according to default rule, after obtaining expansion
Plot area information;
The index value total amount in the plot region after calculating the expansion;
Index value is distributed to multiple plot regions in the plot region to expansion according to default rule.
It should be noted that since the index after the geocoding of part has coverage rate deficiency, it is of the invention thus
This is solved the problems, such as using the distribution method of plot group, expands level-one upwards especially by the basis of Land unit.It is first
First, plot area information is obtained, and the plot region will be expanded into level-one according to default rule, the ground after being expanded
Block area information.For example, if obtain plot area information be street, can expand level-one to commercial circle range, it is therein
Commercial circle range may include multiple plot regions, that is to say, that may include multiple streets.Then, after calculating the expansion
Plot region index value total amount.Index value distributed according still further to default rule more in the plot region to expansion
A plot region.Wherein, the plot region after expansion contains multiple plot regions.In addition, default rule can be technology
Personnel are preset, for example, index value total amount can be averagely divided in multiple plot regions, according to above-described embodiment,
Then the index value total amount of commercial circle can be averagely divided into each street in commercial circle.
According to embodiments of the present invention, the corresponding plot region of coordinate information according to after the optimization carries out efficiency
Index calculates, and obtains calculated result, specifically:
The plot region is divided into first order region and second level region;
The ratio that the commercial power and industry for calculating first order region are paid taxes, obtains input and output index value;
The equal power consumption index of settlement place for calculating second level region is worth by the ratio of residential electricity consumption total amount and block area
To unit electricity consumption density.
It should be noted that the first order region is the region Jiu Cun, second area is old plant area domain.Then, the is calculated
The ratio that the commercial power of level-1 area and industry are paid taxes, obtains input and output index value;Calculate the settlement place in second level region
Equal power consumption index obtains unit electricity consumption density by the ratio of residential electricity consumption total amount and block area.Those skilled in the art can
To calculate different index values according to actual needs, but the scheme that any quantitative estimation method through the invention is calculated is all
It will fall into the scope of the present invention.
According to embodiments of the present invention, further includes:
Obtain first order region and the density of population in second level region, input and output index, unit electricity consumption density index;
By the density of population in the first order region and second level region, input and output index, unit electricity consumption density index
It is normalized;
The unit electricity consumption density criterion value in first order region is calculated multiplied by density of registered inhabitants standard value, then takes radical sign 2
Value;
The input and output index value in second level region is calculated multiplied by working population's density criterion value, then takes the value of radical sign 2.
Third aspect present invention also provides a kind of computer readable storage medium, in the computer readable storage medium
Including a kind of ground utilization quantitative estimation method program based on multi-source data, the ground utilization amount based on multi-source data
When change appraisal procedure program is executed by processor, a kind of land used effect based on multi-source data as described in any one of the above embodiments is realized
The step of rate quantitative estimation method.
Technical solution in order to better illustrate the present invention can be described in detail by following embodiment.
Fig. 3 shows the method flow diagram of quantitative estimation method embodiment of the present invention.
The present invention is on the basis of obtaining the piecemeal in old village and old factory, by geo-coding techniques to residential electricity consumption, industry
Electricity consumption, industry such as pay taxes at the government datas and Tencent's population location data carries out quantification to the service efficiency of two kinds of lands used and comments
Estimate, specifically includes the following steps:
Step 1: evaluation unit obtains, and will study the village area Nei Jiu by fine navigation road network and old factory site block carries out
Piecemeal, the fine cells evaluated;
Step 2: government data geocoding, government affairs geodata supplement are real by the geocoding API of Amap
Coordinate transformation now is carried out to literal address;
Step 3: government data geocoding precision optimizing, there are two addresses letters in each record in government affairs geodata
Breath, Fig. 4 show the schematic diagram data of encoding precision of the present invention optimization, as shown in figure 4, as industry is paid taxes in data comprising receiving
Tax name claims and production and operation address, the two attributes of same record all include the description information in geographical location, simultaneously
Two addresses are subjected to geocoding respectively, the latitude and longitude coordinates and positioning accuracy of two interpretations are obtained, by positioning accuracy height
As final positioning result.
Step 4: government data arest neighbors spatial match, geocoding result are state's survey office coordinate (GCJ02), exist and add
Close Position disturbance causes partial coordinates point there are certain random offset error, and small part point can deviate evaluation unit, uses for reference
Shift point processing method in GPS track data handling procedure is carried out space to offset point and is entangled by arest neighbors spatial match
Just, the function can be realized by Python code program.As shown in figure 5, obtain it is nearest close on plot, nearest is faced
Restitution point of the near-earth block as offset point.
Step 5: there is coverage rate deficiency in the index after the geocoding of part, we are using plot group's thus
Distribution method solves the problems, such as this, expands level-one upwards especially by the basis of Land unit, using committee, village as basic unit
Plot is grouped, the plot group in committee, same village is combined into a plot group, will fall into intramassif index value total
Meter calculates, then entrusts range assignment to plot group by village, to promote the covering of plot index value in committee, village.
Step 6: government affairs index calculates: old factory site block efficiency index chooses input and output from economic benefit angle
Than the ratio (annual commercial power and data of paying taxes) paid taxes by commercial power and industry obtains input and output index value
(member/degree);The plot Jiu Cun efficiency index chooses the equal power consumption index of settlement place (1 year inhabitation electricity consumption number from vacant angle
According to), unit electricity consumption density (degree/square metre/day) is obtained by the ratio of residential electricity consumption total amount and block area.
Step 7: Tencent's density of population index calculates, density data of population (people/square metre) go out from use intensity angle
Hair carries out operation to the service efficiency in old village and old factory.Long period population location data in research area is grabbed by network
(working day and day off), to analyze the resident population working population density (9 points to 6 points at night of morning) in area and live
The density of population (at night 7 points to morning next day 8 points).
Step 8: ground utilization index comprehensive calculates, the comprehensive inefficient index operation of Niche Model, specific formula for calculation
It is as follows:
Old station-service:
The ground utilization of old factory needs to take into account input and output and use intensity, and the inefficient of any index can all use ground
Efficiency decline.
Old village's land used:
The ground utilization in old village needs to take into account vacancy rate and use intensity, and the inefficient of any index can all use ground effect
Rate decline.
It should be noted that influencing to eliminate the dimension of different indexs and extreme value, the village Jiang Jiu and the population of old factory are close
Degree, input and output index, unit electricity consumption density index carry out segmentation normalized, and formula is as follows:
IQR=(Q3-Q1)
xlow_outlier=Q3-1.5×IQR
xup_outlier=Q3+1.5×IQR
Xnorm--- criterion value
xlow_outlier--- it is the interior limit lower limit of index
xup_outlier--- it is the interior limit upper limit of index
IQR --- quartile spacing
Q3--- index is compared with large quantile
Q1--- the smaller quantile of index
Step 9: ground utilization overall target assignment, by ground utilization overall target assignment to plot, as guidance
Plot redevelopment timing standard.The plot Jiu Chang, Jiu Cun is ascending by ground utilization comprehensive index value, and Temporal development is by arriving first
Afterwards.
For example, as shown in fig. 7, being divided based on navigation road network to research plot, 3.17 hectares of plot average area.
By the way that after government data progress geocoding precision optimizing, as a result, data volume of paying taxes is 28243, arrangement amount is 5291, essence
It is 18.73% that degree, which promotes accounting,.The matching of space arest neighbors, matching rate 27.94% is promoted on the basis of step 3, is reduced
The loss of valid data.Data volume of paying taxes is 28243, and matching amount is 7891, and it is 27.94% that matching rate, which promotes accounting,.Index
Block group distribution (corresponding step 5) is evaluated industrial construction land used in survey region by industry data of paying taxes, due to depositing
Mixing situation (part cottage industry workshop etc.) in land used, therefore there are achievement data deficient phenomenas for the industrial plot in part.Such as
The practical coverage rate for data of paying taxes is 62.17%, is promoted by coverage rate after the group's distribution of plot to 92.26%.It pays taxes index
Drop point coverage rate is 62.17%, and plot group's allocation overlay rate is 92.26%, and coverage rate promotes accounting 30.09%.Such as Fig. 7 institute
Show, this is quantitative evaluation result schematic diagram, more intuitive can find out accurately quantizating index.
The present invention is by promoting precision methods, the company for including to government data itself based on multi-address geocoding
Enterprise name and address information carry out geocoding respectively, by comparing positioning accuracy grade, choose high-grade positioning accuracy
Data are as final spatial position.And space arest neighbors matching process is used, may be implemented by arest neighbors matching algorithm
Point data is matched in the planar data of arest neighbors, " drift " phenomenon of data positioning is solved, improves data utilization rate.Separately
Outside, also using plot group's Allocation method, by the way that government affairs Distribution Indexes to plot group, are improved the coverage rate of data.Most
Eventually, the evaluation space implementation method that government data is realized in conjunction with geo-coding techniques is positioned in combination with the real-time population of Tencent
Data correction constructs fine real-time evaluation result to Jiu Cun, old factory's ground utilization comprehensive assessment.
The present invention is by realizing exact evaluation to Jiu Cun, old factory's ground utilization based on geo-coding techniques, and in geography
Two technological improvements have been carried out in cataloged procedure, improve the precision and utilization rate of data, while in the distribution of data target
It is reformed, improves data cover rate.
The development of network positions data also provides preferable technology hand for land present status population real-time investigation at this stage simultaneously
Section provides effective data source for the use intensity assessment in soil.We use urban land ecological niche in soil efficiency evaluation
Model carries out comprehensive operation to ground utilization index and promotes the scientific and comprehensive of evaluation result.
After obtaining fine space orientation index, it is thus necessary to determine that fine evaluation unit passes through fine navigation road network
The village Jiang Jiu and old factory site block carry out piecemeal, can obtain the basic unit of corresponding urban land functional analysis, the soil in unit
Ground homogeney is strong, is suitable for the fine granularity basic unit as analysis urban issues.
In several embodiments provided herein, it should be understood that disclosed device and method can pass through
Other modes are realized.Apparatus embodiments described above are merely indicative, for example, the division of the unit, only
For a kind of logical function partition, there may be another division manner in actual implementation, such as: multiple units or components can combine,
Or it is desirably integrated into another system, or some features can be ignored or not executed.In addition, shown or discussed each composition
Partially mutual coupling or direct-coupling or communication connection can be through some interfaces, equipment or unit it is indirect
Coupling or communication connection, can be electrical, mechanical or other forms.
Above-mentioned unit as illustrated by the separation member, which can be or may not be, to be physically separated, aobvious as unit
The component shown can be or may not be physical unit;Both it can be located in one place, and may be distributed over multiple networks
On unit;Some or all of units can be selected to achieve the purpose of the solution of this embodiment according to the actual needs.
In addition, each functional unit in various embodiments of the present invention can be fully integrated in one processing unit, it can also
To be each unit individually as a unit, can also be integrated in one unit with two or more units;It is above-mentioned
Integrated unit both can take the form of hardware realization, can also realize in the form of hardware adds SFU software functional unit.
Those of ordinary skill in the art will appreciate that: realize that all or part of the steps of above method embodiment can lead to
The relevant hardware of program instruction is crossed to complete, program above-mentioned can store in computer-readable storage medium, the program
When being executed, step including the steps of the foregoing method embodiments is executed;And storage medium above-mentioned include: movable storage device, it is read-only
Memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or
The various media that can store program code such as person's CD.
If alternatively, the above-mentioned integrated unit of the present invention is realized in the form of software function module and as independent production
Product when selling or using, also can store in a computer readable storage medium.Based on this understanding, of the invention
Substantially the part that contributes to existing technology can embody the technical solution of embodiment in the form of software products in other words
Out, which is stored in a storage medium, including some instructions are used so that a computer equipment
(can be personal computer, server or network equipment etc.) execute each embodiment the method for the present invention whole or
Part.And storage medium above-mentioned includes: that movable storage device, ROM, RAM, magnetic or disk etc. are various can store journey
The medium of sequence code.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, appoints
What those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, answer
It is included within the scope of the present invention.Therefore, protection scope of the present invention should be with the scope of protection of the claims
It is quasi-.