// SaveInRes.cpp: implementation of the CSaveInRes class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "SaveInRes.h"
#include "Constant.h"
#include "CtrlProtocolDef.h"
#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

#define AR_CROSS_HOLE_EAMOUNT 48
#define AR_CROSS_HOLE_SURFACE_ELECTRODE_NUM 7
//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
extern HHOOK   hHook;   
extern LRESULT  __stdcall CBTHookProc(long nCode, WPARAM wParam, LPARAM lParam);
extern int g_iUILanguage;
CSaveInRes::CSaveInRes()
{
	m_resFile = NULL;
	m_nFirstPole = 0;
	m_nLandSign = 0;
	m_MediumType = 0;
/*	m_rect = NULL;*/
	m_szDistancePole = _T("");
	m_szPoleStart = _T("1");
	m_methodType = 0;
	m_fCrossHoleParam = 0;
	m_fPoleDistance = 0;
	m_i3DRollDiecetion = 0;
	m_fMinAMDistance = 0;
	m_fAPosX = 0;
	m_fAPosY = 0;
	m_iEAmount = 0;
	m_fInitiDepth = 0;
	m_byLineDirection = 0;
}

CSaveInRes::~CSaveInRes()
{


}

BOOL CSaveInRes::CreateResFile(CString szFilePath)
{
	m_resFile = fopen(szFilePath, "w");
	
	if (m_resFile == NULL)
	{
//		hHook = SetWindowsHookEx(WH_CBT,(HOOKPROC)CBTHookProc,AfxGetInstanceHandle(),NULL);  
//		AfxMessageBox("打开文件失败");
		return FALSE;
	}

	return TRUE;
}

BOOL CSaveInRes::AddDateToFile(CStringArray *szDateArray)
{
	CString szDate;
	szDate.Empty();
	
	CStringArray szDealDateArray;
	szDealDateArray.RemoveAll();
	
	switch(m_MediumType)
	{
	case 1://温纳
		Wenner(szDateArray,&szDealDateArray);
		break;
	case 2://偶极
		dipole(szDateArray,&szDealDateArray);
		break;
	case 3:
		Wenner_Y_2D(szDateArray,&szDealDateArray);
		break;
	case 4:
		Pole_dipole_2d(szDateArray,&szDealDateArray);
		break;
	case 5:
		Reverse_Pole_dipole_2d(szDateArray,&szDealDateArray);
		break;
	case 6:
		A_M_ARRAY_2D(szDateArray,&szDealDateArray);
		break;
	case 7:
		AMN_ARRAY_2D(szDateArray,&szDealDateArray);
		break;
	case 9:
		Dipole_Dipole_2D(szDateArray,&szDealDateArray);
		break;
	case 11://斯伦贝谢
		Schlumberger(szDateArray,&szDealDateArray);
		break;
	case 16://跨孔
		Cross_Hole_2D(szDateArray,&szDealDateArray);
		break;
	case 17:
		CE_Schlumberger(szDateArray,&szDealDateArray);
		break;
	case 22:
		Pole_Pole_3D(szDateArray,&szDealDateArray);
		break;
	case 23:
		Pole_Dipole_3D(szDateArray,&szDealDateArray);
		break;
	case 24:
		Dipole_Dipole_3D(szDateArray,&szDealDateArray);
		break;
	case 25:
		Schlumberger3D(szDateArray,&szDealDateArray);
		break;
	case 26:
		WennerAlfa3D(szDateArray,&szDealDateArray);
		break;
	case 27:
		WennerBeta3D(szDateArray,&szDealDateArray);
		break;
	case 28:
		Gradient_3D(szDateArray,&szDealDateArray);
		break;
	case 29:
		Mid_Gardient_Scan_3D(szDateArray,&szDealDateArray);
		break;
	case 34:
		CE_Res_Marine(szDateArray, &szDealDateArray);
		break;
	case 50://to be do, louis		
	case 39://温-斯装置
		Wenner_Schlumberger_2D(szDateArray,&szDealDateArray);
		break;
	//由于边缘梯度的排列和用户自定义一样，在这里就直接沿用了用户自定义的排列
	case 40:
	case 41:
	case 43:
	case 44:
	case 47:
		Custom_2D(szDateArray,&szDealDateArray);
		break;
	case 45:
	case 46:
		CrossHole_2D(szDateArray, &szDealDateArray);
		break;
	default:
		return FALSE;
	}

	szDate.Empty();
	
	for (int nIndex = 0; nIndex < szDealDateArray.GetSize(); nIndex++)
	{
		szDate += szDealDateArray.GetAt(nIndex);
		if(nIndex != szDealDateArray.GetSize() - 1)
			szDate += " ";
	}

	szDate += "\n";
	fwrite(szDate.GetBuffer(szDate.GetLength()), 1, szDate.GetLength(), m_resFile);
	szDate.Empty();
	fflush(m_resFile);

	return TRUE;	
}

CString CSaveInRes::TransformMediumType(CString szMediumType)
{
	CString szResult;
	int nType;
	
	szResult.Empty();
	szResult = szMediumType;
	
	nType = atoi(szResult);
	szResult.Empty();
	
	switch (nType)//nType是我们的装置类型
	{
	case 1:
		szResult = "1";
		break;
		
	case 2:
		szResult = "4";
		break;
		
	case 3:
		szResult = "5";
		break;	
		
	case 4:
		szResult = "6";
		break;	
	
	case 5:
		szResult = "6";
		break;
	case 6:
		szResult = "2";
		break;
	case 7:
		szResult = "6";
		break;
	case 8:
		szResult = szMediumType;
		break;
	case 9:
		szResult = "3";
		break;
	case 10:
		szResult = szMediumType;
		break;
	case 11:
		szResult = "7";
		break;
	case 12:
		szResult = szMediumType;
		break;
	case 13:
		szResult = szMediumType;
		break;
		
	case 14:
		szResult = szMediumType;
		break;
		
	case 15:
		szResult = szMediumType;
		break;
		
	case 16:
	case 45://AR_CROSS_HOLE_GEOMATIVE
	case 46://AR_CROSS_HOLE_GEOMATIVE_AM
		szResult = "12";
		break;
	case 18:
		break;
	case 22:
		szResult = "2";
		break;
	case 23:
		szResult = "6";
		break;
	case 24:
		szResult = "3";
		break;
	case 25:
		szResult = "7";
		break;
	case 26:
		szResult = "1";
		break;
	case 27:
		szResult = "4";
		break;
	case 28:
		szResult = "11";
		break;
	case 29:
		szResult = "11";
		break;
	case 50://to be do, louis		
	case 39:
		szResult = "7";
		break;
	case 40:
	case 41:
	case 43:
	case 44:
	case 47:
	case 48://AR_LAND_FILL_ONE
	case 49://AR_LAND_FILL_TWO
		szResult = "11";
		break;
	default:
		return "NULL";
	}

	return szResult;
}


BOOL CSaveInRes::CloseFile()
{	
	switch (m_MediumType)
	{
	case 2:
	case 11:
		WriteZero(6);
		break;
	case 25:
	case 26:
	case 27:
		WriteZero(5);
		break;
	default:
		WriteZero(4);
	}

	if (m_resFile)
	{
		fclose(m_resFile);
	}
	return FALSE;
}

