<?xml version="1.0" encoding="UTF-8"?>

<?xml-stylesheet type="text/xsl" href="/static/oaitohtml.xsl"?>

<!--
<?xml-stylesheet type="text/xsl" href="/oaitohtml.xsl"?>
-->

<OAI-PMH xmlns="http://www.openarchives.org/OAI/2.0/" xsi:schemaLocation="http://www.openarchives.org/OAI/2.0/ http://www.openarchives.org/OAI/2.0/OAI-PMH.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
    <responseDate>2026-10-11T14:26:10Z</responseDate>
    <request verb="GetRecord" metadataPrefix="oai_dc" identifier="10.57760/sciencedb.space.03783" >https://www.scidb.cn/oai</request>
<GetRecord>
    <record>
    <header >
    <identifier>10.57760/sciencedb.space.03783</identifier>
    <datestamp>2026-07-13T17:15:02Z</datestamp>
</header>
    <metadata>
        
<oai_dc:dc xmlns:oai_dc="http://www.openarchives.org/OAI/2.0/oai_dc/" xmlns:dc="http://purl.org/dc/elements/1.1/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.openarchives.org/OAI/2.0/oai_dc/ http://www.openarchives.org/OAI/2.0/oai_dc.xsd">
  <dc:date>2026-07-13</dc:date>
  <dc:title>AC-CS: Azimuth Compression based Compressed Sensing for MMW SAR Imaging</dc:title>
  <dc:identifier>doi:10.57760/sciencedb.space.03783</dc:identifier>
  <dc:language>en</dc:language>
  <dc:description>Overview--------This project implements a Synthetic Aperture Radar (SAR) imaging pipeline thatintegrates azimuth compression with Compressed Sensing (CS) techniques toreconstruct high-quality SAR images from randomly downsampled raw data. Thecode is written in MATLAB and operates on FMCW (Frequency-ModulatedContinuous-Wave) radar data acquired at 77 GHz.Workflow--------The processing pipeline consists of the following major stages:1. Data Loading and Preprocessing. Raw 3D radar data is loaded from NSSC.mat,&amp;nbsp;&amp;nbsp;permuted to the correct dimensional order, and downsampled by a factor of 4&amp;nbsp;&amp;nbsp;along the slow-time (pulse) dimension. Key geometric parameters &amp;mdash; such as&amp;nbsp;&amp;nbsp;the spatial sampling intervals (dx, dy), the round-trip delay offset (tI),&amp;nbsp;&amp;nbsp;and the nominal range (z0) &amp;mdash; are defined at this stage.2. Range Compression. A Range-FFT is applied to the raw data along the&amp;nbsp;&amp;nbsp;fast-time axis. The range bin corresponding to the target distance z0 is&amp;nbsp;&amp;nbsp;extracted via index calculation based on the chirp slope K, sampling period&amp;nbsp;&amp;nbsp;Ts, and the round-trip delay, thereby performing range focusing.3. Random Downsampling. To simulate sparse or incomplete data acquisition, a&amp;nbsp;&amp;nbsp;subset of rows (vertical/cross-range direction) is randomly discarded &amp;mdash;&amp;nbsp;&amp;nbsp;only 30% of the original rows are retained. This mimics scenarios where&amp;nbsp;&amp;nbsp;entire scan lines are missing, such as in sub-Nyquist sampling or sensor&amp;nbsp;&amp;nbsp;failures.4. Azimuth Compression and Data Recovery via ADMM. The core of the algorithm&amp;nbsp;&amp;nbsp;recovers the missing cross-range samples using the Alternating Direction&amp;nbsp;&amp;nbsp;Method of Multipliers (ADMM). The observation model combines a DFT basis,&amp;nbsp;&amp;nbsp;an azimuth matched-filter operator derived from the Doppler rate Ka, and a&amp;nbsp;&amp;nbsp;row-selection matrix Phi that encodes which rows are observed. The&amp;nbsp;&amp;nbsp;optimization problem minimizes a composite objective with an L1 sparsity&amp;nbsp;&amp;nbsp;penalty and a total-variation (first-order difference) regularization term.&amp;nbsp;&amp;nbsp;The ADMM solver iteratively updates the primal variable x, two auxiliary&amp;nbsp;&amp;nbsp;variables z1 and z2 via soft-thresholding, and the corresponding dual&amp;nbsp;&amp;nbsp;variables u1, u2. The linear system in the x-update is solved efficiently&amp;nbsp;&amp;nbsp;using a precomputed Cholesky factorization.5. SAR Image Formation. The recovered full-aperture data is transformed into&amp;nbsp;&amp;nbsp;the spatial frequency (k-space) domain via a 2D FFT, multiplied by a&amp;nbsp;&amp;nbsp;matched phase factor that accounts for spherical wavefront propagation, and&amp;nbsp;&amp;nbsp;then transformed back to the spatial domain via a 2D IFFT to produce the&amp;nbsp;&amp;nbsp;final SAR image.Files-----&amp;nbsp;demo_AC_CS.m&amp;nbsp;&amp;nbsp;&amp;nbsp;Main script implementing the full SAR imaging pipeline with CS-based&amp;nbsp;&amp;nbsp;&amp;nbsp;data recovery.&amp;nbsp;first_diff_matrix.m&amp;nbsp;&amp;nbsp;&amp;nbsp;Utility function that constructs an (n-1) x n first-order difference&amp;nbsp;&amp;nbsp;&amp;nbsp;matrix used for total-variation regularization in the ADMM solver.&amp;nbsp;NSSC.mat&amp;nbsp;&amp;nbsp;&amp;nbsp;Input data file containing the raw 3D FMCW radar data (sarData).Key Parameters--------------&amp;nbsp;f0&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;77 GHz&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;Start frequency (after ADC offset&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;correction)&amp;nbsp;K&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;63.343 THz/s&amp;nbsp;&amp;nbsp;Chirp slope&amp;nbsp;fS&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;9.121 MHz&amp;nbsp;&amp;nbsp;&amp;nbsp;ADC sampling rate&amp;nbsp;z0&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;245 mm&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;Nominal target range&amp;nbsp;nFFT_time&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;512&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;FFT points for range dimension&amp;nbsp;nFFT_space&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;1024&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;FFT points for spatial dimensions&amp;nbsp;lambda1&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;0.018&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;L1 sparsity regularization weight&amp;nbsp;lambda2&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;0.1&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;Total-variation regularization weight&amp;nbsp;rho1, rho2&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;1e-5&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;ADMM penalty parameters&amp;nbsp;Max iterations&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;800&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;ADMM maximum iterationsHelper Function: first_diff_matrix(n)--------------------------------------Constructs an (n-1) x n first-order difference matrix D such that for a vectorv in R^n, the product D * v yields the forward differences v(i+1) - v(i) fori = 1, ..., n-1. This matrix is used in the total-variation regularizationterm ||D x||_1 within the ADMM framework to promote piecewise smoothness inthe recovered signal.</dc:description>
  <dc:subject>Compressive sensing (CS); millimeter-wave (MMW) radar; near-field imaging; synthetic aperture radar (SAR)</dc:subject>
  <dc:creator>Ruixiang Fang</dc:creator>
  <dc:creator>Xiaojin Shi</dc:creator>
  <dc:creator>Yunhua Zhang</dc:creator>
  <dc:rights>PUBLIC</dc:rights>
  <dc:rights>https://api.github.com/licenses/mit</dc:rights>
  <dc:type>dataset</dc:type>
  <dc:relation>http://www.doi.org/10.1109/lsens.2026.3712770</dc:relation>
  <dc:publisher>Science Data Bank</dc:publisher>
</oai_dc:dc>

    </metadata>
</record>
</GetRecord>
</OAI-PMH>