// XF_FLT_OMPDlg.cpp : 实现文件
//

#include "stdafx.h"
#include "XF_FLT_OMP.h"
#include "XF_FLT_OMPDlg.h"
#include "afxdialogex.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#endif

int ThreadCount = 5;
int Progress1, Progress2;

int iPro = 40;
int ProCount[40+1] = {0};
int iFinishThread;

bool Is45n;

int w, h;
double x, y;
double cell;

int Rw, Rh;
double Rx, Ry;
double Rcell;

double m_Rotate;

int Bw1, Bh1, Bw2, Bh2;

FILE* fp1;    //源FLT
FILE* fp2;    //目标FLT
CCriticalSection critical_section;
CCriticalSection critical_section_progress;
CCriticalSection critical_section_finish;

CString s_Flt2,s_Flt;

// 用于应用程序“关于”菜单项的 CAboutDlg 对话框

class CAboutDlg : public CDialogEx
{
public:
	CAboutDlg();

// 对话框数据
	enum { IDD = IDD_ABOUTBOX };

	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV 支持

// 实现
protected:
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialogEx(CAboutDlg::IDD)
{
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialogEx::DoDataExchange(pDX);
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialogEx)
END_MESSAGE_MAP()


// CXF_FLT_OMPDlg 对话框




CXF_FLT_OMPDlg::CXF_FLT_OMPDlg(CWnd* pParent /*=NULL*/)
	: CDialogEx(CXF_FLT_OMPDlg::IDD, pParent)
	, m_Flt(_T(""))
	, m_Offsets(_T(""))
	, m_Path(_T(""))
	, m_Name(_T(""))
	, m_ThreadCount(20)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CXF_FLT_OMPDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialogEx::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT1, m_Flt);
	DDX_Text(pDX, IDC_EDIT2, m_Offsets);
	DDX_Text(pDX, IDC_EDIT3, m_Path);
	DDX_Text(pDX, IDC_EDIT4, m_Name);
	DDX_Control(pDX, IDC_COMMONDIALOG1, m_CCommondialog1);
	DDX_Control(pDX, IDC_COMMONDIALOG2, m_CCommondialog2);
	DDX_Control(pDX, IDC_PROGRESS1, m_CProgressCtrl1);
	DDX_Text(pDX, IDC_EDIT5, m_ThreadCount);
}

BEGIN_MESSAGE_MAP(CXF_FLT_OMPDlg, CDialogEx)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDOK, &CXF_FLT_OMPDlg::OnBnClickedOk)
	ON_BN_CLICKED(IDCANCEL, &CXF_FLT_OMPDlg::OnBnClickedCancel)
	ON_BN_CLICKED(IDC_BUTTON1, &CXF_FLT_OMPDlg::OnBnClickedButton1)
	ON_BN_CLICKED(IDC_BUTTON2, &CXF_FLT_OMPDlg::OnBnClickedButton2)
	ON_BN_CLICKED(IDC_BUTTON3, &CXF_FLT_OMPDlg::OnBnClickedButton3)
	ON_BN_CLICKED(IDC_BUTTON4, &CXF_FLT_OMPDlg::OnBnClickedButton4)
	ON_WM_TIMER()
END_MESSAGE_MAP()


// CXF_FLT_OMPDlg 消息处理程序

BOOL CXF_FLT_OMPDlg::OnInitDialog()
{
	CDialogEx::OnInitDialog();

	// 将“关于...”菜单项添加到系统菜单中。

	// IDM_ABOUTBOX 必须在系统命令范围内。
	ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
	ASSERT(IDM_ABOUTBOX < 0xF000);

	CMenu* pSysMenu = GetSystemMenu(FALSE);
	if (pSysMenu != NULL)
	{
		BOOL bNameValid;
		CString strAboutMenu;
		bNameValid = strAboutMenu.LoadString(IDS_ABOUTBOX);
		ASSERT(bNameValid);
		if (!strAboutMenu.IsEmpty())
		{
			pSysMenu->AppendMenu(MF_SEPARATOR);
			pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
		}
	}

	// 设置此对话框的图标。当应用程序主窗口不是对话框时，框架将自动
	//  执行此操作
	SetIcon(m_hIcon, TRUE);			// 设置大图标
	SetIcon(m_hIcon, FALSE);		// 设置小图标

	// TODO: 在此添加额外的初始化代码

	return TRUE;  // 除非将焦点设置到控件，否则返回 TRUE
}