void CSaveInRes::WriteZero(int nCount)
{
	for(int nIndex = 0; nIndex < nCount; nIndex++)
	{
		fwrite("0\r\n",1,3,m_resFile);
		fflush(m_resFile);
	}
}

void CSaveInRes::Wenner(CStringArray *szSRCArray, CStringArray *szTARArray)
{
//szSRCArray这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
//Wenner的计算值有三列，计算公式分别：
//第一列:A * 间距
//第二列:(M-A) * 间距
//第三列:R0

	szTARArray->RemoveAll();

//第一列:
	DOUBLE firstRow;
	CString szFirstRow;
	firstRow = (0.5*(atof(szSRCArray->GetAt(2)) + atof(szSRCArray->GetAt(3)))-1) * atof(m_szDistancePole.GetBuffer(0));
	szFirstRow.Format("%.3f",fabs(firstRow));
	szTARArray->Add(szFirstRow);

//第二列:
	DOUBLE secondRow;
	CString szSecondRow;
	secondRow = (atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0)) ) * atof(m_szDistancePole.GetBuffer(0));
	szSecondRow.Format("%.3f",fabs(secondRow));
	szTARArray->Add(szSecondRow);


//第三列
	DOUBLE thirdRow;
	CString szThirdRow = _T("");
	thirdRow = atof(szSRCArray->GetAt(4));
	szThirdRow.Format("%.3f",thirdRow);
	szTARArray->Add(szThirdRow);

}

void CSaveInRes::Schlumberger(CStringArray *szSRCArray, CStringArray *szTARArray)
{
//szSRCArray这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
//Schlumberger的计算值有四列，计算公式分别：
//第一列:MN的中点 * m_szDistancePole
//第二列:(p2 - p1) * m_szDistancePole
//第三列:(p1 - c1) / (p2 - p1)
//第四列:R0

	szTARArray->RemoveAll();
	
	//第一列:
	DOUBLE firstRow;
	DOUBLE fDistancePole;
	CString szFirstRow;

	fDistancePole = atof(m_szDistancePole.GetBuffer(0));
	firstRow = (0.5*(atof(szSRCArray->GetAt(2))+atof(szSRCArray->GetAt(3)))-1) * fDistancePole;
	szFirstRow.Format("%.3f",fabs(firstRow));
	szTARArray->Add(szFirstRow);
	
	//第二列:
	CString szSencondRow;
	DOUBLE fdistance, fSencondRow;
	fdistance = atof(m_szDistancePole);//m_szDistancePole成员变量，保存了电极间距
	fSencondRow = ( atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)) ) * fdistance;
	szSencondRow.Format("%.3f",fabs(fSencondRow));
	szTARArray->Add(szSencondRow);
	
	//第三列
	DOUBLE thirdRow;
	CString szThirdRow;
	//thirdRow = atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0));
	thirdRow = ( atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0)) ) / ( atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)) );
	szThirdRow.Format("%.5f",fabs(thirdRow));
	szTARArray->Add(szThirdRow);
	
	//第四列
	DOUBLE FourthRow;
	CString szFourthRow;
	FourthRow = atof(szSRCArray->GetAt(4));
	szFourthRow.Format("%.3f",FourthRow);
	szTARArray->Add(szFourthRow);
}
//温钠Beta
void CSaveInRes::dipole(CStringArray *szSRCArray, CStringArray *szTARArray)
{
//szSRCArray这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
//Schlumberger的计算值有四列，计算公式分别：
//第一列:C2 * m_szDistancePole，在这里要减掉第一个电极的基准位置（我们从1开始，反演软件从0开始）
//第二列:(C1 - C2) * m_szDistancePole
//第三列:R0
	
	szTARArray->RemoveAll();
	
	//第一列:
	DOUBLE firstRow;
	CString szFirstRow;
	
	firstRow = (0.5*(atof(szSRCArray->GetAt(2)) + atof(szSRCArray->GetAt(0)))-1) * atof(m_szDistancePole.GetBuffer(0));
	szFirstRow.Format("%.3f",fabs(firstRow));
	szTARArray->Add(szFirstRow);
	
	//第二列:
	CString szSencondRow;
	DOUBLE fSencondRow;
	DOUBLE fdistance;
	fdistance = atof(m_szDistancePole);//m_szDistancePole成员变量，保存了电极间距
	fSencondRow = ( atof(szSRCArray->GetAt(0)) - atof(szSRCArray->GetAt(1)) ) * fdistance;
	szSencondRow.Format("%.3f",fabs(fSencondRow));
	szTARArray->Add(szSencondRow);
	
	//第三列
	DOUBLE thirdRow;
	CString szThirdRow;
	thirdRow = atof(szSRCArray->GetAt(4));
	szThirdRow.Format("%.3f",thirdRow);
	szTARArray->Add(szThirdRow);
	
}

BOOL CSaveInRes::WriteLandInform(CStringArray *f_szArray)
{
	CString szNum;
	
	szNum.Format("%d\n", m_nLandSign);
	fwrite(szNum.GetBuffer(0), szNum.GetLength(), 1, m_resFile);

	//写入地型数据标志 和 地型数据点的个数
	szNum.Format("%d\n", f_szArray->GetSize());
	fwrite(szNum.GetBuffer(0), szNum.GetLength(), 1, m_resFile);

	for (int i = 0; i < f_szArray->GetSize(); i++)
	{
		szNum = f_szArray->GetAt(i);
		szNum += "\n";
		fwrite(szNum.GetBuffer(0), szNum.GetLength(), 1, m_resFile);
	}

	szNum.Format("%d\n", m_nFirstPole);

	fwrite(szNum.GetBuffer(0), szNum.GetLength(), 1, m_resFile);


	return TRUE;
}

void CSaveInRes::WennerAlfa3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	/*
	按顺序传入A、B、M、N、R0 5个值
	*/
	int A = 0, B = 1, M = 2, N = 3, R = 4;
//	int nSort[4] = {A, M, N, B};
	int nSort[4] = {A, B, M, N};
	int nId = 0;

	_pos *Pos = NULL;

	CString szTmp = _T("");
	int nTmp = 0;

	Pos = new _pos();

	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 4; i++)
	{
		nId = nSort[i];
		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
		
		pole_to_pos_3d(&tmpRect, atoi(m_szPoleStart.GetBuffer(0)), m_nPoleStepX, m_nPoleStepY,Pos, nTmp,m_byLineDirection);
		
		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}

		szTmp.Empty();
		szTmp.Format("%.2f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.2f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
	}

	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);

	delete Pos;
	Pos = NULL;

}

void CSaveInRes::WennerBeta3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	/*
	按顺序传入A、B、M、N、R0 5个值
	*/

	int A = 0, B = 1, M = 2, N = 3, R = 4;
	int nSort[4] = {A, B, M, N};
	int nId = 0;
	_pos *Pos = NULL;

	CString szTmp = _T("");
	int nTmp = 0;
	int nPoleStart = 0;

	int nDistancePole = 0;

	nPoleStart = atoi(m_szPoleStart.GetBuffer(m_szPoleStart.GetLength()));
	
	Pos = new _pos();

	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 4; i++)
	{
		nId = nSort[i];

		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
	
		pole_to_pos_3d(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY, Pos, nTmp,m_byLineDirection);

		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
	}

	int q = szSRCArray->GetSize();

	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);

	delete Pos;
	Pos = NULL;
}

