Package 'dicomtk'

Description

Writes a valid DICOM Basic Text SR Storage (⁠1.2.840.10008.5.1.4.1.1.88.11⁠) Part 10 file containing a structured report title, text content, and patient /study identifiers.

Usage

create_dicom_sr(
  out_path,
  title = "Radiology Structured Report",
  text,
  patient_name = "ANONYMOUS",
  patient_id = "ANON001",
  study_uid = ""
)

Arguments

Value

Invisibly returns out_path.

Examples

sr_file <- tempfile(fileext = ".dcm")
create_dicom_sr(
  sr_file,
  title = "AI Radiology Summary",
  text = "FINDINGS: Unremarkable brain parenchyma. IMPRESSION: Normal head CT."
)
read_dicom_sr(sr_file)

Description

Standard DICOM tags with Value Representation DA (Date).

Usage

data(date_tags)

Format

A data.frame with columns tag and name.

Description

Standard DICOM tags with Value Representation DT (DateTime).

Usage

data(datetime_tags)

Format

A data.frame with columns tag and name.

Description

Exports the DICOM header and dataset elements to standard DICOMWeb JSON format (PS3.18).

Usage

dcm2json(file, out_path = NULL)

Arguments

Value

Character string of DICOM JSON (or invisibly out_path if written to a file).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
json_txt <- dcm2json(dcm_file)
substr(json_txt, 1, 200)

Description

Extracts the binary PDF document stored in tag ⁠(0042,0011) EncapsulatedDocument⁠.

Usage

dcm2pdf(file, out_path = tempfile(fileext = ".pdf"), opts = "", ...)

Arguments

Value

Character string of the output .pdf file path.

Description

Exports the DICOM header and dataset elements to XML format.

Usage

dcm2xml(file, out_path = NULL, opts = "", ...)

Arguments

Value

Character string of the output .xml file path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
xml_file <- dcm2xml(dcm_file)
head(readLines(xml_file), 10)

Description

Removes or replaces Protected Health Information (PHI) tags according to HIPAA Safe Harbor and the DICOM PS3.15 Basic Application Level Confidentiality Profile. Supports date shifting, private tag removal, UID regeneration, and custom tag overrides.

Usage

dicom_anonymize(
  path,
  out_path = NULL,
  patient_name = "ANONYMOUS",
  patient_id = "ANON001",
  shift_days = NULL,
  erase_private = TRUE,
  regenerate_uids = TRUE,
  extra_remove_tags = character(),
  replace_tags = list()
)

Arguments

Value

Invisibly returns the path to the anonymized DICOM file.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
anon_file <- dicom_anonymize(dcm_file, patient_name = "STUDY^SUBJECT01")
dicom_tag_value(anon_file, c("PatientName", "PatientIdentityRemoved"))

Description

Answers arbitrary questions about a DICOM file's metadata or image content using the configured LLM (or local metadata inspection).

Usage

dicom_ask(path, question, llm = dicom_llm("ollama"), include_image = FALSE)

Arguments

Value

Character string answer.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
dicom_ask(dcm_file, "What is the modality and image dimensions?", llm = dicom_llm("offline"))

Description

Compresses native uncompressed pixel data into DICOM Lossless RLE transfer syntax (⁠1.2.840.10008.1.2.5⁠, per DICOM PS3.5 Annex G).

Usage

dicom_compress(path, out_path = tempfile(fileext = ".dcm"), method = c("rle"))

Arguments

Value

Invisibly returns out_path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
rle_file <- dicom_compress(dcm_file)
file.exists(rle_file)

Description

Creates a "dicom_connection" object holding PACS / DICOM node connection parameters for use with dicom_echo(), dicom_find(), dicom_get(), and dicom_send(). Supports both DIMSE TCP connections and DICOMweb (QIDO-RS / WADO-RS / STOW-RS) HTTP endpoints.

Usage