void CXF_FLT_OMPDlg::OnSysCommand(UINT nID, LPARAM lParam)
{
	if ((nID & 0xFFF0) == IDM_ABOUTBOX)
	{
		CAboutDlg dlgAbout;
		dlgAbout.DoModal();
	}
	else
	{
		CDialogEx::OnSysCommand(nID, lParam);
	}
}

// 如果向对话框添加最小化按钮，则需要下面的代码
//  来绘制该图标。对于使用文档/视图模型的 MFC 应用程序，
//  这将由框架自动完成。

void CXF_FLT_OMPDlg::OnPaint()
{
	if (IsIconic())
	{
		CPaintDC dc(this); // 用于绘制的设备上下文

		SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);

		// 使图标在工作区矩形中居中
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon + 1) / 2;
		int y = (rect.Height() - cyIcon + 1) / 2;

		// 绘制图标
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialogEx::OnPaint();
	}
}

//当用户拖动最小化窗口时系统调用此函数取得光标
//显示。
HCURSOR CXF_FLT_OMPDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}



void CXF_FLT_OMPDlg::OnBnClickedOk()
{
	// TODO: 在此添加控件通知处理程序代码
	CDialogEx::OnOK();
}


void CXF_FLT_OMPDlg::OnBnClickedCancel()
{
	// TODO: 在此添加控件通知处理程序代码
	CDialogEx::OnCancel();
}


void CXF_FLT_OMPDlg::OnBnClickedButton1()
{
	// TODO: 在此添加控件通知处理程序代码
	m_CCommondialog1.put_DialogTitle(CComBSTR("打开FLT"));
	m_CCommondialog1.put_Filter(CComBSTR("FLT文件(*.flt)|*.flt"));
	m_CCommondialog1.put_InitDir(_T("FLT\\"));
	m_CCommondialog1.ShowOpen();

	CComBSTR bstrFilePath;

	bstrFilePath = m_CCommondialog1.get_FileName();

	if (SysStringLen(bstrFilePath)  == 0) 
	{
		MessageBox(_T("不是有效的flt文件"));
		return;
	}
	GetDlgItem(IDC_EDIT1)->SetWindowText(bstrFilePath);
}

// 读offsets文件
void CXF_FLT_OMPDlg::OnBnClickedButton2()
{
	// TODO: 在此添加控件通知处理程序代码
	m_CCommondialog2.put_DialogTitle(CComBSTR("打开offsets"));
	m_CCommondialog2.put_Filter(CComBSTR("offsets文件(*.dat)|*.dat"));
	m_CCommondialog2.put_InitDir(_T("FLT\\"));
	m_CCommondialog2.ShowOpen();

	CComBSTR bstrFilePath;

	bstrFilePath = m_CCommondialog2.get_FileName();

	if (SysStringLen(bstrFilePath)  == 0) 
	{
		MessageBox(_T("不是有效的offsets文件"));
		return;
	}
	GetDlgItem(IDC_EDIT2)->SetWindowText(bstrFilePath);


}

void CXF_FLT_OMPDlg::OnBnClickedButton3()
{
	// TODO: 在此添加控件通知处理程序代码
	LPITEMIDLIST   pidlPath;   
	TCHAR   pszDplName[MAX_PATH];   
	BROWSEINFO   biFile;   
	IMalloc   *pDllMlc;   
	biFile.hwndOwner   =   m_hWnd;   
	biFile.pidlRoot   =   NULL;   
	biFile.pszDisplayName   =   NULL;   
	biFile.lpszTitle   =   _T("请选择要打开的文件夹:");   
	biFile.ulFlags   =   BIF_RETURNFSANCESTORS;   
	biFile.lpfn   =   NULL;   
	biFile.lParam   =   NULL;   
	biFile.iImage   =   0;   
	if(SUCCEEDED(SHGetMalloc(&pDllMlc)))   
	{   
		pidlPath=SHBrowseForFolder(&biFile);   
		if(pidlPath   !=   NULL)   
		{
			SHGetPathFromIDList(pidlPath,pszDplName);    //pszDplName:文件的路径地址 
			//	this->MessageBox(pszDplName);                //显示获取的地址 
			//	GetDlgItem(IDC_EDIT6)->SetWindowText(pszDplName);
			GetDlgItem(IDC_EDIT3)->SetWindowText(pszDplName);
		}
		else 
			return;   
		pDllMlc->Free(NULL);   
		pDllMlc->Release();   
	}
}

CString CXF_FLT_OMPDlg::ReNewAddress(CString str)
{
	CString str1 = str;
	for (int i = 0; i < str1.GetLength(); i++)
	{
		int unicode = (int)str1.GetAt(i);
		if (unicode == 92)
		{
			str1.Insert(i,_T("\\"));
			i++;
		}
	}
	return str1;
}

