Script Isar Image Processing
Script Isar Image Processing
This application macro computes and displays an inverse synthetic aperture radar (ISAR) image from
backscattered radar cross-section (RCS) data.
Overview
The application macro computes and displays an inverse synthetic aperture radar (ISAR) image from
the backscattered radar cross-section (RCS) data over frequency and angle.
Windowing functions, resampling and extracting the local maxima positions to file are supported. The
application macro requires a far field RCS request or far field RCS data imported from a file containing
far field data (.ffe file). Data imported from a .ffe file must be added manually to a 3D view.
Preconditions
Backscattered RCS data must be in the UV plane of the local coordinate system or the XY plane of the
global coordinate system. Equivalently, the backscattered RCS should be requested in the cut.
Sub-ranges of the total angular phi range could be used to change the viewing angle.
It is assumed that the image will be centred at the global origin. Thus the down-range (X) and cross-
range (Y) should be sufficient to include all the structures / targets relative to the global origin.
The extent and resolution of the ISAR image depend on the frequency range, the number of
frequencies, angular range and number of incident angles.
The CADFEKO CF_ISAR_Resolution.lua application macro is used to calculate the required frequency
and angular ranges and associated sampling to create a 3.5 m x 3 m (down-range x cross-range)
ISAR image with 0.05 m resolution in both down-range (X) and cross-range (Y). The application macro
calculates the required start frequency (f1), end frequency (f2) and the number of frequency samples
(nf) for a given centre frequency (f0). It also calculates the required start angle (phi1), end angle (phi2)
and angular increment (dphi). The application macro assumes that the polar to Cartesian interpolation
will be applied to the RCS data.
The XY_plane solution configuration uses the above-mentioned angular and frequency variables to
create a plane wave excitation (XY_v0) in the global coordinate system with the ISAR view direction
along the X axis.
The UV_plane solution configuration uses the same angular and frequency variables in its plane wave
excitation (UV_v0), but the excitation is defined using a local workplane. The workplane is obtained
from the XY workplane, rotating it by 20° around the local N axis and then by 10° around the local V
axis.
Each configuration contains a far field request set to calculate fields in the plane wave incident direction.
Figure 2: Top view of the XY_plane solution configuration (left) and normal view of
the UV_plane solution configuration (right) showing the positions of the three spheres
with the plane wave excitation. The yellow frame dimensions are 3.5 m x 3 m.
• Dirichlet
•
Bartlett
• Blackman
Interpolate to Cartesian grid The 2D FFT requires that the X and Y distances should be
uniformly spaced. The basic ISAR implementation assumes that
the angle range is small enough such that the polar sampling
approximates a Cartesian grid. The small angle range assumption
can be relaxed by transforming the polar grid to Cartesian. This
reduces smearing and shift of the scattering points and increasing
accuracy of the resulting ISAR image.
Image offset Offset distance of the ISAR image above the plane wave
workplane.
Create 3D view Generates a normal view of the partially transparent ISAR image
with RCS data displayed in dB.
Find local maxima Calculates the specified number of local maxima in descending
order in the ISAR image.
Export local maxima to CSV Export the specified number of local maxima to a CSV file. Find
local maxima has to be checked for this to work.
After accepting the input from the main form, the application macro will process the RCS data and
display a dialog to select the viewing angle from the list. The default viewing angle is the average of the
start angle and end angle.
Note: A different viewing angle could reduce the maximum angular data range.
Figure 4: The Angle selection dialog for specifying the viewing angle.
With the viewing angle specified, the application macro will calculate the maximum angular range that
can be used. The user can then select an angular range from the list.
Figure 5: The Angle selection dialog for specifying the angular range.
The ISAR images associated with the XY_plane and UV_plane solution configurations are shown
below. For this example, a resampling factor of 2 was used with default values for the rest of the input
parameters.
The ISAR dataset can also be used to generate 2D graphs, such as the image below showing the RCS
versus X position corresponding to each sphere’s Y coordinate.
Figure 7: RCS versus X position for the Y coordinates of points A1, A2 and A3.
Example 2: Missile
In this example the view direction and cross-range resolution of the ISAR image is changed by using a
subset of the angular data range. Some of the script's post-processing options are demonstrated.
CADFEKO Model
The physical optics (PO) high-frequency approximation is used to obtain the backscattered RCS for
the model Missile_RCS_Vpol_PO.cfx. The RCS data is computed from phi1 = 18° to phi2 = 57° with
(1/3)° increments (dphi) and over a frequency range from 8 GHz (f1) to 12 GHz (f2) with 81 samples
(nf). Assuming a minimum angular range of 24° the view direction angle can be changed between 30°
and 45°.
POSTFEKO Results
The images below show some results obtained by using different post-processing options available in
the script.
The first ISAR image shows the result obtained without any post-processing of the RCS data for a
view angle phi0 = 45° and an angular range of 24°. The down-range and cross-range resolutions are
0.0370 m and 0.0356 m, respectively.
The angular smear away from the image origin can be reduced by transforming the polar data
to Cartesian coordinates using bi-linear interpolation. The ISAR image with a polar to Cartesian
transformation for a view angle phi0 = 45° and angular range of 24° is shown. The down-range and
cross-range resolutions are 0.0382 m and 0.0435 m, respectively.
The image can be processed further by windowing the RCS data. The next ISAR image is the result
of using a polar to Cartesian transformation and a Hamming window for a view angle phi0 = 45°
and angular range of 24°. The down-range and cross-range resolutions are 0.0382 m and 0.0435 m,
respectively.
Figure 11: ISAR image with a polar to Cartesian transformation and using a Hamming window.
The application macro can be used to change the view angle and the angular range. The ISAR image
below is obtained using a polar to Cartesian transformation and a Hamming window for a view angle
phi0 = 37.33° and an angular range of 38.67°. The down-range and cross-range resolutions are
0.0408 m and 0.0265 m, respectively.
Scripts or application macro can be executed from CADFEKO or POSTFEKO (some scripts work in both
applications). It is easy to add a script or application macro to the macro library so that it is readily
available from the ribbon.
• If the script is already open in the script editor, it could be easier to add the script to the
macro library from the script editor interface. In the script editor, click the icon.
2. Under Macro properties, next to the Script location field, click Browse and specify the main
script file that you want to add to the macro library.
3. [Optional] Under Description, add a short comment to describe the script.
4. In the Label field, enter a label for the script.
5. In the Icon drop-down list, select an icon or specify a custom icon to help identify the script in the
macro library.
6. Click Create to add the script to the macro library and to close the dialog.
The script is added to the macro library. As an example, the macro library menu shows the user-
defined “Tile2Windows” and “Phase center calculator” macros.