get header dtype from spec instead of reading a header

TGSAI · Jun 24, 2024 · b209e4e · b209e4e
1 parent 9d49bc8
commit b209e4e
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/mdio/segy/blocked_io.py b/src/mdio/segy/blocked_io.py
@@ -100,7 +100,8 @@ def to_zarr(
         **kwargs,
     )
 
-    header_dtype = segy_file.header[0].dtype
+    # Get header dtype in native order (little-endian 99.9% of the time)
+    header_dtype = segy_file.spec.trace.header.dtype.newbyteorder("=")
     header_array = metadata_root.create_dataset(
         name="_".join([name, "trace_headers"]),
         shape=grid.shape[:-1],  # Same spatial shape as data