harpy.im.rasterize

Contents

harpy.im.rasterize#

harpy.im.rasterize(sdata, shapes_name, output_labels_name, out_shape=None, chunks=None, client=None, scale_factors=None, overwrite=False)#

Given a shapes element in a SpatialData object, corresponding masks are created, and added as a labels element to the SpatialData object.

The index of the shapes element will be used as the label in the resulting labels element (output_labels_name).

Parameters:

  • sdata (SpatialData) – The SpatialData object.

  • shapes_name (str) – The shapes element to be converted to a labels element.

  • output_labels_name (str) – Name of the resulting labels element that will be added to sdata.

  • out_shape (tuple[int, int] | None (default: None)) – Output shape of the resulting labels element (y,x). Will be automatically calculated if set to None via sdata[shapes_name].geometry.total_bounds. If out_shape is not None, with x_min, y_min, x_max, y_max = sdata[shapes_name].geometry.total_bounds, and out_shape[1]<x_max or out_shape[0]<y_max, then shapes with coordinates outside out_shapes will not be in resulting output_labels_name. For shapes_name with large offset (y_min, x_min), we recommend translating the shapes to the origin, and add the offset via a translation (spatialdata.transformations.Translation).

  • chunks (int | None (default: None)) – If provided, creation of the labels element will be done in a chunked manner, with data divided into chunks for efficient computation.

  • client (Client | None (default: None)) – A Dask client. If specified, a copy of sdata[shapes_name] will be scattered across the workers, reducing the size of the task graph. If not specified, Dask will use the default scheduler as configured on your system.

  • client – A Dask Client instance. If specified, a copy of the GeoDataFrame (sdata[shapes_name]) will be scattered across the workers. This reduces the size of the task graph and can improve performance by minimizing data transfer overhead during computation. If not specified, Dask will use the default scheduler as configured on your system (e.g., single-threaded, multithreaded, or a global client if one is running).

  • scale_factors (Sequence[dict[str, int] | int] | None (default: None)) – Scale factors to apply for multiscale.

  • overwrite (bool (default: False)) – If True, overwrites the output_labels_name if it already exists in sdata.

Return type:

SpatialData

Returns:

: An updated SpatialData object with the added labels element.

Raises:

  • ValueError – If the provided shapes_name contains 3D polygons.

  • ValueError – If the provided shapes_name contains Points.

  • ValueError – If 0 is in the index of the shapes_name. As 0 is used as background in the output_labels_name.

  • TypeError – If chunks is not None and not an instance of int.

Contents