void CSaveInRes::Schlumberger3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	/*
	按顺序传入A、B、M、N、R0 5个值
	*/
	int A = 0, B = 1, M = 2, N = 3, R = 4;
//	int nSort[4] = {A, M, N, B};
	int nSort[4] = {A, B, M, N};
	int nId = 0;

	_pos *Pos = NULL;

	CString szTmp = _T("");
	int nTmp = 0;
	Pos = new _pos();

	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 4; i++)
	{
		nId = nSort[i];
		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
		pole_to_pos_3d(&tmpRect, atoi(m_szPoleStart.GetBuffer(0)), m_nPoleStepX, m_nPoleStepY, Pos, nTmp, m_byLineDirection);
		
		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}

		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
	}
	
	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);

	delete Pos;
	Pos = NULL;
}

void CSaveInRes::AddFileHeader(int method, int mediumType, CStringArray *szDataArray)
{
	m_MediumType = mediumType;	//保存在装置类型
	m_methodType = method;		//保存测试方法(电阻率Res,激电IP...)
	
	/*if (mediumType >= 17 && mediumType <= 21)		//1D脚本装置类型
	{
		InitCEFileHeader(szDataArray);
	}
	
	if ((mediumType >= 1 && mediumType <= 16) || (mediumType >= 39 && mediumType != 46))		//2D脚本装置类型
	{
		Init2DFileHeader(szDataArray);
	}

	if ((mediumType > 21 && mediumType < 39) || mediumType == 46)							//3D脚本装置类型
	{
		Init3DFileHeader(szDataArray);
	}*/
}

void CSaveInRes::Init3DFileHeader(CStringArray *szArray)
{
	//3D的文件头
	CString szTmp = _T("");
	int x0, y0, x1,y1;
	
	for (int i = 0; i < szArray->GetSize(); i++)
	{
		switch (i)
		{
 		case 1:
 			szTmp = szArray->GetAt(i);
			x0 = atoi(szTmp.GetBuffer(0));
 			m_rect.x0 = x0;
 			continue;
 		case 2:
 			szTmp = szArray->GetAt(i);
 			y0 = atoi(szTmp.GetBuffer(0));
 			m_rect.y0 = y0;			
			continue;
		case 3:
			szTmp = szArray->GetAt(i);
			x1 = atoi(szTmp.GetBuffer(0));
			m_rect.x1 = x1;	
			//如果是中梯，则X轴网格大小取的是AB距离/X轴电极间距
			if (29 == atoi(szArray->GetAt(7)))
			{
				float fXPoleDis = atof(szArray->GetAt(5));
				int iXGridSize = (0 == fXPoleDis) ? (int)m_fABDistance : (int)(m_fABDistance/fXPoleDis);
				szTmp.Format(_T("%d"), iXGridSize+1);
			}
			else
			szTmp.Format("%d", x1 - x0 + 1);
			//szTmp.Format("%d", abs(x1 - x0) + 2);
			break;
		case 4:
			szTmp = szArray->GetAt(i);
			y1 = atoi(szTmp.GetBuffer(0));
			m_rect.y1 = y1;
			szTmp.Format("%d", y1 - y0 + 1);
		//	szTmp.Format("%d", abs(y1 - y0) + 2);
			break;
		case 5:
			m_szX_DistancePole = szArray->GetAt(i);
			szTmp = m_szX_DistancePole;
			break;
		case 6:
			m_szY_DistancePole = szArray->GetAt(i);
			szTmp = m_szY_DistancePole;
			break;
		case 7:
			szTmp = szArray->GetAt(i);
			szTmp = TransformMediumType(szTmp);

			if (szTmp == "11")//混合装置
			{
		//		szTmp += "\r\n0\r\n0\r\n0";
				szTmp += _T("\r\n0\r\nType of data(0=apparent resistivity,1=resistance)\r\n0");
			}
			break;
		default://0和8
			szTmp = szArray->GetAt(i);
			break;
		}
		
		szTmp += "\r\n";
		fwrite(szTmp.GetBuffer(szTmp.GetLength()), 1, szTmp.GetLength(), m_resFile);
		fflush(m_resFile);
		szTmp.Empty();
	}
}

void CSaveInRes::Init2DFileHeader(CStringArray *szArray)
{
	//第一行:文件名
	//第二行:电极间距
	//第三行:装置类型

	//(混合装置多三行,比如说微分，多三行0)

	//第四行:测点总数
	//第五行:X位置类型(0代表第一个电极位置，1代表中点)
	//第六行:数据标志(0电阻率,1激电)

	CString szDate = _T("");
	CString szMediumType = _T("");

	for (int nIndex = 0; nIndex < szArray->GetSize(); nIndex++)
	{
		
		if (nIndex == 1)
		{
			m_szDistancePole = szArray->GetAt(nIndex);
			m_fPoleDistance = atof(m_szDistancePole);
		}
		
		if (nIndex == 2)
		{
			szMediumType = szArray->GetAt(nIndex);//将装置类型保存起来，后面根据这个变量进行不同的运算。
			szArray->SetAt(nIndex,TransformMediumType(szMediumType));//这里将装置类型转换成res软件的类型定义表
			if (szArray->GetAt(nIndex) == "NULL")
			{
				return;
			}

			if (szMediumType == "3")//微分属于混合排列类型,文件头需要添加三行0
			{
		//		szArray->InsertAt(nIndex + 1, "0");
		//		szArray->InsertAt(nIndex + 2, "0");
		//		szArray->InsertAt(nIndex + 3, "0");
			}
		}

		szDate += szArray->GetAt(nIndex);
		szDate += "\n";
	}

	fwrite(szDate.GetBuffer(szDate.GetLength()), 1, szDate.GetLength(), m_resFile);
	szDate.Empty();
	fflush(m_resFile);
}

BOOL CSaveInRes::AddElectrodeTitle(CStringArray* szArray)
{
	CString strContent;
	for (int nIndex = 0; nIndex < szArray->GetSize(); nIndex++)
	{
		strContent += szArray->GetAt(nIndex);
		//if (nIndex != szArray->GetSize() - 1)
		strContent += "\n";
	}

	//strContent += "\n";
	fwrite(strContent.GetBuffer(strContent.GetLength()), 1, strContent.GetLength(), m_resFile);
	strContent.Empty();
	fflush(m_resFile);
	return TRUE;
	return TRUE;
}

BOOL CSaveInRes::AddElectrodeCoordinatesInfo(CStringArray* szArray)
{
	CString strContent;
	for (int nIndex = 0; nIndex < szArray->GetSize(); nIndex++)
	{
		strContent += szArray->GetAt(nIndex);
		if (nIndex != szArray->GetSize() - 1)
			strContent += ",";
	}

	strContent += "\n";
	fwrite(strContent.GetBuffer(strContent.GetLength()), 1, strContent.GetLength(), m_resFile);
	strContent.Empty();
	fflush(m_resFile);
	return TRUE;
}

void CSaveInRes::InitCEFileHeader(CStringArray *szArray)
{
	switch (m_methodType)
	{
	case Res:
		InitCeResHeader(szArray);
		break;
	case IP:
		InitCeIpHeader(szArray);
		break;
	default:

		break;
	}
	
}

