@@ -33,18 +33,6 @@ def debug_compare_raw_vs_processed(segy_file, trace_index=0):
3333 # Get processed data with transforms
3434 processed_data = segy_file .header [trace_index ]
3535
36- print ("=== Raw vs Processed Comparison ===" )
37- print (f"Raw data shape: { raw_data .shape } " )
38- print (f"Processed data shape: { processed_data .shape } " )
39-
40- if hasattr (raw_data , 'dtype' ) and raw_data .dtype .names :
41- if 'inline_number' in raw_data .dtype .names :
42- print (f"Raw inline_number: { raw_data ['inline_number' ]} " )
43- print (f"Raw inline_number (hex): { raw_data ['inline_number' ]:08x} " )
44- print (f"Processed inline_number: { processed_data ['inline_number' ]} " )
45- print (f"Processed inline_number (hex): { processed_data ['inline_number' ]:08x} " )
46- print (f"Are they equal? { raw_data ['inline_number' ] == processed_data ['inline_number' ]} " )
47-
4836 return raw_data , processed_data
4937
5038
@@ -72,20 +60,6 @@ def __init__(self, segy_file: SegyFile):
7260 self .header_indexer = segy_file .header
7361 self .transform_pipeline = self .header_indexer .transform_pipeline
7462
75- # Debug: Print transform pipeline information
76- import sys
77- print (f"Debug: System endianness: { sys .byteorder } " )
78- print (f"Debug: File endianness: { self .segy_file .spec .endianness } " )
79- print (f"Debug: Transform pipeline has { len (self .transform_pipeline .transforms )} transforms:" )
80- for i , transform in enumerate (self .transform_pipeline .transforms ):
81- print (f" Transform { i } : { type (transform ).__name__ } " )
82- if hasattr (transform , 'target_order' ):
83- print (f" Target order: { transform .target_order } " )
84- if hasattr (transform , 'direction' ):
85- print (f" Direction: { transform .direction } " )
86- if hasattr (transform , 'keys' ):
87- print (f" Keys: { transform .keys } " )
88-
8963 def get_raw_and_transformed (
9064 self , indices : int | list [int ] | np .ndarray | slice
9165 ) -> tuple [NDArray , NDArray ]:
@@ -101,22 +75,9 @@ def get_raw_and_transformed(
10175 # This reads from filesystem and applies transforms
10276 transformed_data = self .header_indexer [indices ]
10377
104- print (f"Debug: Transformed data shape: { transformed_data .shape } " )
105- if hasattr (transformed_data , 'dtype' ) and transformed_data .dtype .names :
106- print (f"Debug: Transformed data dtype names: { transformed_data .dtype .names [:5 ]} ..." ) # First 5 fields
107- if 'inline_number' in transformed_data .dtype .names :
108- print (f"Debug: First transformed inline_number: { transformed_data ['inline_number' ][0 ]} " )
109- print (f"Debug: First transformed inline_number (hex): { transformed_data ['inline_number' ][0 ]:08x} " )
110-
11178 # Now reverse the transforms to get back to raw data
11279 raw_data = self ._reverse_transforms (transformed_data )
11380
114- print (f"Debug: Raw data shape: { raw_data .shape } " )
115- if hasattr (raw_data , 'dtype' ) and raw_data .dtype .names :
116- if 'inline_number' in raw_data .dtype .names :
117- print (f"Debug: First raw inline_number: { raw_data ['inline_number' ][0 ]} " )
118- print (f"Debug: First raw inline_number (hex): { raw_data ['inline_number' ][0 ]:08x} " )
119-
12081 return raw_data , transformed_data
12182
12283 def _reverse_transforms (self , transformed_data : NDArray ) -> NDArray :
@@ -131,16 +92,10 @@ def _reverse_transforms(self, transformed_data: NDArray) -> NDArray:
13192 # Start with the transformed data
13293 raw_data = transformed_data .copy () if hasattr (transformed_data , 'copy' ) else transformed_data
13394
134- print (f"Debug: Starting reversal with { len (self .transform_pipeline .transforms )} transforms" )
13595
13696 # Apply transforms in reverse order with reversed operations
13797 for i , transform in enumerate (reversed (self .transform_pipeline .transforms )):
138- print (f"Debug: Reversing transform { len (self .transform_pipeline .transforms )- 1 - i } : { type (transform ).__name__ } " )
139- if 'inline_number' in raw_data .dtype .names :
140- print (f"Debug: Before reversal - inline_number: { raw_data ['inline_number' ][0 ]:08x} " )
14198 raw_data = self ._reverse_single_transform (raw_data , transform )
142- if 'inline_number' in raw_data .dtype .names :
143- print (f"Debug: After reversal - inline_number: { raw_data ['inline_number' ][0 ]:08x} " )
14499
145100 return raw_data
146101
@@ -161,37 +116,30 @@ def _reverse_single_transform(self, data: NDArray, transform: Transform) -> NDAr
161116 if isinstance (transform , ByteSwapTransform ):
162117 # For byte swap, we need to reverse the endianness conversion
163118 # If the transform was converting to little-endian, we need to convert back to big-endian
164- print (f"Debug: Reversing byte swap (target was: { transform .target_order } )" )
165119
166120 # Get current data endianness
167121 current_endianness = get_endianness (data )
168- print (f"Debug: Current data endianness: { current_endianness } " )
169122
170123 # If transform was converting TO little-endian, we need to convert TO big-endian
171124 if transform .target_order == Endianness .LITTLE :
172125 reverse_target = Endianness .BIG
173126 else :
174127 reverse_target = Endianness .LITTLE
175128
176- print (f"Debug: Reversing to target: { reverse_target } " )
177129 reverse_transform = ByteSwapTransform (reverse_target )
178130 result = reverse_transform .apply (data )
179131
180- if 'inline_number' in data .dtype .names :
181- print (f"Debug: Byte swap reversal - before: { data ['inline_number' ][0 ]:08x} , after: { result ['inline_number' ][0 ]:08x} " )
182132 return result
183133
184134 elif isinstance (transform , IbmFloatTransform ):
185135 # Reverse IBM float conversion by swapping direction
186136 reverse_direction = "to_ibm" if transform .direction == "to_ieee" else "to_ieee"
187- print (f"Debug: Applying IBM float reversal (direction: { transform .direction } -> { reverse_direction } )" )
188137 reverse_transform = IbmFloatTransform (reverse_direction , transform .keys )
189138 return reverse_transform .apply (data )
190139
191140 else :
192141 # For unknown transforms, return data unchanged
193142 # This maintains compatibility if new transforms are added
194- print (f"Warning: Unknown transform type { type (transform ).__name__ } , cannot reverse" )
195143 return data
196144
197145
0 commit comments