[update] update cali table for 9 hstia

python/biopro/device/device_cc2650.py

@@ -662,7 +662,7 @@ class CC2650Device(Device):
         elif device_type == 'EISZeroOne':
             i = 1
             request_times = 0
-            while i < 13:
+            while i <= 24:
                 try:
                     # send
                     code = self._encode_instruction(DeviceInstruction.TYP_CIS, DeviceInstruction.CIS_CALI, i)

python/biopro/devlib/data.py

@@ -1413,122 +1413,198 @@ class EISZeroOneDataDecoder(RecDataDecoder):
             
     @staticmethod
     def _decode_cali_coeff(cali_coeff: bytes) -> Optional[List[Tuple[int, int]]]:
-        if cali_coeff != b'':
-            cis_data_len = 20
-            cali_table = []
-            hsrtia_a = []
-            hsrtia_b = []
-            rolloff = []
-            phase_coeff = []
-            phase_offset = []
-
-            phase_coeff = numpy.zeros([4, 4], dtype = int)
-            phase_offset = numpy.zeros([4, 4], dtype = int)  
-
-            ########################################
-            # phase_coeff
+            #####################################################
+            # phase_coeff/phase_offset/hsrtia_a/hsrtia_b/rolloff
             # [[gain0, g1, g2, g3] ----->最高頻
             #  [gain0, g1, g2, g3] ----->中頻
             #  [gain0, g1, g2, g3] ----->低頻
             #  [gain0, g1, g2, g3] ----->最低頻
             # ]
-            #######################################
+            #####################################################
+        print('cali_coeff=', cali_coeff)
+        if cali_coeff != b'':
+            cali_table = []
+            hsrtia_a = numpy.zeros([4, 8], dtype = int) #hsrtia_a[freq][gain]
+            hsrtia_b = numpy.zeros([4, 8], dtype = numpy.int64) #hsrtia_b[freq][gain]
+            rolloff = numpy.zeros([4, 8], dtype = int) #rolloff[freq][gain]
+            phase_coeff = numpy.zeros([4, 8], dtype = int) #phase_coeff[freq][gain]
+            phase_offset = numpy.zeros([4, 8], dtype = int) #phase_offset[freq][gain]
+            cis_data_len = 20
 
-            #Lv[0] 160k
+            #gain=0
             cis_cali_packet = 1
             index = (cis_cali_packet - 1) * cis_data_len
-            hsrtia_a.append(struct.unpack('>i', cali_coeff[index+1:index+5])[0])
-            hsrtia_b.append(struct.unpack('>q', cali_coeff[index+5:index+13])[0])
-            rolloff.append(struct.unpack('>i', cali_coeff[index+13:index+17])[0])
+            hsrtia_a[0][0] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][0] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][0]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
             cis_cali_packet = 2
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 0
-            phase_coeff[0][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[0][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[1][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[1][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
-            
+            phase_coeff[0][0]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][0] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][0]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][0] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
             cis_cali_packet = 3
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 0
-            phase_coeff[2][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[2][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[3][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[3][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
-            
-            #Lv[1] 20k
+            phase_coeff[2][0]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][0] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][0]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][0] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            #gain=1
             cis_cali_packet = 4
             index = (cis_cali_packet - 1) * cis_data_len
-            hsrtia_a.append(struct.unpack('>i', cali_coeff[index+1:index+5])[0])
-            hsrtia_b.append(struct.unpack('>q', cali_coeff[index+5:index+13])[0])
-            rolloff.append(struct.unpack('>i', cali_coeff[index+13:index+17])[0])
+            hsrtia_a[0][1] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][1] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][1]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
             cis_cali_packet = 5
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 1
-            phase_coeff[0][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[0][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[1][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[1][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
-            
+            phase_coeff[0][1]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][1] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][1]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][1] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
             cis_cali_packet = 6
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 1
-            phase_coeff[2][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[2][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[3][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[3][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+            phase_coeff[2][1]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][1] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][1]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][1] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
-            #Lv[2] 5k
+            #gain=2
             cis_cali_packet = 7
             index = (cis_cali_packet - 1) * cis_data_len
-            hsrtia_a.append(struct.unpack('>i', cali_coeff[index+1:index+5])[0])
-            hsrtia_b.append(struct.unpack('>q', cali_coeff[index+5:index+13])[0])
-            rolloff.append(struct.unpack('>i', cali_coeff[index+13:index+17])[0])
+            hsrtia_a[0][2] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][2] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][2]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
             cis_cali_packet = 8
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 2
-            phase_coeff[0][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[0][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[1][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[1][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
-            
+            phase_coeff[0][2]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][2] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][2]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][2] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
             cis_cali_packet = 9
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 2
-            phase_coeff[2][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[2][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[3][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[3][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+            phase_coeff[2][2]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][2] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][2]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][2] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
-            #Lv[3] 200R
+            #gain=3
             cis_cali_packet = 10
             index = (cis_cali_packet - 1) * cis_data_len
-            hsrtia_a.append(struct.unpack('>i', cali_coeff[index+1:index+5])[0])
-            hsrtia_b.append(struct.unpack('>q', cali_coeff[index+5:index+13])[0])
-            rolloff.append(struct.unpack('>i', cali_coeff[index+13:index+17])[0])
+            hsrtia_a[0][3] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][3] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][3]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
             cis_cali_packet = 11
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 3
-            phase_coeff[0][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[0][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[1][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[1][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
-            
+            phase_coeff[0][3]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][3] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][3]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][3] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
             cis_cali_packet = 12
             index = (cis_cali_packet - 1) * cis_data_len
-            g = 3
-            phase_coeff[2][g] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
-            phase_offset[2][g] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
-            phase_coeff[3][g] = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
-            phase_offset[3][g] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+            phase_coeff[2][3]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][3] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][3]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][3] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
 
