@@ -556,8 +556,6 @@ def read_segy_trace_header(trace_index: int) -> bytes:
556556 flat_mask = trace_mask .ravel ()
557557 flat_raw_headers = raw_headers_data .ravel () # Flatten to 1D array of 240-byte header records
558558
559- # operation = "w"
560-
561559 for grid_idx in range (flat_mask .size ):
562560 if not flat_mask [grid_idx ]:
563561 print (f"Skipping trace { grid_idx } because it is masked" )
@@ -573,58 +571,4 @@ def read_segy_trace_header(trace_index: int) -> bytes:
573571
574572 assert_array_equal (mdio_header_bytes , segy_header_bytes )
575573
576- # Compare byte-by-byte
577- # Write hexdumps to separate files for analysis
578- # def hexdump_to_string(data: bytes, title: str) -> str:
579- # """Create hexdump string."""
580- # lines = [f"{title}", "=" * len(title), ""]
581-
582- # for i in range(0, len(data), 16):
583- # # Address
584- # addr = i
585- # hex_part = ""
586- # ascii_part = ""
587-
588- # # Process 16 bytes at a time
589- # for j in range(16):
590- # if i + j < len(data):
591- # byte_val = data[i + j]
592- # hex_part += f"{byte_val:02x} "
593- # ascii_part += chr(byte_val) if 32 <= byte_val <= 126 else "."
594- # else:
595- # hex_part += " "
596- # ascii_part += " "
597-
598- # lines.append(f"{addr:08x}: {hex_part} |{ascii_part}|")
599-
600- # return "\n".join(lines)
601-
602- # # Generate filenames for this test case
603- # segy_filename = f"segy_headers_{grid_conf.name}.txt"
604- # mdio_filename = f"mdio_headers_{grid_conf.name}.txt"
605-
606- # # Append SEG-Y hexdump to file
607- # with open(segy_filename, operation) as f:
608- # if segy_trace_idx == 0:
609- # f.write("") # Start fresh for first trace
610- # else:
611- # f.write("\n\n") # Add spacing between traces
612- # f.write(hexdump_to_string(segy_header_bytes,
613- # f"SEG-Y Header - {grid_conf.name} Trace {segy_trace_idx} (240 bytes)"))
614-
615- # # Append MDIO hexdump to file
616- # with open(mdio_filename, operation) as f:
617- # if segy_trace_idx == 0:
618- # f.write("") # Start fresh for first trace
619- # else:
620- # f.write("\n\n") # Add spacing between traces
621- # f.write(hexdump_to_string(mdio_header_bytes,
622- # f"MDIO Raw Header - {grid_conf.name} Trace {segy_trace_idx} (240 bytes)"))
623- # operation = 'a'
624-
625- # if segy_trace_idx == 0:
626- # print(f"\nHeader hexdumps being written for {grid_conf.name}:")
627- # print(f" SEG-Y: {segy_filename}")
628- # print(f" MDIO: {mdio_filename}")
629-
630574 segy_trace_idx += 1
0 commit comments