使用 Python 解析 WIS 格式测井曲线文件

    def _parseWis(self, bfirstEnter, curvename = '', enc = 'cp936'):
        ''' 解析wis文件 返回曲线信息（基础信息+数据体）
            只读取需要的曲线数据
         '''
        try:
            fname = os.path.abspath(self.filename)
            self.wellname = os.path.basename(fname).split(os.path.extsep, 1)[0]
            with open(fname, 'rb') as f:
                f.seek(14)
                tmp = f.read(8)
                objcount, blocksize, objoffset = struct.unpack('<hhi', tmp)
                for i in range(objcount):
                    ioffset = objoffset + i * 72
                    f.seek(ioffset)
                    name = f.read(16).decode(enc).split('\x00', 1)[0]
                    tmp = f.read(12)
                    ostate, otype, _, hoff = struct.unpack('<ihhi', tmp)  # hoff为对象数据体的偏移
                    if ostate == 0 and otype == 1:  # ostate == 0 为正常（非抛弃或删除）otype == 1为通道对象
                        f.seek(hoff)
                        chnl = f.read(52)
                        unit, aliasName, aliasUnit, _, _, mixVal, maxVal = struct.unpack('<8s16s16sHHff', chnl)
                        unit = unit.decode(enc).strip('\x00')
                        f.seek(hoff + 54)
                        tmp = f.read(2)
                        dim, = struct.unpack('<h', tmp)  # attention
                        if dim == 1:  # 维数
                            f.read(32)
                            tmp = f.read(22)
                            fstart, fstep, ncount, maxc, vsize, vtype = struct.unpack('<ffiiih',
                                                                                      tmp)