void CSaveInRes::InitCeResHeader(CStringArray *szArray)
{
	CString szTmp = _T("");

	for (int n = 0; n < szArray->GetSize(); n++)
	{
		switch (n)
		{
		case 0:
			szTmp = szArray->GetAt(n);
			break;
		case 1:
			szTmp = "Array Type\r\nSchlumberger" ;
			break;
		case 2:
			szTmp = "Number of data points\r\n" + szArray->GetAt(n);
			szTmp += "\r\nData Type\r\nResistivity";
			break;
		case 3:
			szTmp = "Error in measurements included\r\n" + szArray->GetAt(n);
			
			szTmp += "\r\nData section";
			break;
		}
		szTmp += "\r\n";
		fwrite(szTmp.GetBuffer(szTmp.GetLength()), 1, szTmp.GetLength(), m_resFile);
		fflush(m_resFile);
		szTmp.Empty();				
	}
}

void CSaveInRes::InitCeIpHeader(CStringArray *szArray)
{
	CString szTmp = _T("");
	
	for (int n = 0; n < szArray->GetSize(); n++)
	{
		switch (n)
		{
		case 0:
			szTmp = szArray->GetAt(n);
			break;
		case 1:
			szTmp = "Array Type\r\nSchlumberger" ;
			break;
		case 2:
			szTmp = "Number of data points\r\n" + szArray->GetAt(n);
			szTmp += "\r\nData Type\r\nIP\r\nChargeability\r\nmsec";
			break;
		case 3:
			szTmp += szArray->GetAt(n);
			break;
		case 4:
			szTmp = "Error in measurements included\r\n" + szArray->GetAt(n);
			szTmp += "\r\nData section";
			break;
		}
		szTmp += "\r\n";
		fwrite(szTmp.GetBuffer(szTmp.GetLength()), 1, szTmp.GetLength(), m_resFile);
		fflush(m_resFile);
		szTmp.Empty();				
	}
}

void CSaveInRes::CE_Schlumberger(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	switch (m_methodType)
	{
	case Res:
		CE_Res_Schlumberger(szSRCArray, szTARArray);
		break;
	case IP:
		CE_Res_Schlumberger(szSRCArray, szTARArray);//IP和电阻率数据数据主体一样
		break;
	}
}

void CSaveInRes::CE_Res_Schlumberger(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	CString szTmp = _T("");
	float fTmp = 0;

	for (int n = 0; n < szSRCArray->GetSize(); n++)	
	{
		if (n == 2)
		{
			fTmp = atof(szSRCArray->GetAt(n));
			szTmp.Format("%.4f", fabs(fTmp));
		}
		else
		{
			fTmp = atof(szSRCArray->GetAt(n));
			szTmp.Format("%.3f", fabs(fTmp));
		}
		
		szTARArray->Add(szTmp);
	}
}
void CSaveInRes::CE_Res_Marine(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	CString szTmp = _T("");
	float fTmp = 0;

	for (int n = 0; n < szSRCArray->GetSize(); n++)
	{
		if (n == 2)
		{
			fTmp = atof(szSRCArray->GetAt(n));
			szTmp.Format("%.4f", fabs(fTmp));
		}
		else
		{
			fTmp = atof(szSRCArray->GetAt(n));
			szTmp.Format("%.3f", fabs(fTmp));
		}

		szTARArray->Add(szTmp);
	}
}

void CSaveInRes::AddCEModel(CStringArray *szArray)
{
	CString szTmp = _T("");
	
	szTmp = szArray->GetAt(0);			//取得用户模型选项
	if (szTmp == "Yes")
	{
		//该部分待开发
	}
	else
	{
		szTmp = "User Starting Model Available\r\n";
		szTmp += "No";
	}
	
	szTmp += "\r\n";
	fwrite(szTmp.GetBuffer(szTmp.GetLength()), 1, szTmp.GetLength(), m_resFile);
	szTmp.Empty();
	fflush(m_resFile);
}

void CSaveInRes::Mid_Gardient_Scan_3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	/*
	按顺序传入A、B、M、N、R0 5个值
	*/

	int A = 0, B = 1, M = 2, N = 3, R = 4;
	int nSort[4] = {A, B, M, N};
	int nId = 0;
	_pos *Pos = NULL;

	CString szTmp = _T("");
	int nTmp = 0;
	int nPoleStart = 0;

	int nDistancePole = 0;

	nPoleStart = atoi(m_szPoleStart.GetBuffer(m_szPoleStart.GetLength()));
	
	Pos = new _pos();

	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}
    szTARArray->Add(_T("4"));
	float fXPoleDistance = atof(m_szX_DistancePole.GetBuffer(0));
    float fYPoleDistance = atof(m_szY_DistancePole.GetBuffer(0));
	int iTmp;
	for(int i = 0; i < 4; i++)
	{
		nId = nSort[i];

		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
	
		pole_to_pos_3d(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY,Pos, nTmp,m_byLineDirection);

		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}
	
		//AB是X轴的基准轴，所以此时Y的位置信息为0
		if (i < 2)
		{
			szTmp.Empty();
			szTmp.Format("%.3f",0 + m_fAPosX+i* m_fABDistance);
			szTARArray->Add(szTmp);
			
			szTmp.Empty();
			szTmp.Format("%.3f",0);
		    szTARArray->Add(szTmp);
		}
		else
		{	
			szTmp.Empty();
			szTmp.Format("%.3f",m_fAPosX + m_fMinAMDistance + fXPoleDistance * Pos->x);
			szTARArray->Add(szTmp);
			
			szTmp.Empty();
			szTmp.Format("%.3f",fYPoleDistance * Pos->y - m_fAPosY);
	    	szTARArray->Add(szTmp);
		}

	}

	int q = szSRCArray->GetSize();

	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);

	delete Pos;
	Pos = NULL;
}

void CSaveInRes::Dipole_Dipole_3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	/*
	按顺序传入A、B、M、N、R0 5个值
	*/

	int A = 0, B = 1, M = 2, N = 3, R = 4;
	int nSort[4] = {A, B, M, N};
	int nId = 0;
	_pos *Pos = NULL;

	CString szTmp = _T("");
	int nTmp = 0;
	int nPoleStart = 0;

	int nDistancePole = 0;

	nPoleStart = atoi(m_szPoleStart.GetBuffer(m_szPoleStart.GetLength()));
	
	Pos = new _pos();

	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 4; i++)
	{
		nId = nSort[i];

		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
	
		pole_to_pos_3d(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY,Pos, nTmp,m_byLineDirection);

		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		}

	int q = szSRCArray->GetSize();

	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);

	delete Pos;
	Pos = NULL;
}

void CSaveInRes::Set_3D_LandSign(int f_type)
{
	CString szTmp = _T("");

	szTmp.Format("Topography\r\n%d", f_type);

	szTmp += "\r\n";
	fwrite(szTmp.GetBuffer(szTmp.GetLength()), 1, szTmp.GetLength(), m_resFile);
	szTmp.Empty();
	fflush(m_resFile);


}

void CSaveInRes::Add_3D_LandFile(CStringArray *szSRCArray)
{
	CString szTmp = _T("");

	for (int n = 0; n < szSRCArray->GetSize(); n++)
	{
		szTmp += szSRCArray->GetAt(n);	
		szTmp +=" ";
	}
	
	szTmp += "\r\n";
	fwrite(szTmp.GetBuffer(szTmp.GetLength()), 1, szTmp.GetLength(), m_resFile);
	szTmp.Empty();
	fflush(m_resFile);
}