-            print('hsrtia_a', hsrtia_a)
-            print('hsrtia_b', hsrtia_b)
-            print('rolloff', rolloff)
+            #gain=4
+            cis_cali_packet = 13
+            index = (cis_cali_packet - 1) * cis_data_len
+            hsrtia_a[0][4] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][4] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][4]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 14
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[0][4]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][4] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][4]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][4] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 15
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[2][4]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][4] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][4]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][4] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            #gain=5
+            cis_cali_packet = 16
+            index = (cis_cali_packet - 1) * cis_data_len
+            hsrtia_a[0][5] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][5] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][5]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 17
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[0][5]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][5] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][5]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][5] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 18
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[2][5]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][5] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][5]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][5] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            #gain=6
+            cis_cali_packet = 19
+            index = (cis_cali_packet - 1) * cis_data_len
+            hsrtia_a[0][6] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][6] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][6]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 20
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[0][6]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][6] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][6]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][6] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 21
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[2][6]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][6] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][6]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][6] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            #gain=7
+            cis_cali_packet = 22
+            index = (cis_cali_packet - 1) * cis_data_len
+            hsrtia_a[0][7] = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            hsrtia_b[0][7] = struct.unpack('>q', cali_coeff[index+5:index+13])[0]
+            rolloff[0][7]  = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 23
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[0][7]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[0][7] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[1][7]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[1][7] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            cis_cali_packet = 24
+            index = (cis_cali_packet - 1) * cis_data_len
+            phase_coeff[2][7]  = struct.unpack('>i', cali_coeff[index+1:index+5])[0]
+            phase_offset[2][7] = struct.unpack('>i', cali_coeff[index+5:index+9])[0]
+            phase_coeff[3][7]  = struct.unpack('>i', cali_coeff[index+9:index+13])[0]
+            phase_offset[3][7] = struct.unpack('>i', cali_coeff[index+13:index+17])[0]
+
+            print('hsrtia_a')
+            print(hsrtia_a)
+            print('hsrtia_b')
+            print(hsrtia_b)
+            print('rolloff')
+            print(rolloff)
             print('phase_coeff')
             print(phase_coeff)
             print('phase_offset')
@@ -1558,16 +1634,16 @@ class EISZeroOneDataDecoder(RecDataDecoder):
         if len(data) < 18:
             return None
 
-        ch1 = struct.unpack('>i', data[1+3:5+3])[0]  # unit: 1/1000 nA
-        ch2 = struct.unpack('>i', data[5+3:9+3])[0]  # unit: mV
-        ch3 = struct.unpack('>i', data[9+3:13+3])[0]  # unit: kOm
+        ch1 = struct.unpack('>i', data[1+3:5+3])[0]
+        ch2 = struct.unpack('>i', data[5+3:9+3])[0]
+        ch3 = struct.unpack('>i', data[9+3:13+3])[0]
         time_stamp: float = struct.unpack('<I', data[13+3:17+3])[0]  # unit: ms
         cycle_number = struct.unpack('>H', data[17+3:19+3])[0]
         d19 = data[19+3]
         gain = data[20+3]
