我在上一篇文章中描述了wis二進位制檔案的結構（wis是石油測井軟體Forward中的一種二進位制檔案格式，很多軟體需要轉換為文字才可使用），然後用純C編寫了轉換程式，利用Python編寫了介面呼叫這個程式，在這裡我利用C++重新編寫了這個程式，當然和C編寫的區別不大，主要的區別在記憶體分配、讀取二進位制檔案和寫入文字檔案的函式上，也是為了說明C與C++檔案操作上的一些區別，有興趣的朋友可以和上一遍文章對比一下，技術有限，僅供參考！

以下是C++原始碼：

/* * Copyright (c) 2011, Sun Yunqiang * All rights reserved. * Software: wis2txt_win.cpp * Version: 1.0 * Author: Sun Yunqiang * English Name: Alex Sun * Release Date: 18/03/2011 */ /** ***************************************************************** */ /** This Program is for converting wis format file to txt format file */ /** Command: wis2txt_win.exe args1 args2 */ /** Description: */ /** wis2txt_win.exe: the executable file as command */ /** args1: the first argument that is the wis file name */ /** args2: the second argument that is the txt file name */ /** ***************************************************************** */ #include <windows.h> #include <iostream> #include <iomanip> #include <fstream> #include <string> #include <cstdlib> #include <cctype> #include <ctime> #include <cmath> #define ROW 65536 #define COL 512 /* **************************************** */ /* Below blocks are the structs of wis file */ /* **************************************** */ typedef struct tagWIS_HEAD { WORD MachineType; // 0-PC 1-SUN 2-IBM 3-HP WORD MaxObjectNumber; WORD ObjectNumber; WORD BlockLen; DWORD EntryOffset; DWORD DataOffset; DWORD FileSize; time_t TimeCreate; char Reserved[32]; }WIS_HEAD; typedef struct tagWIS_OBJECT_ENTRY { char Name[16]; long Status; short Attribute; // 1-通道物件 2-表物件 3-流物件 short SubAttribute; // 1-曲線物件 2-波形物件 3-地層測試物件 4-時深時物件 DWORD Position; // 物件資料體從檔案開始處的偏移量 DWORD BlockNum; time_t TimeCreate; time_t TimeWrite; char Reserved[32]; }WIS_OBJECT_ENTRY; typedef struct tagWIS_CHANNEL_DIMENSION { char Name[8]; char Unit[8]; char AliasName[16]; float StartVal; float Delta; DWORD Samples; DWORD MaxSamples; DWORD Size; WORD RepCode; WORD Reserved; }WIS_CHANNEL_DIMENSION; typedef struct tagWIS_CHANNEL { char Unit[8]; char AliasName[16]; char AliasUnit[16]; WORD RepCode; WORD CodeLen; float MinVal; float MaxVal; WORD Reserved; WORD NumOfDimension; WIS_CHANNEL_DIMENSION DimInfo[4]; }WIS_CHANNEL; typedef struct tagWIS_STREAM { DWORD Length; DWORD Offset; }WIS_STREAM; typedef struct tagWIS_TABLE_FIELD { char Name[32]; WORD RepCode; WORD Length; DWORD Reserved; } WIS_TABLE_FIELD; typedef struct tagWIS_TABLE { DWORD RecordCount; DWORD FieldCount; WIS_TABLE_FIELD *pField; }WIS_TABLE; typedef struct tagWIS_TABLE_DEFAULT_FILED { char Name[16]; char Alias[16]; char Unit[8]; char Type[8]; WORD Length; WORD Count; char DefVal[64][12]; } WIS_TABLE_DEFAULT_FIELD; typedef struct tagWIS_DEFAULT_TABLE { char Name[16]; char Alias[16]; char Attrb[8]; DWORD FieldCount; WIS_TABLE_DEFAULT_FIELD *pField; } WIS_DEFAULT_TABLE; /* 曲線的資訊結構體*/ typedef struct tagCURVE_HEAD { char Name[16]; // 曲線的名字 char Unit[8]; // 曲線的單位 float startMD; // 曲線的起始深度 float stopMD; // 曲線的結束深度 int pointNum; // 資料樣點數 int position; // 資料位置 struct tagCURVE_HEAD *next; } CURVE_HEAD; /* 曲線的資料結構體 */ typedef struct tagCURVE_DATA { char Name[16]; // 曲線的名字 char Unit[8]; // 曲線的單位 float Depth[ROW]; // 曲線的深度值 float Value[ROW]; // 曲線的值 struct tagCURVE_DATA *next; } CURVE_DATA; /* ************************************************************* */ /* Below block is the Main program of processing wis format file */ /* ************************************************************* */ int main(int argc, char *argv[]) { std::ifstream wisfile; // 輸入檔案 std::ofstream txtfile; // 輸出檔案 char *infilename = argv[1]; char *outfilename = argv[2]; // char *infilename = "demoa.wis"; // char *outfilename = "demoa.txt"; if(argc != 3 ) { std::cerr<<"ERROR: Command parameters is not defined correctly!"<<std::endl; exit(1); } wisfile.open(infilename, std::ios::binary); if(!wisfile.is_open()) { std::cerr<<"ERROR: Read file terminated!"<<std::endl; exit(1); } txtfile.open(outfilename); if(!txtfile.is_open()) { std::cerr << "ERROR: Write file terminated!" << std::endl; exit(1); } /* Process wis file, Read it and transfer it to text */ /** wis檔案識別符號從檔案偏移0開始，為10個位元組的字元 */ char wisid[10]; wisfile.read(wisid, 10); std::cout<<"檔案識別符號: WIS" << wisid[4] << wisid[5] << wisid[6] << std::endl; txtfile << "WIS" << wisid[4] << wisid[5] << wisid[6] << " to ASCII TEXT by SUNYQ/n"; /** 標頭檔案緊接檔案標識 */ WIS_HEAD *wishead = new WIS_HEAD; wisfile.seekg(10, std::ios::beg); wisfile.read(reinterpret_cast<char *>(wishead), sizeof(WIS_HEAD)); // 物件入口記錄從檔案開始的偏移量: wishead->EntryOffset // 物件資料記錄從檔案開始的偏移量: wishead->DataOffset // 當前記錄的物件總數: wishead->ObjectNumber /** 物件入口，位置由頭結構中EntryOffset引數指定 最多512條曲線 */ wisfile.seekg(wishead->EntryOffset, std::ios::beg); delete wishead; /** curvePosition陣列 儲存曲線資料開始的位置 WIS_OBJECT_ENTRY->Position的值 */ /** curveName陣列 儲存每條曲線的名稱 */ /** curveNum變數 檔案中儲存曲線的數目 */ int curvePosition[512]; int curveName[512][16]; int curveNum = 0; // curveNum is the identifer of the curves number while(1) { WIS_OBJECT_ENTRY *objectEntry; objectEntry = new WIS_OBJECT_ENTRY; wisfile.read(reinterpret_cast<char *>(objectEntry), sizeof(WIS_OBJECT_ENTRY)); if(objectEntry->Attribute == 0) { break; } else if(objectEntry->Attribute == 1) { // 通道型別: objectEntry->Attribute // 曲線名字: objectEntry->Name // 物件資料體絕對偏移量: objectEntry->Position curvePosition[curveNum] = objectEntry->Position; int i=0; for(i=0; i<16; i++) { curveName[curveNum][i] = objectEntry->Name[i]; } curveNum++; } else { continue; } delete objectEntry; } std::cout << infilename << " 檔案的曲線數目: " << curveNum << std::endl; /* ****************************************************************** */ /* Below block read every curve head information and put it in a link */ /* ****************************************************************** */ /** 建立連結串列 data_node 開始讀取資料 */ /** num變數 每條曲線的行數 */ CURVE_HEAD *first_head = NULL; float delta = 0.0; CURVE_HEAD *data_head = NULL; int num = 0; for(num = 0; num < curveNum ; num++) // num相當於第幾條曲線 { data_head = new CURVE_HEAD; data_head->position = curvePosition[num]; wisfile.seekg(curvePosition[num], std::ios::beg); WIS_CHANNEL *channel; channel = new WIS_CHANNEL; wisfile.read(reinterpret_cast<char *>(channel), sizeof(WIS_CHANNEL)); data_head->pointNum = channel->DimInfo[0].MaxSamples; // 深度取樣點的個數 data_head->startMD = channel->DimInfo[0].StartVal; // 起始深度 delta = channel->DimInfo[0].Delta; // 間隔點值 data_head->stopMD = channel->DimInfo[0].StartVal + delta * channel->DimInfo[0].MaxSamples; // 結束深度 int i = 0; for(i=0; i<16; i++) { data_head->Name[i] = curveName[num][i]; } for(i=0; i<8; i++) { data_head->Unit[i] = channel->Unit[i]; } std::cout << "No." << num+1 << "/t"; // 曲線編號 std::cout << data_head->Name << "/t"; // 曲線名 std::cout << data_head->Unit << "/t"; // 曲線單位 std::cout << data_head->pointNum << "/t"; // 取樣點數 std::cout << delta << "/t"; // 取樣點間隔 std::cout << data_head->startMD << "/t"; // 起始深度 std::cout << data_head->stopMD << std::endl; // 結束深度 data_head->next = first_head; first_head = data_head; } std::cout << "Waiting..." << std::endl; /* *********************************************************************** */ /* Below block is for deciding the min and max Depth, and the max line num */ /* *********************************************************************** */ // 從這裡開始判斷所有曲線中最小深度和最大深度值 float minDepth = 0.0; float maxDepth = 0.0; minDepth = first_head->startMD; maxDepth = first_head->stopMD; for(data_head = first_head; data_head!= NULL; data_head=data_head->next) { if(minDepth>data_head->startMD) { minDepth = data_head->startMD; // 確定最小深度值 } if(maxDepth < data_head->stopMD) { maxDepth = data_head->stopMD; // 確定最大深度值 } } // 利用最小最大深度值和間隔點值，確定輸出的取樣點數 float depth[ROW]; int lineNum = 0; lineNum = (maxDepth - minDepth) / delta; // 取樣點個數：lineNum int row=0; for(row=0; row<lineNum; row++) { depth[row] = minDepth + delta * row; } /* ******************************************************** */ /* Below block put every point value into a link as a curve */ /* ******************************************************** */ // 從這裡開始將資料值放入曲線各自對應的連結串列裡 CURVE_DATA *first_node = NULL; CURVE_DATA *data_node = NULL; data_head = first_head; for(num = 0; num < curveNum ; num++) { data_node = new CURVE_DATA; int i = 0; for(i=0; i<16; i++) { data_node->Name[i] = data_head->Name[i]; // 曲線名稱 } for(i=0; i<8; i++) { data_node->Unit[i] = data_head->Unit[i]; // 曲線單位 } int startLine = 0; for(row=0; row<lineNum; row++) { if( abs(depth[row]-data_head->startMD)<= delta && depth[row] >= data_head->startMD ) { startLine = row; break; } } // 起始點: startLine wisfile.seekg(data_head->position+2*sizeof(WIS_CHANNEL)+464, std::ios::beg); // 測試樣點數: data_head->pointNum for(row=0; row<lineNum; row++) { float data[1]; data[0] = -9999.000; if(row>=startLine && (row-startLine) < data_head->pointNum) { wisfile.read(reinterpret_cast<char *>(&data), sizeof(float)); data_node->Depth[row] = data_head->startMD + delta * (row-startLine); data_node->Value[row] = data[0]; } else { data_node->Depth[row] = minDepth + delta * row; data_node->Value[row] = -9999.000; } // 讀取值：data_node->Value[row] } data_node->next = first_node; first_node = data_node; if(data_head->next == NULL) { break; } data_head = data_head->next; } /* **************************************************** */ /* Below block print all point value for every curve */ /* The important thing is how to decide the right depth */ /* **************************************************** */ // 從這裡開始寫入文字 // Print Curve Name txtfile << std::left << std::setw(12) << "DEPTH"; for(data_node=first_node; data_node!=NULL; data_node=data_node->next) { txtfile << std::left << std::setw(12) << data_node->Name; // txtfile << " "; } txtfile << "/n"; // Print Curve Unit txtfile << std::left << std::setw(12) << "M"; for(data_node=first_node; data_node!=NULL; data_node=data_node->next) { txtfile << std::left << std::setw(12) << data_node->Unit; } txtfile << "/n"; // Print curve value for(row=0; row<lineNum; row++) { if(depth[row] > maxDepth) { break; } txtfile << std::left << std::setw(12) << std::setiosflags(std::ios::fixed) << std::setprecision(3) << depth[row]; // txtfile << " "; for(data_node=first_node; data_node!=NULL; data_node=data_node->next) { txtfile << std::left << std::setw(12) << std::setiosflags(std::ios::fixed) << std::setprecision(3) << data_node->Value[row]; } txtfile << "/n"; } // */ // Process Over wisfile.close(); txtfile.close(); std::cout << infilename << " Complete!" << std::endl; // system("PAUSE"); // Program Over return 0; }