{"id":"0e2ef437-88ec-4bd6-84a6-919a8862f529","arxiv_id":"2506.07686","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A fast numerical method computes the minimum detectable ship radar cross-section for spaceborne SAR by averaging detection probability over a lognormal fluctuation model and using binary integration.","lead":"This letter replaces slow Monte Carlo simulations with fast numerical formulas for deciding whether a spaceborne radar can spot small ships at sea, and applies the method to a very-low-orbit satellite example. It gives engineers a quick feasibility check based on the minimum radar cross-section a ship must have to be detected.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Detection window smaller than ship: n=Nps^w overstates binary-integration gain and may make the claimed feasibility numbers optimistic.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the worked example's detection window is smaller than the target ship, so the n used in the binary-integration expression Eq. (17) is too large. The method itself is otherwise standard: Eq. (16) is a correct Bayesian averaging of the Marcum Q-function over the lognormal SNR PDF, and Eqs. (17)–(18) are valid incomplete-beta formulations. The central empirical claim, that the VLEO X- and Ku-band systems can detect the specified ship, depends directly on the resulting minimum detectable sigma0 values being inside the typical range from [1]. Since the geometry inconsistency inflates the binary-integration gain, the claimed feasibility is optimistic as written. The concern is concrete and testable, not a disagreement with consensus, and it does not by itself invalidate the method if the example is corrected. Therefore the reader's CONDITIONAL verdict remains appropriate; no verdict change is needed, but the manuscript should resolve the Nps^w definition and rerun the example before acceptance.","tokens_in":8184,"tokens_out":8414,"duration_ms":109839,"concrete_test":"Recompute the X- and Ku-band curves in Fig. 4 using Section V parameters and the optimization in Eq. (20), but replace n = Nps^w with the true maximum number of ship pixels that can lie inside a 6 m × 6 m TDW for a 12 m × 4 m ship. One robust way is to grid-search TDW placement and ship orientation over the 12 m × 4 m rectangle and count resolution cells (δaz × δgr) whose centers fall inside the window; also compute the simple upper bound n = floor(Lw^2/A_res). Keep Pfa = 1e-14 and re-optimize m. If the resulting min sigma0 at δr = 0.25 m stays within the 2 dB typical upper bound, the feasibility conclusion survives and the issue is a documentation fix; if it exceeds 2 dB, the headline detectability claim fails as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III sizes the TDW as Lw = min{Lship} and Eq. (13) converts it to pixels, while Fig. 3 block 8 uses this window to bound the number of ship pixels available to the m-of-n decision. Section V sets Lship = 12 m, Wship = 4 m, and Lw = 6 m. A 12 m × 4 m ship cannot be contained in a 6 m × 6 m window in any orientation, so the maximum number of ship-associated pixels that can fall in one TDW is not Nps = floor(A_ship/A_res) = floor(48/1.299) = 36; it is at most floor(Lw^2/A_res) = 27 and, accounting for the 4 m width, closer to 18–24 pixels. Nowhere is Nps^w explicitly distinguished from Nps, so Eq. (17) appears to use n = 36. For Pfa = 1e-14 the optimal m remains 2 as n varies, but P_D ≈ C(n,2)p^2, so using n = 36 instead of n ≈ 20 raises the required pixel-level P_d by a factor of roughly 1.8 and shifts the required mean lognormal SNR upward. Points A and B in Fig. 4 sit at 1.07 dB and -0.80 dB, within but near the upper edge of the [1] typical sigma0 range (-1 to 2 dB); a large enough shift would flip the feasibility conclusion.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an analytical method for evaluating the feasibility of ship detection by spaceborne SAR in the early design stage. The method computes the pixel-level probability of detection by numerically integrating Marcum's Q-function over a lognormal SNR distribution (Eq. 16), then converts pixel-level probabilities to ship-level detection and false-alarm probabilities using the regularized incomplete beta function (Eqs. 17-18). It also proposes a rule for selecting the binary-integration threshold m. The method is validated against Monte Carlo simulation in Fig. 2 and is applied to a VLEO SAR example, yielding minimum detectable RCS values of 61.39 m^2 (X band) and 39.89 m^2 (Ku band) at a slant-range resolution of 0.25 m. The paper claims these values fall within the typical ship backscattering coefficient range from the literature, so the target is deemed detectable.","tokens_in":8546,"tokens_out":12422,"duration_ms":142085,"significance":"If the method is sound, it provides a fast, numerically stable alternative to Monte Carlo simulations for early-stage feasibility assessments, and the use of the incomplete beta function avoids the binomial-coefficient overflow that can occur for large pixel counts. The Monte Carlo validation in Fig. 2, with relative errors below 1%, is a genuine strength. The method yields falsifiable predictions (minimum detectable RCS) that can be checked against measured backscatter statistics. However, the worked example contains a geometric inconsistency that affects the numerical feasibility conclusion, and a key quantity Nps^w is never defined explicitly. Correcting these points is necessary before the example results can be relied upon.","major_comments":[{"comment":"Section V sets Lship = 12 m, Wship = 4 m, and Lw = 6 m, but Section III defines Lw = min{Lship}; these statements are inconsistent. If the definition were applied, Lw would be 12 m, not 6 m. If Lw = 6 m is intentional, then Nps^w, the maximum number of ship-associated pixels in one TDW, is at most floor(Lw * Wship / Ares) = floor(6*4/1.299) = 18 for a ship aligned with the window, not Nps = floor(Aship/Ares) = 36. The manuscript never gives an explicit formula for Nps^w, although Eqs. (17)-(18) and the procedure in Fig. 3 use n = Nps^w. Because PD ≈ C(n,2) Pd^2 for small Pd, using n = 36 instead of n ≈ 18 overstates the binary-integration gain by a factor of about sqrt(C(36,2)/C(18,2)) ≈ 2.0 in the required pixel-level detection probability, shifting the required mean sigma0 upward. Since point A in Fig. 4 lies at 1.07 dB, near the upper edge of the typical sigma0 range (-1 to 2 dB) reported in [1], a corrected n could move the X-band result outside that range and reverse the claimed feasibility. Please define Nps^w explicitly, use a consistent value of Lw, and recompute the numerical example.","section":"Section V / Section III"},{"comment":"The detection model assumes that all Nps^w ship-associated pixels fall within a single TDW, but the TDWs are described as non-overlapping windows that compactly tile the image. A ship straddling a window boundary will have its pixels distributed among two or four windows, so no single window will contain all Nps^w pixels. The computed probability is therefore a best-case upper bound. The paper should state this explicitly and, ideally, quantify the sensitivity of the feasibility result to ship position relative to the TDW grid. Without this, the reported PD and RCSmin_ship values are optimistic in a way that is not disclosed to the reader.","section":"Section III / Section IV"}],"minor_comments":[{"comment":"The exponent contains a stray parenthesis: it reads \"e^{-[lnχ−α′)]^2/2β^2}\", which should be \"e^{-(lnχ−α′)^2/(2β^2)}\".","section":"Eq. (7)"},{"comment":"The sentence \"These plots are created at Pd = 10−10\" appears to be a typo; it should probably read \"Pfa = 10−10\" since Pd is the plotted quantity.","section":"Section IV, Fig. 2"},{"comment":"Reference [2] spells the author's name as \"Sallivan\"; the correct spelling is \"Sullivan\".","section":"References"},{"comment":"The quantity Nps^w is listed in Block 8 of Fig. 3 but is never defined in the text; please provide an explicit formula, as the numerical results depend directly on it.","section":"Section III, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"The core analytical derivation (Eqs. 16-18) appears sound and the Monte Carlo validation is convincing. The problem is confined to the worked example and the missing definition of Nps^w, which directly affect the claimed feasibility numbers. Once the example is corrected and the sensitivity to the TDW alignment is discussed, the paper would be suitable for publication in a letters journal. I would not recommend reject because the issue is fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core method here is a clean combination of textbook results: average the Marcum Q function over a lognormal SNR PDF to get pixel-level detection probability, then use the regularized incomplete beta function for m-of-n binary integration at the ship level. That saves you from Monte Carlo and handles large n without binomial coefficient blowup. The validation in Fig. 2 against 10^7 Monte Carlo runs is solid, and the whole computation takes seconds. If you work in SAR ship-detection feasibility, this is genuinely handy.\n\nThe soft spot is the worked example, and it is not minor. The paper defines the detection window side as L_w = min{L_ship}, but then sets L_w = 6 m for a 12 m ship. A 12 m × 4 m ship cannot fit in a 6 m × 6 m window in any orientation. That means the number of ship pixels available for binary integration in one window is not N_ps = floor(A_ship/A_res) = 36; it is at most about 18–24. The paper never explicitly distinguishes N_ps from N_ps^w, and the flow chart seems to feed n = N_ps into Eq. (17). Because P_D ≈ C(n,2) p^2 for small p, using n = 36 instead of n ≈ 20 overstates the integration gain by roughly a factor of 1.8. That shifts the required pixel-level P_d upward, which raises the minimum detectable sigma0. The claimed points A and B sit near the upper edge of the typical sigma0 range from [1]; a large enough shift could flip the feasibility conclusion. So the example's RCS_min_ship values are optimistic as presented.\n\nThere are also minor gaps: no sensitivity analysis for β or K_pw, and no code supplied. But the main issue is the geometry inconsistency. It is fixable—either set L_w = 12 m for the ship, or use an orientation-aware estimate of how many ship pixels can actually fall in a 6 m window.\n\nThe mathematics of Eq. (16) is correct, and the overall framework is defensible once the window size is reconciled with the ship dimensions. This is a paper for radar system engineers doing early smallsat SAR feasibility studies. It deserves peer review, but a referee should send it back for correction of the example and a clear statement of how N_ps^w is computed. I would not cite it in its current form, but I would revisit it after a revision.","headline":"Useful repackaging of standard detection math, but the example's detection window is smaller than the ship, so the feasibility numbers are optimistic until fixed.","tokens_in":8981,"tokens_out":3474,"would_cite":false,"duration_ms":39995,"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":"This paper claims that the feasibility of spaceborne SAR ship detection can be assessed by numerically integrating Marcum's Q function over a lognormal SNR distribution and converting pixel-level to ship-level probabilities via the…","keywords":["spaceborne SAR","ship detection","very low Earth orbit","minimum detectable RCS","binary integration","lognormal backscatter model","Marcum Q function","feasibility analysis"],"falsifier":"Recompute the worked example with a target detection window large enough for the 12 m ship, or with the maximum number of ship pixels that a 6 m window can geometrically contain, and check whether the minimum detectable RCS values of $61.39\\,\\mathrm{m^2}$ and $39.89\\,\\mathrm{m^2}$ still fall inside the typical backscatter range from [1]. An independent check is to compare the numerical integral in Eq. (16) against the paper's own Monte Carlo benchmark across the SNR range and $\\beta=2$, with the maximum relative error staying below the claimed 1%.","tokens_in":7996,"feed_emoji":"🛰️","tokens_out":8861,"duration_ms":94362,"temperature":0.7,"pith_summary":"To decide whether a planned space radar can see a given ship on the ocean, the paper proposes a fast numerical criterion based on the smallest radar cross-section the system could detect at a required detection probability and a required false-alarm probability. The previous approach needed long Monte Carlo simulations and an ad hoc approximation; the new method replaces those with two exact formulas: one numerical integral for the per-pixel detection chance, and one regularized incomplete $\\beta$ function for the ship-level 'at least $m$ pixels above threshold' rule. The paper also gives a rule for choosing the optimal $m$, and demonstrates the procedure on a very-low-Earth-orbit SAR at X and Ku bands, finding minimum detectable ship RCS values of $61.39\\,\\mathrm{m^2}$ and $39.89\\,\\mathrm{m^2}$ at $0.25$ m slant-range resolution. These values sit inside the typical backscatter range reported for the ship size in question, so the example concludes that such a radar could detect the target. If the method holds up, it gives mission designers a fast feasibility screen before committing to hardware.","feed_headline":"Faster formula decides if a space radar can spot ships","feed_subtitle":"Detection odds come from one integral and one beta function, replacing slow Monte Carlo runs.","key_machinery":"The load-bearing object is a pair of distribution identities. First, because the ship-pixel SNR is a linear scaling of the ship backscatter coefficient, and the latter is lognormal, the SNR is also lognormal with the same shape parameter $\\beta$; equation (7) gives its density. Second, the detection probability for a fixed SNR, $Q_1(\\sqrt{2\\chi},\\sqrt{-2\\ln P_{fa}})$, is averaged over that density in (16), turning the fluctuation problem into a one-dimensional numerical integral. Third, the $m$-of-$n$ rule is represented by the regularized incomplete $\\beta$ function $I_x(m,n-m+1)$ in (17), which is numerically stable for large $n$ and whose inverse (18) lets the designer go from ship-level requirements back to pixel-level thresholds. These identities are combined in the Fig. 3 loop that chooses the optimal $m$ and solves the constrained optimization for the minimum detectable RCS.","core_discovery":"The central claim is that the whole two-step detection chain—a constant-false-alarm-rate threshold on SAR image pixels followed by an $m$-of-$n$ binary integration inside a target detection window—can be evaluated without Monte Carlo, using two numerically stable formulas. The pixel-level detection probability is $P_d = \\int_0^\\infty Q_1\\bigl(\\sqrt{2\\chi},\\sqrt{-2\\ln P_{fa}}\\bigr)\\, p_{\\Xi_{\\rm sp}}(\\chi)\\,d\\chi$, where $Q_1$ is the first-order Marcum Q function (the detection probability for a non-fluctuating signal with random phase in Gaussian noise) and $p_{\\Xi_{\\rm sp}}$ is the lognormal density (7) of the ship-pixel SNR. The ship-level detection and false-alarm probabilities are then $P_F^{sw} = I_{P_f}(m, N_{ps}^w - m + 1)$, the regularized incomplete $\\beta$ function, which is computable and invertible even for large $N_{ps}^w$. The paper further claims that the optimal decision rule is the smallest $m\\ge 2$ for which the ship-level false-alarm probability stays below the pixel-level false-alarm probability, while the $1$-of-$n$ rule maximizes $P_D$ but inflates false alarms roughly by a factor $n$. Applied to a VLEO example, the method yields minimum detectable ship RCS values of $61.39\\,\\mathrm{m^2}$ (X band) and $39.89\\,\\mathrm{m^2}$ (Ku band) at $\\delta_r=0.25$ m, which fall within the typical $\\bar{\\Sigma}_{\\rm sp}^0$ range of the reference data.","pith_inferences":["The framework could be inverted to set radar parameters—peak power, bandwidth, or altitude—for a required minimum detectable ship RCS, rather than only vetting a fixed design as presented.","The apparent advantage of finer resolution (more pixels in the window, lower required per-pixel detection probability) carries hidden costs in false-alarm burden and data rate that the RCS criterion in (19) does not price.","The VLEO example's reported numbers should be re-derived with a consistent detection-window geometry; with $L_w=6$ m and $L_{\\rm ship}=12$ m, the window cannot hold all ship-associated pixels, so those RCS minima are optimistic.","The lognormal parameters come from a single set of 58 ship targets; testing the method against other ship types and sea states would show how far the feasibility criterion generalizes."],"forward_implications":["Mission designers can screen a candidate SAR configuration in seconds rather than running Monte Carlo simulations; the paper reports about 3.5 seconds per RCS-versus-resolution plot on a desktop.","The optimal-$m$ rule turns the binary-integration threshold from a heuristic into a computed quantity: the smallest $m$ for which ship-level false alarm stays below pixel-level false alarm.","The equations work in both directions: given required $P_D$ and $P_{FA}$, the inverse of the regularized incomplete beta function yields the pixel-level thresholds that the feasibility optimization needs.","Because only the SNR density in (16) encodes the target model, the method extends to any statistical model of SAR image data and any detection algorithm, as the paper notes."],"supporting_citations":[{"why":"Supplies the lognormal ship-backscatter model and its fitted parameters, plus the typical backscatter range used to judge detectability.","marker":"[1]"},{"why":"Provides the radar equation for average ship-pixel SNR used to scale backscatter into SNR.","marker":"[2]"},{"why":"Justifies the change-of-variables step that turns the lognormal backscatter density into the lognormal SNR density in (7).","marker":"[3]"},{"why":"Gives the Marcum Q-function detection probability for a non-fluctuating signal with random phase in Gaussian noise, and the 1-of-n false-alarm approximation.","marker":"[4]"},{"why":"Introduces the regularized incomplete beta function representation of m-of-n binary integration and its inverse, which the paper exploits for numeric stability.","marker":"[5]"}],"fun_headline_variants":["One formula predicts if space radar can spot ships","Beta function plus an integral decides SAR ship detectability","Analytic method replaces Monte Carlo for ship detection SAR","No simulation needed to check SAR ship detection odds","Fast check: can a satellite radar spot that ship"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation's load-bearing premise is that the lognormal model fitted to the reference dataset describes the target's radar brightness, together with the geometric assumption that the 6 m detection window can hold all the pixels of the 12 m ship—an assumption the example's own parameters do not satisfy consistently.","fun_headline_variants_meta":{"raw":{"variants":["One formula predicts if space radar can spot ships","Beta function plus an integral decides SAR ship detectability","Analytic method replaces Monte Carlo for ship detection SAR","No simulation needed to check SAR ship detection odds","Fast check: can a satellite radar spot that ship"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001115,"raw_usage":{"total_tokens":4642,"prompt_tokens":946,"completion_tokens":3696,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":3622}},"tokens_in":562,"tokens_out":3696,"duration_ms":28594,"temperature":1.0,"reasoning_tokens":3622,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:28:15.857224+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the worked example with a target detection window large enough for the 12 m ship, or with the maximum number of ship pixels that a 6 m window can geometrically contain, and check whether the minimum detectable RCS values of $61.39\\,\\mathrm{m^2}$ and $39.89\\,\\mathrm{m^2}$ still fall inside the typical backscatter range from [1]. An independent check is to compare the numerical integral in Eq. (16) against the paper's own Monte Carlo benchmark across the SNR range and $\\beta=2$, with the maximum relative error staying below the claimed 1%.","supporting_citations":[{"cited_title":"Detecting Ships in the New Zealand Exclusive Economic Zone: Requirements for a Dedicated SmallSat SAR Mission,","cited_arxiv_id":null,"evidence_quote":"Supplies the lognormal ship-backscatter model and its fitted parameters, plus the typical backscatter range used to judge detectability."},{"cited_title":"Synthetic Aperture Radar,","cited_arxiv_id":null,"evidence_quote":"Provides the radar equation for average ship-pixel SNR used to scale backscatter into SNR."},{"cited_title":"Functions of One Random Variable,","cited_arxiv_id":null,"evidence_quote":"Justifies the change-of-variables step that turns the lognormal backscatter density into the lognormal SNR density in (7)."},{"cited_title":"Detection Fundamentals,","cited_arxiv_id":null,"evidence_quote":"Gives the Marcum Q-function detection probability for a non-fluctuating signal with random phase in Gaussian noise, and the 1-of-n false-alarm approximation."},{"cited_title":"Refinements of Binary Integration for Swerling Target Fluctuations,","cited_arxiv_id":null,"evidence_quote":"Introduces the regularized incomplete beta function representation of m-of-n binary integration and its inverse, which the paper exploits for numeric stability."}],"review_version":1}