void CXF_FLT_OMPDlg::DeleteCh(char* ch1, int ch1count, int delCount, int pos)
{
	for (int i = pos; i<ch1count-delCount; i++)
	{
		ch1[i] = ch1[i + delCount];
	}
}

void CXF_FLT_OMPDlg::InesrtCh(char* ch1, int ch1count, char* ch2, int insCount, int pos)
{
	for (int i = ch1count; i>=pos; i-- )
	{
		ch1[i+insCount] = ch1[i];
	}
	for (int i = pos,k = 0; k<insCount;i++,k++)
	{
		ch1[i] = ch2[k];
	}
}

void CXF_FLT_OMPDlg::CreateHdr()
{
	CFile file_1(_T(""+s_Hdr2+""),CFile::modeCreate);
	file_1.Close();

	fstream file1(s_Hdr);
	char ch_box[1000]="";
	int countbox=0;
	char ch;
	while(file1.get(ch))
	{
		ch_box[countbox]=ch;
		//	cout.put(ch);
		countbox++;
	}
	file1.close();

	//修改ch_box  第一行 ncols
	//删除原ncols
	int countbox1 = 0, position1 = 0;
	for (int j=0; j<countbox; j++)
	{
		if (ch_box[j]==32) 
		{	
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]<=57 && ch_box[k]>=48)
				{
					position1 = k;
					break;
				}
			}
			for (int k=j;;k++)
			{
				if (ch_box[k]==10)
				{
					countbox1 = k - position1;
					break;
				}
			}
			break;
		}
	}
	DeleteCh(ch_box,countbox,countbox1,position1);
	countbox = countbox - countbox1;
	//添加新ncols
	char *buffer=new char[30]; 
	sprintf(buffer,"%d ",Rw);
	int insertCount=0;
	for (int i = 0; buffer[i]<=57 && buffer[i]>=48;i++ )
	{
		insertCount = i + 1;
	}
	InesrtCh(ch_box,countbox,buffer,insertCount,position1);
	countbox = countbox + insertCount;

	//修改ch_box  第2行 nrows
	//删除原nrows
	int count10 = 0;
	int start=0;
	for (int j=0; j<countbox; j++)
	{
		if (ch_box[j]==10) 
		{	
			start = j;
			break;
		}
	}
	for (int j=start; j<countbox; j++)
	{
		if (ch_box[j]==32) 
		{	
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]<=57 && ch_box[k]>=48)
				{
					position1 = k;
					break;
				}
			}
			for (int k=j;;k++)
			{
				if (ch_box[k]==10)
				{
					countbox1 = k - position1;
					break;
				}
			}
			break;
		}
	}
	DeleteCh(ch_box,countbox,countbox1,position1);
	countbox = countbox - countbox1;
	//添加新nrows
	sprintf(buffer,"%d ",Rh);
	for (int i = 0; buffer[i]<=57 && buffer[i]>=48;i++ )
	{
		insertCount = i + 1;
	}
	InesrtCh(ch_box,countbox,buffer,insertCount,position1);
	countbox = countbox + insertCount;

	//修改ch_box  第3行 x
	//删除原x
	count10=0;
	for (int j=0; j<countbox; j++)
	{
		if (ch_box[j]==10) 
		{
			count10++;
			if (count10 == 2)
			{
				start = j;
				break;
			}	
		}	
	}
	for (int j=start; j<countbox; j++)
	{
		if (ch_box[j]==32) 
		{	
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]<=57 && ch_box[k]>=48)
				{
					position1 = k;
					break;
				}
			}
			for (int k=j;;k++)
			{
				if (ch_box[k]==10)
				{
					countbox1 = k - position1;
					break;
				}
			}
			break;
		}
	}
	DeleteCh(ch_box,countbox,countbox1,position1);
	countbox = countbox - countbox1;
	//添加新x
	double dRx;
	dRx = Rx + m_Xshift - m_RXoffset;
	sprintf(buffer,"%lf ",dRx);
	for (int i = 0; buffer[i]<=57 && buffer[i]>=48 || buffer[i]==45;i++ )
	{
		insertCount = i + 6;
	}
	InesrtCh(ch_box,countbox,buffer,insertCount,position1);
	countbox = countbox + insertCount;

	//修改ch_box  第4行 y
	//删除原y
	count10=0;
	for (int j=0; j<countbox; j++)
	{
		if (ch_box[j]==10) 
		{
			count10++;
			if (count10 == 3)
			{
				start = j;
				break;
			}	
		}	
	}
	for (int j=start; j<countbox; j++)
	{
		if (ch_box[j]==32) 
		{	
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]<=57 && ch_box[k]>=48)
				{
					position1 = k;
					break;
				}
			}
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]==10)
				{
					countbox1 = k - position1;
					break;
				}
			}
			break;
		}
	}
	DeleteCh(ch_box,countbox,countbox1,position1);
	countbox = countbox - countbox1;
	//添加新y
	double dRy;
    dRy = Ry - m_RYoffset;
	//	double dRy = Ry;
	sprintf(buffer,"%f ",dRy);
	for (int i = 0; buffer[i]<=57 && buffer[i]>=48 || buffer[i]==45;i++ )
	{
		insertCount = i + 6;
	}
	InesrtCh(ch_box,countbox,buffer,insertCount,position1);
	countbox = countbox + insertCount;


	//修改ch_box  第5行 cell
	//删除原cell
	count10=0;
	for (int j=0; j<countbox; j++)
	{
		if (ch_box[j]==10) 
		{
			count10++;
			if (count10 == 4)
			{
				start = j;
				break;
			}	
		}	
	}
	for (int j=start; j<countbox; j++)
	{
		if (ch_box[j]==32) 
		{	
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]<=57 && ch_box[k]>=48)
				{				 
					position1 = k;
					break;
				}
			}
			for (int k=j;k<countbox;k++)
			{
				if (ch_box[k]==10)
				{
					countbox1 = k - position1;
					break;
				}
			}
			break;
		}
	}
	DeleteCh(ch_box,countbox,countbox1,position1);
	countbox = countbox - countbox1;
	//添加新cell
	sprintf(buffer,"%f ",Rcell);
	for (int i = 0; buffer[i]<=57 && buffer[i]>=48;i++ )
	{
		insertCount = i + 6;
	}
	InesrtCh(ch_box,countbox,buffer,insertCount,position1);
	countbox = countbox + insertCount;

	fstream file5(s_Hdr2);
	file5.write(ch_box,countbox);
	file5.close();
}

