proj.h

一个VC写的直接读取MIF文件的GIS系统

#define    PAI          3.1415926535897932
#define    EARTH_A   	6378245.000000  
#define    EARTH_E12    0.006693421623  

class CMapProj : public CObject
{
private:
	double arctg( double dy, double dx);
	double atanhhh(double x);
	double q (double bb);
	double invq(double q0);
	double r(double b);
	double m(double bb);

public:
	void  rhumb_line(double bb1, double ll1, double bb2, double ll2, double *angle, double* length);//求两经纬度点的方位、距离。（单位为"弧度"和"米"）
	void  rhumb_line1(double bb1, double ll1, double bb2, double ll2, double *angle, double* length);//单位"度""和"米"
	void  inverse_rhumb(double b1,double l1,double angl,double len, double* b2,double* l2);//已知一点和方位、距离，求另一点。（单位为"弧度"和"米"）
	void  inverse_rhumb1(double b1,double l1,double angl,double len, double* b2,double* l2);//单位"度"和"米"
};