Figure 1.
The flowchart of the iterative algorithm.
Figure 1.
The flowchart of the iterative algorithm.
Figure 2.
The flowchart of the proposed method.
Figure 2.
The flowchart of the proposed method.
Figure 3.
The current chart of interferometric phase generation based on M4S software.
Figure 3.
The current chart of interferometric phase generation based on M4S software.
Figure 4.
Basic network framework of the pix2pix model.
Figure 4.
Basic network framework of the pix2pix model.
Figure 5.
The input ocean surface current field velocity: (a) sample 1, (b) sample 2, (c) sample 3, (d) sample 4, (e) sample 5, (f) sample 6, (g) sample 7, (h) sample 8, and (i) sample 9.
Figure 5.
The input ocean surface current field velocity: (a) sample 1, (b) sample 2, (c) sample 3, (d) sample 4, (e) sample 5, (f) sample 6, (g) sample 7, (h) sample 8, and (i) sample 9.
Figure 6.
The estimation ocean surface current field velocity: (a) sample 1, (b) sample 2, (c) sample 3, (d) sample 4, (e) sample 5, (f) sample 6, (g) sample 7, (h) sample 8, and (i) sample 9.
Figure 6.
The estimation ocean surface current field velocity: (a) sample 1, (b) sample 2, (c) sample 3, (d) sample 4, (e) sample 5, (f) sample 6, (g) sample 7, (h) sample 8, and (i) sample 9.
Figure 7.
The control group of the input ocean surface current field velocity.
Figure 7.
The control group of the input ocean surface current field velocity.
Figure 8.
The estimation velocity of ocean surface currents under the different epochs from (a–d): (a) 100 epochs, (b) 200 epochs, (c) 300 epochs, and (d) 400 epochs. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different epochs from (e–h): (e) 100 epochs, (f) 200 epochs, (g) 300 epochs, and (h) 400 epochs.
Figure 8.
The estimation velocity of ocean surface currents under the different epochs from (a–d): (a) 100 epochs, (b) 200 epochs, (c) 300 epochs, and (d) 400 epochs. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different epochs from (e–h): (e) 100 epochs, (f) 200 epochs, (g) 300 epochs, and (h) 400 epochs.
Figure 9.
The estimation velocity of ocean surface currents under the different batch sizes from (a–d): (a) the batch size of 1, (b) the batch size of 2, (c) the batch size of 3, and (d) the batch size of 4. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different batch sizes from (e–h): (e) the batch size of 1, (f) the batch size of 2, (g) the batch size of 3, and (h) the batch size of 4.
Figure 9.
The estimation velocity of ocean surface currents under the different batch sizes from (a–d): (a) the batch size of 1, (b) the batch size of 2, (c) the batch size of 3, and (d) the batch size of 4. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different batch sizes from (e–h): (e) the batch size of 1, (f) the batch size of 2, (g) the batch size of 3, and (h) the batch size of 4.
Figure 10.
The estimation velocity of ocean surface currents under the different sample sizes from (a–d): (a) 100 samples, (b) 200 samples, (c) 300 samples, and (d) 400 samples. The scatter density map between the estimation velocity of the ocean surface currents and the control group velocity of the input ocean surface currents under the different sample sizes from (e–h): (e) 100 samples, (f) 200 samples, (g) 300 samples, and (h) 400 samples.
Figure 10.
The estimation velocity of ocean surface currents under the different sample sizes from (a–d): (a) 100 samples, (b) 200 samples, (c) 300 samples, and (d) 400 samples. The scatter density map between the estimation velocity of the ocean surface currents and the control group velocity of the input ocean surface currents under the different sample sizes from (e–h): (e) 100 samples, (f) 200 samples, (g) 300 samples, and (h) 400 samples.
Figure 11.
The estimation velocity of ocean surface currents under the different learning rates from (a–d): (a) 0.0001, (b) 0.0002, (c) 0.0003, and (d) 0.0004. The scatter density map between the estimation velocity of the ocean surface currents and the control group velocity of the input ocean surface currents under the different learning rates from (e–h): (e) 0.0001, (f) 0.0002, (g) 0.0003, and (h) 0.0004.
Figure 11.
The estimation velocity of ocean surface currents under the different learning rates from (a–d): (a) 0.0001, (b) 0.0002, (c) 0.0003, and (d) 0.0004. The scatter density map between the estimation velocity of the ocean surface currents and the control group velocity of the input ocean surface currents under the different learning rates from (e–h): (e) 0.0001, (f) 0.0002, (g) 0.0003, and (h) 0.0004.
Figure 12.
The velocity calculation results based on the direct method: (a) the input velocity, (b) the current velocity calculated directly from the interferometric phase, and (c) the scatter density map of the input and calculated current velocities.
Figure 12.
The velocity calculation results based on the direct method: (a) the input velocity, (b) the current velocity calculated directly from the interferometric phase, and (c) the scatter density map of the input and calculated current velocities.
Figure 13.
The velocity calculation results based on the iterative method: (a) the input velocity, (b) the iterative current velocity, and (c) the scatter density map of the input and calculated current velocities.
Figure 13.
The velocity calculation results based on the iterative method: (a) the input velocity, (b) the iterative current velocity, and (c) the scatter density map of the input and calculated current velocities.
Figure 14.
The velocity calculation results based on the proposed method: (a) the input velocity, (b) the calculated velocity obtained through the trained pix2pix network, and (c) the scatter density map of the input and calculated current velocities.
Figure 14.
The velocity calculation results based on the proposed method: (a) the input velocity, (b) the calculated velocity obtained through the trained pix2pix network, and (c) the scatter density map of the input and calculated current velocities.
Figure 15.
The estimation velocity of ocean surface currents under the different wind velocities from (a–c): (a) 5 m/s, (b) 10 m/s, and (c) 15 m/s. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different wind velocities (d–f): (d) 5 m/s, (e) 10 m/s, and (f) 15 m/s.
Figure 15.
The estimation velocity of ocean surface currents under the different wind velocities from (a–c): (a) 5 m/s, (b) 10 m/s, and (c) 15 m/s. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different wind velocities (d–f): (d) 5 m/s, (e) 10 m/s, and (f) 15 m/s.
Figure 16.
The estimation velocity of ocean surface currents under the different multi-looks from (a–c): (a) 16 looks, (b) 64 looks, and (c) 100 looks. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different multi-looks from (d–f): (d) 16 looks, (e) 64 looks, and (f) 100 looks.
Figure 16.
The estimation velocity of ocean surface currents under the different multi-looks from (a–c): (a) 16 looks, (b) 64 looks, and (c) 100 looks. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different multi-looks from (d–f): (d) 16 looks, (e) 64 looks, and (f) 100 looks.
Figure 17.
The estimation velocity of ocean surface currents under the different incident angles from (a–e): (a) 13°, (b) 28°, (c) 35°, (d) 43°, and (e) 50°. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different incident angles from (f–j): (f) 13°, (g) 28°, (h) 35°, (i) 43°, and (j) 50°.
Figure 17.
The estimation velocity of ocean surface currents under the different incident angles from (a–e): (a) 13°, (b) 28°, (c) 35°, (d) 43°, and (e) 50°. The scatter density map between the estimation velocity of ocean surface currents and the control group velocity of the input ocean surface currents under the different incident angles from (f–j): (f) 13°, (g) 28°, (h) 35°, (i) 43°, and (j) 50°.
Figure 18.
Real data from the two-channel along-track InSAR mode of the Gaofen-3 spaceborne SAR system (size: 13182(azimuth)*15316(range)).
Figure 18.
Real data from the two-channel along-track InSAR mode of the Gaofen-3 spaceborne SAR system (size: 13182(azimuth)*15316(range)).
Figure 19.
The extracted interferometric phase from the Gaofen-3 along-track InSAR data (size: 8192(azimuth)*8192(range)).
Figure 19.
The extracted interferometric phase from the Gaofen-3 along-track InSAR data (size: 8192(azimuth)*8192(range)).
Figure 20.
The extracted interferometric phase from the Gaofen−3 along−track InSAR data after 1024 multi−looks and three times 10*10 element boxcar averaging (size: 256(azimuth)*256(range)).
Figure 20.
The extracted interferometric phase from the Gaofen−3 along−track InSAR data after 1024 multi−looks and three times 10*10 element boxcar averaging (size: 256(azimuth)*256(range)).
Figure 21.
Results of the direct method: (a) the calculated current filed, (b) the simulated interferometric phase, and (c) the scatter density map of simulated and extracted interferometric phases.
Figure 21.
Results of the direct method: (a) the calculated current filed, (b) the simulated interferometric phase, and (c) the scatter density map of simulated and extracted interferometric phases.
Figure 22.
Results of the iterative method: (a) the calculated current filed, (b) the simulated interferometric phase, and (c) the scatter density map of simulated and extracted interferometric phases.
Figure 22.
Results of the iterative method: (a) the calculated current filed, (b) the simulated interferometric phase, and (c) the scatter density map of simulated and extracted interferometric phases.
Figure 23.
Results of the proposed method: (a) the calculated current filed, (b) the simulated interferometric phase, and (c) the scatter density map of simulated and extracted interferometric phases.
Figure 23.
Results of the proposed method: (a) the calculated current filed, (b) the simulated interferometric phase, and (c) the scatter density map of simulated and extracted interferometric phases.
Table 1.
Experimental environment.
Table 1.
Experimental environment.
Setting | Software | Setting | Software |
---|
System | Ubuntu16.04 | Tool | Anaconda3 |
RAM | 16 GB | Programming | Python3.6 |
CPU | Intel i5-9600KF 3.7GHz × 6 | IDE | Eclipse |
GPU | NVIDIA GTX 1080Ti | Framework | Pytorch-GPU |
Auxiliary tools | MATLAB | Others | CUDA9.0 |
Table 2.
System simulation parameters.
Table 2.
System simulation parameters.
System Parameter | Value | System Parameter | Value |
---|
Center Frequency | 5.4 GHz | Polarization Mode | VV |
Effective Baseline | 28 m | Spatial Resolution | 25 m |
Platform Height | 755 km | Platform Speed | 7000 m/s |
Incident Angle | 35° | NESZ | −25 dB |
The Number of Multi-Look | 100 | Wind Speed | 5 m/s–15 m/s |
Table 3.
Statistics of the estimation results.
Table 3.
Statistics of the estimation results.
| RMSE (m/s) | r | Bias (m/s) |
---|
Sample 1 | 0.019 | 0.760 | 0.028 |
Sample 2 | 0.027 | 0.703 | 0.019 |
Sample 3 | 0.017 | 0.818 | 0.015 |
Sample 4 | 0.019 | 0.785 | 0.009 |
Sample 5 | 0.026 | 0.764 | 0.028 |
Sample 6 | 0.027 | 0.724 | 0.019 |
Sample 7 | 0.018 | 0.779 | 0.013 |
Sample 8 | 0.025 | 0.795 | 0.008 |
Sample 9 | 0.020 | 0.782 | 0.014 |
Statistical average | 0.022 | 0.768 | 0.017 |
Table 4.
Statistics of the estimation results under different epochs.
Table 4.
Statistics of the estimation results under different epochs.
Epoch | RMSE (m/s) | r | Bias (m/s) |
---|
100 | 0.022 | 0.670 | 0.015 |
200 | 0.018 | 0.755 | 0.004 |
300 | 0.017 | 0.795 | 0.002 |
400 | 0.019 | 0.763 | 0.005 |
Table 5.
Statistics of the estimation results under different batch sizes.
Table 5.
Statistics of the estimation results under different batch sizes.
Batch Size | RMSE (m/s) | r | Bias (m/s) |
---|
1 | 0.017 | 0.795 | 0.002 |
2 | 0.025 | 0.736 | 0.009 |
3 | 0.026 | 0.693 | 0.010 |
4 | 0.027 | 0.661 | 0.015 |
Table 6.
Statistics of the estimation results under different sample numbers.
Table 6.
Statistics of the estimation results under different sample numbers.
The Number of Samples | RMSE (m/s) | r | Bias (m/s) |
---|
100 | 0.024 | 0.652 | 0.017 |
200 | 0.023 | 0.664 | 0.014 |
300 | 0.021 | 0.670 | 0.009 |
400 | 0.017 | 0.795 | 0.002 |
Table 7.
Statistics of the estimation results under different learning rates.
Table 7.
Statistics of the estimation results under different learning rates.
Learning Rate | RMSE (m/s) | r | Bias (m/s) |
---|
0.0001 | 0.022 | 0.700 | 0.008 |
0.0002 | 0.017 | 0.800 | 0.002 |
0.0003 | 0.021 | 0.736 | 0.006 |
0.0004 | 0.022 | 0.728 | 0.007 |
Table 8.
Statistics of the estimation results under different methods.
Table 8.
Statistics of the estimation results under different methods.
Method | RMSE (m/s) | r | Bias (m/s) |
---|
Direct Method | 0.104 | 0.714 | 0.028 |
Iterative Method | 0.084 | 0.641 | 0.010 |
Proposed Method | 0.017 | 0.795 | 0.002 |
Table 9.
Statistics of the estimation results under different wind velocities.
Table 9.
Statistics of the estimation results under different wind velocities.
Wind Velocity (m/s) | RMSE (m/s) | r | Bias (m/s) |
---|
5 | 0.017 | 0.795 | 0.002 |
10 | 0.020 | 0.701 | 0.003 |
15 | 0.030 | 0.373 | 0.019 |
Table 10.
Statistics of the estimation results under different multi-looks.
Table 10.
Statistics of the estimation results under different multi-looks.
Multi-Look | RMSE (m/s) | r | Bias (m/s) |
---|
16 | 0.020 | 0.728 | 0.002 |
64 | 0.019 | 0.750 | 0.002 |
100 | 0.017 | 0.795 | 0.002 |
Table 11.
Statistics of the estimation results under different incident angles.
Table 11.
Statistics of the estimation results under different incident angles.
Incident Angle (°) | RMSE (m/s) | r | Bias (m/s) |
---|
20 | 0.031 | 0.296 | 0.016 |
28 | 0.020 | 0.704 | 0.003 |
35 | 0.017 | 0.795 | 0.002 |
43 | 0.018 | 0.761 | 0.005 |
50 | 0.026 | 0.508 | 0.014 |
Table 12.
Statistics of the estimation results under different methods.
Table 12.
Statistics of the estimation results under different methods.
Method | RMSE (rad) | r | Bias (rad) |
---|
Direct Method | 0.113 | 0.759 | 0.109 |
Iterative Method | 0.066 | 0.658 | 0.026 |
Proposed Method | 0.042 | 0.791 | 0.007 |