//联剖A
void CSaveInRes::Pole_dipole_2d(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	//szSRCArray这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
	//Pole_dipole_2d的计算值有四列，计算公式分别：
	//第一列:C1 * 间距
	//第二列:(p2 - p1) * 间距
	//第三列: (p1 - c1) / (p2 - p1)
	//第四列:R0
	//所有三极装置的中点都是MN
 	szTARArray->RemoveAll();
	
	//第一列:
	DOUBLE firstRow;
	DOUBLE fDistancePole;
	CString szFirstRow;
	//modified by lsq 20160716
	//将记录点的位置改为M点位置
//	firstRow = (0.5*(atof(szSRCArray->GetAt(2))+atof(szSRCArray->GetAt(3)))-1) * atof(m_szDistancePole.GetBuffer(0));
	//firstRow = (atof(szSRCArray->GetAt(2)) - 1) * atof(m_szDistancePole.GetBuffer(0));
	firstRow = (atof(szSRCArray->GetAt(0)) - 1 )* atof(m_szDistancePole.GetBuffer(0));
	szFirstRow.Format("%.3f",fabs(firstRow));
	szTARArray->Add(szFirstRow);
	
	//第二列:
	CString szSencondRow;
	DOUBLE fSencondRow;
	DOUBLE fdistance;
	fdistance = atof(m_szDistancePole);//m_szDistancePole成员变量，保存了电极间距
	fSencondRow = ( atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)) ) * atof(m_szDistancePole.GetBuffer(0));
	szSencondRow.Format("%.3f",fabs(fSencondRow));
	szTARArray->Add(szSencondRow);
	
	//第三列
	DOUBLE thirdRow;
	CString szThirdRow;
	thirdRow = 1.0*(atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0))) / (atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)));

// 	DOUBLE a = atof(szSRCArray->GetAt(2));
// 	DOUBLE b = atof(szSRCArray->GetAt(0));
// 	DOUBLE c = atof(szSRCArray->GetAt(3));
// 	DOUBLE d = atof(szSRCArray->GetAt(2));
// 	DOUBLE e = (a - b)/ (c - d);

	szThirdRow.Format("%.3f",fabs(thirdRow));
	szTARArray->Add(szThirdRow);
	
	//第四列
	DOUBLE FourthRow;
	CString szFourthRow;
	FourthRow = atof(szSRCArray->GetAt(5)) / atof(szSRCArray->GetAt(6)) * (2* VAL_PI/(1.0/(atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0))) - 1.0/ (atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(0)))));
	szFourthRow.Format("%.3f",FourthRow);
	szTARArray->Add(szFourthRow);

}

//联剖b
void CSaveInRes::Reverse_Pole_dipole_2d(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	//szSRCArray这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
	//Reverse_Pole_dipole_2d的计算值有四列，计算公式分别：
	//第一列:P1 * 间距
	//第二列:(p1 - p2) * 间距 
	//第三列: (c2-p2)/(p2-p1)
	//第四列:R0
	
	szTARArray->RemoveAll();
	
	//第一列:
	DOUBLE firstRow;
	DOUBLE fDistancePole;
	CString szFirstRow;

	//modified by lsq 20160716
	//将记录点的位置改为N点位置
//	firstRow = (0.5*(atof(szSRCArray->GetAt(2))+atof(szSRCArray->GetAt(3)))-1) * atof(m_szDistancePole.GetBuffer(0));
//	firstRow = (atof(szSRCArray->GetAt(3)) - 1) * atof(m_szDistancePole.GetBuffer(0));
	//将记录点的位置改为M点位置  x
	firstRow = (atof(szSRCArray->GetAt(3)) - 1) * atof(m_szDistancePole.GetBuffer(0));
	szFirstRow.Format("%.3f",fabs(firstRow));
	szTARArray->Add(szFirstRow);
	
	//第二列:
	CString szSencondRow;
	DOUBLE fSencondRow;
// 	DOUBLE fdistance;
// 	fdistance = atof(m_szDistancePole);//m_szDistancePole成员变量，保存了电极间距
	fSencondRow = ( atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(3)) ) * atof(m_szDistancePole.GetBuffer(0));
	szSencondRow.Format("%.3f",fabs(fSencondRow));
	szTARArray->Add(szSencondRow);
	
	//第三列,mnb的n为负值
	DOUBLE thirdRow;
	CString szThirdRow;
	thirdRow = (atof(szSRCArray->GetAt(1)) - atof(szSRCArray->GetAt(2))) / (atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(3)));
	szThirdRow.Format("-%.3f",fabs(thirdRow));
	szTARArray->Add(szThirdRow);
	
	//第四列
	DOUBLE FourthRow;
	CString szFourthRow;
	FourthRow = atof(szSRCArray->GetAt(5)) / atof(szSRCArray->GetAt(6)) * (2 * VAL_PI / (1.0 / (atof(szSRCArray->GetAt(1)) - atof(szSRCArray->GetAt(3))) - 1.0 / (atof(szSRCArray->GetAt(1)) - atof(szSRCArray->GetAt(2)))));
	szFourthRow.Format("%.3f",FourthRow);
	szTARArray->Add(szFourthRow);
}

//微分排列2d
void CSaveInRes::Wenner_Y_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	//Wenner_Y_2D这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
	//Reverse_Pole_dipole_2d的计算值有四列，计算公式分别：
	//第一列:记录点的位置，在这里为BM的中点
	//第二列:电极间距，极AM的距离
	//第三列:R0
    //所有温钠装置的中点都是MN的中点
	szTARArray->RemoveAll();

	//电极间距
	DOUBLE fdistance;
	fdistance = atof(m_szDistancePole);//m_szDistancePole成员变量，保存了电极间距

	//第一列:
	DOUBLE firstRow;
	DOUBLE fDistancePole;
	CString szFirstRow;
	
	firstRow = (0.5*(atof(szSRCArray->GetAt(2)) + atof(szSRCArray->GetAt(1))) -1) * fdistance;
	szFirstRow.Format("%.3f",fabs(firstRow) );
	szTARArray->Add(szFirstRow);
	
	//第二列:
	CString szSencondRow;
	DOUBLE fSencondRow;
	fSencondRow = (atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0))) * fdistance;
	szSencondRow.Format("%.3f",fabs(fSencondRow));
	szTARArray->Add(szSencondRow);
	
	//第三列
	DOUBLE thirdRow;
	CString szThirdRow;
	thirdRow = atof(szSRCArray->GetAt(4));
	szThirdRow.Format("%.3f",thirdRow);
	szTARArray->Add(szThirdRow);

}

void CSaveInRes::Pole_Dipole_3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	int A = 0, B = 1, M = 2, N = 3, R = 4;
	int nSort[3] = {A, M, N};
	int nId = 0;
	_pos *Pos = NULL;
	
	CString szTmp = _T("");
	int nTmp = 0;
	int nPoleStart = 0;
	
	int nDistancePole = 0;
	
	nPoleStart = atoi(m_szPoleStart.GetBuffer(m_szPoleStart.GetLength()));
	
	Pos = new _pos();
	
	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 3; i++)
	{
		nId = nSort[i];
		
		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
		
		pole_to_pos_3d(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY, Pos, nTmp,m_byLineDirection);

		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
	}
	
	int q = szSRCArray->GetSize();
	
	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);
	
	delete Pos;
	Pos = NULL;
}

