Traditional methods generate two-dimensional (2D) inverse synthetic aperture radar (ISAR) images by performing inverse Fourier transforming on computed scattering fields over the frequency–aspect domain. Scattering field calculations using the SBR technique require ray tracing, making the generation of 2D ISAR images time-consuming.
3.1.1. Image Formation
A critical formula is derived, incorporating the ray tracing process into the ISAR image formation formula [
41]. This image-domain ray tube integration formula, developed under small-angle imaging conditions, adheres to the equivalence principle of a radar cross-section (RCS) in both monostatic and bistatic configurations. Given the small-angle imaging condition, second-order terms concerning angle
can be neglected, as expressed in [
18].
From (8), the sum of the contributions from all rays to the 2D image is determined, a process referred to as fast imaging based on forward ray tracing. In (8), the imaging plane is considered as the
plane, with incident plane waves propagating in the
direction, and the observation point is in the
plane, within a small angle range around
. Here,
and
represent the respective components where the ray exits the target,
is the total path traveled by the ray, and
is the number of ray tubes in the forward ray tracing process.
is the bandwidth,
is the central wave number, and
is the observation angle range under the bistatic configuration. The response amplitude of each ray tube
in the image domain is represented as follows:
where
represents the ray-integrating surface element on the target and
represents the geometrical optics (GO) field at the exit of each ray tube.
The computation time required by (8) primarily depends on the number of rays and the discrete points in the imaging scene. Let be the number of rays, the number of grids in the range direction, and the number in the azimuth direction, with the total computational load for the entire imaging process given by . When calculating electrically large targets in rough backgrounds, the increase in the number of ray tubes and discrete points on the imaging window leads to a high time cost. To improve efficiency, the decaying nature of function can be exploited for truncation and acceleration of the computation.
The accuracy of solving the scattering field directly impacts the quality of radar images. As illustrated in
Figure 4,
,
, and
represent three different facets on the target. SBR [
42] techniques follow the path
(forward ray tracing) from the transmitter to the receiver. It may miss certain scattering contributions, especially in complex targets where some regions might not be directly illuminated.
In our proposed bidirectional shooting and bouncing rays (BSBRs) [
43], the scattered path
(backward ray tracing) is also considered in addition to the path from the transmitter to the receiver (i.e.,
). Here, point
acts not only as a reflecting point but also as a scattered source. The inductive electric current generated at point
on the facet scatters electromagnetic waves along path
. Since paths
and
are parallel, additional scattered contributions are captured at the same receiver. Capturing scattered energy from complex surfaces more thoroughly provides additional illumination in areas not directly illuminated, enhancing both the accuracy and completeness of the solution.
The incident electric field in
Figure 4 is denoted as
. During forward ray tracing, the electric field leaving the target is obtained from the incident and reflected electric fields as
. The inductive current induced by the incident wave on the target, represented by
, acts as a new source of scattered electric field
. The interaction between facets, donated by
in
Figure 4, involves the re-radiation of the inductive current. According to Huygens’ principle,
is expressed using the dyadic Green function
in free space as follows:
During the backward ray tracing
, the incident electric field is represented as
. The ray intersects at point
with its plane and undergoes reflection. The reflected field of the ray tube during backward ray tracing is given by
where
is the scattering factor on the reflecting facet and
is the reflection coefficient. The induced electric current during the backward tracing is
The backward scattering field formula is
Ray tracing is employed to establish
forward paths and
backward paths. The total electric field resulting from this BSBR process is computed as
It is assumed that electromagnetic waves radiated by the scattered field behave as plane waves. During the current iterations, facets that have already produced effects are recorded. Marked facets are excluded from further interactions to avoid redundant calculations.
The contribution of each ray to the 2D image in forward ray tracing is defined as
and in backward ray tracing as
. The representation of the 2D image based on BSBRs is given by
It should be noted that although (15) appears formally similar to (8), the total number of ray tubes in (15) and (8) are and , respectively. It implies that a more comprehensive illumination of regions is achieved with BSBRs, enhancing the accuracy and detail of the resulting 2D images.
3.1.2. Validation of Simulated Datasets
A fast image simulation is conducted, followed by a series of experiments on both simulated datasets to validate the simulated images and demonstrate the superiority of the proposed LiOSR-SAR. The experiments were performed on a personal computer (PC) equipped with an Intel i5-10400 CPU, 32 GB of RAM, and a GTX 1660 SUPER graphics card. The development of the architecture was carried out using Pytorch 1.7.0 [
44], with GPU acceleration applied exclusively for the training of LiOSR-SAR.
A simulated dataset was constructed using the FIS-SBR method. Another dataset employed the proposed FIM-BSBR technique. For clarity, these two datasets are referred to as the FIS-SBR and FIM-BSBR. Both datasets utilize consistent radar parameters, as detailed in
Table 1. The imaging scene size was set to 25.5 m by 25.5 m. This scene size is appropriate for the typical dimensions of the vessels used in the experiments, particularly small boats under 10 m in length. It ensures the capture of detailed structural features and the distribution of scattering centers, while also aligning with the required range and azimuth resolution.
The simulated dataset included three target categories, ship 1, ship 2, and ship 3. Each category has distinct dimensions. Ship 1 measures 7.9 m in length, 1 m in width, and 2.2 m in height. Ship 2 has a length of 9.2 m, a width of 1.2 m, and a height of 2.16 m. Ship 3 measures 7.5 m in length, 0.9 m in width, and 2.2 m in height. The dataset contains targets of varying sizes and complexities. These differences reflect the variety in structural features and scattering characteristics [
45]. The size and complexity of the targets directly influence scattering center resolution and classification accuracy. To verify the correctness as well as the efficiency of the dataset, a comparison is presented among three image simulation methods, namely the range Doppler algorithm (RDA) [
46], the FIS-SBR method, and the FIS-BSBR technique.
The comparative imaging results for ship 1 using three different methods—the RDA, FIS-SBR, and the proposed FIS-BSBR—are shown in
Figure 5. To simulate SAR imageries with different polarizations, separate simulations are performed for each polarization mode (HH, HV, VH, and VV). In each mode, the scattering characteristics are calculated based on the specific polarization of the transmitted and received signals. For example, the HH mode uses horizontal polarization for both transmission and reception, while the HV mode uses horizontal transmission and vertical reception. All methods accurately pinpoint the positions of scattering centers, as demonstrated by the discernible basic outline of the ship and the distribution of scattering centers in the ISAR images. It indicates that our method aligns well with the comparative methods in terms of numerical results. However, traditional imaging methods struggle with complex geometric shapes and diverse materials, especially in analyzing coupling effects between scattering centers in intricate structure of the ship. These methods are also hindered by high computational demands and reduced efficiency.
While the FIS-SBR efficiently generates radar images for large targets, it overlooks scattering contributions in complex targets, such as those with cavities, and fails to reveal certain coupling effects. In contrast, as shown in
Figure 5c, the FIS-BSBR significantly enhances the detail representation and the description of coupling among scattering centers compared to
Figure 5a,b. This detailed visualization of interactions between scattering centers is crucial for subsequent open-set target recognition and characteristic analysis.
Quantitative results in
Table 2 confirm that the FIS-BSBR takes more time for image formation compared to the FIS-SBR (64 s vs. 12 s). This is due to the incorporation of both forward and backward ray tracing, which enhances image quality by capturing more detailed scattering information. Unlike the RDA, the FIS-BSBR employs convolution based on the
function, which significantly improves data processing efficiency. This approach leverages the convolution theorem and fast Fourier transform (FFT) to simplify the imaging process by conducting calculations in the frequency domain. As a result, the time spent on multiple scans in the range and azimuth directions is greatly reduced. The data formation time for the FIS-BSBR is considerably shortened, and most of the processing time is concentrated in the convolution step of the imaging process.
Table 3 and
Table 4 detail the distribution of samples per category in the FIS-SBR and FIS-BSBR datasets. Ship 1 and ship 2 are categorized for training as known categories, while ship 3 is designated as an unknown category for testing.