@@ -204,7 +204,12 @@ def data(self, data=None, offset=None, size=None, shape=None, as_memoryview=Fals
             raise ValueError("Size is larger than raster.")
 
         # Create ctypes type array generator
-        ctypes_array = GDAL_TO_CTYPES[self.datatype()] * (shape[0] * shape[1])
+        datatype = self.datatype()
+        if datatype == 15:
+            raise NotImplementedError(
+                "16 bit floating point (GDT_Float16) is not yet supported."
+            )
+        ctypes_array = GDAL_TO_CTYPES[datatype] * (shape[0] * shape[1])
 
         if data is None:
             # Set read mode
 
@@ -32,6 +32,8 @@
     12: "GDT_UInt64",  # 64 bit unsigned integer (GDAL 3.5+).
     13: "GDT_Int64",  # 64 bit signed integer (GDAL 3.5+).
     14: "GDT_Int8",  # 8 bit signed integer (GDAL 3.7+).
+    15: "GDT_Float16",  # 16 bit floating point (GDAL 3.11+).
+    16: "GDT_CFloat16",  # Complex Float16 (GDAL 3.11+).
 }
 
 # A list of gdal datatypes that are integers.
@@ -57,6 +59,8 @@
     c_uint64,
     c_int64,
     c_int8,
+    None,  # _Float16 ticket https://github.com/python/cpython/issues/153740
+    None,
 ]
 
 # List of resampling algorithms that can be used to warp a GDALRaster.
 
@@ -1006,3 +1006,26 @@ def test_band_data_replication(self):
                 )
             else:
                 self.assertEqual(band.data(), list(combo[2]))
+
+    def test_16bit_float(self):
+        rstfile = NamedTemporaryFile(suffix=".tif")
+        self.addCleanup(rstfile.close)
+        rst = GDALRaster(
+            {
+                "driver": "GTiff",
+                "name": rstfile.name,
+                "srid": 4326,
+                "width": 255,
+                "height": 255,
+                "nr_of_bands": 1,
+                "papsz_options": dict(NBITS=16),
+            }
+        )
+        if GDAL_VERSION < (3, 11):
+            self.assertEqual(rst.bands[0].datatype(as_string=True), "GDT_Float32")
+        else:
+            with self.assertRaisesMessage(
+                NotImplementedError,
+                "16 bit floating point (GDT_Float16) is not yet supported.",
+            ):
+                rst.bands[0].data()