xtgeo 4.14.0__cp311-cp311-manylinux_2_26_aarch64.manylinux_2_28_aarch64.whl

xtgeo/xyz/_xyz_io.py

@@ -0,0 +1,737 @@
+"""Private import and export routines for XYZ stuff."""
+
+from __future__ import annotations
+
+import contextlib
+from copy import deepcopy
+
+import numpy as np
+import pandas as pd
+
+from xtgeo.common._xyz_enum import _AttrName, _XYZType
+from xtgeo.common.constants import UNDEF, UNDEF_INT
+from xtgeo.common.exceptions import InvalidFileFormatError
+from xtgeo.common.log import null_logger
+from xtgeo.io._file import FileFormat, FileWrapper
+
+logger = null_logger(__name__)
+
+
+def import_xyz(pfile, zname=_AttrName.ZNAME.value):
+    """Simple X Y Z file. All points as Pandas framework."""
+    return {
+        "zname": zname,
+        "xname": _AttrName.XNAME.value,
+        "yname": _AttrName.YNAME.value,
+        "values": pd.read_csv(
+            pfile.file,
+            sep=r"\s+",
+            skiprows=0,
+            header=None,
+            names=[_AttrName.XNAME.value, _AttrName.YNAME.value, zname],
+            dtype=np.float64,
+            na_values=999.00,
+        ),
+    }
+
+
+def _import_table(pfile, xyztype, file_format="csv"):
+    """Simple CSV or Parquet file. Generic reader for Points and Polygons.
+
+    Table files (CSV, ...) can have some variants, e.g.:
+    * 3 columns, X Y Z
+    * 4 columns, X Y Z POLYID  (if polygon)
+    * 4 columns, <index>, X Y Z
+    * 5 columns, <index>, X Y Z POLYID
+    * 5+ columns, <index>, X Y Z POLYID <name> <attr1> <attr2> ...
+    * 4+ columns, <index>, X Y Z <attr1> <attr2> ... (points)
+    """
+
+    # First do a raw read to get the columns and check if it is a valid CSV
+    # file. If not, then raise an error.
+    if file_format == "parquet":
+        try:
+            dataframe = pd.read_parquet(pfile.file)
+        except ValueError as verr:
+            # If the file is not a valid Parquet
+            raise IOError(f"File is not a valid Parquet file: {verr}")
+    else:
+        try:
+            dataframe = pd.read_csv(pfile.file, comment="#")
+        except ValueError as verr:
+            # If the file is not a valid CSV
+            raise IOError(f"File is not a valid CSV file: {verr}")
+
+    columns = dataframe.columns
+    if "Unnamed" in columns[0]:
+        # The first column is an index column which is removed
+        dataframe = dataframe.iloc[:, 1:]
+        columns = dataframe.columns
+    ncol = len(columns)
+
+    if ncol == 3:
+        if xyztype == _XYZType.POLYGONS.value:
+            # If the file is a polygon file, add a POLY_ID column
+            dataframe["POLY_ID"] = 0
+            return {
+                "xname": columns[0],
+                "yname": columns[1],
+                "zname": columns[2],
+                "pname": _AttrName.PNAME.value,
+                "values": dataframe,
+            }
+
+        # points with 3 columns, X Y Z
+        return {
+            "xname": columns[0],
+            "yname": columns[1],
+            "zname": columns[2],
+            "values": dataframe,
+        }
+    if xyztype == _XYZType.POLYGONS.value and ncol == 4:
+        return {
+            "xname": columns[0],
+            "yname": columns[1],
+            "zname": columns[2],
+            "pname": columns[3],
+            "values": dataframe,
+        }
+    if xyztype == _XYZType.POLYGONS.value and ncol > 4:
+        # need to infer the attrs from column 5...
+        attr_names = columns[4:]
+        attrs = {}
+        for attr in attr_names:
+            # Check if the column is numeric
+            if pd.api.types.is_numeric_dtype(dataframe[attr]):
+                # next infer if it is a float or int
+                if pd.api.types.is_float_dtype(dataframe[attr]):
+                    attrs[attr] = "float"
+                elif pd.api.types.is_integer_dtype(dataframe[attr]):
+                    attrs[attr] = "int"
+            else:
+                attrs[attr] = "str"
+
+        return {
+            "xname": columns[0],
+            "yname": columns[1],
+            "zname": columns[2],
+            "pname": columns[3],
+            "attributes": attrs,
+            "values": dataframe,
+        }
+    if xyztype == _XYZType.POINTS.value and ncol > 3:
+        # need to infer the attrs from column 5...
+        attr_names = columns[3:]
+        attrs = {}
+        for attr in attr_names:
+            # Check if the column is numeric
+            if pd.api.types.is_numeric_dtype(dataframe[attr]):
+                # next infer if it is a float or int
+                if pd.api.types.is_float_dtype(dataframe[attr]):
+                    attrs[attr] = "float"
+                elif pd.api.types.is_integer_dtype(dataframe[attr]):
+                    attrs[attr] = "int"
+            else:
+                attrs[attr] = "str"
+
+        return {
+            "xname": columns[0],
+            "yname": columns[1],
+            "zname": columns[2],
+            "attributes": attrs,
+            "values": dataframe,
+        }
+    return None
+    # More than 4 columns, assume that the first three are X Y Z
+
+
+def import_csv_polygons(pfile):
+    return _import_table(pfile, _XYZType.POLYGONS.value, file_format="csv")
+
+
+def import_csv_points(pfile):
+    return _import_table(pfile, _XYZType.POINTS.value, file_format="csv")
+
+
+def import_parquet_polygons(pfile):
+    return _import_table(pfile, _XYZType.POLYGONS.value, file_format="parquet")
+
+
+def import_parquet_points(pfile):
+    return _import_table(pfile, _XYZType.POINTS.value, file_format="parquet")
+
+
+def import_zmap(pfile, zname=_AttrName.ZNAME.value):
+    """The zmap ascii polygon format; not sure about all details."""
+    # ...seems that I just
+    # take in the columns starting from @<blank> line as is.
+    # Potential here to improve...
+
+    #
+    # !
+    # ! File exported from RMS.
+    # !
+    # ! Project:
+    # ! Date:          2017-11-07T17:22:30
+    # !
+    # ! Polygons/points Z-MAP file generated for ''.
+    # ! Coordinate system is ''.
+    # !
+    # !------------------------------------------------------------------
+    # @FREE POINT        , DATA, 80, 1
+    # X (EASTING)        , 1, 1,  1,      1, 20,,    1.0E+30,,,   4, 0
+    # Y (NORTHING)       , 2, 2,  1,     21, 40,,    1.0E+30,,,   4, 0
+    # Z VALUE            , 3, 3,  1,     41, 60,,    1.0E+30,,,   4, 0
+    # SEG I.D.           , 4, 35, 1,     61, 70,,    1.0E+30,,,   0, 0
+    # @
+    #    457357.781250      6782685.500000      1744.463379         0
+    #    457359.343750      6782676.000000      1744.482056         0
+    #    457370.906250      6782606.000000      1744.619507         0
+    #    457370.468750      6782568.500000      1745.868286         0
+
+    xname = "X_UTME"
+    yname = "Y_UTMN"
+    pname = "POLY_ID"
+
+    dtype = {
+        xname: np.float64,
+        yname: np.float64,
+        zname: np.float64,
+        pname: np.int32,
+    }
+
+    df = pd.read_csv(
+        pfile.file,
+        sep=r"\s+",
+        skiprows=16,
+        header=None,
+        names=[xname, yname, zname, pname],
+        dtype=dtype,
+        na_values=1.0e30,
+    )
+
+    return {"xname": xname, "yname": yname, "zname": zname, "values": df}
+
+
+def import_rms_attr(pfile, zname="Z_TVDSS"):
+    """The RMS ascii file Points format with attributes.
+
+    It appears that the the RMS attributes format is supported for Points only,
+    hence Polygons is not admitted.
+
+    Example::
+
+       Discrete  FaultBlock
+       String    FaultTag
+       Float     VerticalSep
+       519427.941  6733887.914  1968.988     6  UNDEF  UNDEF
+       519446.363  6732037.910  1806.782    19  UNDEF  UNDEF
+       519446.379  6732137.910  1795.707    19  UNDEF  UNDEF
+
+    Returns a kwargs list with the following items:
+        xname
+        yname
+        zname
+        values as a valid dataframe
+        attributes
+
+    Important notes from RMS manual and reverse engineering:
+
+    * For discrete numbers use 'Discrete' or 'Integer', not 'Int'
+    * For Discrete/Integer/Float both UNDEF and -999 will mark as undefined
+    * For Discrete/Integer, numbers less than -999 seems to accepted by RMS
+    * For String, use UNDEF only as undefined
+    """
+
+    kwargs = {}
+    _xn = kwargs["xname"] = "X_UTME"
+    _yn = kwargs["yname"] = "Y_UTMN"
+    _zn = kwargs["zname"] = zname
+
+    dtypes = {_xn: np.float64, _yn: np.float64, _zn: np.float64}
+
+    names = list(dtypes.keys())
+    _attrs = {}
+
+    # parse header
+    skiprows = 0
+    with open(pfile.file, "r") as rmsfile:
+        for iline in range(20):
+            fields = rmsfile.readline().split()
+            if len(fields) != 2:
+                skiprows = iline
+                break
+
+            dty, cname = fields
+            dtyx = None
+
+            # note that Pandas treats dtype str as object, cf:
+            # https://stackoverflow.com/questions/34881079
+            if dty == "Discrete":
+                dtyx = "int"
+            elif dty == "String":
+                dtyx = "str"
+            elif dty == "Float":
+                dtyx = "float"
+            elif dty == "Int":
+                dtyx = "int"
+            else:
+                dtyx = "str"
+            names.append(cname)
+            _attrs[cname] = dtyx
+
+    dfr = pd.read_csv(
+        pfile.file,
+        sep=r"\s+",
+        skiprows=skiprows,
+        header=None,
+        names=names,
+        dtype=dtypes,
+    )
+
+    for col in dfr.columns[3:]:
+        if col in _attrs:
+            # avoid FutureWarning by not using .replace() for mixed dtypes
+            if _attrs[col] == "float":
+                s = dfr[col].mask(dfr[col] == "UNDEF", UNDEF)
+                dfr[col] = pd.to_numeric(s, errors="coerce").astype(float)
+            elif _attrs[col] == "int":
+                s = dfr[col].mask(dfr[col] == "UNDEF", UNDEF_INT)
+                s = pd.to_numeric(s, errors="coerce").fillna(UNDEF_INT).astype(int)
+                dfr[col] = s
+        # cast to numerical if possible
+        with contextlib.suppress(ValueError, TypeError):
+            dfr[col] = pd.to_numeric(dfr[col])
+
+    kwargs["values"] = dfr
+    kwargs["attributes"] = _attrs
+
+    return kwargs
+
+
+def to_file(
+    xyz,  # Points | Polygons instance
+    pfile,
+    fformat="xyz",
+    attributes: bool | list[str] = False,
+    pfilter=None,
+    wcolumn=None,
+    hcolumn=None,
+    mdcolumn="M_MDEPTH",
+    **kwargs,
+):
+    """Export XYZ (Points/Polygons) to file.
+
+    Args:
+        pfile (str): Name of file
+        fformat (str): File format xyz/poi/pol/rms_attr/rms_wellpicks/csv/parquet,
+            dependent if Points or Polygons
+        attributes (bool or list): List of extra columns to export (some formats)
+            or True for all attributes present
+        pfilter (dict): Filter on e.g. top name(s) with keys TopName
+            or ZoneName as {'TopName': ['Top1', 'Top2']}. Only for points!
+        wcolumn (str): Name of well column (rms_wellpicks format only)
+        hcolumn (str): Name of horizons column (rms_wellpicks format only)
+        mdcolumn (str): Name of MD column (rms_wellpicks format only)
+
+    Returns:
+        Number of points exported
+
+    Note that the rms_wellpicks will try to output to:
+
+    * HorizonName, WellName, MD  if a MD (mdcolumn) is present,
+    * HorizonName, WellName, X, Y, Z  otherwise
+
+    Raises:
+        KeyError if pfilter is set and key(s) are invalid
+
+    """
+    filter_deprecated = kwargs.get("filter")
+    if filter_deprecated is not None and pfilter is None:
+        pfilter = filter_deprecated
+
+    pfile = FileWrapper(pfile)
+    pfile.check_folder(raiseerror=OSError)
+
+    ncount = 0
+    if xyz.get_dataframe(copy=False) is None:
+        logger.warning("Nothing to export!")
+        return ncount
+
+    if fformat is None or fformat in FileFormat.XYZ.value:
+        # NB! reuse export_rms_attr function, but no attributes
+        # are possible
+        ncount = export_rms_attr(xyz, pfile.name, attributes=False, pfilter=pfilter)
+
+    elif fformat in FileFormat.RMS_ATTR.value:
+        ncount = export_rms_attr(
+            xyz,
+            pfile.name,
+            attributes=attributes,
+            pfilter=pfilter,
+        )
+    elif fformat in FileFormat.CSV.value:
+        ncount = export_table(
+            xyz,
+            pfile.name,
+            attributes=attributes,
+            pfilter=pfilter,
+            file_format="csv",
+        )
+    elif fformat in FileFormat.PARQUET.value:
+        ncount = export_table(
+            xyz,
+            pfile.name,
+            attributes=attributes,
+            pfilter=pfilter,
+            file_format="parquet",
+        )
+    elif fformat == "rms_wellpicks":
+        ncount = export_rms_wpicks(xyz, pfile.name, hcolumn, wcolumn, mdcolumn=mdcolumn)
+    else:
+        extensions = FileFormat.extensions_string(
+            [FileFormat.XYZ, FileFormat.RMS_ATTR, FileFormat.PARQUET, FileFormat.CSV]
+        )
+        raise InvalidFileFormatError(
+            f"File format {fformat} is invalid for type Points or Polygons. "
+            f"Supported formats are {extensions}, 'rms_wellpicks'."
+        )
+
+    if ncount is None:
+        ncount = 0
+
+    if ncount == 0:
+        logger.warning("Nothing to export!")
+
+    return ncount
+
+
+def export_rms_attr(self, pfile, attributes=True, pfilter=None):
+    """Export til RMS attribute, also called RMS extended set.
+
+    If attributes is None, then it will be a simple XYZ file.
+
+    Attributes can be a bool or a list. If True, then use all attributes.
+
+    Filter is on the form {TopName: ['Name1', 'Name2']}
+
+    Returns:
+        The number of values exported. If value is 0; then no file
+        is made.
+    """
+
+    df = self.get_dataframe()
+
+    if not df.index.any():
+        logger.warning("Nothing to export")
+        return 0
+
+    columns = [self._xname, self._yname, self.zname]
+    df.fillna(value=999.0, inplace=True)
+
+    mode = "w"
+
+    transl = {"int": "Discrete", "float": "Float", "str": "String"}
+
+    logger.info("Attributes is %s", attributes)
+
+    # apply pfilter if any
+    if pfilter:
+        for key, val in pfilter.items():
+            if key in df.columns:
+                df = df.loc[df[key].isin(val)]
+            else:
+                raise KeyError(
+                    f"The requested pfilter key {key} was not found in dataframe. "
+                    f"Valid keys are {df.columns}"
+                )
+
+    if self._xyztype == _XYZType.POLYGONS.value:  # a bit weird: TODO fixup
+        if not attributes and self._pname in df.columns:
+            # need to convert the dataframe
+            df = _convert_idbased_xyz(self, df)
+    elif attributes:
+        use_attributes = list(self._attrs.keys())
+        if isinstance(attributes, bool):
+            all_attrs = list(self._attrs.keys())
+            logger.info("Use all attributes: %s", all_attrs)
+
+        elif attributes and isinstance(attributes, list):
+            all_attrs = list(self._attrs.keys())
+            # Ensure that the attributes are in the dataframe
+            for attr in attributes:
+                if attr not in all_attrs:
+                    raise ValueError(
+                        f"Attribute {attr} is not a valid attribute. "
+                        f"Valid attributes are: {all_attrs}"
+                    )
+            use_attributes = deepcopy(attributes)
+        else:
+            raise TypeError(
+                f"Attributes must be a bool or a list, not {type(attributes)}"
+            )
+
+        mode = "a"
+        columns += use_attributes
+        with open(pfile, "w") as fout:
+            for col in use_attributes:
+                if col in df.columns:
+                    fout.write(transl[self._attrs[col]] + " " + col + "\n")
+                    if self._attrs[col] == "int":
+                        df[col] = df[col].replace(UNDEF_INT, "UNDEF")
+                    elif self._attrs[col] == "float":
+                        df[col] = df[col].replace(UNDEF, "UNDEF")
+
+    elif not attributes:
+        df = df[[self._xname, self._yname, self.zname]]
+
+    with open(pfile, mode) as fc:
+        df.to_csv(fc, sep=" ", header=None, columns=columns, index=False)
+
+    return len(df.index)
+
+
+def export_table(self, pfile, file_format="csv", attributes=False, pfilter=None):
+    """Export to CSV or Parquet file.
+
+    Args:
+        pfile: Output file path.
+        file_format: "csv" or "parquet".
+        attributes: bool or list. If True, use all attributes.
+        pfilter: Optional filter dict.
+
+    Returns:
+        Number of rows exported.
+    """
+    df = self.get_dataframe(copy=True)
+
+    if not df.index.any():
+        logger.warning("Nothing to export")
+        return 0
+
+    # Fill NaNs for CSV (optional: you may want to skip for Parquet)
+    if file_format == "csv":
+        df = df.fillna(value=999.0)
+
+    # Apply filter if any (Points only)
+    if self._xyztype == _XYZType.POINTS.value and pfilter:
+        for key, val in pfilter.items():
+            if key in df.columns:
+                df = df.loc[df[key].isin(val)]
+            else:
+                raise KeyError(
+                    f"The requested pfilter key {key} was not found in dataframe. "
+                    f"Valid keys are {df.columns}"
+                )
+
+    # Select columns based on type and attributes
+    if self._xyztype == _XYZType.POLYGONS.value and not attributes:
+        df = df.iloc[:, 0:4]
+    elif self._xyztype == _XYZType.POINTS.value and not attributes:
+        df = df.iloc[:, 0:3]
+    elif attributes:
+        if isinstance(attributes, bool):
+            attributes = list(self._attrs.keys())
+            logger.info("Use all attributes: %s", attributes)
+            print("Attributes are: ", attributes)
+        elif isinstance(attributes, list):
+            logger.info("Use attributes: %s", attributes)
+            all_attrs = list(self._attrs.keys())
+
+            # Ensure that the attributes are in the dataframe
+            for attr in attributes:
+                if attr not in all_attrs:
+                    raise ValueError(
+                        f"Attribute {attr} is not a valid attribute. "
+                        f"Valid attributes are: {all_attrs}"
+                    )
+
+            for attr in all_attrs:
+                if attr not in attributes:
+                    # If the attribute is not in the list, remove it
+                    del df[attr]
+                    logger.debug(
+                        "Attribute %s is not in the export list, removing", attr
+                    )
+
+        else:
+            raise TypeError(
+                f"Attributes must be a bool or a list, not {type(attributes)}"
+            )
+
+        if self._xyztype == _XYZType.POLYGONS.value:
+            # Ensure POLY_ID is included
+            if self._pname not in df.columns:
+                df[self._pname] = 0
+            df = df[[self._xname, self._yname, self._zname, self._pname] + attributes]
+        else:
+            df = df[[self._xname, self._yname, self._zname] + attributes]
+    # Export
+    if file_format == "csv":
+        with open(pfile, "w") as fc:
+            df.to_csv(fc, index=False)
+    elif file_format == "parquet":
+        df.to_parquet(pfile, index=False)
+    else:
+        raise ValueError(f"Unsupported file format: {file_format}")
+
+    return len(df.index)
+
+
+def _convert_idbased_xyz(self, df):
+    """Conversion of format from ID column to 999 flag."""
+
+    # If polygons, there is a 4th column with POLY_ID. This needs
+    # to replaced by adding 999 line instead (for polygons)
+    # prior to XYZ export or when interactions in CXTGEO
+
+    idgroups = df.groupby(self._pname)
+
+    newdf = pd.DataFrame(
+        columns=[self._xname, self._yname, self._zname], dtype="float64"
+    )
+    udef = pd.DataFrame(
+        [[999.0, 999.0, 999.0]], columns=[self._xname, self._yname, self._zname]
+    )
+
+    for _id, gr in idgroups:
+        dfx = gr.drop(self._pname, axis=1)
+        newdf = pd.concat([newdf, dfx, udef], ignore_index=True)
+
+    return newdf
+
+
+def export_rms_wpicks(self, pfile, hcolumn, wcolumn, mdcolumn="M_MDEPTH"):
+    """Export til RMS wellpicks
+
+    If a MD column (mdcolumn) exists, it will use the MD
+
+    Args:
+        pfile (str): File to export to
+        hcolumn (str): Name of horizon/zone column in the point set
+        wcolumn (str): Name of well column in the point set
+        mdcolumn (str): Name of measured depht column (if any)
+    Returns:
+        The number of values exported. If value is 0; then no file
+        is made.
+
+    """
+
+    df = self.get_dataframe()
+
+    if not df.index.any():
+        logger.warning("Nothing to export")
+        return 0
+
+    columns = []
+
+    if hcolumn in df.columns:
+        columns.append(hcolumn)
+    else:
+        raise ValueError(f"Column for horizons/zones <{hcolumn}> not present")
+
+    if wcolumn in df.columns:
+        columns.append(wcolumn)
+    else:
+        raise ValueError(f"Column for wells <{wcolumn}> not present")
+
+    if mdcolumn in df.columns:
+        columns.append(mdcolumn)
+    else:
+        columns += [self._xname, self._yname, self._zname]
+
+    if not df.index.any():
+        logger.warning("Nothing to export")
+        return 0
+
+    with open(pfile, "w") as fc:
+        df.to_csv(fc, sep=" ", header=None, columns=columns, index=False)
+
+    return len(df.index)
+
+
+def _from_list_like(plist, zname, attrs, xyztype) -> pd.DataFrame:
+    """Import Points or Polygons from a list-like input.
+
+    The following 'list-like' inputs are possible:
+
+    * List of tuples [(x1, y1, z1, <id1>), (x2, y2, z2, <id2>), ...].
+    * List of lists  [[x1, y1, z1, <id1>], [x2, y2, z2, <id2>], ...].
+    * List of numpy arrays  [nparr1, nparr2, ...] where nparr1 is first row.
+    * A numpy array with shape [nrow, ncol], where ncol >= 3
+    * An existing pandas dataframe
+
+    Points scenaria:
+    * 3 columns, X Y Z
+    * 4 or more columns: rest columns are attributes
+
+    Polygons scenaria:
+    * 3 columns, X Y Z. Here P column is assigned 0 afterwards
+    * 4 or more columns:
+        - if totnum = lenattrs + 3 then POLY_ID is missing and will be made
+        - if totnum = lenattrs + 4 then assume that 4'th column is POLY_ID
+
+    It is currently not much error checking that lists/tuples are consistent, e.g.
+    if there always is either 3 or 4 elements per tuple, or that 4 number is
+    an integer.
+
+    Args:
+        plist (str): List of tuples, each tuple is length 3 or 4
+        zname (str): Name of third column
+        attrs (dict): Attributes, for Points
+        xyztype (str): POINTS/POLYGONS/...
+
+    Returns:
+        A valid datafram
+
+    Raises:
+        ValueError: If something is wrong with input
+
+    .. versionadded:: 2.16
+    """
+
+    dfr = None
+    if isinstance(plist, list):
+        plist = np.array(plist)
+
+    if isinstance(plist, np.ndarray):
+        logger.info("Process numpy to points")
+        if len(plist) == 0:
+            return pd.DataFrame([], columns=["X_UTME", "Y_UTMN", zname])
+
+        if plist.ndim != 2:
+            raise ValueError("Input numpy array must two-dimensional")
+        totnum = plist.shape[1]
+        lenattrs = len(attrs) if attrs is not None else 0
+        attr_first_col = 3
+        if totnum == 3 + lenattrs:
+            dfr = pd.DataFrame(plist[:, :3], columns=["X_UTME", "Y_UTMN", zname])
+            dfr = dfr.astype(float)
+            if xyztype == _XYZType.POLYGONS.value:
+                # pname column is missing but assign 0 as ID
+                dfr["POLY_ID"] = 0
+
+        elif totnum == 4 + lenattrs and xyztype == _XYZType.POLYGONS.value:
+            dfr = pd.DataFrame(
+                plist[:, :4],
+                columns=["X_UTME", "Y_UTMN", zname, "POLY_ID"],
+            )
+            attr_first_col = 4
+        else:
+            raise ValueError(
+                f"Wrong length detected of row: {totnum}. Are attributes set correct?"
+            )
+        dfr.dropna()
+        dfr = dfr.astype(np.float64)
+        if xyztype == _XYZType.POLYGONS.value:
+            dfr[_AttrName.PNAME.value] = dfr[_AttrName.PNAME.value].astype(np.int32)
+
+        if lenattrs > 0:
+            for enum, (key, dtype) in enumerate(attrs.items()):
+                dfr[key] = plist[:, attr_first_col + enum]
+                dfr = dfr.astype({key: dtype})
+
+    else:
+        raise TypeError("Not possible to make XYZ from given input")
+
+    return dfr