bool CXF_FLT_OMPDlg::CreateFlt()
{

	UpdateData();
	s_Flt = ReNewAddress(m_Flt);
	s_Hdr = s_Flt;
	s_Hdr.Delete(s_Hdr.GetLength()-3,3);
	s_Hdr.Insert(s_Hdr.GetLength(),_T("hdr"));
	s_Flt2 = ReNewAddress(m_Path);
	if (m_Name == _T(""))
	{
		MessageBox(_T("未输入新文件名"));
		return 0;
	}
	s_Flt2.Insert(s_Flt2.GetLength(),_T("\\\\"));
	s_Flt2.Insert(s_Flt2.GetLength(),_T(""+m_Name+""));
	s_Hdr2 = s_Flt2;
	s_Flt2.Insert(s_Flt2.GetLength(),_T(".flt"));
	s_Hdr2.Insert(s_Hdr2.GetLength(),_T(".hdr"));

	//读取HDR信息
	FILE* fp = _wfopen(s_Hdr,L"rt");
	if (fp==NULL)
	{
		MessageBox(_T("错误的FLT文件！"));
		return 0;
	}
	char buf[_MAX_PATH];
	fscanf(fp,"%s %d",buf,&w);
	fscanf(fp,"%s %d",buf,&h);
	fscanf(fp,"%s %lf",buf,&x);
	fscanf(fp,"%s %lf",buf,&y);
	fscanf(fp,"%s %lf",buf,&cell);
	fclose(fp);

//	x -= m_Xoffset;
//	y -= m_Yoffset;
	//计算原矩阵 顶点坐标
	point Vertex[4];
	for (int i=0; i<2; i++)
	{
		for (int j=0; j<2; j++)
		{
			Vertex[i*2+j].x = x + (w-1) * cell *j;
			Vertex[i*2+j].y = y + (h-1) * cell *i;
		}
	}

	//计算变换后矩阵 顶点坐标
	point RVertex[4];
	for (int i=0; i<2; i++)
	{
		for (int j=0; j<2; j++)
		{
			//RVertex[i*2+j].x =Vertex[i*2+j].x*cos(m_Rotate/180*PI)-Vertex[i*2+j].y*sin(m_Rotate/180*PI) + m_Xoffset;
			//RVertex[i*2+j].y =Vertex[i*2+j].x*sin(m_Rotate/180*PI)+Vertex[i*2+j].y*cos(m_Rotate/180*PI) + m_Yoffset;
			RVertex[i*2+j].x =Vertex[i*2+j].x*cos(-m_Rotate/180*PI)-Vertex[i*2+j].y*sin(-m_Rotate/180*PI);
			RVertex[i*2+j].y =Vertex[i*2+j].x*sin(-m_Rotate/180*PI)+Vertex[i*2+j].y*cos(-m_Rotate/180*PI);

		}
	}

	//计算变换后 Rw,Rh  Rx,Ry Rcell
	Rcell = (abs(cos(-m_Rotate/180*PI)) + abs(sin(-m_Rotate/180*PI))) * cell;
	float xmin = RVertex[1].x,ymin = RVertex[1].y;
	float xmax = RVertex[1].x,ymax = RVertex[1].y;
	for (int i=0; i<2; i++)
	{
		for (int j=0; j<2; j++)
		{
			if (RVertex[i*2+j].x< xmin)
			{
				xmin = RVertex[i*2+j].x;
			}
			if (RVertex[i*2+j].y<ymin)
			{
				ymin = RVertex[i*2+j].y;
			}
			if (RVertex[i*2+j].x> xmax)
			{
				xmax = RVertex[i*2+j].x;
			}
			if (RVertex[i*2+j].y>ymax)
			{
				ymax = RVertex[i*2+j].y;
			}
		}
	}
	Rw = (xmax - xmin) / Rcell + 1;
	Rh = (ymax - ymin) / Rcell + 1;
	Rx = xmin;
	Ry = ymin;
	CreateHdr();

// 	//重新计算 简化坐标 令 X = Y = 0
// 	x = y =0;
// 	for (int i=0; i<2; i++)
// 	{
// 		for (int j=0; j<2; j++)
// 		{
// 			Vertex[i*2+j].x = x + (w-1) * cell *j;
// 			Vertex[i*2+j].y = y + (h-1) * cell *i;
// 		}
// 	}
// 
// 	//计算变换后矩阵 顶点坐标
// 	for (int i=0; i<2; i++)
// 	{
// 		for (int j=0; j<2; j++)
// 		{
// 			RVertex[i*2+j].x =Vertex[i*2+j].x*cos(-m_Rotate/180*PI)-Vertex[i*2+j].y*sin(-m_Rotate/180*PI);
// 			RVertex[i*2+j].y =Vertex[i*2+j].x*sin(-m_Rotate/180*PI)+Vertex[i*2+j].y*cos(-m_Rotate/180*PI);
// 		}
// 	}
// 
// 	//计算变换后 
// 	for (int i=0; i<2; i++)
// 	{
// 		for (int j=0; j<2; j++)
// 		{
// 			if (RVertex[i*2+j].x< xmin)
// 			{
// 				xmin = RVertex[i*2+j].x;
// 			}
// 			if (RVertex[i*2+j].y<ymin)
// 			{
// 				ymin = RVertex[i*2+j].y;
// 			}
// 			if (RVertex[i*2+j].x> xmax)
// 			{
// 				xmax = RVertex[i*2+j].x;
// 			}
// 			if (RVertex[i*2+j].y>ymax)
// 			{
// 				ymax = RVertex[i*2+j].y;
// 			}
// 		}
// 	}
// 	Rx = xmin;
// 	Ry = ymin;

	CFile file_2(_T(""+s_Flt2+""),CFile::modeCreate);
	file_2.Close();
	fp2 = _wfopen(s_Flt2,L"wb");
	if (fp2 == NULL)
	{
		MessageBox(_T("未输入新文件路径或名字"));
		return 0;
	}
	float* z = new float[Rw];
	for (int i = 0; i < Rw ; i++)
	{
		z[i] = -9999;
	}
	for (int i = 0; i < Rh ; i++)
	{
		fwrite(z,1,sizeof(float)*Rw,fp2);
	}
	fclose(fp2);
//	MessageBox(_T("OK"));
	return 1;

}