void CSaveInRes::Pole_Pole_3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{	
	int A = 0, B = 1, M = 2, N = 3, R = 4;
	int nSort[2] = {A, M};
	int nId = 0;
	_pos *Pos = NULL;
	
	CString szTmp = _T("");
	int nTmp = 0;
	int nPoleStart = 0;
	
	int nDistancePole = 0;
	
	nPoleStart = atoi(m_szPoleStart.GetBuffer(m_szPoleStart.GetLength()));
	
	Pos = new _pos();
	
	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 2; i++)
	{
		nId = nSort[i];
		
		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
		pole_to_pos_3d(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY, Pos, nTmp, m_byLineDirection);
		//pole_to_pos_by_normal_axis(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY,Pos, nTmp);
		
		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
	}
	
	int q = szSRCArray->GetSize();
	
	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);
	
	delete Pos;
	Pos = NULL;
}

//中梯,混合类型
void CSaveInRes::Gradient_3D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	/*
	按顺序传入A、B、M、N、R0 5个值
	*/

	int A = 0, B = 1, M = 2, N = 3, R = 4;
	int nSort[2] = {/*A, B, */M, N};
	int nId = 0;
	_pos *Pos = NULL;

	CString szTmp = _T("");
	int nTmp = 0;
	int nPoleStart = 0;

	int nDistancePole = 0;

	nPoleStart = atoi(m_szPoleStart.GetBuffer(m_szPoleStart.GetLength()));
	
	Pos = new _pos();

	szTARArray->Add("2");	//混合类型前面要跟点击数量

	_rect tmpRect;
	//如果是按照Y轴合并，由于m_rect已经是旋转过后的矩阵，那么在计算时，还要将其旋转回去进行计算
	//因为pole_to_pos函数是按照电极的长度对于X轴来进行解析的
	if (2 == m_i3DRollDiecetion)
	{
		tmpRect.x0 = m_rect.y0;
		tmpRect.x1 = m_rect.y1;
		tmpRect.y0 = m_rect.x0;
		tmpRect.y1 = m_rect.x1;
	}
	else
	{
		tmpRect.x0 = m_rect.x0;
		tmpRect.x1 = m_rect.x1;
		tmpRect.y0 = m_rect.y0;
		tmpRect.y1 = m_rect.y1;
	}

	int iTmp;
	for(int i = 0; i < 2; i++)
	{
		nId = nSort[i];

		szTmp = szSRCArray->GetAt(nId);
		nTmp = atoi(szTmp.GetBuffer(0));
	
		pole_to_pos_3d(&tmpRect, nPoleStart, m_nPoleStepX, m_nPoleStepY,Pos, nTmp, m_byLineDirection);

		//如果是按照Y轴方向合并，则需要进行逻辑坐标轴的转换
		//此时转换之需要颠倒X轴和Y轴的位置就好
		if (2 == m_i3DRollDiecetion)
		{
			iTmp = Pos->x;
			Pos->x = Pos->y;
			Pos->y = iTmp;
		}
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->x * atof(m_szX_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
		
		szTmp.Empty();
		szTmp.Format("%.3f",(float)(Pos->y * atof(m_szY_DistancePole.GetBuffer(0))));
		szTARArray->Add(szTmp);
	}

	int q = szSRCArray->GetSize();

	szTmp.Empty();
	szTmp = szSRCArray->GetAt(4);
	szTARArray->Add(szTmp);

	delete Pos;
	Pos = NULL;
}

void CSaveInRes::Dipole_Dipole_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	//szSRCArray这五个值得顺序分别是 0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0
	//Pole_dipole_2d的计算值有四列，计算公式分别：
	//第一列:C1 * 间距
	//第二列:(p2 - p1) * 间距
	//第三列:(p1 - c1) / (p2 - p1)
	//第四列:R0
	//在偶极-偶极装置里，中点是C1和P1之间的中点
	
	szTARArray->RemoveAll();
	
	//第一列:
	DOUBLE firstRow;
	DOUBLE fDistancePole;
	CString szFirstRow;
	
	firstRow = (0.5*(atof(szSRCArray->GetAt(0)) + atof(szSRCArray->GetAt(2)))-1) * atof(m_szDistancePole.GetBuffer(0));
	szFirstRow.Format("%.3f",fabs(firstRow));
	szTARArray->Add(szFirstRow);
	
	//第二列:
	CString szSencondRow;
	DOUBLE fSencondRow;
	DOUBLE fdistance;
	fdistance = atof(m_szDistancePole);//m_szDistancePole成员变量，保存了电极间距
	fSencondRow = ( atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)) ) * atof(m_szDistancePole.GetBuffer(0));
	szSencondRow.Format("%.4f",fabs(fSencondRow));
	szTARArray->Add(szSencondRow);
	
	//第三列
	DOUBLE thirdRow;
	CString szThirdRow;
	thirdRow = ( atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0)) ) / ( atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)) );	
	szThirdRow.Format("%.5f",fabs(thirdRow));
	szTARArray->Add(szThirdRow);
	
	//第四列
	DOUBLE FourthRow;
	CString szFourthRow;
	FourthRow = atof(szSRCArray->GetAt(4));
	szFourthRow.Format("%.4f",FourthRow);
	szTARArray->Add(szFourthRow);
}

void CSaveInRes::SetPoleStep(int f_StepX, int f_StepY)
{
	m_nPoleStepX = f_StepX;
	m_nPoleStepY = f_StepY;
}


void CSaveInRes::Cross_Hole_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	//szSRCArray这顺序分别为0:A的位置 1:B的位置 2:M的位置 3:N的位置 4:R0

    //所生成的内容的格式为
	//1:已经使用电极的个数，这这种方法中，值固定为4
	//2,3列:A的X,Z位置 
	//4,5列:B的X,Z位置 
	//6,7列:M的X,Z位置
	//8,9列:N的X,Z位置 
	//10列：R0
	
	float fXPos = 0;
	float fZPos = 0;
    CString strXPos = _T("");
	CString strZPos = _T("");
	CString strR0 = _T("");
	szTARArray->RemoveAll();

	//如果只有一个字符串_T("1")的话，表示此时需要添加跨孔虚拟电极信息
	if ((2 == szSRCArray->GetSize()) && (_T("1") == szSRCArray->GetAt(0)))
	{
		m_iEAmount = atoi(szSRCArray->GetAt(1));
		GetCrossHoleSurfaceElectrode(m_fPoleDistance, m_fCrossHoleParam, m_iEAmount, szTARArray);
		return;
	}
	
	//第1列
	szTARArray->Add(_T("4"));
//	int nMid = (atof(szSRCArray->GetAt(0)) + atof(szSRCArray->GetAt(1)) -1) /2;
    int nMid = m_iEAmount/2;
    int nFirstRealElecIndex = AR_CROSS_HOLE_SURFACE_ELECTRODE_NUM/2;
	float fFirstRealElecXPos = 2*m_fPoleDistance;
    float fABDistance = m_fCrossHoleParam;
	int iAMParam = (atof(szSRCArray->GetAt(0)) < nMid) ? 0 : 1;
	int iBNParam = (atof(szSRCArray->GetAt(0)) < nMid) ? 1 : 0;
 
	//第2.3列:A点的X,Z位置
    //added by lsq20160613 在这里Z得位置多加一个电极间距，根据祝工的意见
	//modified by lsq 20171009 A在做，B在右边，A和M在同一条井内，电极号从小到大
	//比如24电极，那么该井内就是1到12，另外一条则是13到24
	//modifyed by lsq 10171016 根据最新的文档，孔深的位置=初始孔深+(N-1)*电极间距