dicom_connect(
  host = "127.0.0.1",
  port = 104L,
  called_ae = "ANY-SCP",
  calling_ae = "DICOMTK",
  timeout_ms = 5000L,
  dicomweb_url = NULL,
  local_dir = NULL
)

Arguments

Value

An S3 object of class "dicom_connection".

Examples

conn <- dicom_connect("127.0.0.1", port = 4242, called_ae = "ORTHANC")
conn

Description

Decompresses RLE-compressed (or encapsulated) DICOM files into standard Explicit VR Little Endian (⁠1.2.840.10008.1.2.1⁠) uncompressed DICOM format.

Usage

dicom_decompress(path, out_path = tempfile(fileext = ".dcm"))

Arguments

Value

Invisibly returns out_path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
rle_file <- dicom_compress(dcm_file)
dec_file <- dicom_decompress(rle_file)
is_dicom(dec_file)

Description

Renders the DICOM slice with clinical windowing to a PNG, encodes it as base64, and sends it along with anonymized technical metadata and pixel intensity statistics to a vision-capable LLM (e.g., Ollama llava, OpenAI ⁠gpt-4o⁠, or Anthropic claude-sonnet-4-20250514).

Usage

dicom_describe_image(
  path,
  llm = dicom_llm("ollama"),
  frame = 1L,
  window = "auto",
  prompt = NULL
)

Arguments

Value

A character string with the image description and quantitative stats.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
llm <- dicom_llm("offline")
dicom_describe_image(dcm_file, llm = llm)

Description

Performs a DICOM Upper Layer A-ASSOCIATE-RQ handshake followed by a Verification SOP Class (⁠1.2.840.10008.1.1⁠) C-ECHO-RQ over TCP, or verifies a local directory-backed PACS connection.

Usage

dicom_echo(
  x,
  port = 104L,
  called_ae = "ANY-SCP",
  calling_ae = "DICOMTK",
  timeout_ms = 5000L
)

Arguments

Value

A tibble::tibble with columns host, port, called_ae, success, status, message, and elapsed_ms.

Examples

conn <- dicom_connect(local_dir = system.file("extdata", package = "dicomtk"))
dicom_echo(conn)

Description

Queries a DICOM connection or local PACS directory for matching studies, series, or instances, returning a tidy tibble::tibble that can be piped into dplyr::filter() and dicom_get().

Usage

dicom_find(conn, level = c("STUDY", "SERIES", "IMAGE"), ...)

Arguments

Value

A tibble::tibble of matching DICOM records.

Examples

ext_dir <- system.file("extdata", package = "dicomtk")
conn <- dicom_connect(local_dir = ext_dir)
dicom_find(conn)

Description

Copies or downloads matching DICOM files from a "dicom_connection" or query result tibble into output_dir.

Usage

dicom_get(x, output_dir = tempfile("dicom_get_"), ...)

Arguments

Value

Character vector of retrieved file paths in output_dir.

Examples

ext_dir <- system.file("extdata", package = "dicomtk")
conn <- dicom_connect(local_dir = ext_dir)
files <- conn |> dicom_find() |> dicom_get()
length(files)

Description

Returns a tidy tibble describing the bundled DCMTK version, data dictionary edition, number of standard tags, and active C++ backend.

Usage

dicom_info()

Value

A tibble::tibble with DCMTK configuration details.

Examples

dicom_info()

Description

Creates a "dicom_llm" configuration object connecting dicomtk to:

• Local Ollama (provider = "ollama", default ⁠http://localhost:11434⁠)

• OpenAI / OpenAI-compatible APIs (provider = "openai" or "custom", supporting OpenAI, OpenRouter, LM Studio, vLLM, Azure, or local proxy subscriptions)

• Anthropic Claude (provider = "anthropic", supporting native ⁠/v1/messages⁠ as well as OpenAI-compatible endpoints)

• Offline rule-based engine (provider = "offline", for deterministic local metadata & pixel-statistics summaries without external API calls)

Usage

dicom_llm(
  provider = c("ollama", "openai", "anthropic", "custom", "offline"),
  model = NULL,
  base_url = NULL,
  api_key = NULL,
  temperature = 0.2,
  max_tokens = 1024L,
  anonymize = TRUE,
  offline_fallback = TRUE
)

Arguments

Value

An S3 object of class "dicom_llm".

Examples

# Local Ollama configuration
local_llm <- dicom_llm("ollama", model = "llava:13b")

# OpenAI-compatible configuration
openai_llm <- dicom_llm("openai", model = "gpt-4o")

# Anthropic Claude configuration
claude_llm <- dicom_llm("anthropic", model = "claude-sonnet-4-20250514")

Description

Updates, inserts, or deletes specified tags in a DICOM file using the native C++ engine and writes the modified dataset to out_path.

Usage

dicom_modify(
  path,
  tags = list(),
  out_path = path,
  delete_tags = character(),
  erase_private = FALSE,
  regenerate_uids = FALSE
)

Arguments

Value

Invisibly returns out_path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
tmp <- tempfile(fileext = ".dcm")
dicom_modify(dcm_file, tags = list(PatientName = "ANON^PATIENT"), out_path = tmp)
dicom_tag_value(tmp, "PatientName")

Description

Scans a directory of DICOM files and builds a tidy index table organized by Patient, Study, Series, and Instance UIDs (dcmqrdb model).

Usage

dicom_qrdb_index(directory, recursive = TRUE)

Arguments

Value

A tibble::tibble with columns file, PatientID, PatientName, StudyInstanceUID, StudyDate, Modality, SeriesInstanceUID, SOPInstanceUID, Rows, and Columns.

Examples

ext_dir <- system.file("extdata", package = "dicomtk")
dicom_qrdb_index(ext_dir)

Description

Sets up and returns a local directory-backed "dicom_connection" ready to receive, index, and query DICOM files.

Usage

dicom_receive(
  output_dir = tempfile("dicom_scp_"),
  port = 11112L,
  aetitle = "DICOMTK-SCP"
)

Arguments

Value

A "dicom_connection" object configured with local_dir = output_dir.

Examples

scp <- dicom_receive(tempfile("scp_"))
scp

Description

Combines DICOM technical metadata and vision analysis to produce a structured radiology report (EXAMINATION, TECHNIQUE, FINDINGS, IMPRESSION), and optionally saves it as a DICOM Structured Report (.dcm) via create_dicom_sr().

Usage

dicom_report(
  path,
  llm = dicom_llm("ollama"),
  frame = 1L,
  window = "auto",
  save_sr_path = NULL
)

Arguments

Value

A tibble::tibble with columns section and content, and the full formatted report text in attribute "report_text".

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
rpt <- dicom_report(dcm_file, llm = dicom_llm("offline"))
rpt

Description

Transmits one or more DICOM files to a configured "dicom_connection" (or local PACS repository directory).

Usage

dicom_send(files, conn)

Arguments

Value

A tibble::tibble summarizing the transmission status for each file.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
pacs_dir <- tempfile("pacs_store_")
dir.create(pacs_dir)
conn <- dicom_connect(local_dir = pacs_dir)
dicom_send(dcm_file, conn)

Description

Generates a concise clinical and technical natural-language summary of a DICOM file's header metadata (modality, anatomy, acquisition parameters, image geometry). Automatically strips PHI tags when communicating with non-local LLM endpoints.

Usage

dicom_summarize(path, llm = dicom_llm("ollama"))

Arguments

Value

A character string containing the study summary.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
llm <- dicom_llm("offline")
dicom_summarize(dcm_file, llm = llm)

Description

Checks the official upstream DCMTK GitHub mirror (⁠https://github.com/DCMTK/dcmtk⁠) for new commits, dictionary releases, or C++ module header updates, and optionally synchronizes a local dicomtk source checkout by running tools/bundle_dcmtk.R.

Usage

dicom_sync_upstream(pkg_dir = ".", apply = FALSE, repo = "DCMTK/dcmtk")

Arguments

Value

A tibble::tibble summarizing the bundled DCMTK version/commit, the latest upstream commit SHA/date/message from GitHub, and whether the bundle is up to date.

Examples

dicom_sync_upstream(apply = FALSE)

Description

Searches the 5,354-entry DCMTK 3.7.1 DICOM data dictionary by substring or tag ⁠(GGGG,EEEE)⁠. When an llm object (created by dicom_llm()) is supplied, natural-language questions (e.g., "What tag stores tube current?") are also answered using the configured LLM.

Usage

dicom_tag_lookup(query, llm = NULL)

Arguments

Value

A tibble::tibble of matching DICOM dictionary entries (with an optional "llm_explanation" attribute if llm is provided).

Examples

dicom_tag_lookup("WindowCenter")
dicom_tag_lookup("(0028,1050)")
dicom_tag_lookup("Patient")

Description

Retrieves one or more tag values by tag code (e.g., "(0010,0010)") or standard keyword (e.g., "PatientName", "Modality", "Rows").

Usage

dicom_tag_value(path, tags)

Arguments

Value

A named character vector of tag values (NA_character_ if absent).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
dicom_tag_value(dcm_file, c("Modality", "Rows", "Columns", "(0008,0070)"))

Description

If path is provided, reads and returns the DICOM header tags present in path as a tidy tibble. If path is NULL, returns the complete built-in OFFIS DCMTK 3.7.1 DICOM data dictionary (5,354 tags).

Usage

dicom_tags(path = NULL, ...)

Arguments

Value

A tibble::tibble of DICOM tags.

Examples

# View standard DICOM dictionary
dicom_tags()

# Read tags from a DICOM file
dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
dicom_tags(dcm_file)

Description

Attaches TLS certificate, private key, and CA bundle paths to a "dicom_connection" object for secure DICOM-over-TLS communication (RFC 3240 / DICOM PS3.15).

Usage

dicom_tls_config(
  conn,
  cert_file = NULL,
  key_file = NULL,
  ca_file = NULL,
  profile = c("BCP195", "BCP195-ND", "AES")
)

Arguments

Value

Updated "dicom_connection" object with TLS settings enabled.

Examples

conn <- dicom_connect("pacs.hospital.org", port = 2762) |>
  dicom_tls_config(profile = "BCP195")
conn

Description

Converts a single DICOM file, vector of DICOM files, directory of DICOM slices, or "dicom_series" object into a standard NIfTI-1 (.nii or .nii.gz) volumetric file. Uses RNifti when installed or a built-in native NIfTI-1 binary writer with zero external dependencies.

Usage

dicom_to_nifti(path, out_path = NULL)

Arguments

Value

Invisibly returns out_path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
nii_file <- dicom_to_nifti(dcm_file)
file.exists(nii_file)

Description

Renders a DICOM slice with window/level transformation to a PNG image file.

Usage

dicom_to_png(
  path,
  out_path = tempfile(fileext = ".png"),
  frame = 1L,
  window = "auto",
  width = NULL,
  height = NULL
)

Arguments

Value

Invisibly returns out_path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
png_file <- dicom_to_png(dcm_file)
file.exists(png_file)

Description

Checks whether a DICOM dataset contains digital signatures in DigitalSignaturesSequence (FFFA,FFFA) or ⁠MACParametersSequence (4FFE,0001)⁠ and returns a tidy summary tibble.

Usage

dicom_verify_signature(path)

Arguments

Value

A tibble::tibble of signature elements (0 rows if unsigned).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
dicom_verify_signature(dcm_file)

Description

Launches a Shiny web application for interactive exploration of DICOM files and multi-slice series, featuring slice scrolling, CT window/level presets, custom Window Center & Width sliders, colormap toggling, and a searchable DICOM header tag inspector.

Usage

dicom_viewer(path = NULL, launch.browser = interactive(), ...)

Arguments

Value

Invisibly returns the Shiny application object (or runs the app if interactive).

Examples

if (interactive()) {
  dicom_viewer()
}

Description

Applies linear window/level (Window Center and Window Width) transformation to a DICOM pixel matrix or array, normalizing output values to ⁠[0, 1]⁠. Supports standard clinical CT window presets ("brain", "lung", "bone", "abdomen", "mediastinum", "liver", "stroke", "subdural", "angio"), embedded DICOM WindowCenter/WindowWidth tags ("auto"), or custom numeric center and width.

Usage

dicom_window(x, center = NULL, width = NULL, preset = "auto", invert = NULL)

Arguments

Value

A numeric matrix/array of the same dimensions as x with values scaled to ⁠[0, 1]⁠.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
px <- read_dicom_pixels(dcm_file)
win_auto <- dicom_window(px, preset = "auto")
range(win_auto)

Description

Reads or constructs a tidy Modality Worklist (MWL) schedule table containing ScheduledProcedureStepSequence (0040,0100) attributes (Modality, ScheduledStationAETitle, ScheduledProcedureStepStartDate, PatientName, PatientID, AccessionNumber).

Usage

dicom_worklist(path = NULL)

Arguments

Value

A tibble::tibble representing Modality Worklist scheduled items.

Examples

dicom_worklist()

Description

Alias for plot_dicom() designed for composable ggplot2 pipelines.

Usage

ggdicom(x, frame = 1L, window = "auto", ...)

Arguments

Value

A ggplot2::ggplot object.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
ggdicom(dcm_file, palette = "viridis")

Description

Returns TRUE if the bundled DCMTK C++ engine and data dictionary are initialized and ready.

Usage

have_dcmtk()

Value

A logical scalar (TRUE).

Examples

have_dcmtk()

Description

Inspects the file header for the standard 128-byte preamble + "DICM" magic bytes or valid initial DICOM group/element tags.

Usage

is_dicom(path)

Arguments

Value

A logical vector of the same length as path.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
is_dicom(dcm_file)

Description

Backward-compatible wrappers for functions in the original dcmtk R package (dcmdump, dcmdump_full, dcmhd, wide_hdr, burned_annot, dcmodify, dcm_anon, dcm_subtract_date, dcmj2pnm, dcmdjpeg, dcmdrle, dcmcrle, check_dcmtk_cmd), now powered by the native dicomtk C++17 engine without requiring external CLI binaries.

Usage

dcmdump(file, ...)

dcmdump_full(file, ...)

dcmhd(file, ...)

wide_hdr(hdr, ...)

burned_annot(file, ...)

dcmodify(file, tags = list(), outfile = file, ...)

dcm_anon(file, outfile = NULL, ...)

dcm_subtract_date(file, outfile = NULL, shift_days = -365L, ...)

dcmj2pnm(file, outfile = tempfile(fileext = ".png"), opts = "--write-png", ...)

dcmdjpeg(file, outfile = tempfile(fileext = ".dcm"), ...)

dcmdrle(file, outfile = tempfile(fileext = ".dcm"), ...)

dcmcrle(file, outfile = tempfile(fileext = ".dcm"), ...)

check_dcmtk_cmd(cmd = "dcmdump")

Arguments

Description

Queries the local Ollama daemon (GET /api/tags) and returns a tidy tibble::tibble of installed models (name, size_gb, modified_at). Returns an empty tibble if Ollama is not running.

Usage

ollama_list_models(
  base_url = Sys.getenv("OLLAMA_HOST", unset = "http://localhost:11434")
)

Arguments

Value

A tibble::tibble with columns name, size_gb, and modified_at.

Examples

ollama_list_models()

Description

Renders a 2D slice of a DICOM image as a ggplot2::ggplot object with clinical window/level scaling, proper pixel aspect ratio, and optional radiological metadata overlay.

Usage

plot_dicom(
  x,
  frame = 1L,
  window = "auto",
  center = NULL,
  width = NULL,
  palette = c("gray", "viridis", "inferno", "magma", "plasma"),
  annotate = TRUE,
  title = NULL
)

Arguments

Value

A ggplot2::ggplot object.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
p <- plot_dicom(dcm_file)
p

Description

Reads one or more DICOM files (or directories of DICOM files) using the native C++17 / DCMTK parser and returns a tidy tibble::tibble.

Usage

read_dicom(path, wide = FALSE, include_sequences = FALSE, recursive = TRUE)

Arguments

Value

A tibble::tibble of class c("dicom_tbl", "tbl_df", "tbl", "data.frame").

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
hdr <- read_dicom(dcm_file)
hdr

# Wide format (1 row per file)
read_dicom(dcm_file, wide = TRUE)

Description

Extracts functional groups (SharedFunctionalGroupsSequence (5200,9229), PerFrameFunctionalGroupsSequence (5200,9230)) and tractography sequences (TrackSetSequence (0066,0101)) into a tidy tibble::tibble.

Usage

read_dicom_ect(path)

Arguments

Value

A tibble::tibble of functional group and tractography elements.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
read_dicom_ect(dcm_file)

Description

Reads the header elements of one or more DICOM files into a tidy tibble. Compatible with both dicomtk tidy pipelines and legacy dcmtk::read_dicom_header() arguments.

Usage

read_dicom_header(
  file,
  replace_names = FALSE,
  add_opts = NULL,
  verbose = FALSE,
  wide = FALSE,
  ...
)

Arguments

Value

A tibble::tibble containing the parsed DICOM header elements.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
read_dicom_header(dcm_file)

Description

Reads and decodes 8-bit, 16-bit, or 32-bit monochrome or RGB pixel data (including RLE-compressed frames) directly into an R matrix or array, applying Modality LUT (RescaleSlope and RescaleIntercept) when rescale = TRUE.

Usage

read_dicom_pixels(path, frame = 0L, rescale = TRUE)

Arguments

Value

A numeric matrix ⁠[rows, cols]⁠ (for single-frame monochrome) or 3D/4D array ⁠[rows, cols, frames/channels]⁠ with class "dicom_pixels" and DICOM imaging metadata attached as attributes.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
px <- read_dicom_pixels(dcm_file)
dim(px)

Description

Reads a DICOM Parametric Map (PMAP) dataset and applies RealWorldValueSlope and RealWorldValueIntercept transformations.

Usage

read_dicom_pmap(path)

Arguments

Value

A list of class "dicom_pmap" containing values (real-world mapped numeric array) and metadata (tibble of header tags).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
pmap <- read_dicom_pmap(dcm_file)
dim(pmap$values)

Description

Extracts display shutter, VOI LUT, graphic annotation (0x0070), and spatial transformation parameters from a DICOM Presentation State (PR) object.

Usage

read_dicom_pstate(path)

Arguments

Value

A list of class "dicom_pstate" containing annotations (group 0x0070 elements), voi_lut (window/LUT attributes), and header.

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
ps <- read_dicom_pstate(dcm_file)
ps

Description

Reads DICOM Radiation Therapy objects (RTPLAN, RTDOSE, RTSTRUCT, RTIMAGE), extracting group 0x3002, 0x3004, 0x3006, 0x3008, 0x300A, and 0x300C elements and nested sequences into a tidy tibble::tibble, along with dose grids if pixel data is present (RTDOSE).

Usage

read_dicom_rt(path)

Arguments

Value

A list of class "dicom_rt" containing modality, summary (tidy tibble of top-level RT attributes), rt_elements (full hierarchical RT sequence tibble), and dose_grid (3D dose array if Modality == "RTDOSE").

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
rt_info <- read_dicom_rt(dcm_file)
rt_info$modality

Description

Reads a DICOM Segmentation (SEG) dataset, extracting segment descriptions (SegmentSequence (0062,0002)) and the binary/fractional 3D segmentation mask array.

Usage

read_dicom_seg(path)

Arguments

Value

A list of class "dicom_seg" containing segments (tibble of segment metadata) and mask (pixel/voxel mask array).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
seg <- read_dicom_seg(dcm_file)
dim(seg$mask)

Description

Scans a directory (or vector of DICOM file paths), reads slice headers and pixel data, orders slices spatially by InstanceNumber or ImagePositionPatient, and assembles a 3D numeric array ⁠[rows, cols, slices]⁠ alongside a wide metadata tibble.

Usage

read_dicom_series(path, recursive = TRUE, rescale = TRUE)

Arguments

Value

An object of class "dicom_series" (a named list with elements volume (⁠[rows, cols, slices]⁠ numeric array), metadata (wide tibble of slice headers), and dims (integer vector)).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
series <- read_dicom_series(c(dcm_file, dcm_file))
dim(series$volume)

Description

Parses a DICOM Structured Report (dcmsr) or hierarchical DICOM dataset (including RT Plan/Dose/Structure Set, Segmentation, Parametric Map, and Tractography objects) into a tidy tibble::tibble preserving sequence hierarchy (depth and sequence_path).

Usage

read_dicom_sr(path)

Arguments

Value

A tibble::tibble with columns depth, sequence_path, tag, keyword, vr, and value.

Examples

sr_file <- tempfile(fileext = ".dcm")
create_dicom_sr(
  sr_file,
  title = "Chest CT Report",
  text = "No acute cardiopulmonary abnormality."
)
read_dicom_sr(sr_file)

Description

Reads a DICOM Tractography Results (⁠1.2.840.10008.5.1.4.1.1.66.6⁠) dataset, extracting TrackSetSequence (0066,0101) and streamline metadata into a tidy tibble::tibble.

Usage

read_dicom_tract(path)

Arguments

Value

A list of class "dicom_tract" containing track_sets (tibble of group 0x0066 elements) and header (top-level metadata tibble).

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
tr <- read_dicom_tract(dcm_file)
tr

Description

Standard DICOM tags containing Protected Health Information (PHI) targeted for removal or replacement during anonymization.

Usage

data(removal_tags)

Format

A data.frame with columns tag and name.

Description

Visualizes a 3D DICOM series ("dicom_series" or 3D numeric array) using:

• "mpr": 3-plane Orthogonal Multi-Planar Reconstruction (Axial, Coronal, Sagittal slices) via ggplot2.

• "mip": Maximum Intensity Projection across Axial, Coronal, and Sagittal orientations via ggplot2.

• "rgl": Interactive 3D voxel rendering using the rgl package (falls back to "mip" if rgl is not installed).

Usage

render_volume(
  x,
  method = c("mpr", "mip", "rgl"),
  window = "auto",
  threshold = 0.5,
  palette = c("gray", "viridis", "inferno", "magma", "plasma")
)

Arguments

Value

A ggplot2::ggplot object (for "mpr" and "mip") or invisibly the rgl scene ID (for "rgl").

Examples

dcm_file <- system.file("extdata", "example.dcm", package = "dicomtk")
series <- read_dicom_series(c(dcm_file, dcm_file, dcm_file))
render_volume(series, method = "mpr")
render_volume(series, method = "mip")

Description

Creates a standard DICOM Part 10 Explicit VR Little Endian file from an R 2D matrix or 3D array of pixel values and an optional named list or character vector of DICOM metadata tags.

Usage

write_dicom(x, path, metadata = list())

Arguments

Value

Invisibly returns the normalized output path.

Examples

mat <- matrix(seq(0, 255), nrow = 16, ncol = 16)
out <- tempfile(fileext = ".dcm")
write_dicom(mat, out, metadata = list(PatientName = "TEST^PATIENT", Modality = "CT"))
read_dicom(out, wide = TRUE)