TSX: wrong range pixel from longitude calculated by geo2radar

[falkamelung]

Problem

Notes:

• TSX:

  [1] > /work2/05861/tg851601/stampede2/code/rsmas_insar/sources/MintPy/src/mintpy/objects/coord.py(307)geo2radar()
  -> lat_c = (np.nanmax(lat) + np.nanmin(lat)) / 2.
  (Pdb++) rg_step
  1.9792806721121812
  (Pdb++) ut0.range_ground_resolution(self.src_metadata, print_msg=True)
  near   incidence angle : 43.4144 degree
  center incidence angle : 43.5662 degree
  far    incidence angle : 43.7180 degree
  1.9792806721121812

  -Sentinel:

  [1] > /work2/05861/tg851601/stampede2/code/rsmas_insar/sources/MintPy/src/mintpy/objects/coord.py(307)geo2radar()
  -> lat_c = (np.nanmax(lat) + np.nanmin(lat)) / 2.
  (Pdb++) az_step
  123.70146752891817
  (Pdb++) rg_step
  61.86057572897589
  (Pdb++) rg_step = ut0.range_ground_resolution(self.src_metadata, print_msg=True)
  near   incidence angle : 41.7216 degree
  center incidence angle : 43.8661 degree
  far    incidence angle : 46.0107 degree

  

-TSX:

[11] > /work2/05861/tg851601/stampede2/code/rsmas_insar/sources/MintPy/src/mintpy/utils/utils0.py(133)incidence_angle()
-> range_f = range_n+dR*width
(Pdb++) n
[11] > /work2/05861/tg851601/stampede2/code/rsmas_insar/sources/MintPy/src/mintpy/utils/utils0.py(134)incidence_angle()
-> inc_angle_n = (np.pi - np.arccos((r**2 + range_n**2 - (r+H)**2)/(2*r*range_n))) * 180.0/np.pi
(Pdb++) n
[11] > /work2/05861/tg851601/stampede2/code/rsmas_insar/sources/MintPy/src/mintpy/utils/utils0.py(135)incidence_angle()
-> inc_angle_f = (np.pi - np.arccos((r**2 + range_f**2 - (r+H)**2)/(2*r*range_f))) * 180.0/np.pi
(Pdb++) n
[11] > /work2/05861/tg851601/stampede2/code/rsmas_insar/sources/MintPy/src/mintpy/utils/utils0.py(136)incidence_angle()
-> inc_angle_c = (inc_angle_n + inc_angle_f) / 2.0
(Pdb++) inc_angle_n
43.4144332910256
(Pdb++) inc_angle_f
43.718006404952604
(Pdb++) dR
1.3641053359701238
(Pdb++) width
2197

[falkamelung]

(Pdb++) points_lalo
array([[ 25.8761, -80.1219],
       [ 25.8788, -80.1213]])
(Pdb++) coord_minopy.geo2radar(25.8761,-80.1219)[0:2]
(1875, 602)
(Pdb++) coord_minopy.geo2radar(25.8761,-80.1213)[0:2]
(1870, 570)
(Pdb++) coord_minopy.geo2radar(25.8788,-80.1219)[0:2]
(1726, 626)
(Pdb++) coord_minopy.geo2radar(25.8788,-80.1213)[0:2]
(1721, 595)

-> print (yg_minopy1, xg_minopy1, yg_minopy2, xg_minopy2)
(Pdb++) n
1875 602 1721 595