//	fXPos = fFirstRealElecXPos + fABDistance;	
//	fZPos = m_fPoleDistance*(nMid - atof(szSRCArray->GetAt(0)) + 1);
	fXPos = fFirstRealElecXPos;
	fZPos = m_fInitiDepth + m_fPoleDistance*(atof(szSRCArray->GetAt(0))-iAMParam*nMid-1);	
	strXPos.Format(_T("%.2f"),fXPos);
	strZPos.Format(_T("%.2f"),fZPos);
	szTARArray->Add(strXPos);
	szTARArray->Add(strZPos);
	
	//第4,5列:B的X,Z位置 
	//由于A,B处在同一横向，所以Z轴的坐标是相同的
	strXPos.Empty();
	strZPos.Empty();
	fXPos = fFirstRealElecXPos+fABDistance;
	fZPos = m_fInitiDepth+ m_fPoleDistance*(atof(szSRCArray->GetAt(1)) -iBNParam*nMid-1);
	strXPos.Format(_T("%.2f"),fXPos);
	strZPos.Format(_T("%.2f"),fZPos);
	szTARArray->Add(strXPos);
	szTARArray->Add(strZPos);
	
	//第6,7列:M的X,Z位置 
	
// 	fXPos = fFirstRealElecXPos + fABDistance;
// 	fZPos = m_fPoleDistance*(nMid - atof(szSRCArray->GetAt(2)) + 1);
    fXPos = fFirstRealElecXPos;
	fZPos = m_fInitiDepth+m_fPoleDistance*(atof(szSRCArray->GetAt(2))-iAMParam*nMid-1);
	strXPos.Empty();
	strZPos.Empty();
	strXPos.Format(_T("%.2f"),fXPos);
	strZPos.Format(_T("%.2f"),fZPos);
	szTARArray->Add(strXPos);
	szTARArray->Add(strZPos);
	
	//第8,9列:N的X,Z位置 ,N的Z轴和M的轴Z位置是相同的
	fXPos = fFirstRealElecXPos+fABDistance;
	fZPos = m_fInitiDepth + m_fPoleDistance*(atof(szSRCArray->GetAt(3)) -iBNParam*nMid-1);
	strXPos.Empty();
	strZPos.Empty();
	strXPos.Format(_T("%.2f"),fXPos);
	strZPos.Format(_T("%.2f"),fZPos);
	szTARArray->Add(strXPos);
	szTARArray->Add(strZPos);

	//第10列，R0
	strR0.Format(_T("%.4f"),atof(szSRCArray->GetAt(4)));
	szTARArray->Add(strR0);
}

void CSaveInRes::GetCrossHoleSurfaceElectrode(float fPoleDis, float fCrosshole, int nHoleElectrodeNum,CStringArray* strArray)
{
	strArray->RemoveAll();
	int nElectordeCnt = AR_CROSS_HOLE_SURFACE_ELECTRODE_NUM;
	float fXPos = 0;
	CString strInfo = _T("");
	CString strText = _T("Surface electrode\n");
	strArray->Add(strText);
	float fElecPos[2];
	int iVirElectCnt = (int)(fCrosshole/fPoleDis);
    float fVirElecSpace = fCrosshole/iVirElectCnt;
	strInfo.Format(_T("%d\n"),6+iVirElectCnt-1);
	strArray->Add(strInfo);
//	int nFirstRealElectIndex = nElectordeCnt/2;
	//加入前3个电极，2个虚拟电极和一个实际电极，这3个电极之间的距离为电极间距
	int nFirstRealElectIndex = 3;
	int i = 0;
	for (; i < nFirstRealElectIndex; i++)
	{
		fXPos = i*fPoleDis;
		strInfo.Empty();
		strInfo.Format(_T("%.2f, 0.00\n"),fXPos);
		strArray->Add(strInfo);
	}
	if ( fPoleDis <= 0)
	{
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("电极间距必须为正数!"));
		else
			MessageBoxEx(NULL, _T("The pole space must be postive number!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		fPoleDis = 10;
	}
	if (fCrosshole <= fPoleDis)
	{
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("跨孔电极间距必须大于孔内电极间距!")); 
		else
			MessageBoxEx(NULL, _T("The space of crss-hole must greatter than pole space!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		fCrosshole = fPoleDis;
	}

	//添加虚拟电极

	fElecPos[0] = fXPos;
	for (i = 1; i < iVirElectCnt; i++)
	{
		fXPos = fElecPos[0] + i*fVirElecSpace;
		strInfo.Empty();
		strInfo.Format(_T("%.2f, 0.00\n"),fXPos);
		strArray->Add(strInfo);
	}
	//添加第二个实际电极
	fElecPos[1] = fElecPos[0] + fCrosshole;
	strInfo.Empty();
	strInfo.Format(_T("%.2f, 0.00\n"),fElecPos[1]);
	strArray->Add(strInfo);


// 	//添加两孔中间的虚拟电极
// 	fXPos += fCrosshole*fPoleDis/2;
// 	strInfo.Empty();
// 	strInfo.Format(_T("%.1f, 0.0\n"),fXPos);
// 	strArray->Add(strInfo);
	
	//第二个实际电极的位置
//     i = 0;
// 	fXPos = (nFirstRealElectIndex-1) * fPoleDis + fCrosshole * fPoleDis;
// 	fElecPos[1] = fXPos;
// 	int nRemain = nElectordeCnt - nFirstRealElectIndex - 1;
	//添加第二个实际电极右边的另外2个虚拟电极
	for (i = 1; i <= 2; i++)
	{
		fXPos = fElecPos[1] + i*fPoleDis;
		strInfo.Empty();
		strInfo.Format(_T("%.2f, 0.00\n"),fXPos);
		strArray->Add(strInfo);
	}
	
	
	int nHoleElecNumMid = nHoleElectrodeNum / 2;
	i = 0;
	float fZPos = 0;
	for (; i < 2; i++)
	{
		//添加钻孔1的电极信息
		strText.Empty();
		strText.Format(_T("Drill hole %d\n %d\n"),i+1,nHoleElecNumMid);
		strArray->Add(strText);
		
		for (int j = 0; j < nHoleElecNumMid; j++)
		{
			strInfo.Empty();
			//modifyed by lsq 20171016 按照张杰最新的方案，将电极的深度信息也修改成和实际的一样
			strInfo.Format(_T("%.2f, %.2f\n"), fElecPos[i], m_fInitiDepth+j*fPoleDis);
			strArray->Add(strInfo);
		}
	}
	strText = _T("measured points");
	strArray->Add(strText);
}

void CSaveInRes::Wenner_Schlumberger_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
    if (szSRCArray->GetSize()<5)
    {
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("数据格式错误在Wenner_Schlumberger"));
		else
			MessageBoxEx(NULL, _T("Data format error in Wenner_Schlumberger!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		return;
    }
	szTARArray->RemoveAll();
    float fVal = 0;
	CString strVal = _T("");

	//第一列A点的位置
	fVal = (0.5*(atof(szSRCArray->GetAt(2))+atof(szSRCArray->GetAt(3)))-1);
	strVal.Format(_T("%.2f"), m_fPoleDistance*fVal);
	szTARArray->Add(strVal);

	//第二列最小电极间距
	strVal.Empty();
	strVal.Format(_T("%.2f"), m_fPoleDistance*(atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2))));
	szTARArray->Add(strVal);

	//第三列间隔系数=n/x
	fVal = (atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0))) / (atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)));
	strVal.Empty();
	strVal.Format(_T("%.5f"), fVal);
	szTARArray->Add(strVal);

	//第四列RO的值
	fVal = atof(szSRCArray->GetAt(4));
	strVal.Empty();
	strVal.Format(_T("%.2f"),fVal);
	szTARArray->Add(strVal);

}

