CN110728691A - Multi-temporal water sideline-based coastline automatic judgment method - Google Patents
Multi-temporal water sideline-based coastline automatic judgment method Download PDFInfo
- Publication number
- CN110728691A CN110728691A CN201910949127.0A CN201910949127A CN110728691A CN 110728691 A CN110728691 A CN 110728691A CN 201910949127 A CN201910949127 A CN 201910949127A CN 110728691 A CN110728691 A CN 110728691A
- Authority
- CN
- China
- Prior art keywords
- water
- line
- axis
- lines
- coastline
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 8
- 238000012850 discrimination method Methods 0.000 claims description 3
- 208000035126 Facies Diseases 0.000 claims 2
- 238000012937 correction Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- MQTOSJVFKKJCRP-BICOPXKESA-N azithromycin Chemical compound O([C@@H]1[C@@H](C)C(=O)O[C@@H]([C@@]([C@H](O)[C@@H](C)N(C)C[C@H](C)C[C@@](C)(O)[C@H](O[C@H]2[C@@H]([C@H](C[C@@H](C)O2)N(C)C)O)[C@H]1C)(C)O)CC)[C@H]1C[C@@](C)(OC)[C@@H](O)[C@H](C)O1 MQTOSJVFKKJCRP-BICOPXKESA-N 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20152—Watershed segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The invention fully considers the shoreline position determination principle, namely the boundary on the multi-temporal water sideline obtained based on remote sensing images is similar to the average climax line, namely the shoreline, and invents a method for automatically judging the shoreline based on the multi-temporal water sideline, which comprises the following basic steps: preprocessing a remote sensing image and extracting a water line; calculating the coordinates of the intersection point of the two water line lines and segmenting according to the intersection point; calculating a connecting line and a perpendicular line of each section of end point; establishing a plane rectangular coordinate system; comparing the size of the vertical coordinate of the intersection point of the longitudinal axis and the water line, wherein the water line with the large vertical coordinate is an upper boundary; and judging and combining section by section, and circularly realizing the judgment of a plurality of water lines to obtain the coastline. The method provided by the invention is scientific and reasonable and easy to realize, realizes the automatic discrimination of the coastline based on the multi-temporal water line by utilizing the method, and improves the efficiency of the coastline discrimination.
Description
Technical Field
The invention relates to the field of remote sensing image processing, in particular to the application field of coastline extraction, and specifically relates to a coastline automatic discrimination method based on multi-temporal waterplane lines.
Background
The coastline is a sea-land boundary, which means the sea-land boundary when the average high tide and high tide level is reached for many years in China. Under the influence of natural environment and human development, the coastline is always in a changing state, and the accurate mastering of the coastline position, the transition process and the future change trend has very important significance for guiding activities such as coastline cultivation, sailing transportation and the like.
The traditional shoreline surveying and mapping method mainly comprises a real land measurement method and a photogrammetry method. The coastline scope is wide, the change is fast, the ground feature is broken, has increased coastline investigation's complexity, and traditional detection method intensity of labour is big, working cycle is long, inefficiency, is difficult to realize the dynamic monitoring of coastline, receives the restriction of conditions such as geographical environment moreover, leads to investigation region difficult to reach, and the survey and drawing degree of difficulty is big.
The modern remote sensing technology is a comprehensive application technology for earth observation based on physical means, geological analysis and mathematical methods, has strong data acquisition capacity, has the advantages of large range, high time resolution, multiple spectra, multiple time phases, no limitation of weather and geographic environment conditions and the like, and is an effective technical means for coastline monitoring and extraction.
The existing methods for automatically extracting the coastline based on the remote sensing image are various in types, but most of the results extracted by the automatic extraction methods are instantaneous water lines of the satellite transit time, and the water lines are not necessarily real coastlines, namely average climax and climax lines. The probability that the multi-temporal remote sensing image falls in the time window of the climax is higher, and the boundary of the multi-temporal water borderline obtained based on the remote sensing image can be approximated to the average climax, so that the invention provides the method for automatically judging the boundary of the multi-temporal water borderline, and the automatic judgment of the real coastline is realized.
Disclosure of Invention
Technical problem to be solved
The invention provides a coastline automatic discrimination method based on a multi-temporal waterwall, which realizes automatic discrimination of boundaries on the multi-temporal waterwall by fully considering the principle that the multi-temporal waterwall automatically discriminates the coastline.
(II) technical scheme
The invention comprises the following steps:
(1) and carrying out water sideline extraction on the preprocessed remote sensing image to obtain a multi-temporal water sideline.
Preprocessing the remote sensing images of the same area at different periods and respectively extracting the waterside lines, wherein the waterside line corresponding to each image is named as S in sequence1,S2,…,Sn;
(2) And calculating intersection points of the water lines and segmenting. Selecting two water lines S in the water lines extracted in the step (1)1And S2, S1And S2The coordinates of the intersection points are sequentially (x)1,y1),(x2,y2),…,(xm+1,ym+1) Will S1,S2Dividing the intersection point into m sections, each section including two water lines S1m,S2mThe connecting line of each section of end point is named as L in sequence1,L2,…,Lm;
(3) And establishing a plane rectangular coordinate system on a connecting line of each section of end point. Calculating LmSlope k of1mAnd the midpoint ZmCoordinate (x)zm,yzm) Find ZmAnd with LmVertical straight line CmWith ZmAs the origin of coordinates, in LmIs the x-axis with LmPerpendicular line CmIs the y-axis, to L1,L2,…,LmEstablishing a plane rectangular coordinate system for each section of the water line;
(4) and (5) judging the upper boundary of the water boundary. Respectively calculate CmAnd S1m,S2mPoint of intersection J1m,J2mCoordinate of (D), comparison J1m,J2mThe water side line with the large ordinate is the upper boundary of the mth section, and the upper boundaries of all the sections of water side lines are combined into one water side line H1To obtain S1,S2The upper boundary of the waterside line;
(5) and (4) finishing automatic discrimination of the coastline according to the upper boundary of the water sideline obtained in the step (4). H is to be1And S3Judging according to the steps (2) to (4) to obtain the upper boundary H of the water sideline2And repeating the above steps to sequentially judge the n water side lines until H is obtainedn-2And SnWater edge upper boundary Hn-1And realizing automatic discrimination of the coastline.
Further, the preprocessing in the step (1) comprises radiation correction and geometric correction, and an algorithm for extracting the waterside line is a marked watershed algorithm.
Further, in the step (2), m +1 represents S1,S2The number of the intersection points is specified to only take the part between the head and tail intersection points of the two water lines as the water line to be distinguished, and the rest parts are ignored.
Further, in the step (3), LmThe slope is calculated as:
kLm=(ym+1-ym)/(xm+1-xm)
midpoint ZmCoordinate (x)zm,yzm) The calculation formula is as follows:
xzm=(xm+xm+1)/2
yzm=(ym+ym+1)/2
Cmthe equation of the straight line is as follows:
y-yzm=kCm×(x-xzm)
kCm=-1/kLm
wherein k isLm,kCmAre respectively Lm,Cm(x) slope of (C)zm,yzm) Is LmAnd (4) a midpoint coordinate.
Further, in the step (4), the method for determining the positive direction of the coordinate system comprises: with L1To LmThe positive direction of (2) is the positive x-axis direction, along the positive x-axis direction, if the land is on the right side of the x-axis, the positive y-axis direction is on the right side of the x-axis, and if the land is on the left side of the x-axis, the positive y-axis direction is on the left side of the x-axis.
(III) advantageous effects
The advantages of the invention are embodied in that:
although the automatic coastline extraction based on the remote sensing image has been greatly developed and the existing methods are various, most of the automatic extraction methods have the extraction results of the instantaneous water line of the satellite transit time, and the water line is not necessarily the real coastline, namely the average climax. Therefore, a shoreline position determination principle is proposed: the nearest neighbor land boundary of the multi-temporal water line obtained based on the remote sensing image is approximate to the average climax line. In consideration of the algorithm, the method realizes the discrimination of the nearest neighbor land boundary of the multi-temporal water sideline by using the algorithm, realizes automatic discrimination, improves the accuracy of the discrimination result of the coastline, and has important significance for the research of the coastline.
Drawings
FIG. 1 is a flow chart of the steps performed in the present invention,
fig. 2 is a schematic diagram of a water border line discriminating process.
Detailed Description
In order to make the objects, contents, and advantages of the present invention clearer, the following detailed description of the embodiments of the present invention is provided with reference to the accompanying drawings and examples:
referring to fig. 1, the method comprises the following specific steps:
(1) preprocessing the remote sensing images of the same area at different periods and respectively extracting the waterside lines, wherein the waterside line corresponding to each image is named as S in sequence1,S2,…,Sn;
The preprocessing comprises radiation correction and geometric correction, and the algorithm for extracting the water borderline is a marked watershed algorithm.
(2) Referring to fig. 2, two of the edge lines S extracted in step (1) are selected1And S2,S1,S2The coordinates of the intersection points are sequentially (x)1,y1),(x2,y2),…,(xm+1,ym+1) Will S1,S2Dividing the intersection point into m sections, each section including two water lines S1m,S2mThe connecting line of each section of end point is marked as L in sequence1,L2,…,Lm;
Wherein m +1 represents S1,S2The number of the intersection points is specified to only take the part between the head intersection point and the tail intersection point of the two water side lines as the water side line to be distinguished, and the rest parts are ignored.
(3) And establishing a plane rectangular coordinate system on a connecting line of each section of end point. Calculating LmSlope k of1mAnd the midpoint ZmCoordinate (x)zm,yzm) Find ZmAnd with LmVertical straightLine CmWith ZmIs the origin of coordinates, LmIs the x-axis, LmPerpendicular line CmIs the y-axis, to L1,L2,…,LmEstablishing a plane rectangular coordinate system for each section of the water line;
wherein L ismThe slope is calculated as:
kLm=(ym+1-ym)/(xm+1-xm)
midpoint ZmCoordinate (x)zm,yzm) The calculation formula is as follows:
xzm=(xm+xm+1)/2
yzm=(ym+ym+1)/2
Cmthe equation of the straight line is as follows:
y-yzm=kCmX(x-xzm)
kCm=-1/kLm
wherein k isLm,kCmAre respectively Lm,Cm(x) slope of (C)zm,yzm) Is LmAnd (4) a midpoint coordinate.
Wherein L is defined1To LmThe positive direction of (2) is the positive x-axis direction, along the positive x-axis direction, if the land is on the right side of the x-axis, the positive y-axis direction is on the right side of the x-axis, and if the land is on the left side of the x-axis, the positive y-axis direction is on the left side of the x-axis.
(4) And (5) judging the upper boundary of the water boundary. Respectively calculate CmAnd S1m,S2mPoint of intersection J1m,J2mCoordinate of (D), comparison J1m,J2mThe water side line with the large ordinate is the upper boundary of the mth section, and the upper boundaries of all the sections of water side lines are combined into one water side line H1To obtain S1,S2Is located above the water edge.
(5) And (4) finishing automatic discrimination of the coastline according to the upper boundary of the water sideline obtained in the step (4). H is to be1And S3Judging according to the steps (2) to (4) to obtain the upper boundary H of the water sideline2Analogizing in turn for n water linesSequentially distinguishing until H is obtainedn-2And SnWater edge upper boundary Hn-1And realizing automatic discrimination of the coastline.
Claims (3)
1. A coastline automatic discrimination method based on a multi-temporal waterside line is characterized by comprising the following steps:
(1) and carrying out water sideline extraction on the preprocessed remote sensing image to obtain a multi-temporal water sideline.
Preprocessing the remote sensing images of the same area at different periods and respectively extracting the waterside lines, wherein the waterside line corresponding to each image is named as S in sequence1,S2,…,Sn;
(2) And calculating intersection points of the water lines and segmenting. Selecting two water lines S in the water lines extracted in the step (1)1And S2,S1And S2The coordinates of the intersection points are sequentially (x)1,y1),(x2,y2),…,(xm+1,ym+1) Will S1,S2Dividing the intersection point into m sections, each section including two water lines S1m,S2mThe connecting line of each section of end point is named as L in sequence1,L2,…,Lm;
Wherein m +1 represents S1,S2The number of the intersection points is specified to only take the part between the head and tail intersection points of the two water lines as the water line to be distinguished, and the rest parts are ignored.
(3) And establishing a plane rectangular coordinate system on a connecting line of each section of end point. Calculating LmSlope k of1mAnd the midpoint ZmCoordinate (x)zm,yzm) Find ZmAnd with LmVertical straight line CmWith ZmIs the origin of coordinates, LmIs the x-axis, LmPerpendicular line CmIs the y-axis, to L1,L2,,...,LmEstablishing a plane rectangular coordinate system for each section of the water line;
wherein the positive direction of the x-axis is defined as L1To LmIf the land is on the right side of the positive x-axis direction, the positive y-axis direction is definedThe x-axis right side, and if the land is on the x-axis positive side, the y-axis positive side is defined as the x-axis left side.
(4) And (5) judging the upper boundary of the water boundary. Respectively calculate CmAnd S1m,S2mPoint of intersection J1m,J2mCoordinate of (D), comparison J1m,J2mThe water side line with the large ordinate is the upper boundary of the mth section, and the upper boundaries of all the sections of water side lines are combined into one water side line H1To obtain S1,S2The upper boundary of the waterside line;
(5) and (4) finishing automatic discrimination of the coastline according to the upper boundary of the water sideline obtained in the step (4). H is to be1And S3Judging according to the steps (2) to (4) to obtain the upper boundary H of the water sideline2And repeating the above steps to sequentially judge the n water side lines until H is obtainedn-2And SnWater edge upper boundary Hn-1And realizing automatic discrimination of the coastline.
2. The coastline automatic discrimination algorithm based on multi-temporal facies lines of claim 1, wherein: the rectangular coordinate system of the plane in the step (3) is ZmIs the origin of coordinates, LmIs the x-axis, LmPerpendicular line CmIs the y axis, and the positive direction of the x axis is L1To LmThe positive y-axis direction is the right side of the x-axis if the land is on the right side of the positive x-axis direction, and the positive y-axis direction is the left side of the x-axis if the land is on the left side of the positive x-axis direction.
3. The coastline automatic discrimination algorithm based on multi-temporal facies lines of claim 1, wherein: the multi-temporal borderline judging process in the step (5) is to use H1And S3Judging according to the steps (2) to (5) to obtain the upper boundary H of the water edge line2And repeating the above steps to sequentially judge the n water side lines until H is obtainedn-2And SnWater edge upper boundary Hn-1I.e. coastline Hn-1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910949127.0A CN110728691B (en) | 2019-10-08 | 2019-10-08 | Multi-temporal water sideline-based coastline automatic judgment method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910949127.0A CN110728691B (en) | 2019-10-08 | 2019-10-08 | Multi-temporal water sideline-based coastline automatic judgment method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110728691A true CN110728691A (en) | 2020-01-24 |
CN110728691B CN110728691B (en) | 2021-03-23 |
Family
ID=69220795
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910949127.0A Expired - Fee Related CN110728691B (en) | 2019-10-08 | 2019-10-08 | Multi-temporal water sideline-based coastline automatic judgment method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110728691B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111461046A (en) * | 2020-04-10 | 2020-07-28 | 生态环境部卫星环境应用中心 | Automatic coast erosion and siltation identification method and device based on shoreline data |
CN111680424A (en) * | 2020-06-11 | 2020-09-18 | 南京师范大学 | River attack automatic discrimination method based on x-shaped graphic state characteristics |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105374041A (en) * | 2015-11-18 | 2016-03-02 | 国家海洋局第一海洋研究所 | Method of extracting sandy coastline by using multiple periods of remote sensing images |
WO2017149744A1 (en) * | 2016-03-04 | 2017-09-08 | 株式会社日立国際電気 | Water level measurement system and water level measurement method |
KR20180020421A (en) * | 2016-08-18 | 2018-02-28 | 경북대학교 산학협력단 | Method and system for extracting coastline based on a large-scale high-resolution satellite images |
CN109919070A (en) * | 2019-02-28 | 2019-06-21 | 南京师范大学 | A kind of coastline remote sensing projectional technique of profile morphology auto-adapted fitting |
-
2019
- 2019-10-08 CN CN201910949127.0A patent/CN110728691B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105374041A (en) * | 2015-11-18 | 2016-03-02 | 国家海洋局第一海洋研究所 | Method of extracting sandy coastline by using multiple periods of remote sensing images |
WO2017149744A1 (en) * | 2016-03-04 | 2017-09-08 | 株式会社日立国際電気 | Water level measurement system and water level measurement method |
KR20180020421A (en) * | 2016-08-18 | 2018-02-28 | 경북대학교 산학협력단 | Method and system for extracting coastline based on a large-scale high-resolution satellite images |
CN109919070A (en) * | 2019-02-28 | 2019-06-21 | 南京师范大学 | A kind of coastline remote sensing projectional technique of profile morphology auto-adapted fitting |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111461046A (en) * | 2020-04-10 | 2020-07-28 | 生态环境部卫星环境应用中心 | Automatic coast erosion and siltation identification method and device based on shoreline data |
CN111461046B (en) * | 2020-04-10 | 2020-12-25 | 生态环境部卫星环境应用中心 | Automatic coast erosion and siltation identification method and device based on shoreline data |
CN111680424A (en) * | 2020-06-11 | 2020-09-18 | 南京师范大学 | River attack automatic discrimination method based on x-shaped graphic state characteristics |
Also Published As
Publication number | Publication date |
---|---|
CN110728691B (en) | 2021-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111460984B (en) | Global lane line detection method based on key points and gradient equalization loss | |
CN111860336B (en) | High-resolution remote sensing image inclined ship target detection method based on position sensing | |
CN106989683B (en) | A kind of shield tail clearance of shield machine vision measuring method | |
CN107063197B (en) | Reservoir characteristic curve extraction method based on spatial information technology | |
CN109100741A (en) | A kind of object detection method based on 3D laser radar and image data | |
CN111951392B (en) | Method for reconstructing topography above beach withered water level based on time series remote sensing image and water level monitoring data | |
CN102542551B (en) | Automatic change detection technology for floating ice at edges of polar ice sheets | |
CN110378921B (en) | Intelligent identification method for substrate layer boundary of channel based on floating mud rheological property and gray level co-occurrence matrix | |
CN103871072B (en) | Orthography based on project digital elevation model inlays line extraction method | |
CN112766184B (en) | Remote sensing target detection method based on multi-level feature selection convolutional neural network | |
CN110619368B (en) | Planet surface navigation feature imaging matching detection method | |
CN101672915A (en) | High spatial resolution remote sensing image crown outline delineation system and method | |
CN110728691B (en) | Multi-temporal water sideline-based coastline automatic judgment method | |
Bremer et al. | Eigenvalue and graph-based object extraction from mobile laser scanning point clouds | |
CN109919070A (en) | A kind of coastline remote sensing projectional technique of profile morphology auto-adapted fitting | |
CN106530313A (en) | Sea-sky line real-time detection method based on region segmentation | |
James et al. | Measuring 3D coastal change with a digital camera | |
CN113628227A (en) | Coastline change analysis method based on deep learning | |
CN102567726B (en) | Technology for automatically extracting floating ice in polar ice sheet edge regions | |
CN111597930A (en) | Coastline extraction method based on remote sensing cloud platform | |
CN112926468A (en) | Tidal flat elevation automatic extraction method | |
CN115877400A (en) | Tunnel roof support steel belt drilling positioning method based on radar and vision fusion | |
CN110956083A (en) | Bohai sea ice drift remote sensing detection method based on high-resolution four-signal optical satellite | |
CN115293262B (en) | Shoreline function and structure long time sequence detection method | |
Li et al. | Automatic Road Extraction from High-Resolution Remote Sensing Image Based on Bat Model and Mutual Information Matching. |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210323 |