-        finishMode = (d19 & 0x80) >> 7  
-        ch4 = struct.unpack('<i', data[21+3:25+3])[0] # Amp[uV]
-        
+        finishMode = (d19 & 0x80) >> 7
+        ch4 = struct.unpack('<i', data[21+3:25+3])[0]
+
         notify_one = struct.unpack('<i', data[25+3:29+3])[0]
         notify_two = struct.unpack('<i', data[29+3:33+3])[0]
         notify_three = struct.unpack('<i', data[33+3:37+3])[0]
@@ -1589,61 +1665,59 @@ class EISZeroOneDataDecoder(RecDataDecoder):
                 return None
             else:
                 if self.cali_coeff is not None and (self._mode == 0 or self._mode == 5):
-                    hsrtia_a = []
-                    hsrtia_b = []
-                    rolloff = []
                     phase_coeff, phase_offset, hsrtia_a, hsrtia_b, rolloff = self.cali_coeff[0]
 
                 if (self._mode == 0 or self._mode == 5):
-                    img = ch1
-                    real = ch2
-                    freq = ch3
+                    
+                    img = ch1 #img[ohm]
+                    real = ch2 #real[ohm]
+                    freq = ch3 #freq[10mHz]
                     fre_idx = 0
-                    voltage_amp = round(ch4 / 1000) # Amp[mV]
-                    rolloff_cali = rolloff[gain]
+                    voltage_amp = round(ch4 / 1000) #ch4=Amp[uV] #voltage_amp[mV]
+                    rolloff_cali = rolloff[0][gain]
 
                     voltage_mag = math.sqrt(img ** 2 + real ** 2) * (1 + freq ** 2 /  rolloff_cali ** 2 / 1e4)
-                    current = (voltage_mag ** 2 * hsrtia_a[gain] + voltage_mag * hsrtia_b[gain]) / 1e8 #[nA]
+                    current = (voltage_mag ** 2 * hsrtia_a[0][gain] + voltage_mag * hsrtia_b[0][gain]) / 1e8 #current[nA]
                     
                     if (current != 0):
-                        # impedance[mOhm] = voltage_amp[mv] * 1000000 / 1.414213 / current[nA]       #RMS=amp*SQRT(2), SQRT(2)=1.414213
+                        # impedance[mOhm] = voltage_amp[mV] * 1000000 / 1.414213 / current[nA]       #RMS=amp*SQRT(2), SQRT(2)=1.414213
                         impedance = voltage_amp * 707106.78 / current
                     else:
                         impedance = 0
 
                     raw_phase = math.atan2(img , real) * 180 / math.pi
 
-                    if (freq >= 1000000): # 10000 Hz
+                    if (freq >= 1000000): #10000Hz
                         fre_idx = 0
-                    elif (freq >= 10000): # 100 Hz
+                    elif (freq >= 10000): #100Hz
                         fre_idx = 1
-                    elif (freq >= 1000): # 10 Hz
+                    elif (freq >= 1000): #10Hz
                         fre_idx = 2
-                    elif (freq >= 1): # 0.01 Hz
+                    elif (freq >= 1): #0.01Hz
                         fre_idx = 3
                     
                     ideal_raw_phase = phase_coeff[gain][fre_idx] /1e10 * freq + phase_offset[gain][fre_idx] / 1e6
                     phase = raw_phase - ideal_raw_phase
                     phase = phase % 180 if phase % 180<=90 else phase % 180-180
 
