ffmpeg-skill 3.1.4 → 3.1.6

package/scripts/grid.py

@@ -26,15 +26,34 @@ import os
 import sys
 
 from _common import add_common, apply_common, aac_args, cfr_args, default_output, die, emit, X264_PRESETS, fmt_secs, \
-    draws_chroma, drawn_chain, escape_drawtext, ffmpeg_base, info, label_chain, matrix_chain, probe, run, validate_color, video_args
+    as_limited, draws_chroma, drawn_chain, escape_drawtext, ffmpeg_base, info, label_chain, matrix_chain, range_chain, probe, run, validate_color, video_args
 from _common.fonts import font_choice, record_font_choice
 from _common import drawtext_fontfile  # noqa: E402
 from _common import ffmpeg_filters, filter_missing_hint  # noqa: E402
 from _common.emit import note, record_fallback  # noqa: E402
+from _common.fonts import script_font_for_text  # noqa: E402
+from _common.probe import square_pixels  # noqa: E402
+from _common.wrap import text_width_em  # noqa: E402
 
 LABEL_MARGIN = 10
 
 
+def _fit_label(stem: str, cell_width: int, font_size: int) -> str:
+    """The label as drawn: `stem` when it fits the cell (less the margin and the box's border on
+    each side), else its head with an ellipsis. It is right-aligned, so a label wider than the
+    cell used to lose its start off the left edge and show only its tail (#411). The width is the
+    averaged advance table's estimate."""
+    room = (cell_width - 2 * LABEL_MARGIN - 8) / float(max(1, font_size))
+    if text_width_em(stem) <= room:
+        return stem
+    if room <= 1:
+        return stem[:1]  # no room for an ellipsis: keep the start, as a shortened label does (#416 review)
+    head = stem
+    while head and text_width_em(head + "\u2026") > room:
+        head = head[:-1]
+    return (head.rstrip() + "\u2026") if head else stem[:1]
+
+
 def main() -> int:
     ap = argparse.ArgumentParser(description=__doc__, formatter_class=argparse.RawDescriptionHelpFormatter)
     ap.add_argument("inputs", nargs="+", help="video clips, filled into the grid left-to-right, top-to-bottom")
@@ -69,6 +88,9 @@ def main() -> int:
     n = args.cols * args.rows
     if args.cols <= 0 or args.rows <= 0:
         die(f"--cols/--rows must be > 0, got cols={args.cols} rows={args.rows}")
+    if n == 1:
+        # xstack takes two inputs or more: a 1x1 grid used to reach ffmpeg and fail there (exit 222)
+        die("--cols 1 --rows 1 is a single cell; a grid needs at least two (fit.py letterboxes one clip)")
     if len(args.inputs) != n:
         die(f"--cols {args.cols} --rows {args.rows} needs exactly {n} inputs, got {len(args.inputs)}")
     if args.cell_width <= 0 or args.cell_height <= 0:
@@ -95,10 +117,26 @@ def main() -> int:
     durations = [m.get("duration") or 0.0 for m in metas]
     target_duration = max(durations) if args.pad else min(durations)
 
-    font_used = font_choice(args.font)  # the label font, reported as `font_used` when labels are drawn
-    font_file = font_used["file"]
+    # each label is drawn with a font that covers its own script, as every other text-drawing tool
+    # picks one (through 3.1.3 a Japanese or Thai filename came out as tofu, #411); `font_used` is
+    # the first label's choice, and a grid whose labels needed different fonts says so in a note
+    stems = [os.path.splitext(os.path.basename(p))[0] for p in args.inputs]
+    font_explicit = args.font != ap.get_default("font")
+    label_fonts = []
+    for stem in (stems if args.label == "auto" else []):
+        _sc, script_file, script_family = script_font_for_text(stem, font=args.font, font_explicit=font_explicit)
+        label_fonts.append(font_choice(args.font, script_file=script_file, script_family=script_family))
+    font_used = label_fonts[0] if label_fonts else font_choice(args.font)
     if args.label == "auto":
         record_font_choice(font_used)
