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