//开始
void CXF_FLT_OMPDlg::OnBnClickedButton4()
{
	// TODO: 在此添加控件通知处理程序代码

	UpdateData();
	ThreadCount = m_ThreadCount;
	if (m_ThreadCount<1 || m_ThreadCount>40)
	{
		MessageBox(_T("线程输入错误，请输入1--40的整数！"));
		return;
	}
	//return;
	s_Offsets = ReNewAddress(m_Offsets);
	FILE* fp  = _wfopen(s_Offsets, L"rt");
	if (fp == NULL)
	{
		MessageBox(_T("错误的Offsets文件！"));
		return;
	}
	char buf[_MAX_PATH];
	fscanf(fp,"%s",buf);
	fscanf(fp,"%s %lf",buf,&m_Xoffset);
	fscanf(fp,"%s %lf",buf,&m_Yoffset);
	fscanf(fp,"%s %lf",buf,&m_RXoffset);
	fscanf(fp,"%s %lf",buf,&m_RYoffset);
	fscanf(fp,"%s %lf",buf,&m_Xshift);
	fscanf(fp,"%s %lf",buf,&m_Rotate);
	fclose(fp);

	if (m_Rotate==0 || int(m_Rotate)%45 == 0 && m_Rotate == int(m_Rotate))
	{
		Is45n = true;
	}
	else Is45n = false;
	//return ;
	if (!CreateFlt())
	{
		return;
	}
	GetDlgItem(IDC_BUTTON4)->EnableWindow(FALSE);
	UpdateData();
	ThreadCount = m_ThreadCount;
	Bw1 = Rw;
	Bh1 = Rh / ThreadCount;
//	Bh1 = Bh1;
 //	Bw2 = Rw - Bw1*k;
	Progress1 = 100 / ThreadCount;
    Progress2 = 100 - Progress1 * ThreadCount;
	
 	Bh2 = Rh - Bh1*ThreadCount;
	if (Bh2 == 0)
	{
		if (Progress2!=0)
		{
			m_CProgressCtrl1.SetPos(Progress2);
		}
		else
		{
			m_CProgressCtrl1.SetPos(0);
		}
		iFinishThread = 1;
	}
	else
	{ 
		iFinishThread = 0;
		m_CProgressCtrl1.SetPos(0);
	}
	for (int i = 0; i <= iPro; i++)
	{
		ProCount[i] = 0;
	}
	
	int i = 1;
	CWinThread *pWrite1 = AfxBeginThread(Write1,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	CWinThread *pWriteLast= AfxBeginThread(WriteLast,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
 	CWinThread *pWrite2 = AfxBeginThread(Write2,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
 	CWinThread *pWrite3 = AfxBeginThread(Write3,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
 	CWinThread *pWrite4 = AfxBeginThread(Write4,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
 	CWinThread *pWrite5 = AfxBeginThread(Write5,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
 	CWinThread *pWrite6 = AfxBeginThread(Write6,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
 	CWinThread *pWrite7 = AfxBeginThread(Write7,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite8 = AfxBeginThread(Write8,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite9 = AfxBeginThread(Write9,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite10 = AfxBeginThread(Write10,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite11 = AfxBeginThread(Write11,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite12 = AfxBeginThread(Write12,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite13 = AfxBeginThread(Write13,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite14 = AfxBeginThread(Write14,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite15 = AfxBeginThread(Write15,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite16 = AfxBeginThread(Write16,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite17 = AfxBeginThread(Write17,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite18 = AfxBeginThread(Write18,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite19 = AfxBeginThread(Write19,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite20 = AfxBeginThread(Write20,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite21 = AfxBeginThread(Write21,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite22 = AfxBeginThread(Write22,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite23 = AfxBeginThread(Write23,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite24 = AfxBeginThread(Write24,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite25 = AfxBeginThread(Write25,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite26 = AfxBeginThread(Write26,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite27 = AfxBeginThread(Write27,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite28 = AfxBeginThread(Write28,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite29 = AfxBeginThread(Write29,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite30 = AfxBeginThread(Write30,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite31 = AfxBeginThread(Write31,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite32 = AfxBeginThread(Write32,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite33 = AfxBeginThread(Write33,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite34 = AfxBeginThread(Write34,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite35 = AfxBeginThread(Write35,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite36 = AfxBeginThread(Write36,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite37 = AfxBeginThread(Write37,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL); 
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite38 = AfxBeginThread(Write38,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite39 = AfxBeginThread(Write39,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	if (i == ThreadCount) return; i++;
	CWinThread *pWrite40 = AfxBeginThread(Write40,&m_CProgressCtrl1,THREAD_PRIORITY_NORMAL, 0, NULL);
	
}

UINT Write1(LPVOID lpParam)
{
	int iThread = 1;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write2(LPVOID lpParam)
{
	int iThread = 2;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write3(LPVOID lpParam)
{
	int iThread = 3;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write4(LPVOID lpParam)
{
	int iThread = 4;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write5(LPVOID lpParam)
{
	int iThread = 5;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write6(LPVOID lpParam)
{
	int iThread = 6;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write7(LPVOID lpParam)
{
	int iThread = 7;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write8(LPVOID lpParam)
{
	int iThread = 8;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write9(LPVOID lpParam)
{
	int iThread = 9;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write10(LPVOID lpParam)
{
	int iThread = 10;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write11(LPVOID lpParam)
{
	int iThread = 11;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write12(LPVOID lpParam)
{
	int iThread = 12;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write13(LPVOID lpParam)
{
	int iThread = 13;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write14(LPVOID lpParam)
{
	int iThread = 14;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write15(LPVOID lpParam)
{
	int iThread = 15;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write16(LPVOID lpParam)
{
	int iThread = 16;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write17(LPVOID lpParam)
{
	int iThread = 17;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write18(LPVOID lpParam)
{
	int iThread = 18;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write19(LPVOID lpParam)
{
	int iThread = 19;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write20(LPVOID lpParam)
{
	int iThread = 20;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write21(LPVOID lpParam)
{
	int iThread = 21;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write22(LPVOID lpParam)
{
	int iThread = 22;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write23(LPVOID lpParam)
{
	int iThread = 23;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write24(LPVOID lpParam)
{
	int iThread = 24;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write25(LPVOID lpParam)
{
	int iThread = 25;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write26(LPVOID lpParam)
{
	int iThread = 26;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write27(LPVOID lpParam)
{
	int iThread = 27;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write28(LPVOID lpParam)
{
	int iThread = 28;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write29(LPVOID lpParam)
{
	int iThread = 29;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write30(LPVOID lpParam)
{
	int iThread = 30;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write31(LPVOID lpParam)
{
	int iThread = 31;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write32(LPVOID lpParam)
{
	int iThread = 32;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write33(LPVOID lpParam)
{
	int iThread = 33;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write34(LPVOID lpParam)
{
	int iThread = 34;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write35(LPVOID lpParam)
{
	int iThread = 35;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write36(LPVOID lpParam)
{
	int iThread = 36;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write37(LPVOID lpParam)
{
	int iThread = 37;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write38(LPVOID lpParam)
{
	int iThread = 38;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write39(LPVOID lpParam)
{
	int iThread = 39;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT Write40(LPVOID lpParam)
{
	int iThread = 40;
	if (iThread > ThreadCount)  return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY; //对应源数据 左下角零点的位置
	int numX, numY; //对应源数据 左下角零点的位置
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;       //进度条

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;
	for (int row = Bh1 -1 ; row >= 0  ; row--)
	{
		yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
			offX += xplus;
			offY += yplus;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}

			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*((Bh1 * iThread - row -1)*Rw),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress1 * (Bh1 - row) / Bh1 > ProCount[iThread])
		{
			ProCount[iThread] = Progress1 * (Bh1 - row) / Bh1;
			critical_section_progress.Lock();
			int  nTemp = CPC->GetPos();
			nTemp++;
			CPC->SetPos(nTemp);
			critical_section_progress.Unlock();
		}

	}
	delete z2;
	delete Fxy;
	delete z;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

UINT WriteLast(LPVOID lpParam)
{
	if (Bh2 == 0) return 0;
	FILE* fp0 = _wfopen(s_Flt,L"rb");	
	float* z = new float[1];
	double xformer, yformer, yafter;//, xafter;  //变换前后坐标
	double offX, offY;
	int numX, numY;
	double Bili_X,Bili_Y;  //相对位置 0-1

	float* Fxy = new float[4];
	float* z2 = new float[Bw1];
	CProgressCtrl* CPC = (CProgressCtrl*)lpParam;

	double cosR = cos(-m_Rotate/180*PI);
	double sinR = sin(-m_Rotate/180*PI);
	double xplus = Rcell * cosR / cell;
	double yplus = Rcell * (-sinR) / cell;

	for (int row = Bh2 -1 ; row >= 0  ; row--)
	{
	//	yafter = Ry + (row + Rh - Bh1*iThread) * Rcell;
		yafter = Ry + row * Rcell;
		xformer = cosR * (Rx) + sinR * (yafter);	
		yformer = -sinR * (Rx) + cosR * (yafter);
		offX = (xformer - x) / cell -xplus;
		offY = (yformer - y) / cell -yplus;
		for (int col = 0; col < Bw1  ; col++)
		{
// 			xafter = Rx + col * Rcell;
// 			xformer = cos(-m_Rotate/180.0*PI) * (xafter) + sin(-m_Rotate/180.0*PI) * (yafter);
// 			yformer = -sin(-m_Rotate/180.0*PI) * (xafter) + cos(-m_Rotate/180.0*PI) * (yafter);
// 
// 			offX = (xformer - x) / cell;
// 			offY = (yformer - y) / cell;
			offX += xplus;
			offY += yplus;
			//numX = numY =0;
			if ( offX >=0 && offX <= w-1 && offY>=0 && offY<=h-1)
			{
				if (Is45n)
				{
					numX = (int)(offX+0.5);
					numY = (int)(offY+0.5);
					fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
					fread(z,1,sizeof(float),fp0);
					z2[col] = z[0];
				} 
				else
				{
					numX = (int)offX;
					numY = (int)offY;
					Bili_X = offX - numX;
					Bili_Y = offY - numY;
					if (Bili_X == 0 && Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						z2[col] = Fxy[0];
						continue;
					}
					if (Bili_X == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*1,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*1,fp0);
						Fxy[1]=Fxy[3]=0;
					}
					else if (Bili_Y == 0)
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						Fxy[2]=Fxy[3]=0;
					}
					else
					{
						fseek(fp0,sizeof(float)*((h - numY -1) * w + numX),SEEK_SET);
						fread(Fxy,1,sizeof(float)*2,fp0);
						fseek(fp0,sizeof(float)*((h - numY -1 - 1) * w + numX),SEEK_SET);
						fread(Fxy+2,1,sizeof(float)*2,fp0);
					}
					if (Fxy[0]!=-9999 && Fxy[1]!=-9999 && Fxy[2]!=-9999 &&Fxy[3]!=-9999)
					{
						z2[col] = Fxy[0] * (1 - Bili_X) *(1 - Bili_Y)
							+ Fxy[1] * (Bili_X) *(1 - Bili_Y)
							+ Fxy[2] * (1 - Bili_X) *(Bili_Y)
							+ Fxy[3] * (Bili_X) *(Bili_Y);
					}
					else z2[col] = -9999;
				}
			}
			else 
			{
				z2[col] = -9999;
			}
		}
		critical_section.Lock();
		fp2 = _wfopen(s_Flt2,L"rb+");
		fseek(fp2,sizeof(float)*(((Rh - row -1)*Rw)),SEEK_SET);
		fwrite(z2,sizeof(float),Bw1,fp2);
		fclose(fp2);
		critical_section.Unlock();

		if (Progress2 * (Bh2 - row) / Bh2 > ProCount[0])
		{
			if (Bh2 == 1)
			{
				critical_section_progress.Lock();
				int  nTemp = CPC->GetPos();
				nTemp += Progress2;
				CPC->SetPos(nTemp);
				critical_section_progress.Unlock();
			}
			else if (Progress2>Bh2)
			{
				int p = Progress2 / (Bh2 -1);
				int q = Progress2 - p * (Bh2-1);
				ProCount[0]++;
				critical_section_progress.Lock();
				int  nTemp = CPC->GetPos();
				if (q == 0)
				{
					nTemp += p;
				}
				else
				{
					if (ProCount[0] != Bh2)
					{
						nTemp += p;
					}
					else nTemp += q;
				}
				CPC->SetPos(nTemp);
				critical_section_progress.Unlock();
			}
			else{
				ProCount[0] = Progress2 * (Bh2 - row) / Bh2;
				critical_section_progress.Lock();
				int  nTemp = CPC->GetPos();
				nTemp++;
				CPC->SetPos(nTemp);
				critical_section_progress.Unlock();
			}
		}
		//		int nTemp = 100 * (Bh1 - row) / Bh1 ;
		//		CPC->SetPos(nTemp);
	}
	delete z;
	delete Fxy;
	delete z2;
	fclose(fp0);
	critical_section_finish.Lock();
	iFinishThread++;
	if (iFinishThread == ( ThreadCount +1 ))
	{
		CXF_FLT_OMPDlg *pMMM = (CXF_FLT_OMPDlg*)AfxGetMainWnd();
		pMMM->SetTimer(1,300,NULL);
	}
	critical_section_finish.Unlock();
	return 0;
}

void CXF_FLT_OMPDlg::OnTimer(UINT_PTR nIDEvent)
{
	// TODO: 在此添加消息处理程序代码和/或调用默认值
	if (nIDEvent == 1)
	{
		KillTimer(1);
		GetDlgItem(IDC_BUTTON4)->EnableWindow(TRUE);
		MessageBox(_T("完成"));	
	}
	CDialogEx::OnTimer(nIDEvent);
}