CN109376638A - A method for calculating land coverage based on remote sensing images and geographic information system - Google Patents

A method for calculating land coverage based on remote sensing images and geographic information system Download PDF

Info

Publication number
CN109376638A
CN109376638A CN201811197720.6A CN201811197720A CN109376638A CN 109376638 A CN109376638 A CN 109376638A CN 201811197720 A CN201811197720 A CN 201811197720A CN 109376638 A CN109376638 A CN 109376638A
Authority
CN
China
Prior art keywords
remote sensing
cultural
area
land used
land
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201811197720.6A
Other languages
Chinese (zh)
Other versions
CN109376638B (en
Inventor
王树声
朱旭东
徐玉倩
拜晓蒙
张瑶
高元
李小龙
王凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Architecture and Technology
Original Assignee
Xian University of Architecture and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Architecture and Technology filed Critical Xian University of Architecture and Technology
Priority to CN201811197720.6A priority Critical patent/CN109376638B/en
Publication of CN109376638A publication Critical patent/CN109376638A/en
Application granted granted Critical
Publication of CN109376638B publication Critical patent/CN109376638B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/176Urban or other man-made structures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/25Determination of region of interest [ROI] or a volume of interest [VOI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20104Interactive definition of region of interest [ROI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种基于遥感图像和地理信息系统的文地率计算方法,具体步骤包括:在GIS矢量地图中导入文化用地POI数据;将GIS矢量地图与遥感地图配准,识别文化用地Pi周围主干道、次干道所形成的封闭区域;在封闭区域内导入非文化用地类型POI,将图像分割为若干个区域,选择包含Pi的纹理分割区域作为文化用地区域,记为Ci;对得到的每个文化用地区域Ci进行面积计算,记为Si;计算所有文化用地面积和S:求取所有文化用地面积和S与城市建设总用地面积A的比值R:R=S/A;R为计算求得的文地率。本发明通过识别文化用地兴趣点(Point of Interest,POI),基于遥感图像和GIS进行文地率计算,计算结果较精确,计算效率较高。

The invention discloses a method for calculating land ratio based on remote sensing images and a geographic information system. The specific steps include: importing cultural land POI data into a GIS vector map; registering the GIS vector map and the remote sensing map to identify the cultural land P i The enclosed area formed by the surrounding main roads and secondary roads; the non-cultural land type POI is imported into the enclosed area, the image is divided into several areas, and the texture segmentation area containing P i is selected as the cultural land area, denoted as C i ; Calculate the area of each cultural land area C i obtained, denoted as S i ; calculate the area and S of all cultural land: Calculate the ratio of all cultural land area and S to the total urban construction land area A: R: R=S/A; R is the cultural land ratio obtained by calculation. By identifying the cultural land interest point (POI), the invention calculates the cultural land rate based on remote sensing images and GIS, the calculation result is more accurate, and the calculation efficiency is high.

Description

A kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System
Technical field
It is the invention belongs to urban and rural planning field, in particular to a kind of based on remote sensing images and GIS-Geographic Information System (GIS) Literary ground rate calculation method.
Background technique
GIS-Geographic Information System (Geographic Information System, abbreviation GIS) be it is a kind of acquisition, storage, pipe Reason is analyzed, the computer technology system of display and applied geography information, is the general skill of analysis and processing large quantity geography data Art;It mainly includes that data input shows and exports with editor, data management, data manipulation and data;Its conduct acquisition, Important tool, technology and the subject of processing, management and analysis geographical spatial data, get the attention.
Currently, urban culture gradually declines under Rapid-Urbanization background.City is carrying Urban Spirit and culture literaryly It is multi-functional to have both Public Culture Service and cultural continuity innovation and development etc. for urban land.Text ground systems organization in city aims to solve the problem that Cultural continuity and development problem in Modern Urban Planning construction.Literary ground rate is the core index of literary ground systems organization, is characterization The important material space attribute of Urban Spirit and culture needs a kind of literary ground rate calculation method at present.
Summary of the invention
A kind of the purpose of the present invention is to provide text based on remote sensing images and GIS-Geographic Information System rate calculation method, with Solve above-mentioned technical problem.The present invention is by identifying cultural land used point of interest (Point of Interest, POI), base In remote sensing images and GIS, rate is calculated with carrying out text, and calculated result is more accurate, and computational efficiency is higher.
In order to achieve the above objectives, the invention adopts the following technical scheme:
A kind of text ground rate calculation method, specific steps based on remote sensing images and GIS-Geographic Information System include:
Step 1, it treats zoning and carries out the town site area recognition based on index image, obtain urban construction Total land area A;
Step 2, cultural land used POI data is imported in the GIS map vector in region to be calculated, cultural land used is denoted as Pi, i Indicate that the number of cultural land used, value are 1 to n;
Step 3, the distant of the GIS map vector in the region to be calculated of cultural land used POI data and region to be calculated will be imported Feel map registration, identifies cultural land used PiMajor trunk roads and major trunk roads, subsidiary road and subsidiary road and master within the scope of surrounding 100m Arterial highway and subsidiary road are formed by closed area;
Step 4, m non-culture land-use style POI are imported in the closed area that step 3 obtains is denoted as Qm, used with culture Ground PiCentered on, point based on m non-cultural land-use style POI carries out image using the watershed algorithm based on texture gradient Segmentation, divides the image into m+1 region, and selection includes PiTexture Segmentation region as cultural land used region, be denoted as Ci
Step 5, i is successively traversed from 1 to n, to each of obtaining cultural land used region CiAreal calculation is carried out, S is denoted asi
Step 6, the literate land area of institute and S are calculated:SiFor each cultural land used CiArea;
Step 7, seek literate land area and S and the total land area A of urban construction ratio R: R=S/A;R is With calculating the text acquired rate.
Further, it in step 1, when carrying out town site area recognition, selects and normalizes exposed index conduct The indicator of town site is removed the noise information in non-constructive land region by Threshold segmentation, mask process, obtained Urban Construction Land_use information.
Further, in step 1, the specific steps of town site area recognition include:
(1) calculate and normalize exposed index: the wave band that spectrum gray average is directed between 20-120 range is remembered For wave band 1, the wave band that spectrum gray average is directed between 20-180 range is denoted as wave band 2;
Normalize the calculation formula of exposed index are as follows: normalize exposed index=(wave band 2- wave band 1)/(wave band 1+ wave band 2);
(2) given threshold;It is 1 to the pixel assignment for being higher than threshold value in remote sensing image, the pixel assignment lower than threshold value is 0;Being assigned a value of 1 is town site area.
Further, the threshold value in step (2) is set as 140.
Further, further includes: step (3);(3) interference is carried out to the image after step (2) assignment to handle.
Further, cultural land used be divided into cultural spirit mark land used, commemorate land used, religion land used, cultural heritage land used, Cultural facility land used;POI data includes: classification, title and coordinate.
Further, step 3, will import the GIS map vector in the region to be calculated of cultural land used POI data with it is to be calculated The remote sensing map autoregistration in region passes through remote sensing images automatic identification culture land used PiMajor trunk roads within the scope of surrounding 100m with Major trunk roads, subsidiary road and subsidiary road and major trunk roads and subsidiary road are formed by closed area;Major trunk roads are had a lot of social connections for 45m-55m, Subsidiary road is had a lot of social connections for 40m-50m.
Further, the specific steps of map vector and remote sensing map progress autoregistration include:
(1) select Ground Control Points, reselection are selected in the register of main menu in ENVI Image to map carries out geometric precision correction to polar plot and remote sensing figure;
(2) in ENVI by correct after image based on, use Function/Overlays/ in master image window Polar plot is added on the image after correcting by Vector Layers, obtains registration image.
Further, pass through remote sensing images automatic identification culture land used PiMajor trunk roads and trunk within the scope of surrounding 100m The specific steps that road, subsidiary road and subsidiary road and major trunk roads are formed by closed area with subsidiary road include:
(1) block of pixels for reading predetermined position on Raster by the Read method of IRaster interface in GIS, passes through The GetVal method of block of pixels obtains the pixel value in predetermined position, determines cultural land used P on remote sensing imagesiBelong in 100 meters of surrounding In the pixel of road;
(2) directional information for belonging to the pixel of road is calculated;
(3) expansive working is carried out to the pixel for belonging to road;
(4) road in remote sensing images is identified by road area and directional information, obtains major trunk roads and subsidiary road institute The closed area of formation.
Further, step 4 is when the watershed algorithm based on texture gradient carries out image segmentation, specific algorithm step packet It includes:
(1) including PiClosed area in all POI points carry out initialization tag image, generate image M;
(2) M is traversed, edge enhancing is carried out;
(3) continuous gradient is calculated, obtained point is stored in corresponding height above sea level queue;
(4) pixel in height above sea level queue is traversed, m watershed line is found according to simulative immersion process, m item divides water Closed area is divided into m+1 region by ridge line;Complete the segmentation to closed area.
Compared with prior art, the invention has the following advantages:
The calculation method of to text rate proposed by the present invention is effective, and calculating step is accurate and convenient, has relatively strong Novelty.After the present invention realizes that cultural land used region area automation calculates by optimization area computation method in GIS, base Originally it is calculated as each cultural land area summation, realizes simply, is not required to by complicated equation, calculating speed is very fast, very convenient;This Convenient for realizing using remote sensing and GIS, road Identification and image segmentation are the basic content of remote sensing images identification, area surface for invention Product calculating is the basic function of GIS, and only need to optimize region area algorithm can be obtained.The present invention is by identifying cultural land used interest Point (Point of Interest, POI) carries out text ground rate based on remote sensing images and GIS and calculates, and calculated result is more accurate, calculates Efficiency is higher
Detailed description of the invention
Fig. 1 is a kind of process signal of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System of the invention Block diagram;
Fig. 2 is Xi'an culture land used POI distribution schematic diagram in the embodiment of the present invention;
Fig. 3 is clock tower neighbouring road extraction schematic diagram in Xi'an in the embodiment of the present invention;
Fig. 4 is clock tower areal calculation schematic diagram in Xi'an in the embodiment of the present invention.
Specific embodiment
Invention is further described in detail in the following with reference to the drawings and specific embodiments.
Referring to Fig. 1, a kind of literary ground rate calculation method based on remote sensing images and GIS of the invention, comprising the following steps:
Step 1, it treats zoning and carries out the town site area recognition based on index image, obtain urban construction Total land area A;When carrying out town site area recognition, selects and normalize exposed index as town site Indicator is removed the noise information in non-constructive land region by Threshold segmentation, mask process, obtains more accurately cities and towns Construction land information.
Specific town site area recognition method are as follows:
1, it calculates and normalizes exposed index, the wave band that spectrum gray average is directed between 20-120 range is denoted as Wave band 1, the wave band that spectrum gray average is directed between 20-180 range are denoted as wave band 2.Then normalize exposed index are as follows: wave 1/ wave band 1+ wave band 2 of section 2- wave band.
2, given threshold 140.
It 3, is 1 to the pixel assignment for being higher than threshold value in remote sensing image, the pixel assignment lower than threshold value is 0.
4, interference is carried out to index image to handle.
Step 2, cultural land used POI data is imported in the GIS map vector in region to be calculated, cultural land used is denoted as Pi, Wherein, the value of i is 1 to n;The GIS map vector in region to be calculated can be downloaded from network, such as by cultural land used POI data It imports on the map vector of Xi'an;Cultural land used include: 1. cultural spirit mark land used, i.e., country, area, urban development Cheng Zhong has great cultural significance and spirit mark value, and the natural landscape with the approval of height social mentality and humane mark Know shared land used;2. commemorating land used, that is, carry out cultural memorial and educational activities place land used;3. religion land used, that is, open The lands used such as land used, such as temple, mausoleum, sight, Buddhist nunnery, church needed for opening up religious activities;4. cultural heritage land used has history, skill Art, scientific value and without other lands used using function;5. cultural facility land used meets the daily Public Culture life of resident Land used shared by all kinds of cultural facilities of demand.POI data includes: classification, title and coordinate.
Step 3, the distant of the GIS map vector in the region to be calculated of cultural land used POI data and region to be calculated will be imported Feel map autoregistration, remote sensing images automatic identification culture land used PiMajor trunk roads, subsidiary road are formed by within the scope of neighbouring 100m Closed area.
Referring to Fig. 2, according to the pertinent regulations of national " urban planning quota index temporary provisions ", road can in step 3 It is divided into level Four.Wherein, major trunk roads are had a lot of social connections for 45-55m, and subsidiary road is had a lot of social connections for 40-50m.Major trunk roads are level-one road, secondary dry Road is secondary road.Map vector and remote sensing map carry out autoregistration, need remote sensing images and vector are opened in ENVI Figure, is registrated using original image space with the remote sensing images of the corresponding dot pair trial zone on normed space.Specific steps are as follows: 1. selecting select Ground Control Points in the register of main menu in ENVI, then select Image to map Geometric precision correction is carried out to polar plot and remote sensing figure.2. in ENVI by correct after image based on, used in master image window Polar plot is added to the image after correcting up by Function/Overlays/Vector Layers.By taking Xi'an as an example, match Image after standard is as shown in Fig. 2.
In step 3, automatic identification is carried out to remote sensing images by remote sensing images roads recognition method, specific steps include: (1) block of pixels for reading designated position on Raster by the Read method of IRaster interface in GIS, then passes through pixel The GetVal method of block obtains the pixel value of designated position, so that it is determined that cultural land used P on remote sensing imagesiBelong in neighbouring 100 meters In the pixel of road.(2) directional information for belonging to the pixel of road is directly calculated.(3) pixel for belonging to road is clicked through Row expansive working.(4) road in remote sensing images is identified by road area and directional information.
Referring to Fig. 3,1. determine the pixel for belonging to road on remote sensing images near clock tower by taking clock tower as an example.2. determining Belong to the directional information of the pixel of road.3. the pixel for pair belonging to road carries out expansive working.4. according to expansive working knot Fruit determines road.Final road Identification result is as shown in Fig. 3.
Step 4, the P obtained in step 3iM non-culture land-use style POI are imported in closed area is denoted as Qm, with PiFor in The heart based on the non-cultural land-use style POI of m, carries out image segmentation using the watershed algorithm based on texture gradient, by image It is divided into m+1 region, selection includes PiTexture Segmentation region as cultural land used region, be denoted as Ci
Step 4 is when the watershed algorithm based on texture gradient carries out image segmentation, core algorithm are as follows: 1. include at one PiClosed area in all POI points carry out initialization tag image, generate image M.2. traversing M, edge enhancing is carried out.3. Continuous gradient is calculated, obtained point is stored in corresponding height above sea level queue.4. the pixel in height above sea level queue is traversed, Watershed line is found according to simulative immersion process, closed area is divided into m+1 region by m watershed line.5. pair enclosed area Domain is split.P can be obtained as precisely as possible in closed area using the methodiSurface area.
Referring to Fig. 4, step 5, i is successively traversed from 1 to n, to cultural land used region CiAutomation areal calculation is carried out, is remembered For Si.This step is substantially carried out to cultural land used PiThe automation of area calculates, and by taking clock tower as an example, key step is as follows: 1. call Draw function in GIS, draws polygonal region.Polygon facet Product function this._calArea (polygon) is calculated 2. rewriteeing. It is about 1300 square metres by the way that clock tower area is calculated, as shown in figure 4, it is smaller to calculate error.In addition, wild goose column area is about 2000 square metres, about 1350 square metres of drum tower area.
Step 6, the literate land area of institute and S are calculated:SiFor each cultural land used CiArea.Such as attached drawing 2 It is shown, seek the sum of 109 cultural land areas.
Step 7, seek literate land area and S and the total land area A of urban construction ratio R: R=S/A;R is For with calculating the text acquired rate.
Literary ground rate is the ratio between the cultural land used gross area and planning construction land area within the scope of planning construction land used.This Invention is a kind of method calculated using remote sensing images and GIS as the literary ground rate that main means are unfolded, to solve urban culture land used Assessment, planning and monitoring providing method are supported.The drip irrigation device are as follows: 1, carry out the town site based on index image Area recognition.2, cultural land used point of interest (POI) is imported in GIS-Geographic Information System (GIS);3, map vector and remote sensing images After autoregistration, nearby level-one, subsidiary road are formed by closed area to remote sensing images automatic identification culture land used POI;4, it is sealing Other type POI are imported in closed region, centered on cultural land used POI, based on other types POI, using based on texture ladder The watershed algorithm of degree carries out image segmentation, obtains cultural land used region;5, automation calculates cultural land used region area;6, tire out Add the literate land used region area of institute;7, the cultural land used region gross area and the town site gross area carry out division operation, With obtaining text rate.The present invention is very novel, innovative very strong, with having hewed out text rate calculation method.Calculate step accurately and just It is prompt.Calculating speed of the present invention is fast, can be seen that from step 1 to step 7, realizes culture by optimization area computation method in GIS After the automation calculating of land used region area, basic calculating is each cultural land area summation, realizes simply, is not required to by complexity side Journey.Convenient for realizing using remote sensing and GIS, road Identification and image segmentation are the basic content of remote sensing images identification, area for invention Domain areal calculation is the basic function of GIS, need to only optimize region area algorithm.The assessment of energy objective quantitative of the present invention one The case where urban culture land used, while certain reference role is played to urban culture plan for land.

Claims (10)

1. a kind of text ground rate calculation method based on remote sensing images and GIS-Geographic Information System, which is characterized in that specific steps include:
Step 1, it treats zoning and carries out the town site area recognition based on index image, obtain urban construction and always use Ground area A;
Step 2, cultural land used POI data is imported in the GIS map vector in region to be calculated, cultural land used is denoted as Pi, i expression The number of cultural land used, value are 1 to n;
Step 3, the GIS map vector in the region to be calculated of cultural land used POI data and the remote sensing in region to be calculated will be imported Figure registration identifies cultural land used PiMajor trunk roads and major trunk roads, subsidiary road and subsidiary road and major trunk roads within the scope of surrounding 100m Closed area is formed by with subsidiary road;
Step 4, m non-culture land-use style POI are imported in the closed area that step 3 obtains is denoted as Qm, with cultural land used PiFor Center, point based on the non-cultural land-use style POI of m, carries out image segmentation using the watershed algorithm based on texture gradient, will Image segmentation is m+1 region, and selection includes PiTexture Segmentation region as cultural land used region, be denoted as Ci
Step 5, i is successively traversed from 1 to n, to each of obtaining cultural land used region CiAreal calculation is carried out, S is denoted asi
Step 6, the literate land area of institute and S are calculated:SiFor each cultural land used CiArea;
Step 7, seek literate land area and S and the total land area A of urban construction ratio R: R=S/A;R is to calculate The literary ground rate acquired.
2. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 1, special Sign is, in step 1, when carrying out town site area recognition, selects and normalizes exposed index as urban construction use The indicator on ground removes the noise information in non-constructive land region by Threshold segmentation, mask process, obtains city-building use Ground information.
3. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 1, special Sign is, in step 1, the specific steps of town site area recognition include:
(1) it calculates and normalizes exposed index: the wave band that spectrum gray average is directed between 20-120 range is denoted as wave Section 1, the wave band that spectrum gray average is directed between 20-180 range is denoted as wave band 2;
Normalize the calculation formula of exposed index are as follows: normalize exposed index=(wave band 2- wave band 1)/(wave band 1+ wave band 2);
(2) given threshold;It is 1 to the pixel assignment for being higher than threshold value in remote sensing image, the pixel assignment lower than threshold value is 0;It assigns Value is town site area for 1.
4. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 3, special Sign is that the threshold value in step (2) is set as 140.
5. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 3, special Sign is, further includes: step (3);
(3) interference is carried out to the image after step (2) assignment to handle.
6. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 1, special Sign is that cultural land used is divided into cultural spirit mark land used, commemorates land used, religion land used, cultural heritage land used and cultural facility Land used;
POI data includes: classification, title and coordinate.
7. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 1, special Sign is, step 3, will import the GIS map vector in the region to be calculated of cultural land used POI data and the remote sensing in region to be calculated Map autoregistration passes through remote sensing images automatic identification culture land used PiIt is major trunk roads within the scope of surrounding 100m and major trunk roads, secondary Arterial highway and subsidiary road and major trunk roads and subsidiary road are formed by closed area;Major trunk roads are had a lot of social connections for 45m-55m, and subsidiary road is had a lot of social connections For 40m-50m.
8. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 7, special Sign is that the specific steps that map vector carries out autoregistration with remote sensing map include:
(1) select Ground Control Points, reselection Image are selected in the register of main menu in ENVI To map carries out geometric precision correction to polar plot and remote sensing figure;
(2) in ENVI by correct after image based on, use Function/Overlays/Vector in master image window Polar plot is added on the image after correcting by Layers, obtains registration image.
9. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 7, special Sign is, passes through remote sensing images automatic identification culture land used PiMajor trunk roads within the scope of surrounding 100m and major trunk roads, subsidiary road with The specific steps that subsidiary road and major trunk roads are formed by closed area with subsidiary road include:
(1) block of pixels for reading predetermined position on Raster by the Read method of IRaster interface in GIS, passes through pixel The GetVal method of block obtains the pixel value in predetermined position, determines cultural land used P on remote sensing imagesiBelong in 100 meters of surrounding The pixel on road;
(2) directional information for belonging to the pixel of road is calculated;
(3) expansive working is carried out to the pixel for belonging to road;
(4) road in remote sensing images is identified by road area and directional information, obtains major trunk roads and is formed with subsidiary road Closed area.
10. a kind of literary ground rate calculation method based on remote sensing images and GIS-Geographic Information System according to claim 1, special Sign is that when the watershed algorithm based on texture gradient carries out image segmentation, specific algorithm step includes: step 4
(1) including PiClosed area in all POI points carry out initialization tag image, generate image M;
(2) M is traversed, edge enhancing is carried out;
(3) continuous gradient is calculated, obtained point is stored in corresponding height above sea level queue;
(4) pixel in height above sea level queue is traversed, m watershed line, m watershed line are found according to simulative immersion process Closed area is divided into m+1 region;Complete the segmentation to closed area.
CN201811197720.6A 2018-10-15 2018-10-15 Text-to-ground rate calculation method based on remote sensing image and geographic information system Active CN109376638B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811197720.6A CN109376638B (en) 2018-10-15 2018-10-15 Text-to-ground rate calculation method based on remote sensing image and geographic information system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811197720.6A CN109376638B (en) 2018-10-15 2018-10-15 Text-to-ground rate calculation method based on remote sensing image and geographic information system

Publications (2)

Publication Number Publication Date
CN109376638A true CN109376638A (en) 2019-02-22
CN109376638B CN109376638B (en) 2022-03-04

Family

ID=65400411

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811197720.6A Active CN109376638B (en) 2018-10-15 2018-10-15 Text-to-ground rate calculation method based on remote sensing image and geographic information system

Country Status (1)

Country Link
CN (1) CN109376638B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111063004A (en) * 2019-11-18 2020-04-24 东南大学 Regional road material stock calculation method based on high-resolution remote sensing map and GIS
CN112650822A (en) * 2019-10-09 2021-04-13 浙江宇视科技有限公司 Map access method, map access device, electronic equipment and machine-readable storage medium

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1790052A (en) * 2005-12-19 2006-06-21 武汉大学 Area feature variation detection method based on remote sensing image and GIS data
CN1904875A (en) * 2006-09-07 2007-01-31 上海大学 Normalized distribution image computer generating method for space scattered light on rolling ground
CN101034472A (en) * 2007-03-29 2007-09-12 上海大学 Landform transformation of satellite remote sensing digital image supported by GIS
US20090089078A1 (en) * 2007-09-28 2009-04-02 Great-Circle Technologies, Inc. Bundling of automated work flow
CN101493932A (en) * 2009-03-05 2009-07-29 西安电子科技大学 Watershed texture imaging segmenting method based on morphology Haar small wave texture gradient extraction
CN101614822A (en) * 2009-07-17 2009-12-30 北京大学 A method for detecting road damage based on post-disaster high-resolution remote sensing images
CN102194120A (en) * 2011-01-24 2011-09-21 武汉理工大学 Method for extracting land for teaching by using remote sensing image, relative elevation and geographic ontology
CN103235952A (en) * 2013-04-03 2013-08-07 中国科学院遥感应用研究所 Urban green space measuring method based on high-resolution remote sensing image
CN103903014A (en) * 2014-04-17 2014-07-02 西安煤航卫星数据应用有限公司 Desertification interpretation method based on satellite remote sensing technology
CN103971115A (en) * 2014-05-09 2014-08-06 中国科学院遥感与数字地球研究所 Automatic extraction method for newly-increased construction land image spots in high-resolution remote sensing images based on NDVI and PanTex index
CN104408463A (en) * 2014-10-15 2015-03-11 中国土地勘测规划院 High-resolution construction land graph spot identification method based on PanTex and linear characteristic
CN104616274A (en) * 2015-02-09 2015-05-13 内蒙古科技大学 Algorithm for fusing multi-focusing image based on salient region extraction
CN105912764A (en) * 2016-04-06 2016-08-31 东南大学 Noise zoning model based city floor area ratio optimization method
CN106960436A (en) * 2017-03-20 2017-07-18 哈尔滨工业大学 A kind of remote sensing image process performance appraisal procedure
CN107966116A (en) * 2017-11-20 2018-04-27 苏州市农业科学院 The remote-sensing monitoring method and system of a kind of Monitoring of Paddy Rice Plant Area
CN108596465A (en) * 2018-04-17 2018-09-28 西安建筑科技大学 A kind of urban residence building system carbon energy measuring method

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1790052A (en) * 2005-12-19 2006-06-21 武汉大学 Area feature variation detection method based on remote sensing image and GIS data
CN1904875A (en) * 2006-09-07 2007-01-31 上海大学 Normalized distribution image computer generating method for space scattered light on rolling ground
CN101034472A (en) * 2007-03-29 2007-09-12 上海大学 Landform transformation of satellite remote sensing digital image supported by GIS
US20090089078A1 (en) * 2007-09-28 2009-04-02 Great-Circle Technologies, Inc. Bundling of automated work flow
CN101493932A (en) * 2009-03-05 2009-07-29 西安电子科技大学 Watershed texture imaging segmenting method based on morphology Haar small wave texture gradient extraction
CN101614822A (en) * 2009-07-17 2009-12-30 北京大学 A method for detecting road damage based on post-disaster high-resolution remote sensing images
CN102194120A (en) * 2011-01-24 2011-09-21 武汉理工大学 Method for extracting land for teaching by using remote sensing image, relative elevation and geographic ontology
CN103235952A (en) * 2013-04-03 2013-08-07 中国科学院遥感应用研究所 Urban green space measuring method based on high-resolution remote sensing image
CN103903014A (en) * 2014-04-17 2014-07-02 西安煤航卫星数据应用有限公司 Desertification interpretation method based on satellite remote sensing technology
CN103971115A (en) * 2014-05-09 2014-08-06 中国科学院遥感与数字地球研究所 Automatic extraction method for newly-increased construction land image spots in high-resolution remote sensing images based on NDVI and PanTex index
CN104408463A (en) * 2014-10-15 2015-03-11 中国土地勘测规划院 High-resolution construction land graph spot identification method based on PanTex and linear characteristic
CN104616274A (en) * 2015-02-09 2015-05-13 内蒙古科技大学 Algorithm for fusing multi-focusing image based on salient region extraction
CN105912764A (en) * 2016-04-06 2016-08-31 东南大学 Noise zoning model based city floor area ratio optimization method
CN106960436A (en) * 2017-03-20 2017-07-18 哈尔滨工业大学 A kind of remote sensing image process performance appraisal procedure
CN107966116A (en) * 2017-11-20 2018-04-27 苏州市农业科学院 The remote-sensing monitoring method and system of a kind of Monitoring of Paddy Rice Plant Area
CN108596465A (en) * 2018-04-17 2018-09-28 西安建筑科技大学 A kind of urban residence building system carbon energy measuring method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ASIM ANEES等: ""Detecting beetle infestations in pine forests using MODIS NDVI time-series data"", 《2013 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM - IGARSS》 *
GENGXING ZHAO,等: ""Remote Sensing Image Based Information Extraction for Land Salinized Degradation and Its Evolution"", 《2010 SIXTH INTERNATIONAL CONFERENCE ON NATURAL COMPUTATION (ICNC 2010) 》 *
YU LIJUN等: ""Remote sensing archaeology analysis of Liangzhu culture"", 《2013 DIGITAL HERITAGE INTERNATIONAL CONGRESS (DIGITALHERITAGE)》 *
郑时龄等: ""城市设计的中国智慧"", 《建筑学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112650822A (en) * 2019-10-09 2021-04-13 浙江宇视科技有限公司 Map access method, map access device, electronic equipment and machine-readable storage medium
CN112650822B (en) * 2019-10-09 2023-09-15 浙江宇视科技有限公司 Map access method, map access device, electronic equipment and machine-readable storage medium
CN111063004A (en) * 2019-11-18 2020-04-24 东南大学 Regional road material stock calculation method based on high-resolution remote sensing map and GIS

Also Published As

Publication number Publication date
CN109376638B (en) 2022-03-04

Similar Documents

Publication Publication Date Title
Lin et al. Measuring the relationship between morphological spatial pattern of green space and urban heat island using machine learning methods
Li et al. Quantifying the shade provision of street trees in urban landscape: A case study in Boston, USA, using Google Street View
CN108764193B (en) Urban functional area dividing method fusing POI and remote sensing image
CN110348506B (en) Land use classification method, storage medium and computing device based on remote sensing images
Sui et al. Modeling the dynamics of landscape structure in Asia’s emerging desakota regions: a case study in Shenzhen
CN107145872B (en) A method for obtaining the spatial distribution of desert riparian forest based on GIS buffer analysis
KR102368261B1 (en) Planting location analysis device and method
Dmowska et al. High resolution dasymetric model of US demographics with application to spatial distribution of racial diversity
CN114385712B (en) Country ecological landscape multi-source data space fusion method based on GNSS
CN112070056A (en) Sensitive land use identification method based on object-oriented and deep learning
CN113780307A (en) Method for extracting blue-green space information with maximum regional year
CN109657728B (en) Sample production method and model training method
Shi et al. An improved framework for assessing the impact of different urban development strategies on land cover and ecological quality changes-A case study from Nanjing Jiangbei New Area, China
CN109376638A (en) A method for calculating land coverage based on remote sensing images and geographic information system
CN108287845A (en) A kind of Automatic extraction method for road information and device and hybrid navigation system
Liang et al. Characterizing street trees in three metropolises of central China by using Street View data: From individual trees to landscape mapping
Wang et al. Quantitative relationship between urban green canopy area and urban greening land area
CN116562125B (en) Street security simulation prediction method and system using crowd characteristics and street built environment characteristics in coupling mode
CN116415499B (en) Community comfort simulation prediction method
KR100918008B1 (en) Wetland location search method using wes
CN115393886B (en) A method and storage medium for identifying urban functional areas based on POI and AOI data
CN115880427A (en) Three-dimensional construction method and device of natural environment
See et al. Mapping human impact using crowdsourcing
Fang et al. Data-driven city: an innovative approach to urban area delineation
Ratnaparkhi et al. Classification of land use and land cover using remotely sensed data for Parbhani City, Maharashtra, India

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