Added STAC generator

.github/workflows/python-test.yml

@@ -8,7 +8,7 @@ jobs:
     runs-on: ubuntu-latest
     strategy:
       matrix:
-        python-version: ["3.8", "3.9", "3.10", "3.11", "3.12"]
+        python-version: ["3.9", "3.10", "3.11", "3.12"]
 
     steps:
       - uses: actions/checkout@main

README.md

@@ -4,7 +4,7 @@
 
 # PAINT
 
-[![](https://img.shields.io/badge/Python-3.8+-blue.svg)](https://www.python.org/downloads/)
+[![](https://img.shields.io/badge/Python-3.9+-blue.svg)](https://www.python.org/downloads/)
 [![Code style: black](https://img.shields.io/badge/code%20style-black-000000.svg)](https://github.com/psf/black)
 [![](https://img.shields.io/badge/Contact-artist%40lists.kit.edu-orange?label=Contact)]([email protected])
 ![](./coverage.svg)
@@ -14,7 +14,7 @@
 PAINT is a FAIR database for Concentrating Solar Power plants (CSP).
 
 ## Installation
-We heavily recommend installing the `PAINT` package in a dedicated `Python3.8+` virtual environment. You can
+We heavily recommend installing the `PAINT` package in a dedicated `Python3.9+` virtual environment. You can
 install ``PAINT`` directly from the GitHub repository via:
 ```bash
 pip install git+https://github.com/ARTIST-Association/PAINT

paint/data/binary_extractor.py

@@ -0,0 +1,236 @@
+import json
+import struct
+from pathlib import Path
+from typing import Union
+
+import h5py
+import torch
+
+import paint.util.paint_mappings as mappings
+from paint.util.utils import to_utc_single
+
+
+class BinaryExtractor:
+    """
+    Implement an extractor that extracts data from a binary file and saves it to h5 and json.
+
+    This extractor considers data form a binary file that contains deflectometry data and heliostat properties. The data
+    is extracted and the deflectometry data saved in a h5 format and the heliostat properties as a json.
+
+    Attributes
+    ----------
+    input_path : Path
+        The file path to the binary data file that will be converted.
+    output_path : Path
+        The file path to save the converted h5 file.
+    file_name : str
+        The file name of the converted h5 file.
+    raw_data : bool
+        Whether the raw data or filled data is extracted.
+    heliostat_id : str
+        The heliostat ID of the heliostat considered in the binary file.
+    json_handle : str
+        The file path to save the json containing the heliostat properties data.
+    deflectometry_created_at : str
+        The time stamp for when the deflectometry data was created. Required for properties later.
+    surface_header_name : str
+        The name for the surface header in the binary file.
+    facet_header_name : str
+        The name for the facet header in the binary file.
+    points_on_facet_struct_name : str
+        The name of the point on facet structure in the binary file.
+
+    Methods
+    -------
+    nwu_to_enu()
+        Cast from an NWU to an ENU coordinate system.
+    convert_to_h5_and_extract_properties()
+        Convert binary data to h5 and extract heliostat properties not to be saved in the deflectometry file.
+    """
+
+    def __init__(
+        self,
+        input_path: Union[str, Path],
+        output_path: Union[str, Path],
+        surface_header_name: str,
+        facet_header_name: str,
+        points_on_facet_struct_name: str,
+    ) -> None:
+        """
+        Initialize the extractor.
+
+        Parameters
+        ----------
+        input_path : Union[str, Path]
+            The file path to the binary data file that will be converted.
+        output_path : Union[str, Path]
+            The file path to save the converted h5 deflectometry file.
+        surface_header_name : str
+            The name for the surface header in the binary file.
+        facet_header_name : str
+            The name for the facet header in the binary file.
+        points_on_facet_struct_name : str
+            The name of the point on facet structure in the binary file.
+        """
+        self.input_path = Path(input_path)
+        self.output_path = Path(output_path)
+        name_string = self.input_path.name.split("_")
+        if len(name_string) == 6:
+            file_name = (
+                name_string[1]
+                + "-"
+                + name_string[4]
+                + "-"
+                + str(to_utc_single(name_string[-1].split(".")[0]))
+            )
+            self.raw_data = False
+        else:
+            file_name = (
+                name_string[1] + "-" + str(to_utc_single(name_string[-1].split(".")[0]))
+            )
+            self.raw_data = True
+        self.heliostat_id = name_string[1]
+        self.file_name = file_name + mappings.DEFLECTOMETRY_SUFFIX
+        self.json_handle = name_string[1] + mappings.FACET_PROPERTIES_SUFFIX
+        self.deflectometry_created_at = to_utc_single(name_string[-1].split(".")[0])
+        self.surface_header_name = surface_header_name
+        self.facet_header_name = facet_header_name
+        self.points_on_facet_struct_name = points_on_facet_struct_name
+
+    @staticmethod
+    def nwu_to_enu(nwu_tensor: torch.Tensor) -> torch.Tensor:
+        """
+        Cast the coordinate system from NWU to ENU.
+
+        Parameters
+        ----------
+        nwu_tensor : torch.Tensor
+            The tensor in the NWU coordinate system.
+
+        Returns
+        -------
+        torch.Tensor
+            The converted tensor in the ENU coordinate system.
+        """
+        return torch.tensor(
+            [-nwu_tensor[1], nwu_tensor[0], nwu_tensor[2]], dtype=torch.float
+        )
+
+    def convert_to_h5_and_extract_properties(
+        self,
+    ) -> None:
+        """
+        Extract data from a binary file and save the deflectometry measurements and heliostat properties.
+
+        The binary files we consider, contain both deflectometry measurements and certain heliostat properties, such as
+        the number of facets, the facet translation vectors, and the facet canting vectors. Therefore, the deflectometry
+        measurements are extracted and saved as a h5 file, whilst the heliostat properties are extracted and saved in a
+        json file.
+        """
+        # Create structures for reading binary file correctly.
+        surface_header_struct = struct.Struct(self.surface_header_name)
+        facet_header_struct = struct.Struct(self.facet_header_name)
+        points_on_facet_struct = struct.Struct(self.points_on_facet_struct_name)
+
+        with open(self.input_path, "rb") as file:
+            surface_header_data = surface_header_struct.unpack_from(
+                file.read(surface_header_struct.size)
+            )
+
+            # Calculate the number of facets.
+            n_xy = surface_header_data[5:7]
+            number_of_facets = int(n_xy[0] * n_xy[1])
+
+            # Create empty tensors for storing data.
+            facet_translation_vectors = torch.empty(number_of_facets, 3)
+            canting_e = torch.empty(number_of_facets, 3)
+            canting_n = torch.empty(number_of_facets, 3)
+            surface_points_with_facets = torch.empty(0)
+            surface_normals_with_facets = torch.empty(0)
+            for f in range(number_of_facets):
+                facet_header_data = facet_header_struct.unpack_from(
+                    file.read(facet_header_struct.size)
+                )
+
+                facet_translation_vectors[f] = torch.tensor(
+                    facet_header_data[1:4], dtype=torch.float
+                )
+                canting_n[f] = self.nwu_to_enu(
+                    torch.tensor(
+                        facet_header_data[4:7],
+                        dtype=torch.float,
+                    )
+                )
+                canting_e[f] = self.nwu_to_enu(
+                    torch.tensor(
+                        facet_header_data[7:10],
+                        dtype=torch.float,
+                    )
+                )
+                number_of_points = facet_header_data[10]
+                if f == 0:
+                    surface_points_with_facets = torch.empty(
+                        number_of_facets, number_of_points, 3
+                    )
+                    surface_normals_with_facets = torch.empty(
+                        number_of_facets, number_of_points, 3
+                    )
+
+                points_data = points_on_facet_struct.iter_unpack(
+                    file.read(points_on_facet_struct.size * number_of_points)
+                )
+                for i, point_data in enumerate(points_data):
+                    surface_points_with_facets[f, i, :] = torch.tensor(
+                        point_data[:3], dtype=torch.float
+                    )
+                    surface_normals_with_facets[f, i, :] = torch.tensor(
+                        point_data[3:6], dtype=torch.float
+                    )
+
+        # to maintain consistency, we cast the west direction to east direction
+        canting_e[:, 0] = -canting_e[:, 0]
+
+        # extract deflectometry data and save
+        saved_deflectometry_path = (
+            Path(self.output_path)
+            / self.heliostat_id
+            / mappings.SAVE_DEFLECTOMETRY
+            / self.file_name
+        )
+        saved_deflectometry_path.parent.mkdir(parents=True, exist_ok=True)
+        with h5py.File(saved_deflectometry_path, "w") as file:
+            for i in range(number_of_facets):
+                facet = file.create_group(name=f"{mappings.FACET_KEY}{i+1}")
+                facet.create_dataset(
+                    name=f"{mappings.SURFACE_NORMAL_KEY}",
+                    data=surface_normals_with_facets[i, :, :],
+                )
+                facet.create_dataset(
+                    name=f"{mappings.SURFACE_POINT_KEY}",
+                    data=surface_points_with_facets[i, :, :],
+                )
+
+        # extract facet properties data and save
+        saved_facet_path = (
+            Path(self.output_path)
+            / self.heliostat_id
+            / mappings.SAVE_PROPERTIES
+            / self.json_handle
+        )
+        saved_facet_path.parent.mkdir(parents=True, exist_ok=True)
+        if self.raw_data:
+            with open(saved_facet_path, "w") as handle:
+                properties = {
+                    mappings.NUM_FACETS: number_of_facets,
+                    mappings.FACETS_LIST: [
+                        {
+                            mappings.TRANSLATION_VECTOR: facet_translation_vectors[
+                                i, :
+                            ].tolist(),
+                            mappings.CANTING_E: canting_e[i, :].tolist(),
+                            mappings.CANTING_N: canting_n[i, :].tolist(),
+                        }
+                        for i in range(number_of_facets)
+                    ],
+                }
+                json.dump(properties, handle)