o i5 @s`dZddlmZddlmZmZddlmZmZzddl Z dZ Wn e y+dZ Ynwzddl Z ddl mZdZWn e yEdZYnwdd lmZdd lmZdd lmZdd lmZerldd lmZddlmZddddZdddddZ dHdddddeeBdBddfddZ    dIddddd d!d"e dBd#edBddf d$d%Z!dJd'd(Z"dJd)d*Z#e"Z$ dKdddd&d+e d,e%ddf d-d.Z&dLdddd&d+e ddfd/d0Z'dLdddd&d+e ddfd1d2Z( dMddd3ed4ed5e%ddf d6d7Z)dNd8d9Z*dNd:d;Z+  <  dOddd=e,d"e dBd#edBdd>f d?d@Z-  dHdAe.e dBfdCedDe dBdEe dBde/e f dFdGZ0dS)Pz Pixel data processing functions.)BytesIO)unpack unpack_from) TYPE_CHECKINGcastNTF)ImageCms)get_palette_files) warn_and_log)UID)VR)Iterable)DatasetsRGBXYZLAB)SRGBrr proceduralzrelative colorimetric saturationzabsolute colorimetric)rarr np.ndarraydszDataset | Nonepalettereturnc s|s|std|r_dddddddd d }t|js;zd d d dddddd}||}Wnty:td|dwzddlm}||}|t|d}Wnty^td|dwtd|}t|dddvrptdd|vrxtdtt t |j }|dpd}|d } |d!} t d"| d#g} d$|vrtt|j} tt|j} tt|j}ttdBt|d%d}t| |d&}t d"|d#| fd'd(| | ||fDnd)|vratt|j} tt|j} tt|j}ttdBt|d*d}t|d+r |jjd }n|jd }|dur td,t|jd-d.|}| d&}|d kr/d/nd0}| }| | ||fD]&}|r^t||}|||}tt|||}| t j!|d1q9ntd|d2vrstd3|d4d5d6d7| Dt"fd8d(d dDstd9t j#|j$d1}|| k}||t %| ||<t j&|d|d |d:t j't |j$t| gd1}t(| D] \}}|||d;|f<q|S)` or :dcm:`Supplemental Palette Color LUT Module` is not currently supported. Parameters ---------- arr : numpy.ndarray The pixel data to apply the color palette to. ds : dataset.Dataset, optional Required if `palette` is not supplied. A :class:`~pydicom.dataset.Dataset` containing a suitable :dcm:`Image Pixel` or :dcm:`Palette Color Lookup Table` Module. palette : str or uid.UID, optional Required if `ds` is not supplied. The name of one of the :dcm:`well-known` color palettes defined by the DICOM Standard. One of: ``'HOT_IRON'``, ``'PET'``, ``'HOT_METAL_BLUE'``, ``'PET_20_STEP'``, ``'SPRING'``, ``'SUMMER'``, ``'FALL'``, ``'WINTER'`` or the corresponding well-known (0008,0018) *SOP Instance UID*. Returns ------- numpy.ndarray The RGB or RGBA pixel data as an array of ``np.uint8`` or ``np.uint16`` values, depending on the 3rd value of (0028,1201) *Red Palette Color Lookup Table Descriptor*. References ---------- * :dcm:`Image Pixel Module` * :dcm:`Supplemental Palette Color LUT Module` * :dcm:`Enhanced Palette Color LUT Module` * :dcm:`Palette Colour LUT Module` * :dcm:`Supplemental Palette Color LUTs ` z$Either 'ds' or 'palette' is requiredz hotiron.dcmzpet.dcmzhotmetalblue.dcmz pet20step.dcmz spring.dcmz summer.dcmzfall.dcmz winter.dcm)1.2.840.10008.1.5.11.2.840.10008.1.5.21.2.840.10008.1.5.31.2.840.10008.1.5.41.2.840.10008.1.5.51.2.840.10008.1.5.61.2.840.10008.1.5.71.2.840.10008.1.5.8rrrrr r!r#r")HOT_IRONPETHOT_METAL_BLUE PET_20_STEPSPRINGSUMMERFALLWINTERzUnknown palette ''r)dcmreadr PixelPresentationN)MIXEDCOLORzgUse of this function with the Supplemental Palette Color Lookup Table Module is not currently supported$RedPaletteColorLookupTableDescriptorz3No suitable Palette Color Lookup Table Module foundrruintz.0fRedPaletteColorLookupTableData AlphaPaletteColorLookupTableDatac3s"|] }|rtj|dVqdS)dtypeN)np frombuffer).0 lut_bytesr7O/mnt/sdb/aimis/docanh/lib/python3.10/site-packages/pydicom/pixels/processing.py s  z"apply_color_lut..'SegmentedRedPaletteColorLookupTableData)SegmentedAlphaPaletteColorLookupTableData file_metaeUnable to determine the endianness of the dataset, please set an appropriate Transfer Syntax UID in ' .file_meta'><BHr7)r6zThe bit depth of the LUT data 'z.1fz2' is invalid (only 8 or 16 bits per entry allowed)cSg|]}t|qSr=)lenr;iir=r=r> z#apply_color_lut..c3s|] }|dkVqdS)rNr=rK) lut_lengthsr=r>r?sz LUT data must be the same lengthout.)) ValueErrorr is_validKeyErrorpydicomr-rrgetattrlistintr1r9r8bytesr4 GreenPaletteColorLookupTableDataBluePaletteColorLookupTableDatarJextendr@)SegmentedGreenPaletteColorLookupTableData(SegmentedBluePaletteColorLookupTableDatahasattrrB_tsyntax_encodingoriginal_encodingAttributeErrortype__name___expand_segmented_lutrappendasarrayallzerosshapeint32clipempty enumerate)rrrdatasetsuidsr-fnamelut_desc nr_entries first_map nominal_depthlutsr_lutg_lutb_luta_lut actual_depthis_little_endian endianness byte_depthfmtseglen_segs_fmtlut_ints clipped_iv mapped_pixelsrQrLlutr=)r8rOr>apply_color_lut*s3                     r transformz%PIL.ImageCms.ImageCmsTransform | Noneintent color_spacecCststd|dur|std|dur|rtd|jdvr'td|j|jddkr6td |j|s?t|||d }|jd krdt|D]\}}tj |}t j ||d d t |||<qH|Stj |}t j ||d d t ||d<|S)a{Apply an `ICC Profile `_ to `arr`, either from the dataset `ds` or an existing Pillow color transformation object `transform`. .. versionadded:: 3.0 .. warning:: This function requires `NumPy `_ and `Pillow `_. Parameters ---------- arr : numpy.ndarray 8-bit RGB pixel data to apply an ICC profile to, shaped as either (rows, columns, samples) or (frames, rows, columns, samples). ds : pydicom.dataset.Dataset, optional Required if `transform` is not supplied, a :class:`~pydicom.dataset.Dataset` containing elements from the :dcm:`ICC Profile` module. transform : PIL.ImageCms.ImageCmsTransform, optional An :class:`~PIL.ImageCms.ImageCmsTransform` instance such as is returned by the :func:`create_icc_transform` function. Required if `ds` is not used and recommended when the same ICC profile is to be re-used on multiple ndarrays. intent : int, optional If `transform` is not supplied this is the rendering intent of the transformation that will be created from `ds`, one of: * ``0``: perceptual * ``1``: relative colorimetric * ``2``: saturation * ``3``: absolute colorimetric If no `intent` is specified then the default rendering intent in the ICC profile will be used. color_space : str, optional If `transform` is not supplied this is the output color space to use for the transformation created from `ds`. If not used then defaults to the (0028,2002) *Color Space* value if its present in `ds`, otherwise defaults to ``"sRGB"``. If used then it will override the (0028,2002) *Color Space* element value (if present). Must be one of the :func:`supported Pillow color spaces `, although in practice only ``"sRGB"`` is likely to be used. Returns ------- np.ndarray The input ndarray with the color transformation applied in-place. z8Pillow is required to apply an ICC profile to an ndarrayNz+Either 'ds' or 'transform' must be suppliedz9Only one of 'ds' and 'transform' should be used, not both)rz-The ndarray must have 3 or 4 dimensions, not rzUInvalid ndarray shape, must be (rows, columns, 3) or (frames, rows, columns, 3), not )rrrrT)inPlace.)HAVE_PIL ImportErrorrRndimrjcreate_icc_transformrnPILImage fromarrayrapplyTransformr9array)rrrrridxframeimr=r=r>apply_icc_profiles<9      rr cCs|drttd|jd}ttt|jdpd}ttt|jd}ttt|jd}d|}|djtjkrqt|d rH|j j d}n|j d}|d ur\t d t |jd d |}||d} t| tt|j} nttt|j} tj| |d} tj|j|jd} ||k} || || | <tj| d|d| d| | Sd|vrd|vr|tjtt|j}|tt|j7}|S)aApply a modality lookup table or rescale operation to `arr`. Parameters ---------- arr : numpy.ndarray The :class:`~numpy.ndarray` to apply the modality LUT or rescale operation to. ds : dataset.Dataset A dataset containing a :dcm:`Modality LUT Module `. Returns ------- numpy.ndarray An array with applied modality LUT or rescale operation. If (0028,3000) *Modality LUT Sequence* is present then returns an array of ``np.uint8`` or ``np.uint16``, depending on the 3rd value of (0028,3002) *LUT Descriptor*. If (0028,1052) *Rescale Intercept* and (0028,1053) *Rescale Slope* are present then returns an array of ``np.float64``. If neither are present then `arr` will be returned unchanged. Notes ----- When *Rescale Slope* and *Rescale Intercept* are used, the output range is from (min. pixel value * Rescale Slope + Rescale Intercept) to (max. pixel value * Rescale Slope + Rescale Intercept), where min. and max. pixel value are determined from (0028,0101) *Bits Stored* and (0028,0103) *Pixel Representation*. References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.11.1 ` * DICOM Standard, Part 4, :dcm:`Annex N.2.1.1 ` ModalityLUTSequencer rr2rrr3LUTDatarBNrCrDrErGr7rP RescaleSlopeRescaleIntercept)getrrWrrX LUTDescriptorr OWr_rBr`rarbrcrdrrYrr9rgrirjr8rlastypefloat64floatrr)rritemrsrtrur8r|r} unpack_fmtunc_datalut_datarrr=r=r>apply_modality_lutOs< &   rc CsTd|vr|jd}|j\}}}|dkrdn|}|dkrdnd}||d}|d}|jtjkr8tj|jdd} ntj|jd |d |d} ||} re| j j sU| } tj | | | d tj | | | d |d }||8}|||d }|d}td| |Sd|vr|j} | dvrtd|jd| dkr||}|S)a-Apply a Presentation LUT to `arr` and return the P-values. Parameters ---------- arr : numpy.ndarray The :class:`~numpy.ndarray` to apply the presentation LUT operation to. ds : dataset.Dataset A dataset containing :dcm:`Presentation LUT Module ` elements. Returns ------- numpy.ndarray If a Presentation LUT Module is present in `ds` then returns an array of P-values, otherwise returns `arr` unchanged. Notes ----- If the dataset the *Pixel Data* originated from contains a Modality LUT and/or VOI LUT then they must be applied before the Presentation LUT. See Also -------- :func:`~pydicom.pixels.processing.apply_modality_lut` :func:`~pydicom.pixels.processing.apply_voi_lut` PresentationLUTSequencerr2r6rHru2r7Nr3rPfloat32ruint16rPresentationLUTShape)IDENTITYINVERSEz1A (2050,0020) 'Presentation LUT Shape' value of 'z' is not supportedr)rrr USr9rgvaluer:rflags writeablecopy left_shift right_shiftrminmaxrrstripupperNotImplementedError) rrrrsrt bit_depthitemsizenr_byteselemr bit_shiftrr=r=r>apply_presentation_luts>         rindex prefer_lutcCsd}|dr$ttd|j|_d|jddd|jdddfv}d|dd|d dfv}|rD|rD|r>t|||St|||S|rLt|||S|rTt|||S|S) aApply a VOI lookup table or windowing operation to `arr`. .. versionchanged:: 2.1 Added the `prefer_lut` keyword parameter Parameters ---------- arr : numpy.ndarray The :class:`~numpy.ndarray` to apply the VOI LUT or windowing operation to. ds : dataset.Dataset A dataset containing a :dcm:`VOI LUT Module`. If (0028,3010) *VOI LUT Sequence* is present then returns an array of ``np.uint8`` or ``np.uint16``, depending on the 3rd value of (0028,3002) *LUT Descriptor*. If (0028,1050) *Window Center* and (0028,1051) *Window Width* are present then returns an array of ``np.float64``. If neither are present then `arr` will be returned unchanged. index : int, optional When the VOI LUT Module contains multiple alternative views, this is the index of the view to return (default ``0``). prefer_lut : bool When the VOI LUT Module contains both *Window Width*/*Window Center* and *VOI LUT Sequence*, if ``True`` (default) then apply the VOI LUT, otherwise apply the windowing operation. Returns ------- numpy.ndarray An array with applied VOI LUT or windowing operation. Notes ----- When the dataset requires a modality LUT or rescale operation as part of the Modality LUT module then that must be applied before any windowing operation. See Also -------- :func:`~pydicom.pixels.processing.apply_modality_lut` :func:`~pydicom.pixels.processing.apply_voi` :func:`~pydicom.pixels.processing.apply_windowing` References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.11.2 ` * DICOM Standard, Part 3, :dcm:`Annex C.8.11.3.1.5 ` * DICOM Standard, Part 4, :dcm:`Annex N.2.1.1 ` FVOILUTSequencer Nrrr WindowCenter WindowWidth)rrrWr apply_voiapply_windowing)rrrr valid_voivalid_windowingr=r=r> apply_voi_luts(8       rcCs||ds|St|jtjstdttd|j|}ttt |j }|dp*d}|d}|d}|tt dd vr?d }n|d krFd }nt d |d|dj t jkrt|drb|jjd} n|jd} | durvtdt|jdd| |d} t| tt|j} nttt |j} tj| |d} tj|j|d} ||k}|||| |<tj| d|d| dtd| | S)aApply a VOI lookup table to `arr`. .. versionadded:: 2.1 See :func:`~pydicom.pixels.processing.apply_voi_lut` for applying *Window Width*/*Window Center* as a fallback if no *VOI LUT Sequence* is present. Parameters ---------- arr : numpy.ndarray The :class:`~numpy.ndarray` to apply the VOI LUT to. ds : dataset.Dataset A dataset containing a :dcm:`VOI LUT Module`. If (0028,3010) *VOI LUT Sequence* is present then returns an array of ``np.uint8`` or ``np.uint16``, depending on the 3rd value of (0028,3002) *LUT Descriptor*, otherwise `arr` will be returned unchanged. index : int, optional When the VOI LUT Module contains multiple alternative views, this is the index of the view to return (default ``0``). Returns ------- numpy.ndarray An array with applied VOI LUT. See Also -------- :func:`~pydicom.pixels.processing.apply_modality_lut` :func:`~pydicom.pixels.processing.apply_windowing` :func:`~pydicom.pixels.processing.apply_voi_lut` References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.11.2 ` * DICOM Standard, Part 3, :dcm:`Annex C.8.11.3.1.5 ` * DICOM Standard, Part 4, :dcm:`Annex N.2.1.1 ` rzDApplying a VOI LUT on a float input array may give incorrect resultsr rr2rr rr6uint8r,z%' bits per LUT entry is not supportedrrBNrCrDrErGr7rPr)rr9 issubdtyper8integerr rrWrrXrrangerr rr_rBr`rarbrcrdrrYrrgrirjrl)rrrrlut_descriptorrsrtrur8r|rrrrrr=r=r>rGsJ *    rcCsd|vr d|vr |S|jdvrtdttt|dd}|d}|jdkr1ttt|j |n|j }tt|}|d}|jdkrLttt|j |n|j }tt|}tt |j |_ | dr}d }ttd |j d }ttt |jd } d | d} n |jd krd }d |j d} nd |j d }d |j dd} | d d } | dd } | d ur| d urtt|j|_tt|j|_||j|j}| |j|j} | |} |d}|dvr1|dkr|dkrtd|d8}|d8}n|d krtd|||d k}|||d k}t||}|||<| ||<|r/||||d| |||<|S|dkrR|d kr?td| dtd||||}|Std|d)aApply a windowing operation to `arr`. .. versionadded:: 2.1 Parameters ---------- arr : numpy.ndarray The :class:`~numpy.ndarray` to apply the windowing operation to. ds : dataset.Dataset A dataset containing a :dcm:`VOI LUT Module`. If (0028,1050) *Window Center* and (0028,1051) *Window Width* are present then returns an array of ``np.float64``, otherwise `arr` will be returned unchanged. index : int, optional When the VOI LUT Module contains multiple alternative views, this is the index of the view to return (default ``0``). Returns ------- numpy.ndarray An array with applied windowing operation. Notes ----- When the dataset requires a modality LUT or rescale operation as part of the Modality LUT module then that must be applied before any windowing operation. See Also -------- :func:`~pydicom.pixels.processing.apply_modality_lut` :func:`~pydicom.pixels.processing.apply_voi` :func:`~pydicom.pixels.processing.apply_voi_lut` References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.11.2 ` * DICOM Standard, Part 3, :dcm:`Annex C.8.11.3.1.5 ` * DICOM Standard, Part 4, :dcm:`Annex N.2.1.1 ` rr) MONOCHROME1 MONOCHROME2zWhen performing a windowing operation only 'MONOCHROME1' and 'MONOCHROME2' are allowed for (0028,0004) Photometric InterpretationVOILUTFunctionLINEARrrrr rrNrr)r LINEAR_EXACTzbThe (0028,1051) Window Width must be greater than or equal to 1 for a 'LINEAR' windowing operation?z\The (0028,1051) Window Width must be greater than 0 for a 'LINEAR_EXACT' windowing operationSIGMOIDzWThe (0028,1051) Window Width must be greater than 0 for a 'SIGMOID' windowing operationz0Unsupported (0028,1056) VOI LUT Function value 'r,)PhotometricInterpretationrRrstrrVrVMrWrrrX BitsStoredrrrPixelRepresentationrrrr9 logical_andanyexp)rrrvoi_funcrcenterwidthy_minrry_maxslope intercepty_rangebelowabovebetweenr=r=r>rsz, $ $           "rcurrentdesired per_framec Cs|jtdkrtd|jdddd}||td ||td |ttd d }z|||}WntyBtd |d |dwt|jdkr]|r]t |D] \}}||||<qP|S||S)aConvert the image(s) in `arr` from one color space to another. .. versionchanged:: 2.2 Added `per_frame` keyword parameter. Parameters ---------- arr : numpy.ndarray The image(s) as :class:`numpy.ndarray` with :attr:`~numpy.ndarray.shape` (frames, rows, columns, 3) or (rows, columns, 3) and a 'uint8' :class:`~numpy.dtype` (unsigned 8-bit). current : str The current color space, should be a valid value for (0028,0004) *Photometric Interpretation*. One of ``'RGB'``, ``'YBR_FULL'``, ``'YBR_FULL_422'``. desired : str The desired color space, should be a valid value for (0028,0004) *Photometric Interpretation*. One of ``'RGB'``, ``'YBR_FULL'``, ``'YBR_FULL_422'``. per_frame : bool, optional If ``True`` and the input array contains multiple frames then process each frame individually and update `arr` in-place to reduce memory usage. Default ``False``. Returns ------- numpy.ndarray The image(s) converted to the desired color space. If `per_frame` is ``False`` (the default) then a new :class:`~numpy.ndarray` will be returned, otherwise `arr` will be updated in-place. References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.7.6.3.1.2` * ISO/IEC 10918-5:2012 (`ITU T.871 `_), Section 7 u1zInvalid ndarray.dtype 'z>' for color space conversion, must be 'uint8' or an equivalentrrrcSs|S)Nr=)rr=r=r> _no_changeesz'convert_color_space.._no_change) YBR_FULL_422YBR_FULLRGB)rrr)rrrzConversion from z to z is not supported.rNrrrr) r8r9rR_convert_YBR_FULL_to_RGB_convert_RGB_to_YBR_FULLrTrrJrjrn) rrrrr _converters converterrrr=r=r>convert_color_space3s<,   rcCsp|j}tjgdgdgdgtjd}tj||tjd}|gd7}tj||dtj|dd|d||S) aReturn an ndarray converted from RGB to YBR_FULL color space. Parameters ---------- arr : numpy.ndarray An ndarray of an 8-bit per channel images in RGB color space. Returns ------- numpy.ndarray The array in YBR_FULL color space. References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.7.6.3.1.2` * ISO/IEC 10918-5:2012 (`ITU T.871 `_), Section 7 )gA`"?g$@#ſr)gbX9?g_n3տg{Fڿ)gv/?rgRдr7)r`@rrPr)r8r9rgrmatmulfloorrlr)r orig_dtype rgb_to_ybrr=r=r>rs  rcCs|j}tjgdgdgdgtjd}|tj}|gd8}tj|||d|d7}tj||dtj|dd |d||S) aReturn an ndarray converted from YBR_FULL to RGB color space. Parameters ---------- arr : numpy.ndarray An ndarray of an 8-bit per channel images in YBR_FULL color space. Returns ------- numpy.ndarray The array in RGB color space. References ---------- * DICOM Standard, Part 3, :dcm:`Annex C.7.6.3.1.2` * ISO/IEC 10918-5:2012, Section 7 )?rr)gk3Tֿg'1Z?)g;On?gbCW4rr7)rrrPrrr)r8r9rgrrrrrl)rr ybr_to_rgbr=r=r>rs  r icc_profilezPIL.ImageCms.ImageCmsTransformcCststd|dur|std|r|rtdt|d|}|s$td|r(|nt|dd}|rL|tvrL|sBd |d }t|d |d }t|t|}tt|}|r]|nt |}|t vrntd |dtj t ||ddt |dS)amReturn a Pillow color transformation object from either the dataset `ds` or an ICC profile `icc_profile`. .. versionadded:: 3.0 .. warning:: This function requires `NumPy `_ and `Pillow `_. Parameters ---------- ds : pydicom.dataset.Dataset, optional Required if `icc_profile` is not supplied, a :class:`~pydicom.dataset.Dataset` containing elements from the :dcm:`ICC Profile` module. icc_profile : bytes, optional Required if `ds` is not supplied, an ICC profile encoded as :class:`bytes`. intent : int, optional The rendering intent of the transformation, one of: * ``0``: perceptual * ``1``: relative colorimetric * ``2``: saturation * ``3``: absolute colorimetric If no `intent` is specified then the default rendering intent in the ICC profile will be used. color_space : str, optional The output color space to use for the transformation created from `ds`. If not used then defaults to the (0028,2002) *Color Space* value if its present in `ds`, otherwise defaults to ``"sRGB"``. If used then it will override the (0028,2002) *Color Space* element value (if present). Must be one of the :func:`supported Pillow color spaces`, although in practice only ``"sRGB"`` is likely to be used. Returns ------- PIL.ImageCms.ImageCmsTransform A color transformation object that can be used with :func:`apply_icc_profile`. z:Pillow is required to create a color transformation objectNz-Either 'ds' or 'icc_profile' must be suppliedz;Only one of 'ds' and 'icc_profile' should be used, not both ICCProfilez6No (0028,2000) 'ICC Profile' element was found in 'ds' ColorSpacerz%The (0028,2002) 'Color Space' value 'z`' is not supported by Pillow, please use the 'color_space' argument to specify a supported valuez!Unsupported 'color_space' value 'z!', must be 'sRGB', 'LAB' or 'XYZ'zInvalid 'intent' value 'z', must be 0, 1, 2 or 3r) outputProfile inputProfileinModeoutModerenderingIntent)rrrRrVr_CMS_COLOR_SPACESrImageCmsProfilergetDefaultIntent _CMS_INTENTSbuildTransform createProfileIntent)rrrrcsmsgprofiler=r=r>rs<.     rdata.r nr_segments last_valuecsd|vrgdngd}g}dd}dtkr}d}d7|dkr<|||7n|dkr|rG|d} n |rL|} ntd} d7| | kre|| g|ny| | |} tt| | | |} |d d | Dn\|dkr|std d |vrfd d |D} | dd>| dBd>| dd>| dBB}d7ndd>B}d7|t|d|||dntd|d|d7}||kr|Sdtks|S)aReturn a list containing the expanded lookup table data. Parameters ---------- data : tuple of int The decoded segmented palette lookup table data. May be padded by a trailing null. fmt : str The format of the data, should contain `'B'` for 8-bit, `'H'` for 16-bit, `'<'` for little endian and `'>'` for big endian. nr_segments : int, optional Expand at most `nr_segments` from the data. Should be used when the opcode is ``2`` (indirect). If used then `last_value` should also be used. last_value : int, optional The previous value in the expanded lookup table. Should be used when the opcode is ``2`` (indirect). If used then `nr_segments` should also be used. Returns ------- list of int The reconstructed lookup table data. References ---------- * DICOM Standard, Part 3, Annex C.7.9 <)rrrr)rrrrrrrrzdError expanding a segmented palette color lookup table: the first segment cannot be a linear segmentcSrIr=)rXr;vvr=r=r>rMrNz)_expand_segmented_lut..zgError expanding a segmented palette color lookup table: the first segment cannot be an indirect segmentrFcsg|]}|qSr=r=rroffsetr=r>rMsr6rHrrNzHError expanding a segmented palette lookup table: unknown segment type 'r,)rJr\rRr9aroundlinspacere)rrrr indirect_iir segments_readopcodelengthy0y1stepvalsrL byte_offsetr=r r>re?s`%    , "s              E c VL Se  V ,. `