-                    imag_after_cal = round(impedance * math.sin(phase * math.pi / 180))
-                    real_after_cal = round(impedance * math.cos(phase * math.pi / 180))
+                    imag_after_cal = impedance * math.sin(phase * math.pi / 180)
+                    real_after_cal = impedance * math.cos(phase * math.pi / 180)
 
                 if self._show_data:
                     if (self._mode == 0 or self._mode == 5):
                         print('|', '{:10}'.format(time_stamp),
                               '|', '{:5}'.format(delta),
-                              '|', '{:5}'.format(ch1),                            #raw_img
-                              '|', '{:5}'.format(ch2),                            #raw_real
-                              '|', '{:8}'.format(ch3 * 10), '[mHz]',              #Frequency [mHz]
-                              '|', '{:5}'.format(cycle_number),                   #cycle
-                              '|', '{:5}'.format(round(imag_after_cal)), '[Ohm]', #Z_imag [Ohm]
-                              '|', '{:5}'.format(round(real_after_cal)), '[Ohm]', #Z_real [Ohm]
-                              '|', '{:5}'.format(round(impedance)), '[Ohm]',      #Impedance [Ohm]
-                              '|', '{:5}'.format(round(phase*1000)), '[mdegree]', #Phase [millidegree]
-                              '|', '{:5}'.format(round(current)), '[nA]',         #Current [nA]
-                              '|', '{:1}'.format(gain),                           #gain
-                              '|', '{:1}'.format(finishMode),                     #finishMode
+                              '|', '{:5}'.format(img),
+                              '|', '{:5}'.format(real),
+                              '|', '{:9}'.format(freq*10), '[mHz]',
+                              '|', '{:5}'.format(cycle_number),
+                              '|', '{:5}'.format(round(imag_after_cal)), '[Ohm]', #Z_imag[Ohm]
+                              '|', '{:5}'.format(round(real_after_cal)), '[Ohm]', #Z_real[Ohm]
+                              '|', '{:5}'.format(round(impedance)), '[mOhm]',
+                              '|', '{:5}'.format(round(phase*1000)), '[mdegree]',
+                              '|', '{:10}'.format(round(current)), '[nA]',
+                              '|', '{:1}'.format(gain),
+                              '|', '{:1}'.format(finishMode),
                               '@', str(self.device), '|', flush = True)
 
                         print('|', '{:10}'.format(time_stamp),
@@ -1651,7 +1725,7 @@ class EISZeroOneDataDecoder(RecDataDecoder):
                               '|', '{:5}'.format(notify_one),
                               '|', '{:5}'.format(notify_two),
                               '|', '{:5}'.format(notify_three),
-                              '|', '{:5}'.format(voltage_amp),                    #amp[mV]
+                              '|', '{:5}'.format(voltage_amp), '[mV]',
                               '|', flush = True)
                         pass
                     else:
@@ -1673,28 +1747,30 @@ class EISZeroOneDataDecoder(RecDataDecoder):
             self._mode_stop = 0
 
         ret = RecordingData(self.device, int(time_stamp * 1000 / 2), 0)
-        if self._mode == 0 or self._mode == 5:
-            ret.append_data(0, ch1)                     #raw_img
-            ret.append_data(1, ch2)                     #raw_real
-            ret.append_data(2, ch3 * 10)                #Frequency [mHz]
-            ret.append_data(3, cycle_number)            #cycle
-            ret.append_data(4, imag_after_cal)          #Z_imag [Ohm]
-            ret.append_data(5, real_after_cal)          #Z_real [Ohm]
-            ret.append_data(6, round(impedance))        #Impedance [Ohm]
-            ret.append_data(7, round(phase*1000))       #Phase [millidegree]
-            ret.append_data(8, round(current))          #Current [nA]
-            ret.append_data(9, gain)                    #gain
+        if self._mode == 0 or self._mode == 5:          #EIS/CF Mode
+            ret.append_data(0, img)
+            ret.append_data(1, real)
+            ret.append_data(2, freq * 10)               #[mHz]
+            ret.append_data(3, cycle_number)
+            ret.append_data(4, round(imag_after_cal))   #Z_imag [Ohm]
+            ret.append_data(5, round(real_after_cal))   #Z_real [Ohm]
+            ret.append_data(6, round(impedance))        #[mOhm]
+            ret.append_data(7, round(phase*1000))       #[millidegree]
+            ret.append_data(8, round(current))          #[nA]
+            ret.append_data(9, gain)
             
             #debug data
             ret.append_data(10, notify_one)
             ret.append_data(11, notify_two)
             ret.append_data(12, notify_three)
-            ret.append_data(13, voltage_amp)            #amp[mV]
+            ret.append_data(13, voltage_amp)            #mV
+
+
 
         else:                                           #CV Mode
-            ret.append_data(0, ch1)                     #Iin     [nA]
-            ret.append_data(1, ch2)                     #Vset    [nV]
-            ret.append_data(2, ch3)                     #Vout    [nV]
+            ret.append_data(0, ch1)                     #Iin     [nA]?
+            ret.append_data(1, ch2)                     #Vset    [nV]?
+            ret.append_data(2, ch3)                     #Vout    [nV]?
             ret.append_data(3, cycle_number)
 
         if cycle_number != self._cycle_number: