Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

dem_error w/ phase velocity: rm all-zero column in design matrix #905

Merged
merged 2 commits into from
Nov 12, 2022
Merged

dem_error w/ phase velocity: rm all-zero column in design matrix #905

Changes from 1 commit
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
29d986a
dem_error w/ phase velocity: rm all-zero column in design matrix
yunjunz Nov 9, 2022
c30adb4
more comments
yunjunz Nov 12, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
more comments
  • Loading branch information
@yunjunz
yunjunz committed Nov 12, 2022
commit c30adb45045d1d80b6e1f6339a46679efe8af1f5
6 changes: 4 additions & 2 deletions src/mintpy/dem_error.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -249,9 +249,11 @@ def estimate_dem_error(ts0, G0, tbase, date_flag=None, phase_velocity=False,
tbase_diff = np.diff(tbase[date_flag], axis=0).reshape(-1,1)
ts = np.diff(ts, axis=0) / np.repeat(tbase_diff, ts.shape[1], axis=1)
G = np.diff(G, axis=0) / np.repeat(tbase_diff, G.shape[1], axis=1)
# remove the all-zero column in G
# remove the all-zero column in G and all-one column in G0
# as the unknown constant term of the polynomial func is not estimated,
# resulting in a larger time-series residual ts_res, thus, not suitable
# for RMS based noise evaluation with timeseries_rms.py
G = np.hstack((G[:,:1], G[:,2:]))
# remove the all-one column in G0 because it is not estimated
G0 = np.hstack((G0[:,:1], G0[:,2:]))

# Inverse using L-2 norm to get unknown parameters X
Expand Down