Wrap GMT_Read_Data

pygmt/clib/session.py

@@ -1068,6 +1068,49 @@ def put_matrix(self, dataset, matrix, pad=0):
         if status != 0:
             raise GMTCLibError(f"Failed to put matrix of type {matrix.dtype}.")
 
+    def read_data(self, family: str, method, geometry, mode, wesn, input, data):
+        """
+        Read a data file into a GMT data container.
+
+        Wraps ``GMT_Read_Data`` but only allows reading from a file.
+
+        Parameters
+        ----------
+        family
+            A valid GMT data family name (e.g., ``'GMT_IS_DATASET'``). See the
+            ``FAMILIES`` attribute for valid names.
+
+        """
+        c_read_data = self.get_libgmt_func(
+            "GMT_Read_Data",
+            argtypes=[
+                ctp.c_void_p,
+                ctp.c_uint,
+                ctp.c_uint,
+                ctp.c_uint,
+                ctp.c_uint,
+                ctp.POINTER(ctp.c_double),
+                ctp.c_char_p,
+                ctp.c_void_p,
+            ],
+            restype=ctp.c_void_p,
+        )
+        family_int = self._parse_constant(family, valid=FAMILIES, valid_modifiers=VIAS)
+        geometry_int = self._parse_constant(geometry, valid=GEOMETRIES)
+
+        pointer = c_read_data(
+            self.session_pointer,
+            family_int,
+            self["GMT_IS_FILE"],
+            geometry_int,
+            self[mode],
+            None,
+            # sequence_to_ctypes_array(wesn, ctp.c_double, 6),
+            input.encode(),
+            data,
+        )
+        return pointer
+
     def write_data(self, family, geometry, mode, wesn, output, data):
         """
         Write a GMT data container to a file.