void CSaveInRes::A_M_ARRAY_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
    if (szSRCArray->GetSize()<5)
    {
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("数据格式错误在A_M_array"));
		else
			MessageBoxEx(NULL, _T("Data format error in A_M_array!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		return;
    }
	szTARArray->RemoveAll();
    float fVal = 0;
	CString strVal = _T("");
	
	//第一列AM中点的位置
	fVal = 0.5*(atof(szSRCArray->GetAt(0))+atof(szSRCArray->GetAt(2)))-1;
	strVal.Format(_T("%.2f"), fVal*m_fPoleDistance);
	szTARArray->Add(strVal);
	
	//第二列AM之间的电极间距
	strVal.Empty();
	strVal.Format(_T("%.2f"), m_fPoleDistance*(atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0))));
	szTARArray->Add(strVal);
	
	//第三列RO的值
	fVal = atof(szSRCArray->GetAt(4));
	strVal.Empty();
	strVal.Format(_T("%.2f"),fVal);
	szTARArray->Add(strVal);
	
}
void CSaveInRes::SetMidGardientInfo(float fAPosX, float fAPosY, float fABDistance, float fMinAMDistance)
{
	m_fAPosX = fAPosX;
	m_fAPosY = fAPosY;
	m_fABDistance = fABDistance;
	m_fMinAMDistance = fMinAMDistance;

}

void CSaveInRes::AMN_ARRAY_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
    if (szSRCArray->GetSize()<5)
    {
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("数据格式错误在A_M_array"));
		else
			MessageBoxEx(NULL, _T("Data format error in A_M_array!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		return;
    }
	szTARArray->RemoveAll();
    float fVal = 0;
	CString strVal = _T("");
	
	//第一列A点的位置
	fVal = atof(szSRCArray->GetAt(2)) - 1;
	strVal.Format(_T("%.2f"), m_fPoleDistance*fVal);
	szTARArray->Add(strVal);
	
	//第二列AM之间的电极间距
	strVal.Empty();
	strVal.Format(_T("%.2f"), m_fPoleDistance*(atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2))));
	szTARArray->Add(strVal);

	//第三列AM/MN
	fVal =  (atof(szSRCArray->GetAt(2)) - atof(szSRCArray->GetAt(0))) / (atof(szSRCArray->GetAt(3)) - atof(szSRCArray->GetAt(2)));
	strVal.Empty();
	strVal.Format(_T("%.2f"), fVal);
	szTARArray->Add(strVal);

	//第三列RO的值
	fVal = atof(szSRCArray->GetAt(4));
	strVal.Empty();
	strVal.Format(_T("%.2f"),fVal);
	szTARArray->Add(strVal);
	
}

void CSaveInRes::Custom_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
    if (szSRCArray->GetSize()<5)
    {
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("数据格式错误在Custom_array"));
		else
			MessageBoxEx(NULL, _T("Data format error in Custom_array!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		return;
    }
	szTARArray->RemoveAll();
    float fVal = 0;
	CString strZVal = _T("0.000");
	int iRealPoleCnt = 0;
	for (int i = 0; i < 4; i++)
	{
		if (atof(szSRCArray->GetAt(i)) > 0)
		{
			iRealPoleCnt++;
		}
	}
	CString strVal = _T("");
	strVal.Format(_T("%d"),iRealPoleCnt);
	//modified by lsq 20160718
	//第1列电极的个数，在这里修改为实际电极个数
	szTARArray->Add(strVal);
	
	
	//第2,3列A点的位置和Z位置
	fVal = atof(szSRCArray->GetAt(0));
	if (fVal > 0)
	{
		strVal.Empty();
		strVal.Format(_T("%.3f"), m_fPoleDistance*fVal);
		szTARArray->Add(strVal);
		szTARArray->Add(strZVal);
	}
	
	
	//第4,5列B点的位置和Z位置
	fVal = atof(szSRCArray->GetAt(1));
	if (fVal > 0)
	{
		strVal.Empty();
		strVal.Format(_T("%.3f"), m_fPoleDistance*fVal);
		szTARArray->Add(strVal);
		szTARArray->Add(strZVal);
	}
	
	
	//第6,7列M点的位置和Z位置
	fVal = atof(szSRCArray->GetAt(2));
	if (fVal > 0)
	{
		strVal.Empty();
		strVal.Format(_T("%.3f"), m_fPoleDistance*fVal);
		szTARArray->Add(strVal);
		szTARArray->Add(strZVal);
	}
	
	//第8,9列N点的位置和Z位置
	fVal = atof(szSRCArray->GetAt(3));
	if (fVal > 0)
	{
		strVal.Empty();
		strVal.Format(_T("%.3f"), m_fPoleDistance*fVal);
		szTARArray->Add(strVal);
		szTARArray->Add(strZVal);
	}
	
	//第10列RO的值
	fVal = atof(szSRCArray->GetAt(4));
	strVal.Empty();
	strVal.Format(_T("%.3f"),fVal);
	szTARArray->Add(strVal);
}

void CSaveInRes::CrossHole_2D(CStringArray *szSRCArray, CStringArray *szTARArray)
{
	if (szSRCArray->GetSize() < 5)
	{
		if (LANG_ZHCN == g_iUILanguage)
			AfxMessageBox(_T("数据格式错误在Crosshole_array"));
		else
			MessageBoxEx(NULL, _T("Data format error in Crosshole_array!"), STRING_MESSAGEBOXEX_TITLE, MB_OK, MAKELANGID(LANG_ENGLISH, SUBLANG_ENGLISH_US));
		return;
	}
	szTARArray->RemoveAll();
	float fVal = 0;
	CString strZVal = _T("0.000");
	int iRealPoleCnt = 0;
	for (int i = 0; i < 4; i++)
	{
		if (atof(szSRCArray->GetAt(i)) > 0)
		{
			iRealPoleCnt++;
		}
	}
	CString strVal = _T("");
	strVal.Format(_T("%d"), iRealPoleCnt);
	//modified by lsq 20160718
	//第1列电极的个数，在这里修改为实际电极个数
	szTARArray->Add(strVal);


	//第2,3列A点的位置和Z位置
	szTARArray->Add(szSRCArray->GetAt(0));

	//第4,5列B点的位置和Z位置
	szTARArray->Add(szSRCArray->GetAt(1));


	//第6,7列M点的位置和Z位置
	szTARArray->Add(szSRCArray->GetAt(2));

	//第8,9列N点的位置和Z位置
	szTARArray->Add(szSRCArray->GetAt(3));

	//第10列RO的值
	fVal = atof(szSRCArray->GetAt(4));
	strVal.Empty();
	strVal.Format(_T("%.3f"), fVal);
	szTARArray->Add(strVal);
}