大地测量 平面坐标转换

实现大地测量课程中，两个平面坐标之间的相互转换的c程序

#include 
#include 

#define PI 3.1415926

typedef struct Data
{
	double x;
	double y;
} d1[4]d2[4];

void ReadDat（Data d1[4]Data d2[4]）
{
	
	FILE *fp1*fp2;
	fp1=fopen（“old.txt““r“）;	
	fp2=fopen（“new.txt““r“）;
	double txtypxpy;
	int cnt1=0cnt2=0;
	do
	{
		fscanf（fp1“%lf%lf“&tx&ty）;
		d1[cnt1].x=tx;
		d1[cnt1].y=ty;
		cnt1++;
	}while（cnt1<4）;

	do
	{
		fscanf（fp2“%lf%lf“&px&py）;
		d2[cnt2].x=px;
		d2[cnt2].y=py;
		cnt2++;
	}while（cnt2<4）;	
	fclose（fp1）;
	fclose（fp2）;	
}

//矩阵转置
template
void Transpose（T1 *mat1T2 *mat2int aint b）
{  
	int ij;
	for（i=0;i		for（j=0;j			mat2[j][i]=mat1[i][j];
		return;
}
//矩阵相乘
template
void Multiply（T1 *mat1T2 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    .CA....     41984  2011-11-21 00:04  平面坐标转换\Debug\vc60.idb

    .CA....     53248  2011-11-21 00:01  平面坐标转换\Debug\vc60.pdb

    .CA....    188526  2011-11-21 00:01  平面坐标转换\Debug\坐标转换.exe

    .CA....    199136  2011-11-21 00:01  平面坐标转换\Debug\坐标转换.ilk

    .CA....     11756  2011-11-21 00:04  平面坐标转换\Debug\坐标转换.obj

    .CA....    221452  2011-11-21 00:01  平面坐标转换\Debug\坐标转换.pch

    .CA....    484352  2011-11-21 00:01  平面坐标转换\Debug\坐标转换.pdb

    .CA....       126  2011-11-20 19:11  平面坐标转换\new.txt

    .CA....       134  2011-11-20 19:11  平面坐标转换\old.txt

    .CA....      2805  2011-11-21 00:01  平面坐标转换\坐标转换.cpp

    .CA....      3427  2011-11-20 19:11  平面坐标转换\坐标转换.dsp

    .CA....       524  2011-11-20 19:11  平面坐标转换\坐标转换.dsw

    .CA....     50176  2011-11-24 15:59  平面坐标转换\坐标转换.ncb

    .CA....     48640  2011-11-24 15:59  平面坐标转换\坐标转换.opt

    .CA....       757  2011-11-21 00:04  平面坐标转换\坐标转换.plg

     目录          0  2011-11-21 00:09  平面坐标转换\Debug

     目录          0  2011-11-24 15:59  平面坐标转换

----------- ---------  ---------- -----  ----

              1307043                    17