CN107239756A - Density of population analysis system based on high score satellite remote sensing date combination type of ground objects - Google Patents
Density of population analysis system based on high score satellite remote sensing date combination type of ground objects Download PDFInfo
- Publication number
- CN107239756A CN107239756A CN201710370420.2A CN201710370420A CN107239756A CN 107239756 A CN107239756 A CN 107239756A CN 201710370420 A CN201710370420 A CN 201710370420A CN 107239756 A CN107239756 A CN 107239756A
- Authority
- CN
- China
- Prior art keywords
- population
- density
- ground objects
- grid
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/50—Context or environment of the image
- G06V20/52—Surveillance or monitoring of activities, e.g. for recognising suspicious objects
- G06V20/53—Recognition of crowd images, e.g. recognition of crowd congestion
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of density of population analysis system based on high score satellite remote sensing date combination type of ground objects, including:Remotely-sensed data acquisition module, type of ground objects determining module, population coefficient determination module, grid partition module and density of population determining module.Compared to prior art, the density of population analysis system based on high score satellite remote sensing date combination type of ground objects that the application is provided, by the way that target area is divided into multiple grid, then the density of population in each grid is calculated respectively in units of grid, it is more accurate compared to prior art so as to calculate population dispersal more specific in target area.On the other hand, the application can accurately determine the type of ground objects composition in target area based on remotely-sensed data, based on above-mentioned type of ground objects composition, the corresponding population coefficient of different types of ground objects can be accurately determined, so as to ensure that the density of population finally calculated has the higher degree of accuracy.
Description
Technical field
The present invention relates to density of population analysis technical field, and in particular to one kind is based on high score satellite remote sensing date in combination
The density of population analysis system of species type.
Background technology
The density of population is the population lived on unit area soil, and it is the finger for representing the dense degree of regional population
Mark, level of economic development and urban construction level that can be for weighing a region etc., or country is grand with place
See regulation and control, urban development planning provide data supporting, in addition, accurately density of population distributed data contribute to enterprises and institutions and
Entrepreneur its make the decision-makings such as rational addressing, industrial pattern.
At present, the density of population be usually using administrative division as unit of account, such as the population of a certain counties and cities divided by
Area is the density of population as the counties and cities, and precision is very poor, and population is unknowable in the specific distribution situation of some in this county city
's.
To sum up, at present in the urgent need to a kind of higher density of population analysis system of precision.
The content of the invention
For defect of the prior art, the present invention provides a kind of based on high score satellite remote sensing date combination type of ground objects
Density of population analysis system, is macro adjustments and controls, the urban development rule in country and place to improve the precision of density of population calculating
Draw and data supporting is provided, and addressing for enterprises and institutions and entrepreneur and industrial pattern provide data and supported.
A kind of density of population analysis system based on high score satellite remote sensing date combination type of ground objects that the present invention is provided, bag
Include:Remotely-sensed data acquisition module, type of ground objects determining module, population coefficient determination module, grid partition module and the density of population
Determining module;Wherein,
The remotely-sensed data acquisition module, for obtaining the corresponding remotely-sensed data in target area;
The type of ground objects determining module, the type of ground objects group for determining the target area according to the remotely-sensed data
Into;
The population coefficient determination module, for determining the corresponding population system of difference type of ground objects in the target area
Number, the population coefficient is the ratio of the size of population and unit area;
The grid partition module, for target area to be divided into multiple grid;
The density of population determining module, for according to the type of ground objects of each grid composition and the population system
Number, calculates the density of population of each grid, to determine the population dispersal of the target area.
Optionally, the type of ground objects determining module, including:
Radar data type of ground objects determining unit, for special to the reflection of radar signal and scattering based on different types of ground objects
Property, the radar remote sensing data obtained according to the remotely-sensed data acquisition module determine that the type of ground objects of the target area is constituted.
Optionally, the type of ground objects determining module, including:
Multispectral data type of ground objects determining unit, for based on different types of ground objects to different-waveband spectral reflectivity
Difference, the Multi-spectral Remote Sensing Data obtained according to the remotely-sensed data acquisition module determines the type of ground objects group of the target area
Into.
Optionally, the population coefficient determination module, including:
Computing unit is returned, for the sample data according to the region of clear and definite type of ground objects composition and the size of population,
The corresponding population coefficient of each type of ground objects is calculated using regression algorithm.
Optionally, the density of population determining module, including:
Density of population determining unit, the density of population for calculating each grid according to following mathematical algorithm:
Wherein,The corresponding density of population of i-th of grid is represented, j numbers for different types of ground objects, ajRepresent jth kind
The corresponding population coefficient of type of ground objects, XjFor accounting of the area in the grid of jth kind type of ground objects, n represents the grid
The quantity of middle type of ground objects.
Optionally, the density of population analysis system based on high score satellite remote sensing date combination type of ground objects, in addition to:
First density of population optimization module, for the corresponding relation based on nighttime light intensity and the density of population, according to night
Between the density of population of each grid that calculates the density of population determining unit of light remotely-sensed data optimize, with excellent
Change the population dispersal of the target area.
Optionally, first density of population optimization module, including:
First density of population optimizes unit, for being carried out according to following mathematical algorithm to the density of population of grid each described
Optimization:
Wherein, PiThe corresponding density of population of i-th of grid obtained after optimization is represented,Represent that the density of population is determined
Unit calculates the corresponding density of population of i-th of grid obtained;LjThe corresponding intensity of light of j-th of grid is represented,Represent described
The average intensity of light of target area;PlRepresent the size of population that unit light is represented;S is regulation coefficient.
Optionally, the type of ground objects includes Urban Construction Land_use and urban residents' land used;
The density of population computing system, in addition to:
Second density of population optimization module, for the corresponding relation based on urban heat island strength and the density of population, according to red
The density of population for each grid that outer remotely-sensed data is calculated the density of population determining unit is optimized, to optimize
State the population dispersal of target area.
Optionally, second density of population optimization module, including:
Second density of population optimizes unit, for being carried out according to following mathematical algorithm to the density of population of grid each described
Optimization:
Wherein, PiThe corresponding density of population of i-th of grid obtained after optimization is represented,Represent that i-th of grid is corresponding
Belong to the density of population of Urban Construction Land_use,Represent the corresponding density of population for belonging to urban residents' land used of i-th of grid;Ii
The corresponding urban heat island strength of i-th of grid is represented,Represent average urban heat island strength;PIRepresent unit urban heat island strength
The size of population of representative;WiuThe corresponding population coefficient of Urban Construction Land_use, WirRepresent the corresponding population system of urban residents' land used
Number;AiuRepresent accounting of the area of Urban Construction Land_use in i-th of grid, AirRepresent the area of urban residents' land used i-th
Accounting in individual grid.
Optionally, the density of population analysis system based on high score satellite remote sensing date combination type of ground objects, in addition to:
Density of population distribution map generation module, will be each described for the mapping relations according to the density of population and different colours
Color corresponding with the grid density of population is filled in the corresponding position of grid, is distributed with the density of population for drawing the target area
Figure.
As shown from the above technical solution, what the present invention was provided is a kind of based on high score satellite remote sensing date combination type of ground objects
Density of population analysis system, including:Remotely-sensed data acquisition module, type of ground objects determining module, population coefficient determination module, grid
Division module and density of population determining module.Compared to prior art, it is described based on high score satellite remote sensing number that the application is provided
According to the density of population analysis system for combining type of ground objects, by the way that target area is divided into multiple grid, then using grid to be single
Position calculates the density of population in each grid respectively, so as to calculate density of population distribution more specific in target area
Situation is more accurate compared to prior art.On the other hand, the application, which is based on remotely-sensed data, can accurately determine target
Type of ground objects composition in region, based on above-mentioned type of ground objects composition, can accurately determine different type of ground objects correspondences
Population coefficient, so as to ensure that the density of population that finally calculates has the higher degree of accuracy.To sum up, can be more based on the application
Accurately, the population dispersal in target area is accurately determined, so as to be the macro adjustments and controls in country and place, city hair
Exhibition planning provides data supporting, and addressing for enterprises and institutions and entrepreneur and industrial pattern provide data and supported.
Brief description of the drawings
, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical scheme of the prior art
The accompanying drawing used required in embodiment or description of the prior art is briefly described.
Fig. 1 shows that one kind that first embodiment of the invention is provided is based on high score satellite remote sensing date combination type of ground objects
Density of population analysis system schematic diagram;
Fig. 2 shows a kind of schematic diagram of remotely-sensed data acquisition module;
Fig. 3 shows schematic diagram of each type of ground objects to the reflectivity of different-waveband spectrum;
Fig. 4 shows some region of population dispersal design sketch provided in an embodiment of the present invention;
Fig. 5 shows that the density of population that a certain region provided in an embodiment of the present invention optimizes through urban heat island strength is distributed feelings
Condition design sketch.
Embodiment
The embodiment of technical solution of the present invention is described in detail below in conjunction with accompanying drawing.Following examples are only used for
Clearly illustrate technical scheme, therefore be intended only as example, and the protection of the present invention can not be limited with this
Scope.
It should be noted that unless otherwise indicated, technical term or scientific terminology used in this application should be this hair
The ordinary meaning that bright one of ordinary skill in the art are understood.
The present invention provides a kind of density of population analysis system based on high score satellite remote sensing date combination type of ground objects.Below
Embodiments of the present invention are described with reference to the accompanying drawings.
Fig. 1 shows that one kind that first embodiment of the invention is provided is based on high score satellite remote sensing date combination type of ground objects
Density of population analysis system schematic diagram.As shown in figure 1, one kind that first embodiment of the invention is provided is distant based on high score satellite
The density of population analysis system of sense data combination type of ground objects includes:
Remotely-sensed data acquisition module 1, type of ground objects determining module 2, population coefficient determination module 3, the and of grid partition module 4
Density of population determining module 5;Wherein,
The remotely-sensed data acquisition module 1, for obtaining the corresponding remotely-sensed data in target area;
The type of ground objects determining module 2, the type of ground objects for determining the target area according to the remotely-sensed data
Composition;
The population coefficient determination module 3, for determining the corresponding population system of difference type of ground objects in the target area
Number, the population coefficient is the ratio of the size of population and unit area;
The grid partition module 4, for target area to be divided into multiple grid;
The density of population determining module 5, for according to the type of ground objects of each grid composition and the population system
Number, calculates the density of population of each grid, to determine the population dispersal of the target area.
Wherein, the type of ground objects is the classification of the different demarcation according to mulching material, can be clever according to the actual requirements
It is living to divide, for example, blue top building, red top building, cement top building, bare area, lake, river, farmland can be divided into
With forest land etc.;What the application was calculated is the density of population, because population distribution is in settlement place, accordingly it is also possible to which type of ground objects is straight
Connect and be divided into settlement place and non-resident;And consider cities and towns and rural difference, Urban Construction Land_use, rural area can also be divided into and occupied
Civilianly with non-resident land used;Etc..It is the change embodiment of the application above, those skilled in the art can make accordingly
Go out numerous variations embodiment, it all should be within the protection domain of the application.
With the development of remote sensing technology and high-resolution data acquisition technique, the resolution ratio more and more higher of remotely-sensed data, number
Increasingly enriched according to type, therefore, the type of ground objects on ground is made a distinction using high-definition remote sensing data and has been possibly realized
And the accuracy more and more higher of identification, based on this, the grid partition can also be obtained smaller, in terms of realizing the embodiment of the present invention
Calculate the purpose of more accurate, the accurate density of population.
Because the remote sensing mode that different remote sensing satellite is used is different, the remotely-sensed data of collection may be also different, for example I
No. 3 satellites of high score of state's transmitting and the Radarsat-2 satellites of Canada's transmitting are using synthetic aperture radar collection remote sensing number
According to its data mode is radar data, and No. 5 satellites of high score and Landsat series of satellites use full spectral coverage imager etc.
Instrument gathers remotely-sensed data, and its data mode is multispectral data, and above radar remote sensing data and Multi-spectral Remote Sensing Data all may be used
For being carried out to target area in the determination of type of ground objects, the embodiment of the present invention, the remotely-sensed data acquisition module 1 can root
Corresponding remotely-sensed data is selected according to the operation principle of the type of ground objects determining module 2 of configuration.
Wherein, the type of ground objects that can more accurately determine target area using high-quality multispectral data is constituted,
But because the weather such as sexual intercourse mist snow can influence the accuracy of multispectral data, therefore, in one embodiment that the application is provided,
Using according to weather condition selection, using by the way of different remotely-sensed datas, clear sky area uses Multi-spectral Remote Sensing Data, cloudy
Area of heavy rainfull then uses radar remote sensing data, so as to ensure the accuracy subsequently calculated to greatest extent since data source, has
Body embodiment is refer to Fig. 2, and it illustrates a kind of schematic diagram of remotely-sensed data acquisition module 1, the remotely-sensed data is obtained
Module 1 includes:Weather judging unit 11, radar data acquiring unit 12 and multispectral data acquiring unit 13;
The weather judging unit 11 is used for according to the weather conditions of target area selection triggering radar data acquiring unit
12 obtain radar remote sensing data or the acquisition Multi-spectral Remote Sensing Data of triggering multispectral data acquiring unit 13;Specifically, can be
Multispectral data acquiring unit 13 is triggered under the conditions of bright day gas and obtains Multi-spectral Remote Sensing Data, in day gas bars such as sexual intercourse mist snow
Radar data acquiring unit 12 is triggered under part and obtains radar remote sensing data;
The radar data acquiring unit 12 is used to obtain the target area under the triggering of the weather judging unit 11
The radar remote sensing data in domain;
The multispectral data acquiring unit 13 is used to obtain the target under the triggering of the weather judging unit 11
The Multi-spectral Remote Sensing Data in region.
Wherein, the weather judging unit 11, can be advance according to weather statistics result to the basis for estimation of weather conditions
Each department and the database of common weather conditions corresponding relation are set, and the weather judging unit 11 is as needed in real time from described
Transferred in database;It can also be the weather records for transferring the target area scheduled date, determine that target area is specified according to record
The weather conditions on date;It is the change embodiment of the application, within the protection domain of the application.
According to the difference of remotely-sensed data, the type of ground objects determining module 2 also carries out carrying for data in different ways
Take, for example, to radar remote sensing data, due to the difference of the layout of building, material, structure and surrounding environment, in SAR image
Building distribution in different textural characteristics, such as city is presented in (i.e. radar remote sensing data), and neatly, spacing is larger between building,
It is mostly the neat high-rise building of flat-top, there is good reflection rate mostly using material, strong luminance area is shown as on image,
And the road between building, coarse vegetation such as lawn etc., due to surface scattering, dark areas is shown as, therefore, city dweller
Light and dark texture is shown as on image, similitude is smaller;Be distributed to urban residents it is relatively at random, without evident regularity,
And the region such as road on image not substantially, therefore irregular speck shape is presented, similitude is larger.
Therefore, in one embodiment that the application is provided, the type of ground objects determining module 2, including:Radar data
Thing type determining units, for the reflection based on different types of ground objects to radar signal and scattering properties, according to the remote sensing number
The radar remote sensing data obtained according to acquisition module 1 determine the type of ground objects composition of the target area.
Specifically, so that the type of ground objects is settlement place and non-resident ground as an example, the radar data type of ground objects is determined
Unit, can be theoretical based on variogram, in analysis middle high-resolution SAR image on the basis of settlement place textural characteristics, adopts
With the Threshold based on iteration P parametric methods, weights are assigned to for the pixel point that meets threshold range, with increase settlement place with
The variogram on non-resident ground is poor, so as to extract settlement place, can not only ensure higher verification and measurement ratio, can also significantly reduce
False alarm rate, the above-mentioned technology from radar remote sensing extracting data settlement place is prior art, and the present embodiment is repeated no more, due to ground
The division of species type is different, and specific embodiment is also not quite similar, and those skilled in the art can carry out change implementation accordingly
Constituted with the type of ground objects for determining target area.
And for Multi-spectral Remote Sensing Data, the method that although prior art discloses the extracting section water surface, building, but hair
A person of good sense has found that its extraction accuracy, accuracy are unsatisfactory in the application, therefore, present applicant proposes more accurate, accuracy more
High mode, in one embodiment that the application is provided, the type of ground objects determining module 2, including:
Multispectral data type of ground objects determining unit, for based on different types of ground objects to different-waveband spectral reflectivity
Difference, the Multi-spectral Remote Sensing Data obtained according to the remotely-sensed data acquisition module 1 determines the type of ground objects of the target area
Composition.
Specifically, in one embodiment that the application is provided, the multispectral data type of ground objects determining unit, bag
Include:
Terrain classification subelement, for type of ground objects to be divided into blue top building, red top building, cement top building
Thing, bare area, lake, river, farmland and forest land;Wherein, blue top building, red top building, cement top building belong to settlement place,
For the main determination object of the embodiment of the present invention, other types of ground objects can be using universal formulation to be non-resident, due to non-resident
The size of population can consider to be zero, therefore it may only be necessary to be extracted according to the remotely-sensed data, settlement place can (other regions be straight
Connect and be defined as non-residently);
Atural object determines scheduling subelement, right for according to division result of the terrain classification subelement to type of ground objects
Belong to the type of ground objects of settlement place, call following extracting index to build subelement, exponential quantity computation subunit and binaryzation respectively
Processing subelement extracts the corresponding region of the atural object from the target area, so that it is determined that the ground species of the target area
Type is constituted;
Extracting index builds subelement, for according to type of ground objects to be extracted and other types of ground objects to different-waveband light
The difference of spectrum reflectivity builds the Objects extraction index that can make a distinction the type of ground objects to be extracted and other atural objects;
Exponential quantity computation subunit, for calculating the corresponding Objects extraction index of each pixel in the remotely-sensed data
Exponential quantity;
Binary conversion treatment subelement, for the exponential quantity of the Objects extraction index of each pixel to be carried out at binaryzation
Reason, and the remotely-sensed data is split according to binaryzation result, extract the corresponding region of type of ground objects to be extracted.
In the above-described embodiments, the terrain classification subelement according to different types of atural object to the anti-of different-waveband spectrum
The type of ground objects that the difference and settlement place of rate are included is penetrated, it is more careful, accurate that type of ground objects is divided into blue top building
(predominantly factory's canopy of enterprise), red top building (predominantly Hong Ding houses, are partly factory of enterprise canopy), cement top building are (main
Will be for town dweller area, road etc.), bare area, lake (man-made lake, reservoir etc.), river, farmland (vegetation) and forest land,
Because the corresponding population coefficient of different types of ground objects is different, so careful division is favorably improved the density of population finally calculated
Accuracy.
The extracting index builds subelement by relatively more each type of ground objects to the reflectivity of different-waveband spectrum, Jin Ergen
The difference of different-waveband spectral reflectivity is built according to type of ground objects to be extracted and other types of ground objects can will be to be extracted
The Objects extraction index that type of ground objects makes a distinction with other atural objects, refer to Fig. 3, it illustrates each type of ground objects to different ripples
In the schematic diagram of Duan Guangpu reflectivity, figure, wave band 2 represents blue wave band, and wave band 3 represents green light band, and wave band 4 represents feux rouges
Wave band, as seen from the figure, indigo plant top building blue wave band reflectivity apparently higher than the reflectivity of green light band, and other atural objects
Type then remain basically stable either green light band reflectivity be higher than blue wave band reflectivity, so, if calculate blue wave band
Reflectivity subtract the reflectivity of green light band, building corresponding numerical value in indigo plant top is larger positive number, and other types of ground objects
Corresponding numerical value is then negative or the positive number close to zero, accordingly can extract indigo plant top building;Using same reason
By, it is red top building and cement top building (including exposed soil) red spectral band reflectivity apparently higher than green light band reflection
Rate, and other types of ground objects are then reflectivity of the reflectivity higher than red spectral band of green light band, so, if calculating feux rouges ripple
The reflectivity of section subtracts the reflectivity of green light band, and the corresponding numerical value of red top building and cement top building (including exposed soil) is
Larger positive number, and the corresponding numerical value of other types of ground objects is then negative, accordingly can build red top building and cement top
Thing (including exposed soil) is extracted.Wherein it is possible to after blue top building is extracted, then the reflectivity based on blue wave band and green
The reflectivity of optical band extracts exposed soil, and so deduction exposed soil can extract more accurate red top building and cement top is built
Build thing.
It can ignore because exposed soil real area is less, in order to simplify calculating, the embodiment of the present invention uses and includes exposed soil
Objects extraction method is illustrative, and those skilled in the art can change implementation on the basis of the above description, enter one
Step is extracted and deducted after exposed soil, and to extract the corresponding region of more accurate type of ground objects, it is also in the protection model of the application
Within enclosing.
By taking the Objects extraction containing exposed soil as an example, the extracting index builds subelement and passes through above-mentioned calculating, you can according to each
Type of ground objects builds can make a distinction atural object to be extracted and other atural objects to the difference of different-waveband spectral reflectivity
Objects extraction index, if for example, type of ground objects to be extracted is blue top building, the extracting index builds subelement, can
With the difference according to corresponding first reflection differences of blue top building the first reflection differences corresponding with other types of ground objects, build
Below for the Objects extraction index of blue top building, wherein, first reflection differences refer to the anti-of blue wave band spectrum
Penetrate rate and the difference of the reflectivity to green light band spectrum:
In formula, NDBIB2-B3Represent the Objects extraction index for blue top building, OLI2Represent to blue wave band spectrum
Reflectivity, OLI3Represent the reflectivity to green light band spectrum.
And for example, if type of ground objects to be extracted for it is red top building and cement top building (be not easy to distinguish, can be in the lump
Extract), then the extracting index builds subelement, can be according to corresponding second reflection of red top building and cement top building
The difference of rate difference the second reflection differences corresponding with other types of ground objects, builds and is built below for red top building and cement top
The Objects extraction index of thing, wherein, second reflection differences refer to the reflectivity of red spectral band spectrum and to green light band
The difference of the reflectivity of spectrum:
In formula, NDBIB4-B3Represent the Objects extraction index for red top building and cement top building, OLI4Expression pair
The reflectivity of red spectral band spectrum, OLI3Represent the reflectivity to green light band spectrum.
Using the specific Objects extraction index of above-mentioned two, it can further amplify atural object to be extracted and other types of ground objects
The difference of the corresponding index, so as to help accurately to come out Objects extraction to be extracted in subsequent treatment, specific real
Shi Shi, can also subtract an adjusting parameter to above-mentioned formula, by larger positive number and less positive number be adjusted to positive number with
Negative, to reduce the noise produced during follow-up binary conversion treatment or error.
Accordingly, the exponential quantity of the Objects extraction index of each pixel is being carried out two by the binary conversion treatment subelement
Value is handled, and after being split according to binaryzation result to the remotely-sensed data, you can according to described for blue top building
The binary conversion treatment result of the exponential quantity of the Objects extraction index of thing extracts blue top building, is built according to described for red top
The binary conversion treatment result for building the exponential quantity of the Objects extraction index of thing and cement top building extracts red top building and water
Mud top building.
Wherein, the population coefficient refers to specify the size of population in type of ground objects unit area, can be existing by statistics
There are more accurate data to determine, in order to ensure the accuracy of the population coefficient calculating, the implementation provided in the present invention
In example, the population coefficient determination module 3, including:
Computing unit is returned, for the sample data according to the region of clear and definite type of ground objects composition and the size of population,
The corresponding population coefficient of each type of ground objects is calculated using regression algorithm.
For example, the recurrence computing unit can obtain the sample data in multiple regions similar to target area, each
Sample data includes clear and definite type of ground objects composition situation (real area for including each type of ground objects) and population in sample areas
Sum, based on above-mentioned sample data, can using the population coefficient of each type of ground objects as independent variable, using population as dependent variable,
Regression model is set up, sample data is then inputted into the regression model, determines that each type of ground objects is corresponding by data fitting
Population coefficient.The corresponding population coefficient of more accurate each type of ground objects can be obtained using above-mentioned regression algorithm, so as to help
The more accurately density of population is obtained in final calculate.
The grid partition module 4, for target area to be divided into multiple grid, the division of the grid can basis
The height of actual demand and the remotely-sensed data resolution ratio is flexibly set, such as target area can be divided into multiple ten meters
Grid, hundred meters of grid or km grid etc., it is within the protection domain of the application.
Situation and population coefficient are constituted according to the type of ground objects of the target area, the density of population determining module 5 is
Can calculate each grid of determination (grid be target area is divided obtained by, the type of ground objects composition situation of target area is determined,
The type of ground objects composition situation of corresponding each grid is also determined that) the density of population, circular can be according to reality
Demand is flexibly set, in one embodiment that the application is provided, the density of population determining module 5, including:
Density of population determining unit, the density of population for calculating each grid according to following mathematical algorithm:
Wherein,The corresponding density of population of i-th of grid is represented, j numbers for different types of ground objects, ajRepresent jth kind
The corresponding population coefficient of type of ground objects, XjFor accounting of the area in the grid of jth kind type of ground objects, n represents the grid
The quantity of middle type of ground objects.
Due to grid be target area is divided obtained by, the density of population of each grid is determined, then target area
The population dispersal (i.e. population spatial distribution) in domain is also determined that.
Illustrated based on above example, first embodiment of the invention by target area by being divided into multiple grid, then
The density of population in each grid is calculated respectively in units of grid, so as to calculate people more specific in target area
Mouth Density Distribution situation is more accurate compared to prior art.On the other hand, can accurately it be determined based on remotely-sensed data
Type of ground objects composition in target area, based on above-mentioned type of ground objects composition, can accurately determine different types of ground objects
Corresponding population coefficient, so as to ensure that the density of population finally calculated has the higher degree of accuracy.To sum up, based on the present invention first
Embodiment can more accurate, accurately determine the population dispersal in target area, so as to be country and place
Macro adjustments and controls, urban development planning provide data supporting, and are that enterprises and institutions and the addressing of entrepreneur and industrial pattern are carried
Supported for data.
In order to more intuitively show the population dispersal, in one embodiment that the application is provided, institute
The density of population analysis system based on high score satellite remote sensing date combination type of ground objects is stated, in addition to:
Density of population distribution map generation module, will be each described for the mapping relations according to the density of population and different colours
Color corresponding with the grid density of population is filled in the corresponding position of grid, is distributed with the density of population for drawing the target area
Figure.
As the change embodiment of above-described embodiment, cromogram can be replaced to characterize the people of target area using gray-scale map
Mouth density profile, as shown in figure 4, it is some region of population dispersal effect provided in an embodiment of the present invention
In figure, figure, color is whiter to represent that the density of population is bigger, as seen from the figure, compared to existing simple and crude use administrative division meter
The mode of population dispersal is calculated and characterized, using mode provided in an embodiment of the present invention, can more accurately be determined
The population dispersal of target area.
It is considered that type of ground objects is the key factor for reflecting population distribution, but only analyze the density of population with type of ground objects
During the distribution situation of distribution situation, it is difficult to distinguish the density of population difference between identical type of ground objects.And study and show, night lamp
There is height correlation with the density of population in light data.Therefore, it is described to be based on high score satellite in one embodiment that the application is provided
The density of population analysis system of remotely-sensed data combination type of ground objects, in addition to:
First density of population optimization module, for the corresponding relation based on nighttime light intensity and the density of population, according to night
Between the density of population of each grid that calculates the density of population determining unit of light remotely-sensed data optimize, with excellent
Change the population dispersal of the target area.
Wherein, the nighttime light intensity can be obtained from the corresponding night lights remotely-sensed data in the target area, and
Night lights remotely-sensed data can by with stare full-color camera or stare multispectral camera remote sensing satellite gather obtain, example
No. 4 satellites of high score as China launches can gather night lights remotely-sensed data at night, according to the night of the target area of collection
Between light remotely-sensed data and the grid partition to target area, you can it is determined that each corresponding nighttime light intensity of grid and mesh
The average intensity of light in region is marked, the population for each grid that can be calculated accordingly the density of population determining unit is close
Degree is optimized.
Specifically, in one embodiment that the application is provided, first density of population optimization module, including:
First density of population optimizes unit, for being carried out according to following mathematical algorithm to the density of population of grid each described
Optimization:
Wherein, PiThe corresponding density of population of i-th of grid obtained after optimization is represented,Represent that the density of population is determined
Unit calculates the corresponding density of population of i-th of grid obtained;LjThe corresponding intensity of light of j-th of grid is represented,Represent described
The average intensity of light of target area;PlRepresent the size of population that unit light is represented;S is regulation coefficient.
Those skilled in the art can carry out various reasonable change, tool to specific mathematical algorithm based on above-described embodiment explanation
Body is repeated no more, and it all should be within the protection domain of the application.
The density of population of the grid is optimized by using nighttime light data, can be by by night lights number
Characterized according to by the density of population difference between identical type of ground objects, so that the density of population calculated is more accurate.
It is easily understood that the density of population of each grid for calculating the density of population determining unit is entered
The situation of row optimization, the density of population distribution map generation module can utilize the density of population of each grid after optimization
Draw the density of population distribution map of the target area.
It can effectively reflect the density of population difference in cities and towns and rural area in view of urban heat island strength, and actually cities and towns
Differed greatly with the density of population in rural area, for not only including rural area and including the region in cities and towns again, if type of ground objects not area
Point cities and towns and rural area, error is there may be using the corresponding population coefficient of the counted each type of ground objects of above-mentioned regression algorithm, and then
Cause the population dispersal finally determined there is also error, therefore, in one embodiment that the application is provided, may be used also
To be optimized using urban heat island strength to the density of population of each grid, the type of ground objects is used including city-building
Ground, urban residents' land used and non-residential areas;The density of population computing system, in addition to:
Second density of population optimization module, for the corresponding relation based on urban heat island strength and the density of population, according to red
The density of population for each grid that outer remotely-sensed data is calculated the density of population determining unit is optimized, to optimize
State the population dispersal of target area.
Wherein, the urban heat island strength is calculated according to the corresponding infrared remote sensing data in the target area and determined, and red
Outer remotely-sensed data can be gathered by the remote sensing satellite with full spectral coverage imager and obtained, and the high score 5 of such as China's transmitting is defended
Star carries full spectral coverage imager, can obtain that earth's surface 20m is visible near-infrared and 40m in long infrared remote sensing data, according to collection
Infrared remote sensing data, you can each pixel correspondence in the infrared remote sensing data is extrapolated by land surface temperature inversion technique
Surface temperature, and urban heat island strength is defined as the temperature gap in downtown area temperature and suburb, for characterizing due to city
Urban area temperature caused by city's structure is higher than the degree of suburb temperature, therefore can calculate institute according to the surface temperature
The corresponding urban heat island strength of each pixel in infrared remote sensing data is stated, so that it is determined that the corresponding urban heat island strength of each grid
The average urban heat island strength of the urban heat island strength of all pixels in grid (can be average) and target area, according to
The density of population of this each grid that can be calculated the density of population determining unit is optimized.
Specifically, in one embodiment that the application is provided, second density of population optimization module, including:
Second density of population optimizes unit, for being carried out according to following mathematical algorithm to the density of population of grid each described
Optimization:
Wherein, PiThe corresponding density of population of i-th of grid obtained after optimization is represented,Represent that i-th of grid is corresponding
Belong to the density of population of Urban Construction Land_use,Represent the corresponding density of population for belonging to urban residents' land used of i-th of grid;Ii
The corresponding urban heat island strength of i-th of grid is represented,Represent average urban heat island strength;PIRepresent unit urban heat island strength
The size of population of representative;WiuThe corresponding population coefficient of Urban Construction Land_use, WirRepresent the corresponding population system of urban residents' land used
Number;AiuRepresent accounting of the area of Urban Construction Land_use in i-th of grid, AirRepresent the area of urban residents' land used i-th
Accounting in individual grid.
Wherein, land surface temperature inversion technique is used to calculate surface temperature to be existing ripe based on infrared remote sensing data
Technology, here is omitted, is applied to the whole implementation scheme constituted in the embodiment of the present invention, the protection in the application
Within the scope of.The corresponding population coefficient of the Urban Construction Land_use and the corresponding population coefficient of urban residents' land used may be referred to
The explanation change for stating embodiment calculates determination by the population coefficient determination module 3.
It is easily understood that the density of population of each grid for calculating the density of population determining unit is entered
The situation of row optimization, the density of population distribution map generation module can utilize the density of population of each grid after optimization
Draw the density of population distribution map of the target area.Fig. 5 is refer to, a certain region provided in an embodiment of the present invention is through Urban Thermal
In the population dispersal design sketch of island strength optimization, figure, the more white corresponding density of population of color is higher, as seen from the figure,
After the strength optimization of heat island city, the density of population distribution map of acquisition is more accurate, accurate.
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means to combine specific features, structure, material or the spy that the embodiment or example are described
Point is contained at least one embodiment of the present invention or example.In this manual, to the schematic representation of above-mentioned term not
Identical embodiment or example must be directed to.Moreover, specific features, structure, material or the feature of description can be with office
Combined in an appropriate manner in one or more embodiments or example.In addition, in the case of not conflicting, the skill of this area
Art personnel can be tied the not be the same as Example or the feature of example and non-be the same as Example or example described in this specification
Close and combine.
It should be noted that the flow chart and block diagram in accompanying drawing show according to the present invention multiple embodiments system,
Architectural framework in the cards, function and the operation of method and computer program product.At this point, in flow chart or block diagram
Each square frame can represent a part for a module, program segment or code, the part bag of the module, program segment or code
Containing one or more executable instructions for being used to realize defined logic function.It should also be noted that in some realities as replacement
In existing, the function of being marked in square frame can also be with different from the order marked in accompanying drawing generation.For example, two continuous sides
Frame can essentially be performed substantially in parallel, and they can also be performed in the opposite order sometimes, and this is according to involved function
It is fixed.It is also noted that the group of each square frame in block diagram and/or flow chart and the square frame in block diagram and/or flow chart
Close, can be realized with the special hardware based system of defined function or action is performed, or specialized hardware can be used
Combination with computer instruction is realized.
The analysis of the density of population based on the high score satellite remote sensing date combination type of ground objects system that the embodiment of the present invention is provided
System can be computer program product, including store the computer-readable recording medium of program code, described program code bag
The instruction included can be used for performing the method described in previous methods embodiment, implements and can be found in embodiment of the method, herein not
Repeat again.
It is apparent to those skilled in the art that, for convenience and simplicity of description, the system of foregoing description,
The specific work process of system and unit, may be referred to the corresponding process in preceding method embodiment, will not be repeated here.
, can be with several embodiments provided herein, it should be understood that disclosed system, system and method
Realize by another way.System embodiment described above is only schematical, for example, the division of the unit,
It is only a kind of division of logic function, there can be other dividing mode when actually realizing, in another example, multiple units or component can
To combine or be desirably integrated into another system, or some features can be ignored, or not perform.It is another, it is shown or beg for
The coupling each other of opinion or direct-coupling or communication connection can be by some communication interfaces, system or unit it is indirect
Coupling is communicated to connect, and can be electrical, machinery or other forms.
The unit illustrated as separating component can be or may not be it is physically separate, it is aobvious as unit
The part shown can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be selected to realize the mesh of this embodiment scheme according to the actual needs
's.
In addition, each functional unit in each embodiment of the invention can be integrated in a processing unit, can also
That unit is individually physically present, can also two or more units it is integrated in a unit.
If the function is realized using in the form of SFU software functional unit and is used as independent production marketing or in use, can be with
It is stored in a computer read/write memory medium.Understood based on such, technical scheme is substantially in other words
The part contributed to prior art or the part of the technical scheme can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, including some instructions are to cause a computer equipment (can be individual
People's computer, server, or network equipment etc.) perform all or part of step of each of the invention embodiment methods described.
And foregoing storage medium includes:USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.
Finally it should be noted that:Various embodiments above is merely illustrative of the technical solution of the present invention, rather than its limitations;To the greatest extent
The present invention is described in detail with reference to foregoing embodiments for pipe, it will be understood by those within the art that:Its according to
The technical scheme described in foregoing embodiments can so be modified, or which part or all technical characteristic are entered
Row equivalent substitution;And these modifications or replacement, the essence of appropriate technical solution is departed from various embodiments of the present invention technology
The scope of scheme, it all should cover among the claim of the present invention and the scope of specification.
Claims (10)
1. a kind of density of population analysis system based on high score satellite remote sensing date combination type of ground objects, it is characterised in that including:
Remotely-sensed data acquisition module, type of ground objects determining module, population coefficient determination module, grid partition module and the density of population are determined
Module;Wherein,
The remotely-sensed data acquisition module, for obtaining the corresponding remotely-sensed data in target area;
The type of ground objects determining module, for determining that the type of ground objects of the target area is constituted according to the remotely-sensed data;
The population coefficient determination module, for determining the corresponding population coefficient of difference type of ground objects, institute in the target area
State ratio of the population coefficient for the size of population and unit area;
The grid partition module, for target area to be divided into multiple grid;
The density of population determining module, for according to the type of ground objects of each grid composition and the population coefficient, meter
The density of population of each grid is calculated, to determine the population dispersal of the target area.
2. the density of population analysis system according to claim 1 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, the type of ground objects determining module, including:
Radar data type of ground objects determining unit, for the reflection based on different types of ground objects to radar signal and scattering properties,
The radar remote sensing data obtained according to the remotely-sensed data acquisition module determine the type of ground objects composition of the target area.
3. the density of population analysis system according to claim 1 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, the type of ground objects determining module, including:
Multispectral data type of ground objects determining unit, for based on difference of the different types of ground objects to different-waveband spectral reflectivity
Different, the Multi-spectral Remote Sensing Data obtained according to the remotely-sensed data acquisition module determines the type of ground objects group of the target area
Into.
4. the density of population analysis system according to claim 1 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, the population coefficient determination module, including:
Computing unit is returned, for the sample data according to the region of clear and definite type of ground objects composition and the size of population, is used
Regression algorithm calculates the corresponding population coefficient of each type of ground objects.
5. the density of population analysis system according to claim 1 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, the density of population determining module, including:
Density of population determining unit, the density of population for calculating each grid according to following mathematical algorithm:
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
Wherein,The corresponding density of population of i-th of grid is represented, j numbers for different types of ground objects, ajRepresent jth kind ground species
The corresponding population coefficient of type, XjFor accounting of the area in the grid of jth kind type of ground objects, n represents atural object in the grid
The quantity of type.
6. the density of population analysis system according to claim 5 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, also including:
First density of population optimization module, for the corresponding relation based on nighttime light intensity and the density of population, according to night lamp
The density of population for each grid that light remotely-sensed data is calculated the density of population determining unit is optimized, to optimize
State the population dispersal of target area.
7. the density of population analysis system according to claim 6 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, first density of population optimization module, including:
First density of population optimizes unit, excellent for being carried out according to following mathematical algorithm to the density of population of grid each described
Change:
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>i</mi>
<mo>-</mo>
<mover>
<mi>L</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>*</mo>
<msub>
<mi>P</mi>
<mi>l</mi>
</msub>
<mo>*</mo>
<mi>s</mi>
</mrow>
Wherein, PiThe corresponding density of population of i-th of grid obtained after optimization is represented,Represent the density of population determining unit
Calculate the corresponding density of population of i-th of grid obtained;LjThe corresponding intensity of light of j-th of grid is represented,Represent the target
The average intensity of light in region;PlRepresent the size of population that unit light is represented;S is regulation coefficient.
8. the density of population analysis system according to claim 5 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, the type of ground objects includes Urban Construction Land_use and urban residents' land used;
The density of population computing system, in addition to:
Second density of population optimization module, for the corresponding relation based on urban heat island strength and the density of population, according to infrared distant
The density of population for each grid that sense data are calculated the density of population determining unit is optimized, to optimize the mesh
Mark the population dispersal in region.
9. the density of population analysis system according to claim 8 based on high score satellite remote sensing date combination type of ground objects,
Characterized in that, second density of population optimization module, including:
Second density of population optimizes unit, excellent for being carried out according to following mathematical algorithm to the density of population of grid each described
Change:
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>+</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
<mo>*</mo>
<msub>
<mi>P</mi>
<mi>I</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mi>u</mi>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>u</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
Wherein, PiThe corresponding density of population of i-th of grid obtained after optimization is represented,Represent that i-th of grid is corresponding and belong to city
The density of population of town construction land,Represent the corresponding density of population for belonging to urban residents' land used of i-th of grid;IiRepresent the
The corresponding urban heat island strength of i grid,Represent average urban heat island strength;PIRepresent what unit urban heat island strength was represented
The size of population;WiuThe corresponding population coefficient of Urban Construction Land_use, WirRepresent the corresponding population coefficient of urban residents' land used;AiuTable
Show accounting of the area of Urban Construction Land_use in i-th of grid, AirRepresent the area of urban residents' land used in i-th of grid
Accounting.
10. the density of population based on high score satellite remote sensing date combination type of ground objects according to claim any one of 1-9
Analysis system, it is characterised in that also include:
Density of population distribution map generation module, for the mapping relations according to the density of population and different colours, by each grid
Color corresponding with the grid density of population is filled in corresponding position, to draw the density of population distribution map of the target area.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710370420.2A CN107239756A (en) | 2017-05-23 | 2017-05-23 | Density of population analysis system based on high score satellite remote sensing date combination type of ground objects |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710370420.2A CN107239756A (en) | 2017-05-23 | 2017-05-23 | Density of population analysis system based on high score satellite remote sensing date combination type of ground objects |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107239756A true CN107239756A (en) | 2017-10-10 |
Family
ID=59985126
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710370420.2A Pending CN107239756A (en) | 2017-05-23 | 2017-05-23 | Density of population analysis system based on high score satellite remote sensing date combination type of ground objects |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107239756A (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109543592A (en) * | 2018-11-19 | 2019-03-29 | 北京英视睿达科技有限公司 | The method and device of atmosphere pollution hot spot grid is determined based on remote sensing atural object |
CN110309381A (en) * | 2019-07-04 | 2019-10-08 | 北京光启元数字科技有限公司 | A kind of display methods of demographic data, device and equipment |
CN110390277A (en) * | 2019-07-04 | 2019-10-29 | 中科卫星应用德清研究院 | Complex Underlying Surface identifying water boy method and black and odorous water prediction technique |
CN110570470A (en) * | 2019-09-06 | 2019-12-13 | 南京大学 | ghost community identification and housing vacancy rate estimation method based on multi-source remote sensing data |
CN110807724A (en) * | 2019-10-26 | 2020-02-18 | 福建省伟志地理信息科学研究院 | Population flow monitoring system and method based on satellite remote sensing |
CN112215059A (en) * | 2020-08-26 | 2021-01-12 | 厦门大学 | Urban village identification and population estimation method and system based on deep learning and computer readable storage medium |
CN112818747A (en) * | 2020-12-31 | 2021-05-18 | 上海应用技术大学 | Urban characteristic neighborhood population density estimation method and system method based on spatial big data |
CN113849694A (en) * | 2021-09-22 | 2021-12-28 | 上海妙一生物科技有限公司 | Analysis and device of medicine registration data |
CN114064831A (en) * | 2021-10-27 | 2022-02-18 | 北京融信数联科技有限公司 | Data rasterization display method and system and storage medium |
CN116108124A (en) * | 2023-04-13 | 2023-05-12 | 中国铁塔股份有限公司 | Land type similarity determination method, system, equipment and medium based on GIS |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794164A (en) * | 2015-03-26 | 2015-07-22 | 华南理工大学 | Method for recognizing settlement parking spaces meeting social parking requirement on basis of open source data |
CN106455058A (en) * | 2016-12-02 | 2017-02-22 | 中国联合网络通信集团有限公司 | Method and device for determining population distribution condition |
CN106650618A (en) * | 2016-11-15 | 2017-05-10 | 中山大学 | Random forest model-based population data spatialization method |
-
2017
- 2017-05-23 CN CN201710370420.2A patent/CN107239756A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794164A (en) * | 2015-03-26 | 2015-07-22 | 华南理工大学 | Method for recognizing settlement parking spaces meeting social parking requirement on basis of open source data |
CN106650618A (en) * | 2016-11-15 | 2017-05-10 | 中山大学 | Random forest model-based population data spatialization method |
CN106455058A (en) * | 2016-12-02 | 2017-02-22 | 中国联合网络通信集团有限公司 | Method and device for determining population distribution condition |
Non-Patent Citations (10)
Title |
---|
唐小龙: "高分辨率城市人口密度模拟——以重庆市北碚城区为例", 《中国优秀硕士学位论文全文数据库》 * |
张珣等: "一种尺度效应指数修正的格网示意地图制图算法", 《武汉大学学报》 * |
徐建刚等: "城市居住人口密度估算模型的研究", 《环境遥感》 * |
杜培军等: "《遥感科学与进展》", 28 February 2007, 中国矿业大学出版社 * |
毛莹莹: "城市人口数据空间化研究——以福州市中心城区为例", 《中国优秀硕士学位论文全文数据库》 * |
满其霞: "激光雷达和高光谱数据融合的城市土地利用分类方法研究", 《中国博士学位论文全文数据库》 * |
范文义等: "《3S理论与技术》", 31 March 2003, 东北林业大学出版社 * |
袁飞等: "《环境变化与水文过程模拟》", 31 December 2012, 河海大学出版社 * |
贺华翔等: "松辽流域人口信息空间分布规律研究", 《中国人口•资源与环境》 * |
鲁楠等: "顾及城乡差异的大区域人口密度估算——以山东省为例", 《测绘学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109543592A (en) * | 2018-11-19 | 2019-03-29 | 北京英视睿达科技有限公司 | The method and device of atmosphere pollution hot spot grid is determined based on remote sensing atural object |
CN110309381A (en) * | 2019-07-04 | 2019-10-08 | 北京光启元数字科技有限公司 | A kind of display methods of demographic data, device and equipment |
CN110390277A (en) * | 2019-07-04 | 2019-10-29 | 中科卫星应用德清研究院 | Complex Underlying Surface identifying water boy method and black and odorous water prediction technique |
CN110570470A (en) * | 2019-09-06 | 2019-12-13 | 南京大学 | ghost community identification and housing vacancy rate estimation method based on multi-source remote sensing data |
CN110807724A (en) * | 2019-10-26 | 2020-02-18 | 福建省伟志地理信息科学研究院 | Population flow monitoring system and method based on satellite remote sensing |
CN110807724B (en) * | 2019-10-26 | 2022-09-13 | 福建省伟志地理信息科学研究院 | Population flow monitoring system and method based on satellite remote sensing |
CN112215059A (en) * | 2020-08-26 | 2021-01-12 | 厦门大学 | Urban village identification and population estimation method and system based on deep learning and computer readable storage medium |
CN112215059B (en) * | 2020-08-26 | 2023-10-27 | 厦门大学 | Deep learning-based urban village identification and population estimation method, system and computer-readable storage medium |
CN112818747A (en) * | 2020-12-31 | 2021-05-18 | 上海应用技术大学 | Urban characteristic neighborhood population density estimation method and system method based on spatial big data |
CN113849694A (en) * | 2021-09-22 | 2021-12-28 | 上海妙一生物科技有限公司 | Analysis and device of medicine registration data |
CN114064831A (en) * | 2021-10-27 | 2022-02-18 | 北京融信数联科技有限公司 | Data rasterization display method and system and storage medium |
CN116108124A (en) * | 2023-04-13 | 2023-05-12 | 中国铁塔股份有限公司 | Land type similarity determination method, system, equipment and medium based on GIS |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107239756A (en) | Density of population analysis system based on high score satellite remote sensing date combination type of ground objects | |
Deng et al. | Land use/land cover classification using time series Landsat 8 images in a heavily urbanized area | |
Lu et al. | Regional mapping of human settlements in southeastern China with multisensor remotely sensed data | |
Pu et al. | Spectral mixture analysis for mapping abundance of urban surface components from the Terra/ASTER data | |
CN107229913A (en) | Density of population analysis system based on high score satellite remote sensing date combination building height | |
Gao et al. | Using the random forest model and validated MODIS with the field spectrometer measurement promote the accuracy of estimating aboveground biomass and coverage of alpine grasslands on the Qinghai-Tibetan Plateau | |
CN107966116B (en) | A kind of remote-sensing monitoring method and system of Monitoring of Paddy Rice Plant Area | |
CN105608293A (en) | Forest aboveground biomass inversion method and system fused with spectrum and texture features | |
Yu et al. | Investigating impacts of urban morphology on spatio-temporal variations of solar radiation with airborne LIDAR data and a solar flux model: a case study of downtown Houston | |
CN103942420A (en) | Rapid estimation method for solar energy in construction size | |
CN104142142A (en) | Method for estimating global vegetation coverage | |
Chen et al. | Future “local climate zone” spatial change simulation in Greater Bay Area under the shared socioeconomic pathways and ecological control line | |
CN104408258A (en) | Large-scale soil organic carbon spatial distribution simulation method involving environmental factors | |
Han et al. | Extraction of multilayer vegetation coverage using airborne LiDAR discrete points with intensity information in urban areas: A case study in Nanjing City, China | |
CN110990511B (en) | Local climate zone classification method considering two-dimensional and three-dimensional forms of city | |
CN107491758A (en) | Yangtze river basin water body information and its space encoding method | |
Raduła et al. | Effect of environmental gradients, habitat continuity and spatial structure on vascular plant species richness in semi-natural grasslands | |
Hwang et al. | Expanding vegetated areas by human activities and strengthening vegetation growth concurrently explain the greening of Seoul | |
CN107239886A (en) | GDP density Analysis Systems based on high score satellite remote sensing date | |
CN107121681A (en) | Residential area extraction system based on high score satellite remote sensing date | |
Kang et al. | Construction of multidimensional features to identify tea plantations using multisource remote sensing data: A case study of Hangzhou city, China | |
CN108596425A (en) | A kind of polynucleated city development boundary demarcation method | |
CN116665073A (en) | Corn yield remote sensing estimation method based on multi-source data | |
Tokula et al. | The impact of urban expansion on agricultural land and crop output in Ankpa, Kogi State, Nigeria | |
Sun et al. | Evaluating the street greening with the multiview data fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171010 |