ska_sdp_instrumental_calibration.numpy_processors.rotation_matrix module

ska_sdp_instrumental_calibration.numpy_processors.rotation_matrix.generate_rotation_matrices(rm, frequency, output_dtype=<class 'numpy.complex64'>)[source]

Generate station rotation matrices based on Rotation Measure (RM) values.

This function calculates the Faraday rotation angle for each station and frequency channel and constructs the corresponding 2x2 rotation matrices. These matrices effectively rotate the polarization plane of the incident radiation.

Parameters:
  • rm (numpy.ndarray) -- A 1D array containing the Rotation Measure values for each station in rad/m^2. Shape: (n_stations,).

  • frequency (numpy.ndarray) -- A 1D array of frequency channels in Hz. Shape: (n_channels,).

  • output_dtype (data-type, optional) -- The desired data type of the output rotation matrix. Default is np.complex64.

Returns:

A 4D array containing the rotation matrices. Shape: (n_stations, n_channels, 2, 2).

Return type:

numpy.ndarray

Notes

The rotation angle :math: phi is calculated as :math: phi = RM . lambda^2.