+        if len({f.get("file") for f in label_fonts}) > 1:
+            note("cell labels drawn with more than one font, by script: "
+                 + "; ".join(sorted({str(f.get("family") or f.get("file")) for f in label_fonts})))
+    labels = [_fit_label(stem, args.cell_width, args.font_size) for stem in stems]
+    cut = [stem for stem, lab in zip(stems, labels) if lab != stem]
+    if args.label == "auto" and cut:
+        note(f"{len(cut)} cell label(s) wider than the {args.cell_width} px cell at --font-size "
+             f"{args.font_size} were shortened from the end with an ellipsis")
     output = args.output or default_output(args.inputs[0], "grid")
 
     cmd = ffmpeg_base()
@@ -117,6 +155,8 @@ def main() -> int:
     convert = [matrix_chain(m, None) for m in metas]
     if any(convert) or chroma:
         convert = [c or label_chain(m, None) for c, m in zip(convert, metas)]
+    # xstack mixes samples: a full-range cell beside a limited one is taken to limited range (#410)
+    convert = [c + range_chain(m, None) for c, m in zip(convert, metas)]
     parts = []
     for i, p in enumerate(args.inputs):
         drawn = [  # (with this colour, with black)
@@ -129,17 +169,19 @@ def main() -> int:
         elif not args.pad:
             drawn.append((f"trim=duration={target_duration:.3f}",) * 2)
         if args.label == "auto":
-            stem = os.path.splitext(os.path.basename(p))[0]
+            font_file = label_fonts[i]["file"]
             font_opt = drawtext_fontfile(font_file) if font_file else f"font='{escape_drawtext(args.font)}'"
-            drawn.append((f"drawtext=text='{escape_drawtext(stem)}':{font_opt}:fontsize={args.font_size}:"
+            drawn.append((f"drawtext=text='{escape_drawtext(labels[i])}':{font_opt}:fontsize={args.font_size}:"
                           f"fontcolor={args.font_color}:x=w-tw-{LABEL_MARGIN}:y=h-th-{LABEL_MARGIN}:"
                           f"box=1:boxcolor=black@0.5:boxborderw=4", "null"))
         gap_r = args.gap if i % args.cols < args.cols - 1 else 0
         gap_b = args.gap if i // args.cols < args.rows - 1 else 0
         if chroma and (gap_r or gap_b):
             drawn.append((f"pad=iw+{gap_r}:ih+{gap_b}:0:0:color={args.background}", f"pad=iw+{gap_r}:ih+{gap_b}:0:0:color=black"))
-        chain = f"fps={args.fps:g},scale={args.cell_width}:{args.cell_height}:force_original_aspect_ratio=decrease,"
-        chain += drawn_chain(",".join(c for c, _n in drawn), metas[i], colours, ",".join(n for _c, n in drawn),
+        # a non-square-pixel clip is squared at its display aspect before it is fitted (#411)
+        sq = (square_pixels(metas[i]["video"]) + ",") if square_pixels(metas[i]["video"]) else ""
+        chain = f"fps={args.fps:g},{sq}scale={args.cell_width}:{args.cell_height}:force_original_aspect_ratio=decrease,"
+        chain += drawn_chain(",".join(c for c, _n in drawn), as_limited(metas[i]), colours, ",".join(n for _c, n in drawn),
                              line_meta=None, tag=f"c{i}", labelled=True)
         parts.append(f"[{i}:v]{convert[i]}{chain}[v{i}]")
 

package/scripts/insert.py

@@ -27,7 +27,10 @@ import argparse
 import math
 import sys
 
+from _common import still_input_args  # noqa: E402
 from _common import add_common, apply_common, default_output, die, emit, ffmpeg_base, info, parse_time, probe, run, video_args, X264_PRESETS, time_arg, fmt_secs
+from _common.emit import note
+from _common.probe import display_size, square_pixels, square_pixels_note
 
 
 def even(n: float) -> int:
@@ -68,7 +71,11 @@ def main() -> int:
     if (meta.get("duration") or 0) > 0.5 or meta.get("audio"):
         die(f"{args.input} is a video, not a still image; insert.py animates a still (Ken Burns). "
             f"For a clip use broll.py (cutaway) or cut.py/join.py")
-    sw, sh = meta["video"]["width"], meta["video"]["height"]
+    # the displayed picture: ffmpeg turns an EXIF-rotated still upright on decode (probe reads
+    # the Orientation off the decoded frame), and non-square pixels are squared (#411)
+    sw, sh = display_size(meta["video"])
+    if square_pixels(meta["video"]):
+        note(square_pixels_note(meta["video"]))
     ratio = sw / sh
 
     if args.width and args.height:
@@ -114,8 +121,10 @@ def main() -> int:
             f"fps={args.fps:g}",
         ]
 
+    if square_pixels(meta["video"]):
+        vf.insert(0, square_pixels(meta["video"]))  # the scale and crop above size the squared picture
     output = args.output or default_output(args.input, "insert", "mp4")
-    cmd = ffmpeg_base() + ["-loop", "1", "-i", args.input, "-t", f"{target:.3f}", "-vf", ",".join(vf)]
+    cmd = ffmpeg_base() + ["-loop", "1", *still_input_args(args.input), "-i", args.input, "-t", f"{target:.3f}", "-vf", ",".join(vf)]
     # the SDR line, whatever the image is; its RGB (PNG) or BT.601 (JPEG) pixels are converted to
     # the BT.709 that line is tagged with
     cmd += video_args(None, args.crf, args.preset, source=meta)

package/scripts/join.py

@@ -42,8 +42,10 @@ import sys
 from typing import Any, Dict, List, Optional, Tuple
 
 from _common.probe import SILENT_PEAK_DBFS, signatures_match, stream_signature
-from _common import STATE, AUDIO_MEDIA_EXT, concat_list_line, guard_sidecar, place_output, refuse, sibling_temp, draws_chroma, drawn_chain, label_chain, matrix_chain, video_args, aac_args, add_common, apply_common, audio_codec_for, default_output, die, dry_run_input_pending, emit, ffmpeg_base, info, measured_level_dbfs, is_audio_output, probe, require_tool, run, validate_color, X264_PRESETS, fmt_secs
-from _common.emit import warn
+from _common import STATE, AUDIO_MEDIA_EXT, concat_list_line, guard_sidecar, place_output, refuse, sibling_temp, as_limited, draws_chroma, drawn_chain, label_chain, matrix_chain, range_chain, video_args, aac_args, add_common, apply_common, audio_codec_for, default_output, die, dry_run_input_pending, emit, ffmpeg_base, info, measured_level_dbfs, is_audio_output, probe, require_tool, run, validate_color, X264_PRESETS, fmt_secs
+from _common.emit import note, warn
+from _common.decision import av_clip_length, av_trim, stream_starts
+from _common.probe import display_size, square_pixels, square_pixels_note
 
 TRANSITIONS = ["fade", "dissolve", "wipeleft", "wiperight", "wipeup", "wipedown", "slideleft", "slideright",
                "circleopen", "circleclose", "fadeblack", "fadewhite", "smoothleft", "smoothright", "radial", "none"]
@@ -103,16 +105,11 @@ def overlap_report(durs: List[float], expected: Optional[float]) -> Dict[str, An
 
 def clip_length(meta: Dict[str, Any], fps: float) -> float:
     """The one length a clip holds in a join, for its picture and its sound alike: the
-    video stream's, or the audio stream's when the sound runs more than a frame past the picture
-    (narration is never cut; an AAC tail adds no frame). The container duration (audio priming,
-    a longer track) is only the fallback when a stream length is unknown (a pending input)."""
-    v = (meta.get("video") or {}).get("duration")
-    a = (meta.get("audio") or {}).get("duration")
-    if not v:
-        return meta.get("duration") or 0.0
-    if a and fps and a > v + 1.0 / fps:
-        return a
-    return v
+    video stream's, or up to where the sound ends when it runs more than a frame past the picture
+    (narration is never cut; an AAC tail adds no frame), counted from the first frame. The container
+    duration (audio priming, a longer track) is only the fallback when a stream length is unknown
+    (a pending input). Shared with reverse.py as av_clip_length."""
+    return av_clip_length(meta, fps)
 
 
 def find_duplicates(paths: List[str]) -> List[Dict[str, Any]]:
@@ -655,9 +652,8 @@ def main() -> int:
     if args.copy:
         return join_copy(args, metas)
     first = metas[0]["video"]
-    fw, fh = first["width"], first["height"]
-    if first.get("rotation") in (90, -90, 270, -270):
-        fw, fh = fh, fw
+    # the first clip's displayed frame: rotation applied, non-square pixels squared (#411)
+    fw, fh = display_size(first)
     if args.width and args.height:
         w, h = args.width, args.height
     elif args.width:
@@ -721,13 +717,22 @@ def main() -> int:
     chroma_pad = args.fit != "crop" and draws_chroma([args.pad_color], line)
     if any(convert) or chroma_pad:
         convert = [c or label_chain(m, line) for c, m in zip(convert, metas)]
+    # concat/xfade mix samples: a full-range clip (VP9/AV1 full range, 10-bit pc) beside a limited
+    # one is taken to limited range first, or it was shown with the other's range (#410)
+    convert = [c + range_chain(m, line) for c, m in zip(convert, metas)]
+
+    for p, m in zip(args.inputs, metas):
+        if square_pixels(m.get("video")):
+            note(square_pixels_note(m["video"], p))
 
     def geo(i: int) -> str:
+        # every clip is squared at its display aspect before it is fitted (#411)
+        sq = (square_pixels(metas[i].get("video")) + ",") if square_pixels(metas[i].get("video")) else ""
         if args.fit == "crop":
-            return f"scale={w}:{h}:force_original_aspect_ratio=increase,crop={w}:{h}"
+            return f"{sq}scale={w}:{h}:force_original_aspect_ratio=increase,crop={w}:{h}"
         pad = f"pad={w}:{h}:(ow-iw)/2:(oh-ih)/2:color="
-        return (f"scale={w}:{h}:force_original_aspect_ratio=decrease,"
-                + drawn_chain(pad + args.pad_color, metas[i], [args.pad_color], pad + "black", line_meta=line,
+        return (f"{sq}scale={w}:{h}:force_original_aspect_ratio=decrease,"
+                + drawn_chain(pad + args.pad_color, as_limited(metas[i]), [args.pad_color], pad + "black", line_meta=line,
                               tag=f"pad{i}", labelled=True))
     # the audio-only join keeps the widest layout; the video join used to force stereo and
     # silently dropped the centre/LFE of 5.1 material. A pending input's stub has no measured
@@ -741,7 +746,12 @@ def main() -> int:
         info(f"channel counts differ ({', '.join(str(c) for c in chans)}); every clip becomes {layout}")
     for i in range(n):
         parts.append(f"[{i}:v]{convert[i]}{geo(i)},setsar=1,setpts=PTS-STARTPTS,fps={fps:g},format={pixfmt},settb=AVTB[v{i}]")
-        parts.append(f"[{audio_src[i]}]aformat=sample_rates=48000:channel_layouts={layout},asetpts=PTS-STARTPTS[a{i}]")
+        # the sound is shifted by the picture's start, not to its own (#409): #338 started each
+        # clip's picture at 0, and a per-stream PTS-STARTPTS then moved a sound that starts 0.3 s
+        # after its picture 0.3 s early. Sound before the first frame is dropped and a late sound
+        # is padded, as cut.py's re-cut does; a generated silent track starts at 0 already.
+        shift = (av_trim(stream_starts(metas[i])[0], None)[1] + ",") if metas[i].get("audio") else ""
+        parts.append(f"[{audio_src[i]}]{shift}aformat=sample_rates=48000:channel_layouts={layout},asetpts=PTS-STARTPTS[a{i}]")
 
     # pad the shorter stream to the clip's length: the audio with silence, the picture by holding
     # its last frame, then cut both there -- each clip boundary is then the same for both chains

package/scripts/look.py

@@ -17,10 +17,13 @@ from pathlib import Path
 from typing import List, Optional
 
 from _platforms import PLATFORMS, PLATFORM_CHOICES, resolve as resolve_platform
+from _common import still_input_args  # noqa: E402
 from _common import STATE, require_tool, add_common, apply_common, default_font_file, die, emit, escape_drawtext, ffmpeg_base, info, parse_time, probe, run, run_analysis, time_arg
 from _common import drawtext_fontfile  # noqa: E402
 from _common import ffmpeg_filters, filter_missing_hint  # noqa: E402
 from _common.emit import note, record_fallback  # noqa: E402
+from _common.color import bt709_label_filter  # noqa: E402
+from _common.probe import display_size, square_pixels, square_pixels_note  # noqa: E402
 
 FONT = "fontcolor=white:fontsize=h/18:box=1:boxcolor=black@0.55:boxborderw=6:x=8:y=8"
 
@@ -40,6 +43,12 @@ def timecode_filter(font_prefix: str, offset: float = 0.0) -> str:
     return f"drawtext=text='%{{pts\\:hms}}':{font_prefix}{FONT}"
 
 
+def _half_height(half: int, size: "tuple") -> int:
+    """The even height a `half`-wide picture of displayed `size` scales to (16:9 when unknown)."""
+    w, h = size
+    return max(2, int(round(half * (h / w if w and h else 9 / 16) / 2)) * 2)
+
+
 def _png_size(png_path: str) -> "Optional[tuple]":
     # run_analysis(): the same wall-clock ceiling and #234 UTF-8 decoding as every other
     # ffprobe/ffmpeg measurement in this codebase, check=False since a probe failure here is
@@ -188,10 +197,7 @@ def main() -> int:
     if args.safe:
         z = PLATFORMS[args.safe]["safe"]
         frame = PLATFORMS[args.safe]["frame"]
-        src = meta["video"]
-        dw, dh = src.get("width") or 0, src.get("height") or 0
-        if src.get("rotation") in (90, -90, 270, -270):
-            dw, dh = dh, dw
+        dw, dh = display_size(meta["video"])
         if dw and dh and abs(dw / dh - frame["w"] / frame["h"]) > 0.02:
             info(f"--safe {args.safe}: this source is {dw}x{dh}, not {args.safe}'s "
                  f"{frame['w']}x{frame['h']} -- the zones are drawn as fractions of the frame you gave, "
@@ -220,11 +226,25 @@ def main() -> int:
         tc = "," + tm.rstrip(",") + tc
         info("HDR source: frames are tone-mapped to SDR for display")
     outputs: List[str] = []
+    # An untagged SDR source is BT.709 (what players and every encoder line here take it as), but the
+    # conversion to the PNG's RGB reads "unknown" as BT.601: C03028 came out B31F28 (#410). Its
+    # frames are labelled first (tags only), so that conversion reads them as BT.709.
+    lab = bt709_label_filter(meta)
+    lab = lab + "," if lab else ""
+    # every picture is squared at its display aspect before it is scaled (#411): a DV or HDV
+    # frame's stored size is not the picture a player shows
+    square = square_pixels(meta["video"])
+    sq = square + "," if square else ""
+    if square:
+        note(square_pixels_note(meta["video"]))
 
     if args.compare:
         if not args.at:
             die("--compare needs --at TIME")
-        probe(args.compare)
+        cmeta = probe(args.compare)
+        lab2 = bt709_label_filter(cmeta)
+        lab2 = lab2 + "," if lab2 else ""
+        csq = (square_pixels(cmeta.get("video")) + ",") if square_pixels(cmeta.get("video")) else ""
         for idx, t in enumerate(args.at):
             sec = time_arg(t, "--at", meta["video"].get("fps") if meta.get("video") else None)
             if args.output and len(args.at) == 1:
@@ -237,10 +257,17 @@ def main() -> int:
             half = args.width // 2
             stamp = "" if args.no_timecode else f",drawtext=text='{escape_drawtext(fmt_hms(sec))}':{font_prefix}{FONT}"
             tcs = tc.replace("," + timecode_filter(font_prefix), "") + stamp
-            fc = (f"[0:v]scale={half}:-2{tcs}[a];[1:v]scale={half}:-2{tcs}[b];"
-                  f"[a][b]scale2ref=w=iw:h=ih[a2][b2];[a2][b2]hstack=inputs=2[out]")
-            cmd = ffmpeg_base() + ["-ss", f"{sec:.3f}", "-i", args.input, "-ss", f"{sec:.3f}", "-i", args.compare,
-                                   "-filter_complex", fc, "-map", "[out]", "-frames:v", "1", out]
+            # each picture keeps its own aspect inside a half-width box as tall as the taller of the
+            # two (letterboxed, not stretched): through 3.1.3 scale2ref forced the first to the
+            # second's size, so a 4:3 beside a 16:9 came out squashed (#411)
+            box_h = max(_half_height(half, display_size(meta["video"])),
+                        _half_height(half, display_size(cmeta.get("video"))))
+            fit = (f"scale={half}:{box_h}:force_original_aspect_ratio=decrease,"
+                   f"pad={half}:{box_h}:(ow-iw)/2:(oh-ih)/2:color=0x202020,setsar=1")
+            fc = (f"[0:v]{lab}{sq}{fit}{tcs}[a];[1:v]{lab2}{csq}{fit}{tcs}[b];"
+                  f"[a][b]hstack=inputs=2[out]")
+            cmd = ffmpeg_base() + ["-ss", f"{sec:.3f}", *still_input_args(args.input), "-i", args.input, "-ss", f"{sec:.3f}", *still_input_args(args.compare), "-i", args.compare,
+                                   "-filter_complex", fc, "-map", "[out]", "-frames:v", "1", "-update", "1", out]
             run(cmd)
             outputs.append(out)
     elif args.at:
@@ -254,7 +281,7 @@ def main() -> int:
                 # several frames, or -o given as a stem/prefix without an image extension
                 out = os.path.join(outdir, f"{args.output and Path(args.output).stem or stem}_{sec:.3f}s.png")
             stamp = "" if args.no_timecode else f",drawtext=text='{escape_drawtext(fmt_hms(sec))}':{font_prefix}{FONT}"
-            cmd = ffmpeg_base() + ["-ss", f"{sec:.3f}", "-i", args.input, "-vf", f"scale={args.width}:-2{tc.replace(',' + timecode_filter(font_prefix), '')}{stamp}", "-frames:v", "1", out]
+            cmd = ffmpeg_base() + ["-ss", f"{sec:.3f}", *still_input_args(args.input), "-i", args.input, "-vf", f"{lab}{sq}scale={args.width}:-2{tc.replace(',' + timecode_filter(font_prefix), '')}{stamp}", "-frames:v", "1", "-update", "1", out]
             run(cmd)
             outputs.append(out)
     else:
@@ -277,11 +304,14 @@ def main() -> int:
         # sample at the middle of each slice so the first/last tiles are not black lead-in/out frames
         # with only a handful of frames a mid-slice seek skips past them all: start at 0 instead
         seek = 0.0 if frames_avail < 4 else step / 2
-        # the seek resets pts to 0, so the burned-in timecode gets the seek added back
-        tc_sheet = tc.replace(timecode_filter(font_prefix), timecode_filter(font_prefix, seek)) if seek else tc
-        vf = (f"select='isnan(prev_selected_t)+gte(t-prev_selected_t\\,{step * 0.98:.6f})',scale={tile_w}:-2{tc_sheet},"
+        # -copyts -start_at_zero keeps each frame's own timestamp through the seek, counted from
+        # the file's start, so the stamp is the frame's pts as it is. What an input -ss alone
+        # resets pts to differs by container: an MPEG-TS (which starts at 1.4 s) came out
+        # labelled one seek late when the seek was added back (#411)
+        tc_sheet = tc
+        vf = (f"select='isnan(prev_selected_t)+gte(t-prev_selected_t\\,{step * 0.98:.6f})',{lab}{sq}scale={tile_w}:-2{tc_sheet},"
               f"tile={cols}x{rows}:padding=2:margin=2:color=0x202020")
-        cmd = ffmpeg_base() + ["-ss", f"{seek:.6f}", "-i", args.input, "-vf", vf, "-frames:v", "1", out]
+        cmd = ffmpeg_base() + ["-copyts", "-start_at_zero", "-ss", f"{seek:.6f}", *still_input_args(args.input), "-i", args.input, "-vf", vf, "-frames:v", "1", "-update", "1", out]
         run(cmd)
         outputs.append(out)
         info(f"contact sheet: {n} frames every {step:.2f}s")

package/scripts/loop.py

@@ -21,11 +21,15 @@ Examples:
   python3 loop.py wave.mp4 --boomerang --times 4
 """
 import argparse
+import json
 import math
+import os
 import sys
+import tempfile
 
-from _common import add_common, aac_args, apply_common, cfr_args, default_output, die, emit, ffmpeg_base, info, parse_time, probe, run, video_args, X264_PRESETS, time_arg, fmt_secs
+from _common import concat_list_line, add_common, aac_args, apply_common, cfr_args, default_output, die, emit, ffmpeg_base, info, parse_time, probe, refuse_output_is_input, require_tool, run, video_args, X264_PRESETS, time_arg, fmt_secs
 from _common.emit import warn
+from _common.decision import av_trim
 
 
 # The loop filter's size is capped at 32767 frames, so one boomerang cycle (2N-2 frames) must fit.
@@ -95,6 +99,21 @@ def _nominal_rate(v) -> str:
     return f"{counted:g}" if 0 < counted <= MAX_PLAUSIBLE_FPS else "30"
 
 
+def concat_first_frame(listfile: str) -> float:
+    """Where the first frame sits on the clock ffmpeg reads a concat list on: the list's first video
+    start minus its start_time. The concat demuxer counts an MP4's hidden AAC priming packet into its
+    start (-0.021 s) where the MP4 demuxer does not, so ffmpeg moves both streams 21 ms later, and a
+    frame-rate time base then rounded the picture's move to 33 ms (12 ms off). 0.0 when unmeasured."""
+    proc = run([require_tool("ffprobe"), "-v", "error", "-f", "concat", "-safe", "0", "-i", listfile,
+                "-select_streams", "v:0", "-show_entries", "stream=start_time:format=start_time", "-of", "json"],
+               quiet=True, check=False)
+    try:
+        doc = json.loads(proc.stdout or "")
+        return float(doc["streams"][0]["start_time"]) - float(doc["format"]["start_time"])
+    except (ValueError, KeyError, IndexError, TypeError):
+        return 0.0
+
+
 def main() -> int:
     ap = argparse.ArgumentParser(description=__doc__, formatter_class=argparse.RawDescriptionHelpFormatter)
     ap.add_argument("input")
@@ -119,6 +138,8 @@ def main() -> int:
     has_audio = bool(meta.get("audio")) and not args.boomerang
     notes = []
     output = args.output or default_output(args.input, "loop")
+    # the concat list is the only -i, so run() never sees the source next to -o (#416 review)
+    refuse_output_is_input(output, args.input)
 
     if args.times is not None:
         if args.times < (1 if args.boomerang else 2):
@@ -137,10 +158,24 @@ def main() -> int:
         if meta.get("audio"):
             notes.append("audio dropped: a boomerang plays half its sound backwards -- add a bed with audio.py")
     else:
-        # -stream_loop repeats the whole input read (video and audio together) at the demuxer level
-        # -- exact and lossless-in-intent for a re-encode target, unlike a filter-graph loop that
-        # would need separate video/audio filters kept in lockstep by hand.
-        cmd = ffmpeg_base() + ["-stream_loop", str(stream_loop), "-i", args.input]
+        # The concat demuxer reads the clip once per repeat and starts each repeat at the previous
+        # one's start plus the file's duration (#409): both streams keep their offset from the
+        # file's start in every repeat. -stream_loop took the repeat's length from its streams'
+        # spans without the last video frame, so a sound that starts 0.3 s late shortened every
+        # repeat by a frame and slid against the picture (346 frames for 360, 68 ms off).
+        # Both streams start at the first frame, shifted by one constant (av_trim), and the gap a
+        # late sound leaves at each repeat's start is filled with silence, at any length.
+        workdir = tempfile.TemporaryDirectory(prefix="ffskill_loop_")
+        listfile = os.path.join(workdir.name, "loop.txt")
+        with open(listfile, "w", encoding="utf-8") as fh:
+            fh.write((concat_list_line(os.path.abspath(args.input)) + "\n") * (stream_loop + 1))
+        vtrim, atrim = av_trim(concat_first_frame(listfile), None, fill_gaps=True)
+        cmd = ffmpeg_base() + ["-f", "concat", "-safe", "0", "-i", listfile]
+        if has_audio:
+            cmd += ["-filter_complex", f"[0:v:0]{vtrim}[v];[0:a:0]{atrim}[a]",
+                    "-map", "[v]", "-map", "[a]"]
+        else:
+            cmd += ["-map", "0:v:0", "-vf", vtrim]
     if target is not None:
         cmd += ["-t", f"{target:.3f}"]
     cmd += video_args(meta, args.crf, args.preset)
@@ -150,7 +185,11 @@ def main() -> int:
     else:
         cmd += ["-an"]
     cmd.append(output)
-    run(cmd)
+    try:
+        run(cmd)
+    finally:
+        if not args.boomerang:
+            workdir.cleanup()
 
     result = probe(output, role="output")
     v = result["video"]

package/scripts/multicam.py

@@ -31,8 +31,10 @@ import math
 import sys
 from typing import Dict, List, Tuple
 
-from _common import guard_sidecar, label_chain, matrix_chain, video_args, aac_args, add_common, apply_common, default_output, die, emit, ffmpeg_base, info, time_arg, probe, run, x264_args, X264_PRESETS, fmt_secs, decode_pcm_mono, rms_envelope, STATE
-from _common.emit import record_fallback, warn
+from _common import guard_sidecar, label_chain, matrix_chain, range_chain, video_args, aac_args, add_common, apply_common, default_output, die, emit, ffmpeg_base, info, time_arg, probe, run, x264_args, X264_PRESETS, fmt_secs, decode_pcm_mono, rms_envelope, STATE
+from _common.emit import note, record_fallback, warn
+from _common.decision import av_trim
+from _common.probe import display_size, square_pixels, square_pixels_note
 from sync import measure_offset
 
 ENERGY_WINDOW_S = 0.25  # --switch energy: loudness measured over this window on the reference timeline
@@ -263,14 +265,16 @@ def main() -> int:
         info(f"wrote {args.edl}")
 
     v0 = metas[0]["video"]
-    w, h = args.width or v0["width"], args.height or v0["height"]
-    if v0.get("rotation") in (90, -90, 270, -270) and not (args.width or args.height):
-        w, h = h, w
+    # camera 0's displayed frame: rotation applied, non-square pixels squared (#411)
+    w, h = args.width or display_size(v0)[0], args.height or display_size(v0)[1]
     fps = args.fps or v0.get("fps") or 30.0
     fps = round(fps) if abs(fps - round(fps)) < 0.02 else fps
     pixfmt = "yuv420p10le" if v0.get("bt2020_or_hdr") else "yuv420p"
     geo = f"scale={w}:{h}:force_original_aspect_ratio=decrease,pad={w}:{h}:(ow-iw)/2:(oh-ih)/2,setsar=1,fps={fps:g},format={pixfmt}"
 
+    for p, m in zip(args.inputs, metas):
+        if square_pixels(m.get("video")):
+            note(square_pixels_note(m["video"], p))
     cmd = ffmpeg_base()
     for p in args.inputs:
         cmd += ["-i", p]
@@ -283,6 +287,8 @@ def main() -> int:
     convert = [matrix_chain(m, metas[0]) for m in metas]
     if any(convert):
         convert = [c or label_chain(m, metas[0]) for c, m in zip(convert, metas)]
+    # concat mixes samples: a full-range camera beside a limited one is taken to limited range (#410)
+    convert = [c + range_chain(m, metas[0]) for c, m in zip(convert, metas)]
     # what is actually rendered: a range whose camera has not started yet plays camera 0, and the
     # result says so (`cuts` names the camera on screen, diagnostics.fallbacks the one asked for)
     rendered: List[Tuple[float, float, int]] = []
@@ -296,7 +302,8 @@ def main() -> int:
                             requested_camera=c, camera=0, starts_at=round(offsets[c], 3))
             c, src_s, src_e = 0, s, e
         rendered.append((s, e, c))
-        parts.append(f"[{c}:v]trim=start={src_s:.4f}:end={src_e:.4f},setpts=PTS-STARTPTS,{convert[c]}{geo}[v{i}]")
+        sq = (square_pixels(metas[c]["video"]) + ",") if square_pixels(metas[c]["video"]) else ""
+        parts.append(f"[{c}:v]trim=start={src_s:.4f}:end={src_e:.4f},setpts=PTS-STARTPTS,{convert[c]}{sq}{geo}[v{i}]")
         labels.append(f"[v{i}]")
         project_clips.append({"src": args.inputs[c], "in": round(src_s, 3), "out": round(src_e, 3)})
     parts.append("".join(labels) + f"concat=n={len(filled)}:v=1:a=0[vout]")
@@ -306,7 +313,9 @@ def main() -> int:
     # applied before asetrate/aresample rescale that axis -- otherwise the trim lands at the wrong
     # point once the stream's timebase has already been stretched/compressed by the drift ratio
     # (mirrors sync.py, which seeks with -ss, an input-level operation, before its drift_af filters).
-    afx = [f"atrim=start={a_start:.4f}", "asetpts=PTS-STARTPTS"]
+    # shifted by a_start, not to its own first sample: a sound that starts after its file's start
+    # keeps that offset (padded with silence), as the offset was measured on the file's clock (#409)
+    afx = [av_trim(a_start, None)[1]]
     if abs(ratios[a] - 1.0) > 1e-7:
         sr = metas[a]["audio"].get("sample_rate") or 48000
         afx += [f"asetrate={sr * ratios[a]:.6f}", f"aresample={sr}"]