when sampling is not given, it's calculated based on spataxis and taxis

AmirMardan · AmirMardan · commit f8cc8e90c655 · 2024-12-13T11:13:47.000-05:00
diff --git a/pylops/waveeqprocessing/seismicinterpolation.py b/pylops/waveeqprocessing/seismicinterpolation.py
@@ -272,8 +272,8 @@ def SeismicInterpolation(
                     )
                 else:
                     sampling = (
-                        np.abs(spataxis[1] - spataxis[1]),
-                        np.abs(taxis[1] - taxis[1]),
+                        np.abs(spataxis[1] - spataxis[0]),
+                        np.abs(taxis[1] - taxis[0]),
                     )
             Pop = FFT2D(dims=dims, nffts=nffts, sampling=sampling)
             Pop = Pop.H

Original file line number	Diff line number	Diff line change
`@@ -272,8 +272,8 @@ def SeismicInterpolation(`
`272`	`272`	`)`
`273`	`273`	`else:`
`274`	`274`	`sampling = (`
`275`		`- np.abs(spataxis[1] - spataxis[1]),`
`276`		`- np.abs(taxis[1] - taxis[1]),`
	`275`	`+ np.abs(spataxis[1] - spataxis[0]),`
	`276`	`+ np.abs(taxis[1] - taxis[0]),`
`277`	`277`	`)`
`278`	`278`	`Pop = FFT2D(dims=dims, nffts=nffts, sampling=sampling)`
`279`	`279`	`Pop = Pop.H`