{"id":"dccd2837-82bc-4026-a2e6-f44fc00c86a8","arxiv_id":"2506.08620","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"An extended Spherical Geometry Algorithm with instantaneous Doppler centroid removal and matched filtering processes stripmap and TOPS spaceborne SAR data without azimuth aliasing.","lead":"The paper extends a spotlight-only synthetic aperture radar imaging algorithm (SGA) to stripmap and TOPS modes by adding a Doppler centroid removal step and a matched-filter azimuth compression. This could make the spherical-geometry algorithm, which handles curved orbit and Earth shape, usable for wider-swath imaging modes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The deramp rate k_t is fixed to a scene-center value without resolving whether the SGA's azimuth-frequency scale is the orbit radius or the target slant range, so the central alias-free- and well-focused claim rests on an unproven range-invariance assumption.","rationale":"The reader's weakest-assumption correctly identifies the single-k_t deramp as a fragile point, and I agree that the paper does not analyze its space-variance. My concern is slightly more specific: the paper's Eq. (11) derives the azimuth frequency using an R_e from the spherical geometry of Fig. 1, while Section III.C defines k_t using R_scene, with no explicit identification of these two quantities. This ambiguity is load-bearing because the central claim is that the extended SGA is alias-free and well-focused. If R_e is the orbit radius, then the derived k_t is wrong by an order of magnitude and the entire explanation of why the original SGA aliases is questionable; if R_e is the scene range, then k_t genuinely varies across the swath, and the single-reference-range deramp and matched filter may produce range-dependent scaling or residual phase error. Either way, the paper's argument is incomplete. The real-data figures are suggestive but do not settle the question: they show no PSLR/IRW measurements, no geolocation check, and no comparison with a standard spaceborne processor, so a visual inspection cannot distinguish a correct algorithm from one with a spatially varying distortion. A point-target simulation in the true spherical geometry is the natural and decisive check, and it should be a condition of acceptance. Because the concern is addressable by additional analysis and simulation rather than by a demonstrated contradiction in the presented results, the reader's CONDITIONAL verdict remains appropriate.","tokens_in":16204,"tokens_out":28281,"duration_ms":367793,"concrete_test":"Run a point-target simulation using the full spherical geometry underlying Eqs. (2)-(8), with targets placed at near, mid, and far slant range across the Table I stripmap and TOPS swaths. Process the simulated raw data with the proposed algorithm using a single mid-swath k_t, and separately using a per-target k_t(R_t). For each target, measure the azimuth peak position against the true x_t, the azimuth impulse-response width, and the peak sidelobe ratio. If the edge-target position shift exceeds one resolution cell or PSLR degrades by more than 1 dB relative to the mid-swath target, the single-k_t assumption is load-bearing and a range-dependent deramp correction is required. Also check whether the original SGA aliases in the same simulation; this will disambiguate whether the correct R_e is the orbit radius or the scene slant range.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III.C defines a single deramp rate k_t for all scatterers: k_t = 2v^2/(lambda*R_scene) for stripmap and k_t = 2v^2/lambda*(1/R_scene + 1/R_centre) for TOPS, introduced near Eq. (21), then uses that same rate in the matched filter Eq. (30). However, the azimuth frequency obtained after SGA range preprocessing in Eq. (11) is f_a = (2v/(lambda*R_e))*x_t, where R_e comes from the spherical geometry of Fig. 1. The paper never states that R_e equals the scene-center slant range R_scene. If R_e is the radar-to-Earth-center distance (roughly 6900 km at the stated 532 km altitude), the azimuth frequency span for the 30 km stripmap scene is only about 1.2 kHz, below the 3900 Hz PRF, so the original SGA should not alias as shown in Fig. 9(a). If R_e is instead the scene slant range of 597 km, the span is about 13.6 kHz and the aliasing analysis applies, but then the effective R_e varies with slant range across the 20 km stripmap and 36 km TOPS swaths, so k_t is range-dependent. The paper uses a single reference-range k_t and does not analyze the resulting mismatch. Depending on which interpretation is correct, this ambiguity either contradicts the motivating aliasing observation or leaves a space-variant deramp/matched-filter error at swath edges. The real-data validation in Figs. 9-11 contains no quantitative focus or geolocation metrics and no comparison against a reference processor, so it cannot discriminate between these cases.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the Spherical Geometry Algorithm (SGA), originally developed for spotlight SAR, to stripmap and TOPS modes. The authors first analyze why the original SGA fails in these modes: after SGA range preprocessing the azimuth signal becomes a set of single-frequency bursts whose instantaneous Doppler centroid sweeps linearly in time, so the total azimuth bandwidth can exceed the PRF, causing aliasing in azimuth resampling and in the final azimuth IFFT. The proposed extended SGA adds an instantaneous Doppler centroid removal step (a range-frequency-dependent deramp) before the azimuth resampling, and replaces the final spectral analysis with a matched filter. The authors validate the algorithm with real data from the Chaohu-1 satellite in both stripmap and TOPS modes, showing images produced by the classical SGA and by the proposed method.","tokens_in":16556,"tokens_out":22394,"duration_ms":256976,"significance":"If the proposed algorithm is correct, it is a useful practical extension of the SGA to two important wide-swath modes, and the two modifications (deramp before azimuth resampling and matched filtering instead of IFFT) are conceptually interesting. The paper's central idea is independent of the prior SGA publication in that it introduces new phase compensation and a new azimuth compression scheme rather than merely rearranging known steps. The real-data demonstration shows a clear qualitative improvement over the unmodified SGA. However, the validity of the central claim is not fully established: a key geometric parameter (R_e in Eq. 11) is left undefined, the range invariance of the deramp rate is assumed without analysis, and the experimental validation lacks quantitative focus and geolocation metrics as well as comparison with an independent reference processor.","major_comments":[{"comment":"The parameter R_e in Eq. (11) is never defined. If R_e is the Earth-center-to-radar distance (about 6,900 km for the quoted orbit), the azimuth frequency span after SGA range preprocessing for the 30 km stripmap scene is roughly 1.2 kHz, well below the 3.9 kHz PRF, which would contradict the aliasing behavior shown in Fig. 9(a). If R_e is instead the scene-center slant range (597 km), the span is about 13.6 kHz and the aliasing analysis applies, but then the relationship between R_e and R_scene must be stated explicitly. Moreover, in that interpretation R_e effectively varies with slant range across the 20 km stripmap and 36 km TOPS swaths, so the single deramp rate k_t = 2v^2/(λ R_scene) (or the TOPS counterpart) is range-dependent; the paper does not analyze the resulting space-variant deramp/matched-filter error at the swath edges, which is load-bearing for the alias-free and well-focused claim.","section":"Section III.C, Eq. (11) and Eq. (21)"},{"comment":"The experimental validation is entirely qualitative. The paper reports no impulse-response measurements (azimuth/range resolution, peak sidelobe ratio, integrated sidelobe ratio), no geolocation accuracy, and no comparison with a reference processor such as backprojection or a standard Range-Doppler algorithm. The central claim that scatterers are 'well focused' and 'accurately focused at their true positions' in Figs. 9 to 11 cannot be verified from the displayed images alone, especially since the TOPS image in Fig. 11(b) does not include an enlarged point-target view.","section":"Section IV.A and IV.B"},{"comment":"The derivation explicitly restricts itself to the basic SGA formulation that neglects non-coplanar effects induced by Earth's rotation, and states that the extension to the enhanced formulation is direct. However, the experimental data are collected by a LEO satellite at roughly 532 km altitude, where orbital curvature and Earth rotation are not negligible over the quoted coherent times. The paper does not provide the extended derivation of the new deramp and matched-filter steps under the non-coplanar geometry, nor does it quantify the residual phase errors for the experimental scenario; this is a load-bearing omission because the real-data validation is performed in the full geometry while the algorithm analysis is carried out in the simplified geometry.","section":"Section II, first paragraph, and Section III"},{"comment":"The derivation of the modified deramp function and the claim that it avoids new coupling after the Keystone azimuth resampling are not fully transparent; the equations as printed are ambiguous. The authors should provide a clean step-by-step derivation showing that after substituting the Keystone transform, the quadratic phase term is independent of range frequency, and should explicitly justify that the LFM signal after deramping has a total bandwidth not exceeding the PRF, so that matched filtering can be performed without additional oversampling.","section":"Section III.C, Eq. (24) and Eq. (25)"}],"minor_comments":[{"comment":"There are several grammatical errors, e.g., 'can't directly applied' in the abstract, and 'an groundbreaking approach' in Section I; these should be corrected.","section":"Abstract and Section I"},{"comment":"The sentence describing the signal frequency after range processing is garbled: 'the signal frequency is /t s tf k x v= −' should be written as a clear formula, e.g., f = -k_s x_t/v, with k_s defined in terms of the geometry parameters.","section":"Section III.C, text near Eq. (21)"},{"comment":"The parameter 'Pulse Repeat Frequency' should be 'Pulse Repetition Frequency'.","section":"Table I"},{"comment":"The sentence 'the data acquiring time is 0.8s, which corresponding to a synthetic aperture length of 6135m' has a subject-verb agreement error and should also state the platform velocity used for the conversion.","section":"Section IV.B"},{"comment":"The claim that matched filtering 'improve computational efficiency' is not quantified; a brief complexity comparison with the oversampling approach or with the original SGA would be helpful.","section":"Section III.C, matched filtering discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a plausible extension of the SGA, but the central claim is currently under-supported because of the undefined R_e and the assumed range invariance of k_t, and because the experimental validation lacks quantitative metrics. The equations are also difficult to read due to numerous transcription errors. I recommend major revision with a request to clarify the geometry, provide a range-dependence analysis or a justification of range invariance, and add quantitative validation. The reliance on the author's earlier SGA paper [20] is acceptable for an extension, though the self-citation pattern should not obscure the novelty of the new deramp and matched-filter steps."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Xinhua Mao and colleagues have extended their spherical geometry algorithm (SGA) from spotlight to stripmap and TOPS modes. The two changes are straightforward: add an instantaneous Doppler centroid removal before the azimuth resampling, and replace the final spectral analysis with a matched filter. That combination is new in the SGA context, and the real data from Chaohu-1 show that it fixes the aliasing that plagues the unmodified SGA. The TOPS result in particular looks much cleaner. Credit where due: the problem identification is correct—after SGA range preprocessing, the azimuth signal bandwidth can exceed the PRF in stripmap/TOPS—and the modified deramp of Eq. (24) is a sensible way to avoid a new coupling term in the Keystone transform. The paper is readable and the flow is logical.\n\nThe soft spots are mostly about things left unsaid. First, Eq. (11) drops in R_e without a definition. The stress-test note worries this could be the orbital radius, which would kill the aliasing claim. It doesn't hold up: if R_e were 6900 km, the original SGA would not alias at 3900 Hz PRF, but Fig. 9(a) clearly shows aliasing. So R_e has to be the scene slant range. That said, the paper should define it properly.\n\nSecond, the real concern: the deramp rate k_t is a single value computed at scene center. For the 20 km and 36 km swaths, the slant range changes by tens of km, which changes the true Doppler rate and the azimuth frequency scaling. The paper does not analyze the resulting defocus at swath edges, and the figures are qualitative. There are no impulse response measurements, no sidelobe numbers, no geolocation errors, and no comparison against a standard processor like Range-Doppler or chirp scaling. So while the images look good, the range-invariance assumption is unverified. This is not necessarily fatal—the real data may simply be in a regime where the variation is tolerable—but it needs to be addressed.\n\nThe citation pattern is fine; this is a direct extension of the author's own SGA, and the earlier paper is cited. The math is mostly coherent, though the presentation in places is rough (some symbols are dropped in the OCR, and the \"endomorphism property\" in the abstract is never explained).\n\nWho should read this: practitioners working on spaceborne SAR image formation who care about handling curved orbits and spherical Earth without giving up stripmap/TOPS operation. It is not a fundamental advance, but it is a useful incremental one.\n\nYes, it deserves peer review. A serious referee should ask for quantitative validation and an analysis of the range-dependence of k_t. With those additions, this could be a solid journal paper.","headline":"A plausible extension of the SGA to stripmap/TOPS with real-data support, but the neglected range-dependence of the deramp rate and missing quantitative metrics keep it from being a slam dunk.","tokens_in":17060,"tokens_out":8637,"would_cite":true,"duration_ms":99909,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The spherical geometry algorithm, originally a spotlight SAR method, is extended to stripmap and TOPS modes by removing the instantaneous Doppler centroid before azimuth resampling and by compressing the azimuth signal with matched…","keywords":["synthetic aperture radar","spherical geometry algorithm","stripmap mode","TOPS mode","spaceborne SAR","azimuth resampling","Doppler centroid removal","matched filtering"],"falsifier":"Take a stripmap or TOPS data set and compute $k_t$ from the orbit at the near and far edges of the swath. If the residual quadratic phase after applying the reference-range deramp exceeds about $\\pi/4$ over the aperture for edge scatterers, the matched filter will broaden those targets; measuring their azimuth impulse-response width or contrast and showing it grows with range offset would falsify the single-rate assumption.","tokens_in":16007,"feed_emoji":"🛰️","tokens_out":6744,"duration_ms":80235,"temperature":0.7,"pith_summary":"The paper claims that the spherical geometry algorithm (SGA), a SAR image-formation method originally built for spotlight mode, can be made to work for stripmap (fixed-beam) and TOPS (beam-swept wide-swath) modes without giving up its curved-orbit, spherical-Earth advantages. In those modes, the azimuth signal occupies a total bandwidth much larger than the pulse repetition frequency, so the original algorithm's azimuth resampling and final Fourier transform alias scatterers at the scene edges. The authors add an instantaneous Doppler-centroid removal step before azimuth resampling and replace the final spectral-analysis step with a matched filter. Real stripmap and TOPS data then focus cleanly across the demonstrated 20 km to 36 km swaths, where the original SGA produces aliased and defocused edge targets.","feed_headline":"Spherical-geometry SAR now handles stripmap and TOPS modes","feed_subtitle":"Adding Doppler-centroid removal and matched filtering keeps scene edges sharp without oversampling.","key_machinery":"The load-bearing object is the instantaneous Doppler centroid and its removal before the Keystone azimuth resampling. The paper multiplies the data by $\\exp(-j\\pi k_t t_a^2 (f_c+f_r)/f_c)$, with the range-frequency factor chosen so that the deramp commutes with the resampling operation and no new two-dimensional coupling appears. It then deliberately leaves the resulting quadratic azimuth phase in place and compresses it with a matched filter whose reference is $\\exp(-j\\pi f_a^2/k_t)$ after the azimuth FFT. This combination keeps the total azimuth bandwidth inside the pulse repetition frequency during interpolation while avoiding the image-domain aliasing that would follow if the quadratic term were removed before a final IFFT.","core_discovery":"On the paper's own terms, the central discovery is that the failure of the original SGA outside spotlight mode has a single cause and a two-part cure. The cause is that, after SGA's range preprocessing, the azimuth signal of each scatterer becomes a single-frequency tone, and in stripmap and TOPS modes the total bandwidth of these tones across the scene exceeds the pulse repetition frequency. The cure is to remove the time-varying Doppler centroid before azimuth resampling, and then to compress the deliberately retained quadratic phase by matched filtering instead of by an azimuth FFT followed by peak detection. The paper derives the compensation rate $k_t = 2v^2/(\\lambda R_{\\text{scene}})$ for stripmap and $k_t = 2v^2/\\lambda\\,(1/R_{\\text{scene}} + 1/R_{\\text{centre}})$ for TOPS, adjusts the deramp so that the Keystone resampling introduces no new range–azimuth coupling, and shows on measured satellite data that edge scatterers that alias under the original SGA are fully focused by the extended version.","pith_inferences":["Implicit in the paper, but not tested, is how $k_t$ varies with slant range across very wide swaths; its single reference-range value may need to become range-dependent for swaths well beyond the 20 km to 36 km demonstrated.","The same Doppler-centroid-removal-plus-matched-filter recipe should apply to sliding spotlight and squinted stripmap geometries, where the beam steering lies between the modes analyzed here.","A quantitative metric, such as azimuth impulse-response width or image contrast at the swath edges before and after the fix, would sharpen the visual comparison the paper provides."],"forward_implications":["Stripmap and TOPS data can be processed by the same geometric imaging chain that handles spotlight mode, preserving the curved-orbit and spherical-Earth corrections.","Scene-edge scatterers, which alias and defocus under the original SGA, focus correctly without increasing the pulse repetition frequency or oversampling the data.","The matched-filter azimuth compression avoids enlarging the data volume, so the extended algorithm is computationally cheaper than an oversampling-based fix.","The same two modifications should carry over to the non-coplanar, Earth-rotation version of SGA, which the paper notes is a direct extension.","The extended SGA broadens the applicability of the spherical geometry algorithm to the standard operational modes of modern spaceborne SAR systems."],"supporting_citations":[{"why":"Supplies the original spherical geometry algorithm for spotlight SAR that the paper extends to stripmap and TOPS modes.","marker":"[20]"},{"why":"Provides the spotlight SAR polar-format processing framework whose azimuth resampling and final spectral analysis are the steps being modified.","marker":"[4]"},{"why":"Supplies the spotlight-mode SAR signal processing formulation behind the original algorithm's azimuth compression.","marker":"[5]"}],"fun_headline_variants":["SGA extended to stripmap and TOPS without aliasing","Doppler-centroid removal extends SGA to stripmap and TOPS","Matched filtering replaces FFT in SGA for stripmap and TOPS","SGA broadened: stripmap and TOPS focused with matched filtering"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fix assumes that the azimuth signal's center frequency drifts at one fixed rate for all targets across the whole swath. If that rate changes with distance from the satellite, the correction is slightly wrong at the swath edges and those targets will defocus.","fun_headline_variants_meta":{"raw":{"variants":["SGA extended to stripmap and TOPS without aliasing","Doppler-centroid removal extends SGA to stripmap and TOPS","Matched filtering replaces FFT in SGA for stripmap and TOPS","SGA broadened: stripmap and TOPS focused with matched filtering"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000928,"raw_usage":{"total_tokens":3989,"prompt_tokens":973,"completion_tokens":3016,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":2934}},"tokens_in":589,"tokens_out":3016,"duration_ms":25316,"temperature":1.0,"reasoning_tokens":2934,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:06:59.464620+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a stripmap or TOPS data set and compute $k_t$ from the orbit at the near and far edges of the swath. If the residual quadratic phase after applying the reference-range deramp exceeds about $\\pi/4$ over the aperture for edge scatterers, the matched filter will broaden those targets; measuring their azimuth impulse-response width or contrast and showing it grows with range offset would falsify the single-rate assumption.","supporting_citations":[{"cited_title":"Spherical geometry algorithm for spaceborne synthetic aperture radar imaging,","cited_arxiv_id":null,"evidence_quote":"Supplies the original spherical geometry algorithm for spotlight SAR that the paper extends to stripmap and TOPS modes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the spotlight SAR polar-format processing framework whose azimuth resampling and final spectral analysis are the steps being modified."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the spotlight-mode SAR signal processing formulation behind the original algorithm's azimuth compression."}],"review_version":1}