def get_G_P_L(self):
# 增加G矢量计算
# 第一组经纬度（成像仪）的ECEF坐标系下的值
G_pos = np.zeros(np.append(self.Lons.shape, 3))
high = np.zeros_like(self.Lons)
G_pos[:, :, 0], G_pos[:, :, 1], G_pos[:, :, 2] = pb_space.LLA2ECEF(
self.Lons, self.Lats, high)
self.G_pos = G_pos
def get_G_P_L(self):
# 增加G矢量计算
# 第一组经纬度（成像仪）的ECEF坐标系下的值
G_pos = np.zeros(np.append(self.Lons.shape, 3))
high = np.zeros_like(self.Lons)
G_pos[:, :, 0], G_pos[:, :, 1], G_pos[:, :, 2] = pb_space.LLA2ECEF(
self.Lons, self.Lats, high)
self.G_pos = G_pos
# compute CrIS LOS Vector1 局地球面坐标系 RAE--->东-北-天坐标系ENU
cris_east, cris_north, cris_up = pb_space.RAE2ENU(
self.satAzimuth, self.satZenith, self.satRange)
print cris_east.shape
# compute CrIS LOS Vector2 东-北-天坐标系ENU--->地球中心地球固定坐标系ECEF
L_pos = np.zeros(np.append(self.Lons.shape, 3))
L_pos[:, :, 0], L_pos[:, :, 1], L_pos[:, :, 2] = \
pb_space.ENU2ECEF(
cris_east, cris_north, cris_up, self.Lons, self.Lats)
L_pos = L_pos * -1.0
self.L_pos = L_pos
self.P_pos = G_pos - L_pos