CN110440722A - A kind of building index construction method suitable for the infrared data without in - Google Patents
A kind of building index construction method suitable for the infrared data without in Download PDFInfo
- Publication number
- CN110440722A CN110440722A CN201910766452.3A CN201910766452A CN110440722A CN 110440722 A CN110440722 A CN 110440722A CN 201910766452 A CN201910766452 A CN 201910766452A CN 110440722 A CN110440722 A CN 110440722A
- Authority
- CN
- China
- Prior art keywords
- infrared
- building
- data
- information
- index
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/28—Measuring arrangements characterised by the use of optical techniques for measuring areas
- G01B11/285—Measuring arrangements characterised by the use of optical techniques for measuring areas using photoelectric detection means
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
Abstract
The present invention relates to a kind of building index construction methods suitable for the infrared data without in, and first based on the data that same phase has middle infrared information, the spectral information of construction area, infrared-near-infrared relational model in foundation are read using random point;Followed by infrared-near-infrared relational model in established, common building exponential model is modified, establishes and does not need the new building exponential model that middle infrared information participates in;Then using the data of the infrared information without in as object, the building index of the extrinsic information data without in is calculated based on revised new building exponential model;Step-length is finally set, the architecture information extraction accuracy under different threshold values is calculated, optimal threshold is determined according to this, if architecture information extraction accuracy is up to standard, then it represents that the modified building index of institute has validity.The present invention can extract building site in the spectral region for only arriving near-infrared, thus make these data using more flexible, extensive.
Description
Technical field
The present invention relates to remote sensing technology field, especially a kind of building index construction method suitable for the infrared data without in.
Background technique
With the continuous quickening of urbanization process, the form of land use had occurred significantly in a very long time in past
Variation.Wherein, building site is as a kind of important land use pattern, gradually instead of with the lands used class such as vegetation, water body
Natural landscape based on type.Just because of the presence of construction area, the natural environment and ground mulching type on city and its periphery
Occur from natural terrain to be converted into low-density building, be converted into that high density building area etc. is a series of to be changed from low-density building
Become, so that environment ability to bear weakens;The appearance of large area building, has blocked the infitration process of ground water, has been to China
Ecological environment to global urban all brings negative effect.The increase of construction area can not only reflect the development and expansion in city
, while city surface structure will be caused to change, and then change regional climate, influence urban ecological system service function.Therefore,
The area of accurate estimation building area has a significance for the development plan, Ecology Environment Evaluation and resource management etc. in city.
Tradition is unable to satisfy Development of Urbanization with rapid changepl. never-ending changes and improvements to the monitoring method of construction area in land use already.
As the deep and remote sensing technology of global change research due is in the development in China, change with satellite image monitoring building site
Means also emerge one after another, and this greatly reduces the testing costs of building site.Meanwhile the features such as remote sensing large scale, multidate
Also strong approach is provided for the continuous observation of building site.Therefore, how to be constructed with remote sensing image and reasonably build index,
So as to more rapidly, calculate to a nicety future city expansion trend, also become the important topic of remote sensing application.
It is currently used extract building site index, generally to use in infrared or even far infrared data.However, for
For many target in hyperspectral remotely sensed image, spectral region generally only arrives near-infrared;In addition to this, China's high score series of satellites, can
For the image that masses are downloaded, spectral region is also to arrive near-infrared.Moreover, non-for many ground monitored on the spot
Imaging Hyperspectral Data and unmanned plane high-spectral data etc., wave-length coverage generally arrives near-infrared and stops.These remote sensing
Although data spectral region is relatively narrow, they have the features such as spatial resolution, spectral resolution is higher.If because light
The limitation of spectral limit and the advantages of these data cannot be made full use of, undoubtedly a kind of huge waste.
Summary of the invention
In view of this, the purpose of the present invention is to propose to a kind of building index construction method suitable for the infrared data without in, it can
To extract building site in the spectral region for only arriving near-infrared, thus make these satellites using more flexible, extensive.
The present invention is realized using following scheme: a kind of building index construction method suitable for the infrared data without in, specifically:
First based on the data that same phase has middle infrared information, believed using the spectrum that random point reads construction area
It ceases, infrared-near-infrared relational model in foundation;
Followed by infrared-near-infrared relational model in established, common building exponential model is modified, is established
The new building exponential model that middle infrared information participates in is not needed;
Then it using the data of the infrared information without in as object, is calculated based on revised new building exponential model and is believed without China and foreign countries
Cease the building index of data;
Step-length is finally set, the architecture information extraction accuracy under different threshold values is calculated, determines optimal threshold according to this, if building
Information extraction precision is up to standard, then it represents that the modified building index of institute has validity.
Further, described based on the data that same phase has middle infrared information, construction area is read using random point
Spectral information, infrared-near-infrared relational model in foundation specifically:
Image data is pre-processed, N number of check post is chosen on the building area in image, counts these check posts
Upper near-infrared and in infrared reflectivity, be fitted regression equation, obtain in infrared-near-infrared relational model are as follows:
Y=1.160x-0.015;
In formula, y is dependent variable, representative simulation go out in infrared reflectivity;X is independent variable, represents the near-infrared of input
Reflectivity.
Further, described to utilize infrared-near-infrared relational model in established, common building exponential model is carried out
Amendment, establish do not need new building exponential model that middle infrared information participates in specifically includes the following steps:
Step S1: IBI index is established:
In formula, NDBI is building index, and SAVI is vegetation index, and MNDWI is water body index;Wherein, NDBI and MNDWI are equal
The infrared participation operation in;
Step S2: infrared in foundation-near-infrared relational model is brought into the formula of step S1 foundation, is not needed
The new building exponential model that middle infrared information participates in.
Finally, setting step-length, calculates the architecture information extraction accuracy under different threshold values, optimal threshold is determined according to this, if building
It is up to standard to build information extraction precision, then it represents that the modified building index of institute has validity.By choosing random point, square is obscured in building
Battle array, verifying extract the precision of building and the validity of the above method.
Compared with prior art, the invention has the following beneficial effects: the present invention can be in the spectral region for only arriving near-infrared
Interior extraction building site, thus make these satellites using more flexible, extensive.Through testing, the present invention can quickly, preparatively
The architecture information in remote sensing image is extracted, and suitable for the data such as the unmanned plane image of external information, ground actual measurement EO-1 hyperion without.
Using this characteristic, by its application range expand to other use in infrared index, also can yet be regarded as and a kind of do well
Method.
Detailed description of the invention
Fig. 1 is that the Landsat 8/OLI threshold value of the embodiment of the present invention chooses result schematic diagram.
Fig. 2 is that the GF-1/WFV threshold value of the embodiment of the present invention chooses result schematic diagram.
Fig. 3 is that the PROBA/CHRIS threshold value of the embodiment of the present invention chooses result schematic diagram.
Fig. 4 is experimental result schematic diagram of the embodiment of the present invention in three kinds of different images.Wherein, (a) is Landsat
8/OLI extract as a result, (b) be GF-1/WFV extract as a result, (c) be PROBA/CHRIS extraction result.
Fig. 5 is IBI and the present embodiment extracts part Architecture contrast schematic diagram.Wherein, (a) is that Landsat 8/OLI is original
Image is (b) that IBI is extracted as a result, (c) being this method extraction result.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and embodiments.
It is noted that described further below be all exemplary, it is intended to provide further instruction to the application.Unless another
It indicates, all technical and scientific terms used herein has usual with the application person of an ordinary skill in the technical field
The identical meanings of understanding.
It should be noted that term used herein above is merely to describe specific embodiment, and be not intended to restricted root
According to the illustrative embodiments of the application.As used herein, unless the context clearly indicates otherwise, otherwise singular
Also it is intended to include plural form, additionally, it should be understood that, when in the present specification using term "comprising" and/or " packet
Include " when, indicate existing characteristics, step, operation, device, component and/or their combination.
Present embodiments provide a kind of building index construction method suitable for the infrared data without in, the experiment number of the present embodiment
According to for multispectral data Landsat 8/OLI, GF-1/WFV and high-spectral data PROBA/CHRIS.Data are not limited to be mentioned
The embodiment of confession.Method specifically:
First based on the data that same phase has middle infrared information, believed using the spectrum that random point reads construction area
It ceases, infrared-near-infrared relational model in foundation;
Followed by infrared-near-infrared relational model in established, common building exponential model is modified, is established
The new building exponential model that middle infrared information participates in is not needed;
Then it using the data of the infrared information without in as object, is calculated based on revised new building exponential model and is believed without China and foreign countries
Cease the building index of data.
In the present embodiment, described based on the data that same phase has middle infrared information, it is read and is built using random point
The spectral information in region, infrared-near-infrared relational model in foundation specifically:
(including: atmospheric correction, geometric correction, cutting), Zhi Houyong are pre-processed to Landsat 8/OLI image data
The function of random point is generated in ArcGIS software, and is aided with visual interpretation, and 300 verifyings are chosen on the building area in image
Point counts the reflectivity that near-infrared is infrared on these check posts, imported into SPSS, is fitted regression equation, as a result table
It is bright, near-infrared and in it is infrared between linear relationship is presented, therefore infrared-near-infrared relational model in obtaining are as follows:
Y=1.160x-0.015;
In formula, y is dependent variable, representative simulation go out in infrared reflectivity;X is independent variable, represents the near-infrared of input
Reflectivity.Through examining, fitting precision has reached extremely significant horizontal (P < 0.01).Therefore pass through above formula i.e. available near-infrared generation
Infrared carry out exponent arithmetic in replacing.
In the present embodiment, described to utilize infrared-near-infrared relational model in established, to common building exponential model
Be modified, establish do not need new building exponential model that middle infrared information participates in specifically includes the following steps:
Step S1: IBI index is established:
In formula, NDBI is building index, and SAVI is vegetation index, and MNDWI is water body index;By testing for many cases
Card, the index can preferably extract earth's surface building.But NDBI and MNDWI index require in infrared participation operation;
Step S2: infrared in foundation-near-infrared relational model is brought into the formula of step S1 foundation, is not needed
The new building exponential model that middle infrared information participates in.
Finally, setting step-length, calculates the architecture information extraction accuracy under different threshold values, optimal threshold is determined according to this, if building
It is up to standard to build information extraction precision, then it represents that the modified building index of institute has validity.By choosing random point, square is obscured in building
Battle array, verifying extract the precision of building and the validity of the above method.
In the present embodiment, for after calculating in aforementioned manners as a result, drawing out curve with step length selected threshold
Figure, it is apparent that the position of optimal threshold, that is, select this threshold value to extract building site.Fig. 1 to Fig. 3 is three kinds of differences
The result of image selected threshold.
According to above method selected threshold, it can be clearly seen that: (1) for different images, the range of optimal threshold is different;
(2) for different images, the selection section of step-length is different;(3) for different images, optimum precision and best Kappa coefficient
Threshold value is not also identical.Since response of the partial image to the index is more sensitive, if equally choosing identical step-length with other images,
Although the variation of its extraction accuracy and Kappa coefficient, which is still presented, first rises the variation tendency declined afterwards, its variation range mistake
Greatly, it is unfavorable for statistical research.Therefore choose lesser step-length, so as to it is more intuitive, meticulously reflect its variation tendency, also easily facilitate
Choose optimal threshold.From the point of view of above-mentioned statistical result, three kinds of image optimum threshold values are respectively: 0.034,0.026,0.6856,
Its extraction accuracy is respectively: 96%, 92.67%, 94.33%, Kappa coefficient be respectively: 0.9188,0.8477,0.8863,
It is as shown in Figure 4 to extract result difference.Landsat 8/OLI image is directly extracted with IBI, then the finger with the present embodiment building
Number extracts the amplification of result some areas and is compared, as figure 5 illustrates.
As can be seen that the method for the present embodiment can be extracted more completely under the premise of index chooses reasonable threshold value
Building site out.Using this characteristic, by its application range expand to other use in infrared index, also can yet be regarded as one
The good method of kind.
The above described is only a preferred embodiment of the present invention, being not that the invention has other forms of limitations, appoint
What those skilled in the art changed or be modified as possibly also with the technology contents of the disclosure above equivalent variations etc.
Imitate embodiment.But without departing from the technical solutions of the present invention, according to the technical essence of the invention to above embodiments institute
Any simple modification, equivalent variations and the remodeling made, still fall within the protection scope of technical solution of the present invention.
Claims (3)
1. a kind of building index construction method suitable for the infrared data without in, which is characterized in that
First based on the data that same phase has middle infrared information, the spectral information of construction area is read using random point, is built
Infrared-near-infrared relational model in vertical;
Followed by infrared-near-infrared relational model in established, common building exponential model is modified, foundation is not required to
The new building exponential model for wanting middle infrared information to participate in;
Then it using the data of the infrared information without in as object, is calculated based on revised new building exponential model without China and foreign countries' Information Number
According to building index.
2. a kind of building index construction method suitable for the infrared data without according to claim 1, which is characterized in that institute
It states based on the data that same phase has middle infrared information, it is red in foundation using the spectral information of random point reading construction area
Outside-near-infrared relational model, specifically:
Image data is pre-processed, N number of check post is chosen on the building area in image, is counted close on these check posts
It is infrared and in infrared reflectivity, be fitted regression equation, obtain in infrared-near-infrared relational model are as follows:
Y=1.160x-0.015;
In formula, y is dependent variable, representative simulation go out in infrared reflectivity;X is independent variable, represents the anti-of the near-infrared of input
Penetrate rate.
3. a kind of building index construction method suitable for the infrared data without according to claim 1, which is characterized in that institute
It states and utilizes infrared-near-infrared relational model in established, common building exponential model is modified, establish red in not needing
External information participate in new building exponential model specifically includes the following steps:
Step S1: IBI index is established:
In formula, NDBI is building index, and SAVI is vegetation index, and MNDWI is water body index;Wherein, NDBI and MNDWI are all made of
In infrared participation operation;
Step S2: by infrared in foundation-near-infrared relational model bring into step S1 foundation formula in, do not needed in it is red
The new building exponential model that external information participates in.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910766452.3A CN110440722B (en) | 2019-08-20 | 2019-08-20 | Construction index construction method suitable for medium infrared-free data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910766452.3A CN110440722B (en) | 2019-08-20 | 2019-08-20 | Construction index construction method suitable for medium infrared-free data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110440722A true CN110440722A (en) | 2019-11-12 |
CN110440722B CN110440722B (en) | 2021-07-13 |
Family
ID=68436488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910766452.3A Active CN110440722B (en) | 2019-08-20 | 2019-08-20 | Construction index construction method suitable for medium infrared-free data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110440722B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708307A (en) * | 2012-06-26 | 2012-10-03 | 上海大学 | Vegetation index construction method applied to city |
CN103268587A (en) * | 2013-05-06 | 2013-08-28 | 重庆市勘测院 | Method for obtaining urban building land information by means of imitated building land index |
US20130330481A1 (en) * | 2005-10-05 | 2013-12-12 | Certainteed Corporation | Thin Films With High Near-Infrared Reflectivity Deposited on Building Materials |
CN108647704A (en) * | 2018-04-20 | 2018-10-12 | 北京英视睿达科技有限公司 | The information acquisition method and device merged with related coefficient based on NDBI |
KR101914061B1 (en) * | 2018-03-29 | 2018-11-02 | 대한민국 | The method for analyzing heat island characteristic by Satellite |
CN108896185A (en) * | 2018-05-14 | 2018-11-27 | 河海大学 | Remote Sensing temperature space NO emissions reduction method based on normalization desert index |
CN109598273A (en) * | 2018-12-03 | 2019-04-09 | 中国矿业大学 | A kind of city entity boundary recognition methods of fusion surface temperature and building index |
CN109801304A (en) * | 2019-01-11 | 2019-05-24 | 重庆邮电大学 | A kind of building index construction and construction area rapid extracting method inhibiting bare area |
-
2019
- 2019-08-20 CN CN201910766452.3A patent/CN110440722B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130330481A1 (en) * | 2005-10-05 | 2013-12-12 | Certainteed Corporation | Thin Films With High Near-Infrared Reflectivity Deposited on Building Materials |
CN102708307A (en) * | 2012-06-26 | 2012-10-03 | 上海大学 | Vegetation index construction method applied to city |
CN103268587A (en) * | 2013-05-06 | 2013-08-28 | 重庆市勘测院 | Method for obtaining urban building land information by means of imitated building land index |
KR101914061B1 (en) * | 2018-03-29 | 2018-11-02 | 대한민국 | The method for analyzing heat island characteristic by Satellite |
CN108647704A (en) * | 2018-04-20 | 2018-10-12 | 北京英视睿达科技有限公司 | The information acquisition method and device merged with related coefficient based on NDBI |
CN108896185A (en) * | 2018-05-14 | 2018-11-27 | 河海大学 | Remote Sensing temperature space NO emissions reduction method based on normalization desert index |
CN109598273A (en) * | 2018-12-03 | 2019-04-09 | 中国矿业大学 | A kind of city entity boundary recognition methods of fusion surface temperature and building index |
CN109801304A (en) * | 2019-01-11 | 2019-05-24 | 重庆邮电大学 | A kind of building index construction and construction area rapid extracting method inhibiting bare area |
Non-Patent Citations (2)
Title |
---|
NARAYAN KAYET EL.: "Urban heat island explored by co-relationship between land surface temperature vs multiple vegetation indices", 《SPATIAL INFORMATION RESEARCH》 * |
马红: "一种基于遥感指数的城市建筑用地信息提取新方法", 《城市勘测》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110440722B (en) | 2021-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100390565C (en) | Method for measuring forest by unmanned aerial vehicle aerial photography remote sensing | |
CN107063197B (en) | Reservoir characteristic curve extraction method based on spatial information technology | |
CN107389036A (en) | A kind of large spatial scale vegetation coverage computational methods of combination unmanned plane image | |
CN104063718B (en) | The method with selection remotely-sensed data and sorting algorithm in area reckoning is recognized in crop | |
CN110136194A (en) | Snow Cover measuring method based on satellite-borne multispectral remotely-sensed data | |
Zhu et al. | Foliage clumping index over China's landmass retrieved from the MODIS BRDF parameters product | |
CN105913017A (en) | Corresponding period double high resolution remote sensing image-based forest biomass estimation method | |
CN111046613B (en) | Optimal river channel calculation method based on path tracking and river network extraction method based on multi-temporal remote sensing image | |
CN114091613B (en) | Forest biomass estimation method based on high-score joint networking data | |
CN105004320A (en) | Gaofen satellite data land surface vegetation coverage retrieval method and system | |
CN103424405B (en) | Drought monitoring method based on HJ-1A/1B CCD data | |
CN108984803A (en) | A kind of method and system of crop yield spatialization | |
CN107479065A (en) | A kind of three-dimensional structure of forest gap method for measurement based on laser radar | |
CN114595975A (en) | Unmanned aerial vehicle remote sensing pine wood nematode disease monitoring method based on deep learning model | |
CN108896021B (en) | Method for extracting artificial forest stand structure parameters based on aerial photogrammetry point cloud | |
CN105426881B (en) | Mountain background thermal field model constrained underground heat source daytime remote sensing detection locating method | |
CN110988909A (en) | TLS-based vegetation coverage determination method for sandy land vegetation in alpine and fragile areas | |
CN108195736A (en) | A kind of method of three-dimensional laser point cloud extraction Vegetation canopy clearance rate | |
CN106485718A (en) | One kind overdoes slash recognition methodss and device | |
CN110189043A (en) | It is a kind of based on high score satellite remote sensing date using land resource analysis system | |
CN104463971A (en) | Method for creating green-degree spatial arrangement curve used for evaluating urban landscaping three-dimensional layout | |
Hu et al. | Unmanned aerial vehicle (UAV) remote sensing estimation of wheat chlorophyll in subsidence area of coal mine with high phreatic level | |
CN105894006A (en) | Space-time probability model rice remote sensing recognition method | |
Menon | Remote sensing applications in agriculture and forestry | |
CN105403201B (en) | A kind of remote sensing images atmospheric path radiation acquisition methods based on pixel analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |