CN114170524A - Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion - Google Patents

Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion Download PDF

Info

Publication number
CN114170524A
CN114170524A CN202111262990.2A CN202111262990A CN114170524A CN 114170524 A CN114170524 A CN 114170524A CN 202111262990 A CN202111262990 A CN 202111262990A CN 114170524 A CN114170524 A CN 114170524A
Authority
CN
China
Prior art keywords
time
tidal
goci
pixel
wfv
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111262990.2A
Other languages
Chinese (zh)
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.)
East China Normal University
Original Assignee
East China Normal University
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 East China Normal University filed Critical East China Normal University
Priority to CN202111262990.2A priority Critical patent/CN114170524A/en
Publication of CN114170524A publication Critical patent/CN114170524A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a single-day tidal flat digital terrain construction method for high time-space satellite data fusion, which fuses GF-1/WFV satellite data and geostationary satellite GOCI data by using a STARFM time-space data fusion method to generate hourly images with the spatial resolution of 16 m; the tidal water level at the image transit time is obtained by combining the hydrodynamic model, the daily terrain elevation of the tidal flat is obtained, and the remote sensing dynamic monitoring of the tidal flat terrain is realized. Through the comprehensive application of high-resolution satellite observation and dynamic simulation, the difficulty of dynamic monitoring of tidal flat terrain is solved, and a new method for dynamic monitoring of silt tidal flat terrain is established.

Description

Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion
Technical Field
The invention belongs to the technical field of remote sensing construction of tidal flat terrain Digital Elevation models (Digital Elevation models DEMs), and particularly relates to a remote sensing method for generating single-day tidal flat terrain Digital elevations by fusing high-time-high spatial resolution satellite data.
Background
The tidal flat is a shoal with exposed intertidal zones, contains precious wetland biological resources and land resources, is influenced by human engineering, sea level rise and the like, and silts and erosion of the shoal are converted along with the tidal flat, so that the terrain elevation is changed rapidly, and certain threats are caused to the tidal flat wetland resources and the land resources. Therefore, the method has important significance for timely and real-timely grasping the rapid change of the tidal flat terrain elevation. In general, conventional topographic surveying techniques (e.g., leveling, total station, etc.) can provide very accurate tidal flat topographic elevation measurements, but the survey space coverage and work efficiency are very low, and a lot of manpower and financial resources are consumed, and furthermore, tidal flat wetland, mudflat, mud and marsh environments are very difficult to measure under the influence of tidal water fluctuation. The satellite remote sensing technology has the advantages of space-time continuous and large-area synchronous observation, short measurement period, low cost, non-contact type and the like, and realizes the tidal flat terrain measurement by constructing the DEM through the instantaneous water line and the synchronous water level information of the satellite. Generally, optical remote sensing images are richer in data and easier to obtain compared with radar images and stereopair images. Therefore, it is currently most widespread and practical to construct tidal flat digital terrain elevation by using optical remote sensing image data.
Due to tidal power factors, coastal zone tidal beaches are in dynamic change at any moment, but for the remote sensing technology, the spatial distribution range of the coastal zone silt tidal beaches is not large, and a satellite image with high resolution is required. The number of time-frequency low high-resolution satellite images is extremely limited, and in addition, due to the influence of cloud covers in coastal zone areas, the DEM is difficult to construct under the condition that a time sequence series image data set with high-low tide level changes is met in a short time, most scholars adopt different data sources to achieve the condition, the shortest time is shortened to within 1 year, within 0.5 year and the like, but the rapid dynamic change of the wetland in the intertidal zone cannot be really solved.
Disclosure of Invention
The invention aims to provide a method for constructing a single-day tidal flat terrain DEM (digital elevation model) by data space-time fusion of hourly revisit frequency and low spatial resolution (hundred meter level) remote sensing images of a static orbit satellite and multi-day revisit frequency and high spatial resolution (ten meter level) remote sensing images of an polar orbit satellite. The method has the advantages that a single-day tidal flat terrain digital elevation model is constructed, dynamic and rapid change monitoring of tidal flat terrain elevation is achieved, and the major requirements of real-time dynamic monitoring of tidal flat wetland resources and sustainable development of land resources in national economy and ecological civilization construction are met.
The specific technical scheme for realizing the purpose of the invention is as follows:
a single-day tidal flat digital terrain construction method for high-time-space satellite data fusion is characterized by comprising the following steps: the method comprises the following specific steps:
step 1, fusing high-space-time resolution images, and generating hourly images with the spatial resolution of 16m by adopting a space-time adaptive reflection fusion model STARFM;
step 2, acquiring a tide level by using a tide power model TPXO8, and acquiring tide data by using a tide power model TPXO 8;
step 3, based on the water line extraction of the BP neural network, accurately classifying the muddy tidal wetland and the water body by adopting a BP neural network supervision classification method, and further extracting the water line;
and 4, constructing the DTM, and constructing the tidal flat terrain elevation by combining the tidal level and the water line data.
The step 1 specifically comprises:
step 1.1, firstly, introducing information of adjacent and similar pixels in a window, firstly resampling a GOCI image after geometric correction and geographic registration to a GF-1/WFV image with 16m resolution by adopting a STARFM algorithm, and setting pixels (x)i,yj) Type of ground coverage and system error over time t0And image date tkIf not, then:
GF(xi,yj,t0)=G(xi,yj,t0)+GF(xi,yj,tk)-G(xi,yj,tk) (1)
in the formula: GF and G represent GF-1/WFV and GOCI pictures, respectively, (x)i,yj) Is the pixel position, t0To predict time, tkTime of an existing image pair;
step 1.2, predicting the central pixel value by utilizing a sliding window and a weight function
The STARFM algorithm assumes that the errors of the GF-1/WFV and GOCI imagery are constant for the same ground target, and that the errors are systematic in a short time, relying only on the characteristics of the individual pixels; therefore, the error of each pixel is estimated by using a pair of GF-1/WFV-GOCI images at a certain time, and then the high spatial resolution GF-1/WFV images are predicted by using the GOCI images at other times;
the STARFM algorithm first resamples the geometrically corrected and geographically registered GOCI image to a GF-1/WFV image of 16m resolution, sets the pixels (x)i,yj) Type of ground coverage and system error over time t0And image date tkIf not, then:
GF(xi,yj,t0)=G(xi,yj,t0)+GF(xi,yj,tk)-G(xi,yj,tk) (1)
in the formula: GF and G represent GF-1/WFV and GOCI pictures, respectively, (x)i,yj) Is the pixel position, t0To predict time, tkTime of an existing image pair;
however, the GOCI image is usually a mixed pixel under the GF-1/WFV 16m resolution, or within the forecast date, the tidal bank water side line is often changed due to the influence of tide; the predicted time t is calculated using the weighting function in the STARFM algorithm by introducing additional information of neighboring pixels0Surface reflectance of center pixel:
Figure BDA0003326115220000031
in the formula: ω is the size of the search window, (x)ω/2,yω/2) Is the center pixel position of the search window; weight WijkDetermining a contribution of each neighboring pixel to the estimated reflectivity of the central pixel; wijkFrom spectral differences, time differences, central pixel position (x)ω/2,yω/2) And candidate pixel location (x)i,yj) The distance between them is determined;
the spectral difference between GF-1/WFV and GOCI data at a given location is:
Sijk=|GF(xi,yj,tk)-G(xi,yj,tk)| (3)
the difference between the GF-1/WFV observation value and the average value of the adjacent pixels under low spatial resolution is approximately obtained by formula (3); sijkA smaller value of (a) indicates that the pixel has a spectral characteristic closer to the average of its surrounding pixels, i.e., the variation of the high resolution pixel is comparable to the average variation of the surrounding pixels, where the reflectance of the pixel is assigned a higher weight in equation (2); the time difference between the input GOCI data and the predicted GOCI data is:
Tijk=|G(xi,yj,tk)-G(xi,yj,t0)| (4)
Tijkthe variation between prediction and acquisition times is measured, the smaller TijkMeaning the time tkAnd t0The change between is smaller and the pixel should be assigned a higher weight; at the date tkTemporal center pixel position (x)ω/2,yω/2) And candidate pixel location (x)i,yj) The distance between them is:
Figure BDA0003326115220000032
dijkthe spatial distance between the center predicted pixel and the surrounding spectrally similar candidate pixel is assigned a higher weight for its candidate pixel for closer pixels.
The step 2 specifically comprises:
step 2.1, the tidal height of the tidal model TPXO8 is given in the form of complex amplitude, so the partial tide at time t for a single component of frequency f at location x is given by the following equation:
h(t,x)=pu(t,x)·Re[h(x)exp{i[f(t-t0)+V0(t0)+ph(t,x)]}] (6)
wherein V0(t0) Is t0The astronomical parameters at time, pu (t, x) and ph (t, x) are node correction quantities; then the amplitude is | h | and the phase is tan-1(-IM (h)/RE (h)), where IM and RE represent imaginary Number part and Real Number part of complex Number, i x i ═ 1;
and 2.2, obtaining Matlab toolbox codes from the OSU-TPXO Tide Models (https:// www.tpxo.net/home), and inputting the accurate position coordinates and time information into the TPXO8 model to obtain Tide data.
The step 3 specifically includes:
step 3.1, selecting two types of interested areas of the intertidal zone wetland and the water body as training samples, and taking the training samples as an input layer for supervision and classification;
step 3.2, utilizing a BP neural network based on ENVI5.5 (trial version) to supervise and classify the images;
and 3.3, outputting the supervised classification result and the separation degree of the sample, wherein the separation degree of the sample is larger if the difference between the two classification areas is obvious, the separation degree is better, the separation degree is maximum 2.0, and a parting line between the two classification areas is a water line.
The step 4 specifically includes:
step 4.1, constructing a Delaunay (triangulation method) irregular triangulation network, and performing subdivision by using the Delaunay irregular triangulation;
and 4.2, converting TIN (triangulation) into grids, and interpolating the numerical values into a grid DTM (digital terrain elevation) to obtain the digital terrain elevation.
The invention has the beneficial effects that:
compared with the prior art, the method has the advantages that the time interval of the images is shortened to 1h, the spatial resolution is 16m, and the monitoring of the DTM of the tidal flat is facilitated. The tidal flat wetland resource real-time dynamic monitoring is realized, and a foundation and a premise are provided for the protection and treatment of the estuary coastal zone.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a STARFM spatiotemporal fusion flow diagram;
fig. 3 is a TPXO8 model interface diagram;
FIG. 4 is a flow chart of the waterside line extraction;
FIG. 5 is a flow chart of a DTM construction using a water sideline;
FIG. 6 is a DTM verification graph constructed from Lianxing harbor measured tide level and TPXO8 model tide level;
FIG. 7 is a Chongming east beach DTM and difference plot across a one-year scale;
FIG. 8 is a comparison of water lines at 500m and 16m spatial resolution;
FIG. 9 is a DTMLIDAnd annual DTMON-fusion (DTM)LID-DTMnon-fusion).
Detailed Description
Referring to fig. 1, the present invention performs data fusion of GF-1/WFV (spatial resolution: 16 m; revisit frequency: 4 days) and GOCI (spatial resolution: 500 m; revisit frequency: 1 hour) using a spatio-temporal adaptive reflection fusion model (STARFM), and generates an hourly video with a spatial resolution of 16 m. Extracting a water line by using a BP (Back propagation) neural network supervision and classification method, and constructing a digital ground model (DTM) of the tidal flat of the Yangtze river estuary by combining with the tide level data acquired by the TPXO8 tide model; the method specifically comprises the following steps:
step 1: high spatial-temporal resolution image fusion
In order to obtain high-precision intertidal wetland terrain elevation by using a water line extraction method, a satellite image is required to have higher spatial resolution so as to meet the precision requirement of DTM. Meanwhile, for the water line extraction, the revisit period of the satellite images is required to be short enough to capture the rapid change dynamics of the intertidal zone wetland elevation. Aiming at the problem of insufficient space-time resolution of remote sensing images, the STARFM fusion method is used for fusing GF-1 data with high spatial resolution and GOCI data. The model considers the difference of space and the difference of time in the fusion process, and is one of the most widely applied space-time fusion models at present. Compared with other fusion models, STARFM is more suitable for regions dominated by spatio-temporal changes, and is very suitable for the study of dynamically changing tidal flat DTM. Therefore, STARFM is adopted to fuse GF-1 and GOCI data, and the space-time resolution of the remote sensing image is further improved.
STARFM uses a sliding window and a weight function to predict the central pixel value of the window by introducing information of neighboring pixels and similar pixels in the window. The algorithm assumes that the errors of the GF-1 and GOCI shots are invariant when observing the same ground target. The error depends on the characteristics of the individual pixels and is systematic over a short time interval. Therefore, if one image pair of GF-1 and GOCI can be obtained as a reference at a certain time, an error of each pixel in the image can be estimated. And then predicting the fine resolution image by using the GOCI images at other time instants.
Step 1.1: the STARFM algorithm first resamples the geometry-corrected and geo-registered GOCI image onto a GF-1 image of 16m resolution, assuming pixels (x)i,yj) Type of ground coverage and system error over time t0And image date tkIf not, then:
GF(xi,yj,t0)=G(xi,yj,t0)+GF(xi,yj,tk)-G(xi,yj,tk) (1)
in the formula: GF and G represent GF-1/WFV and GOCI pictures, respectively, (x)i,yj) Is the pixel position, t0To predict time, tkThe existing image is compared to time.
Step 1.2: because GOCI images are usually mixed pixels at GF-116 m resolution, or because of tides and the like within the forecast date, tidal shoal water lines are often changed. Thus STARFM uses a sliding window and a weighting function to predict the value of the center pel within the sliding window by introducing information of neighboring and similar pels within the window. By introducing additional information on neighboring pixels, the STARFM algorithm calculates the surface reflectance of the pixel in the center of the predicted date using a weighting function:
Figure BDA0003326115220000061
in the formula: ω is the size of the search window, (x)ω/2,yω/2) Is the center pixel position of the search window. Weight WijkThe contribution of each neighboring pixel to the estimated reflectivity of the central pixel is determined. WijkFrom spectral differences, time differences, central pixel position (x)ω/2,yω/2) And candidate pixel location (x)i,yj) The distance between them, etc.
The spectral difference between GF-1 and GOCI data at a given location is:
Sijk=|GF(xi,yj,tk)-G(xi,yj,tk)| (3)
equation (3) approximately measures the difference between the GF-1 observation and the average of the neighboring pixels at coarse resolution. SijkA smaller value of (a) indicates that the pixel has a spectral characteristic closer to the average of the surrounding pixels, i.e. the variation of the high resolution pixel should be comparable to the average variation of the surrounding pixels. Therefore, the reflectance of the pixel should be assigned a higher weight in equation (2). Input GOCI dataAnd the predicted GOCI data is:
Tijk=|G(xi,yj,tk)-G(xi,yj,t0)| (4)
Tijkthe variation between prediction and acquisition time is measured. Smaller TijkMeaning times t and t0The change between is smaller and the pixel should be assigned a higher weight. At the date tkTemporal center pixel position (x)ω/2,yω/2) And candidate pixel location (x)i,yj) The distance between them is:
Figure BDA0003326115220000062
dijkthe spatial distance between the center predicted pixel and the surrounding spectrally similar candidate pixel is predicted. Spatial similarity is generally better for closer pixels, and therefore, closer candidates should be assigned higher weights. A flow chart of the process is shown in fig. 2. Step 2: TPXO8 tidal power model acquisition tidal level
Step 2.1: the tidal height of TPXO8 is given in complex amplitude, so the partial tide at time t for a single component of frequency f at location x is given by the following equation:
h(t,x)=pu(t,x)·Re[h(x)exp{i[f(t-t0)+V0(t0)+ph(t,x)1}] (6)
wherein V0(t0) Is t0The astronomical parameters at time, pu (t, x) and ph (t, x) are the node correction quantities. Then the amplitude is | h | and the phase is tan-1(-IM (h)/RE (h)) IM and RE represent the imaginary Number part and the Real Number part of the complex Number i ═ 1, respectively.
Step 2.2: matlab kit code was obtained from OSU-TPXO Tide Models (https:// www.tpxo.net/home). The precise position coordinates and time information are input into the TPXO8 model to obtain tidal level data. The input, simulation and output interfaces are shown in fig. 3.
Step 3, water line extraction and DTM construction based on BP neural network
Considering that a data source is a passive non-stereo optical remote sensing image, a water line detection method is adopted to construct the DTM. In order to automatically extract the water sideline, a supervision and classification method based on a BP neural network in ENVI5.5 (trial edition) is adopted, which consists of hidden layers among an input layer, an output layer and an input-output layer, and respectively corresponds to a selected intertidal wetland and water body supervision sample, a supervision and classification result (intertidal wetland and water body) and a neural network supervision and classification process.
Step 3.1: firstly, two types of interesting areas of the intertidal zone wetland and the water body are selected as training samples.
Step 3.2: then, images were supervised classified using a BP neural network based on ENVI5.5 (trial version). The results show that the error is gradually reduced and the separation degree of the final sample is more than 1.9. Typically, if the difference between the two classification regions is significant, the maximum separability of the sample is 2.0. In the silt tidal flat area, the spectrum difference between the intertidal zone wetland and the water body in the visible light range is not obvious, so that if the separability of the sample is more than 1.9, the differentiation is qualified. The intertidal wetland and the water body are classified, and the boundary between the intertidal wetland and the water body is used as a water boundary line.
And 4, step 4: construction of DTM
Based on arcgis10.1, DTM is constructed with a water line (fig. 5). The DTM is generated by the TIN method of constructing a Delaunay irregular triangulation network.
Step 4.1: firstly, data are generated by utilizing Delaunay irregular triangulation, and the specific operation is realized in ArcGIS 10.1. First, in the 3D analysis extension module of arcgis10.1, TIN is generated by selecting Create TIN.
Step 4.2, interpolating the data into a grid DTM, selecting a TIN to rate in a 3D analysis expansion module of arcgis10.1, converting the TIN into a grid to obtain corresponding DTM data, and an operation flow chart is shown in fig. 5.
Examples
The embodiment realizes dynamic monitoring of tidal flat terrain through a space-time fusion model, a remote sensing image and a tidal model. And performing data fusion on GF-1/WFV and GOCI by using a STARFM space-time adaptive reflection fusion model to generate an hourly image with the spatial resolution of 16 meters. A water line is extracted by using a BP (Back propagation) neural network supervision and classification method, and a Yangtze river mouth tidal flat DTM is constructed by combining tide level data acquired by a TPXO8 tide model.
a. Elevation verification of the Yangtze river mouth Chongming northern beach: with DTMLiD(contemporaneous unmanned aerial vehicle-mounted LiDAR data) as a ground truth value, and the DTM constructed by the geostationary satellite is verified. Satellite DTM data is based on 2019 annual hourly time series satellite fusion images of 8, 17, 8 and 9 months, in order to verify DTM without actually measured tide level datamsat accuracy using the data of the station of the sea level (DTM) of FIG. 6assat) and the model tide level Data (DTM) of FIG. 6bTPXO8msat) are respectively assigned as the waterside line elevations to construct satellite DTMs, wherein fig. 6c and 6d are statistical deviation values corresponding to fig. 6a and 6b, and the root mean square errors thereof are 0.22m and 0.16m, respectively. Indicating that the DTM accuracy built using the tide level modeled with TPXO8 is acceptable, it is reasonable to assign a water line using the tide level data of the TPXO8 model in areas without tide stations.
b. Annual change of the landforms of the Yangtze river mouth Chongming east beach: the annual change of the Chongming east beach DTM is researched by the method, 6 time-by-time fusion images of each day of 12-18 days in 2018 and 12-13 days in 2019 are generated by a STARFM space-time fusion method, and two groups of day DTM (shown in FIGS. 7a and 7b) corresponding to 2018 and 2019 are respectively constructed. It was found that there were some areas of greater variation (north in fig. 7 c) in chongming east beach in 2018 and 2019, with significant overall variation, with an average elevation variation of 0.21m and a maximum of 0.98 m. This means that the period of the time-series images should be as short as possible, which is compatible with the method proposed by the present invention.
c. The effect of spatial and temporal resolution on tidal flat terrain elevation: to quantify the effect of spatial resolution on DTM accuracy, high spatial resolution video (GF-1/WFV video, 16m) and low spatial resolution video (GOCI video, 500m) were selected for horizontal line extraction, respectively, for the same region. The maximum horizontal distance between the water line lines of the two data extractions was found to be 372.9m, with an average distance of 134.6 m. If the slope of the tidal flat is assumed to be 0.15 °, the error in estimating the vertical elevation is about 0.35m, and therefore the spatial resolution of the satellite images has a significant impact on the DTM accuracy.
In order to investigate the influence of the time scale of time series image acquisition on the DTM (DTMON-fusion), 6 GF-1/WFV images are used as an experimental area in Chongming North beach, and the DTM (DTMON-fusion) is constructed in a time scale of one year (24 months and 24 days in 2019 to 13 months and 12 months in 2019). The root mean square error is verified to be 0.27m and is larger than 0.16m of daily DTM, which shows that the error of the DTMON-fusion contains the offset caused by the change of the terrain elevation per se within one year.
In fig. 8, (a) is a position diagram, (b) is a water line at a spatial resolution of 500m, and (c) is a water line at a spatial resolution of 16 m. As can be seen from fig. 8(b), the extraction of the water boundary in the GOCI image is affected by the mixed pixels due to the low resolution, resulting in a large deviation of the position of the water boundary, compared to the water boundary extracted in the GF-1/WFV image which is less affected. The maximum horizontal distance between the two waterlines is 372.9m, and the average distance between the two straight lines is 134.6m (the area between the two straight lines divided by the straight line distance between the two ends of the straight line) calculated by ArcGISI 10.1. In the local area, the gradient of the intertidal zone wetland is 0.15 degrees, the error in the horizontal distance is 134.6m, the error in the vertical elevation can be calculated to be about 0.35m, which shows that the spatial resolution of the satellite image has an important influence on the DTM precision, and the higher the spatial resolution of the satellite image is, the higher the DTM construction precision and reliability are.
FIG. 9 also illustrates that in Chongming North beach, the time span also has some effect on the tidal beach DTM. DTMLIDAnd annual DTMON-fusion (i.e., DTM)LIDDTMnon-fusion) and statistical analysis of pixel-by-pixel differences are shown in fig. 9. In FIG. 9(b), 12% of the pixels are distributed in the range of-0.10 m-0.10m, 76% are distributed in the range of-0.30 m-0.30m, and RMSE is 0.27 m. Comparing the error statistics of fig. 6(c) and fig. 9(b), the absolute error of 24% of the pixels in DTMnon-fusion exceeds 0.3m, and the absolute error of 14% of the pixels in dtmffusion exceeds 0.3m, which indicates that the error of DTMnon-fusion includes the offset caused by the elevation change of the terrain itself within one year.

Claims (5)

1.一种高时-空卫星数据融合的单日潮滩数字地形构建方法,其特征在于,该方法包括以下具体步骤:1. a single-day tidal flat digital terrain construction method of high time-space satellite data fusion, is characterized in that, this method comprises the following concrete steps: 步骤1,高时空分辨率影像融合,采用时空自适应反射融合模型STARFM生成空间分辨率为16m的逐小时影像;Step 1, high spatiotemporal resolution image fusion, using the spatiotemporal adaptive reflection fusion model STARFM to generate hourly images with a spatial resolution of 16m; 步骤2,潮汐动力模型TPXO8获取潮位,使用潮汐动力模型TPXO8获取潮汐数据;Step 2, the tidal dynamic model TPXO8 obtains the tidal level, and the tidal dynamic model TPXO8 is used to obtain the tidal data; 步骤3,基于BP神经网络的水边线提取,采用BP神经网络监督分类方法对淤泥质潮汐湿地和水体进行准确分类,进而提取水边线;Step 3, based on the extraction of the water edge of the BP neural network, using the BP neural network supervised classification method to accurately classify the muddy tidal wetlands and water bodies, and then extract the water edge; 步骤4,DTM的构建,结合潮位与水边线数据构建潮滩地形高程。Step 4, DTM construction, combined with tide level and waterline data to construct tidal flat terrain elevation. 2.根据权利要求1所述一种高时-空卫星数据融合的单日潮滩数字地形构建方法,其特征在于,所述步骤1具体包括:2. a kind of single-day tidal flat digital terrain construction method of high time-space satellite data fusion according to claim 1, is characterized in that, described step 1 specifically comprises: 步骤1.1,首先引入窗口中相邻和相似像素的信息,采用STARFM算法首先将几何校正与地理配准之后的GOCI影像重采样到16m分辨率的GF-1/WFV影像上,设像素(xi,yj)上的地面覆盖类型和系统误差在预测时间t0和影像日期tk上不变,则:Step 1.1, first introduce the information of adjacent and similar pixels in the window, and use the STARFM algorithm to first resample the GOCI image after geometric correction and georeferencing to the GF-1/WFV image with a resolution of 16m. , y j ) on the type of ground cover and systematic errors are constant at prediction time t 0 and image date t k , then: GF(xi,yj,t0)=G(xi,yj,t0)+GF(xi,yj,tk)-G(xi,yj,tk) (1)GF( xi , y j , t 0 )=G( xi , y j , t 0 )+GF( xi , y j , t k )-G( xi , y j , t k ) (1) 式中:GF和G分别代表GF-1/WFV与GOCI影像,(xi,yj)为像素位置,t0为预测时间,tk为已有影像对的时间;where GF and G represent GF-1/WFV and GOCI images, respectively, (x i , y j ) is the pixel position, t 0 is the prediction time, and t k is the time of the existing image pair; 步骤1.2,利用滑动窗口和权值函数对中心像素值进行预测Step 1.2, use the sliding window and weight function to predict the center pixel value STARFM算法设对于相同的地面目标,GF-1/WFV和GOCI影像的误差是恒定的,且所述误差在短时间内是系统性的,仅依赖于单个像素的特征;所以,利用某时间的一对GF-1/WFV-GOCI影像来估计各个像素的误差,然后利用其他时刻的GOCI影像预测其高空间分辨率即GF-1/WFV的影像;The STARFM algorithm assumes that for the same ground target, the error of GF-1/WFV and GOCI images is constant, and the error is systematic in a short time and only depends on the characteristics of a single pixel; therefore, using a certain time A pair of GF-1/WFV-GOCI images is used to estimate the error of each pixel, and then the GOCI images at other times are used to predict the high spatial resolution images of GF-1/WFV; STARFM算法首先将几何校正与地理配准之后的GOCI影像重采样到16m分辨率的GF-1/WFV影像上,设像素(xi,yj)上的地面覆盖类型和系统误差在预测时间t0和影像日期tk上不变,则:The STARFM algorithm first resamples the GOCI image after geometric correction and geo-referencing to the GF-1/WFV image with a resolution of 16m. Set the ground coverage type and system error on the pixel (x i , y j ) at the prediction time t 0 and the image date t k unchanged, then: GF(xi,yj,t0)=G(xi,yj,t0)+GF(xi,yj,tk)-G(xi,yj,tk) (1)GF( xi , y j , t 0 )=G( xi , y j , t 0 )+GF( xi , y j , t k )-G( xi , y j , t k ) (1) 式中:GF和G分别代表GF-1/WFV与GOCI影像,(xi,yj)为像素位置,t0为预测时间,tk为已有影像对的时间;where GF and G represent GF-1/WFV and GOCI images, respectively, (x i , y j ) is the pixel position, t 0 is the prediction time, and t k is the time of the existing image pair; 但GOCI影像在GF-1/WFV 16m分辨率下通常是混合像元,或者在预测日期内,潮滩水边线因潮汐影响往往发生了改变;则通过引入相邻像素的附加信息,利用STARFM算法中的加权函数计算预测时间t0中心像素的地表反射率:However, GOCI images are usually mixed pixels at GF-1/WFV 16m resolution, or within the forecast date, the tidal flat waterline is often changed due to tidal effects; by introducing additional information of adjacent pixels, the STARFM algorithm is used. The weighting function in calculates the surface reflectance of the pixel at the center of predicted time t 0 :
Figure FDA0003326115210000021
Figure FDA0003326115210000021
式中:ω是搜索窗口的大小,(xω/2,yω/2)是该搜索窗口的中心像素位置;权重Wijk确定每个相邻像素对中心像素的估计反射率的贡献;Wijk由光谱差异、时间差、中心像素位置(xω/2,yω/2)和候选像素位置(xi,yj)之间的距离决定;where ω is the size of the search window, (x ω/2 , y ω/2 ) is the center pixel position of the search window; the weight W ijk determines the contribution of each adjacent pixel to the estimated reflectivity of the center pixel; W ijk is determined by the spectral difference, the time difference, the distance between the central pixel position (x ω/2 , y ω/2 ) and the candidate pixel position ( xi , y j ); 在给定位置处,GF-1/WFV和GOCI数据之间的光谱差异是:The spectral differences between the GF-1/WFV and GOCI data at a given location are: Sijk=|GF(xi,yj,tk)-G(xi,yj,tk)| (3)S ijk =|GF( xi , y j , t k )-G( xi , y j , t k )| (3) 公式(3)近似得到了GF-1/WFV观测值和低空间分辨率下相邻像素的平均值之间的差异;Sijk的较小值表明该像素具有更接近于其周围像素平均值的光谱特征,即高分辨率像素的变化与周围像素的平均变化相当,此时,在公式(2)中为像素的反射率分配更高的权重;输入的GOCI数据和预测的GOCI数据之间的时间差是:Equation (3) approximates the difference between the GF-1/WFV observations and the average value of adjacent pixels at low spatial resolution; a smaller value of S ijk indicates that the pixel has a value closer to the average value of its surrounding pixels. Spectral features, i.e. the variation of high-resolution pixels is comparable to the average variation of surrounding pixels, at this time, a higher weight is assigned to the reflectance of the pixel in Eq. (2); the difference between the input GOCI data and the predicted GOCI data is The time difference is: Tijk=|G(xi,yj,tk)-G(xi,yj,t0)| (4)T ijk =|G( xi , y j , t k )-G( xi , y j , t 0 )| (4) Tijk衡量了预测和采集时间之间的变化,较小的Tijk意味着时间tk和t0之间的变化较小,应为像素分配更高的权重;在日期tk时中心像素位置(xω/2,yω/2)和候选像素位置(xi,yj)之间的距离为:T ijk measures the change between prediction and acquisition time, a smaller T ijk means less change between time t k and t 0 , pixels should be assigned higher weights; center pixel position at date t k The distance between (x ω/2 , y ω/2 ) and the candidate pixel position ( xi , y j ) is:
Figure FDA0003326115210000022
Figure FDA0003326115210000022
dijk为中心预测像素和周围光谱相似的候选像素之间的空间距离,对于较近的像素来说,为其候选像素指定更高的权重。d ijk is the spatial distance between the central predicted pixel and the surrounding spectrally similar candidate pixels, and for closer pixels, assign a higher weight to its candidate pixels.
3.根据权利要求1所述一种高时-空卫星数据融合的单日潮滩数字地形构建方法,其特征在于,所述步骤2具体包括:3. a kind of single-day tidal flat digital terrain construction method of high time-space satellite data fusion according to claim 1, is characterized in that, described step 2 specifically comprises: 步骤2.1,TPXO8潮汐模型的潮汐高度以复振幅的形式给出,因此,位置x处频率f的单个分量在时间t处的分潮由以下公式给出:In step 2.1, the tidal heights of the TPXO8 tidal model are given as complex amplitudes, so the subtidal division of a single component of frequency f at location x at time t is given by: h(t,x)=pu(t,x)·Re[h(x)exp{i[f(t-t0)+V0(t0)+ph(t,x)]}] (6)h(t,x)=pu(t,x)·Re[h(x)exp{i[f(tt 0 )+V 0 (t 0 )+ph(t,x)]}] (6) 其中V0(t0)是t0时刻的天文参数,pu(t,x)和ph(t,x)是节点改正量;则振幅为|h|,相位为tan-1(-Im(h)/Re(h)),IM与RE分别代表复数的虚数即Imagic Number部分和实数即RealNumber部分,i*i=-1;where V 0 (t 0 ) is the astronomical parameter at time t 0 , and pu(t, x) and ph(t, x) are the node corrections; then the amplitude is |h|, and the phase is tan -1 (-Im(h )/Re(h)), IM and RE represent the imaginary number of the complex number, namely the Imagic Number part and the real number, namely the RealNumber part, i*i=-1; 步骤2.2,从OSU-TPXO Tide Models(https://www,tpxo.net/home)得到Matlab工具箱代码,将精确的位置坐标和时间信息输入TPXO8模型,即得到潮位数据。Step 2.2, get the Matlab toolbox code from OSU-TPXO Tide Models (https://www, tpxo.net/home), and input the precise position coordinates and time information into the TPXO8 model to get the tide level data. 4.根据权利要求1所述一种高时-空卫星数据融合的单日潮滩数字地形构建方法,其特征在于,所述步骤3具体包括:4. the single-day digital terrain construction method of a kind of high time-space satellite data fusion according to claim 1, is characterized in that, described step 3 specifically comprises: 步骤3.1,选取潮间带湿地和水体两类感兴趣区域作为训练样本,由此训练样本作为监督分类的输入层;Step 3.1, select two types of interest areas of intertidal wetlands and water bodies as training samples, so that the training samples are used as the input layer of supervised classification; 步骤3.2,利用基于ENVI5.5(试用版)的BP神经网络对影像进行监督分类;Step 3.2, use the BP neural network based on ENVI5.5 (trial version) to supervise and classify the images; 步骤3.3,输出监督分类结果及样品的分离度,通常两个分类区域之间的差异明显则样本的可分度为越大,并且表示分离程度越好,分离度最大为2.0,两个分类区域之间的分割线即为水边线。Step 3.3, output the supervised classification results and the degree of separation of the samples. Usually, if the difference between the two classification areas is obvious, the greater the degree of divisibility of the sample, and the better the degree of separation, the maximum degree of separation is 2.0, and the two classification areas The dividing line between them is the water edge. 5.根据权利要求1所述一种高时-空卫星数据融合的单日潮滩数字地形构建方法,其特征在于,所述步骤4具体包括:5. the single-day digital terrain construction method of a kind of high time-space satellite data fusion according to claim 1, is characterized in that, described step 4 specifically comprises: 步骤4.1,构造Delaunay不规则三角网,利用Delaunay不规则三角进行剖分;Step 4.1, construct Delaunay irregular triangulation, and use Delaunay irregular triangulation to divide; 步骤4.2,将TIN转换为栅格,把数值插值到网格DTM中,即得到数字地形高程。Step 4.2, convert the TIN to grid, and interpolate the value into the grid DTM to obtain the digital terrain elevation.
CN202111262990.2A 2021-10-28 2021-10-28 Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion Pending CN114170524A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111262990.2A CN114170524A (en) 2021-10-28 2021-10-28 Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111262990.2A CN114170524A (en) 2021-10-28 2021-10-28 Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion

Publications (1)

Publication Number Publication Date
CN114170524A true CN114170524A (en) 2022-03-11

Family

ID=80477571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111262990.2A Pending CN114170524A (en) 2021-10-28 2021-10-28 Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion

Country Status (1)

Country Link
CN (1) CN114170524A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115840761A (en) * 2023-03-01 2023-03-24 自然资源部第三地理信息制图院 Satellite image pixel value modification method, system, equipment and medium
CN117330040A (en) * 2023-12-01 2024-01-02 南通大地测绘有限公司 Tidal flat topography mapping method and system based on unmanned water mapping ship

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110376623A (en) * 2019-07-08 2019-10-25 中国空间技术研究院 Spaceborne GNSS-R specular reflection point oceanic tide amendment localization method and system

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110376623A (en) * 2019-07-08 2019-10-25 中国空间技术研究院 Spaceborne GNSS-R specular reflection point oceanic tide amendment localization method and system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
FENG GAO等: "On the blending of the Landsat and MODIS surface reflectance: predicting daily Landsat surface reflectance", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》, vol. 44, no. 8, 24 July 2006 (2006-07-24), pages 2, XP055650426, DOI: 10.1109/TGRS.2006.872081 *
WENLI GAO等: "Monitoring terrain elevation of intertidal wetlands by utilising the spatial-temporal fusion of multi-source satellite data: A case study in the Yangtze (Changjiang) Estuary", 《GEOMORPHOLOGY》, vol. 383, 5 March 2021 (2021-03-05), pages 1 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115840761A (en) * 2023-03-01 2023-03-24 自然资源部第三地理信息制图院 Satellite image pixel value modification method, system, equipment and medium
CN115840761B (en) * 2023-03-01 2023-04-14 自然资源部第三地理信息制图院 Satellite image pixel value modification method, system, equipment and medium
CN117330040A (en) * 2023-12-01 2024-01-02 南通大地测绘有限公司 Tidal flat topography mapping method and system based on unmanned water mapping ship
CN117330040B (en) * 2023-12-01 2024-02-02 南通大地测绘有限公司 Tidal flat topography mapping method and system based on unmanned water mapping ship

Similar Documents

Publication Publication Date Title
Kang et al. Topographic mapping on large-scale tidal flats with an iterative approach on the waterline method
Li et al. Effects of DEM sources on hydrologic applications
Schumann et al. The utility of spaceborne radar to render flood inundation maps based on multialgorithm ensembles
Jin et al. Downscaling AMSR-2 soil moisture data with geographically weighted area-to-area regression kriging
CN102819023B (en) Method and system of landslide recognition of complicated geological background area based on LiDAR
Zhang et al. Automatic waterline extraction and topographic mapping of tidal flats from SAR images based on deep learning
CN112926468B (en) Tidal flat elevation automatic extraction method
Al-Zubieri et al. Shoreline evolution between Al Lith and Ras Mahāsin on the red sea coast, Saudi Arabia using GIS and DSAS techniques
Liu et al. Seasonal pattern of tidal-flat topography along the Jiangsu middle coast, China, using HJ-1 optical images
CN117423002B (en) Small-scale tidal flat image processing method based on DEM
CN114170524A (en) Single-day tidal flat digital terrain construction method based on high-time-space satellite data fusion
Papasaika et al. Fusion of digital elevation models using sparse representations
CN112285808B (en) A downscaling method for APHRODITE precipitation data
Gao et al. Monitoring terrain elevation of intertidal wetlands by utilising the spatial-temporal fusion of multi-source satellite data: A case study in the Yangtze (Changjiang) Estuary
Yılmaz et al. Determination of the appropriate zone on dam surface for floating photovoltaic system installation using RS and GISc technologies
CN103761375A (en) Method for calibrating tidal flat DEM (digital elevation model) on basis of remote sensing water lines
Mandal et al. An early warning system for inundation forecast due to a tropical cyclone along the east coast of India
Rawat et al. Effect of DEM data resolution on low relief region sub-watershed boundaries delineating using of SWAT model and DEM derived from CARTOSAT-1 (IRS-P5), SRTM and ASTER
Liu et al. Reconstructing Three-dimensional Thermohaline Structures for Mesoscale Eddies Using Satellite Observations and Deep Learning
Zhang et al. A study on coastline extraction and its trend based on remote sensing image data mining
Feng et al. A hierarchical network densification approach for reconstruction of historical ice velocity fields in East Antarctica
Kowalewski et al. Sensitivity of the Baltic Sea level prediction to spatial model resolution
Zhang et al. A method to derive tidal flat topography in Nantong, China using MODIS data and tidal levels
Ihsan et al. Comparison of Potential Energy of Solar Radiation in Rooftop Modeling Using Different Building Levels of Detail
CN112836449B (en) Method for calibrating hydrologic model

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