地图学实践心得体会.doc

131309040441) .32) .33) .34) .55) .66) .7mapinfomapinfomapinfo1word3602mapinfo3sshiftssshiftshift4mapinfo5 620122008-137-1windows2e130904044e3lne4mapinfomapinfo-mapinfo,mapinfo567mapinfo89drawsshiftshiftshift. 12 120904009 20136 : mapinfomapinfomapinfo 1234567 mapinfo 89 1mapinfo mapinfomapinfomapinfo2322004ssnap56789 1 2 1231245 1 2 3 4 1. 2. 1231234541234.11-1 11016xxx123123171742004200451234.5. 2004123451mapinfo professionalraster imagesdjt.jpg1-11-1 1-21-2 1-31-3 1-41-4 1-51-5 32mapinfo professional1-61-6 1-71-7 1-81-8 1-91-9 mapinfo17 tg(450?lgr1?lgr2a2)?(u?e,sin?esin?,e?0lgu2?lgu1btg(45?2)sin?0?kcm100r1u1?m0?a2?b2)a2100r2u2?m0?ku?x?s?cos?y?sin?m?n?kru?p?m?n?00:; m0:?n?43000?,?s?36000?e?120000?,?w?112000?0=11600001=3730,2=410015000=1 01a(1?e2)aa2?b22m?,n?,e?223/2221/2(1?esin?)(1?esin?)a2iag-1975a=6378140mb=6356755m2u1,u2;k0xy3. #include <stdio.h>#include <iostream>#include <cmath>using namespace std;const double pi=3.1415926;const double a1=6378140;const double b1=6356755;double get_e()return(sqrt(a1*a1-b1*b1)/(a1*a1);double get_n(double m)return(a/sqrt(1-get_e()*get_e()*sin(m*pi/180)*sin(m*pi/180);double get_r(double m) return(get_n(m)*cos(m*pi/180);double get_u(double m) /ureturn(tan(45*pi/180+m*pi/360)/pow(1-e*sin(m*pi/360)/(1+e*sin(m*pi/360),e/2);double get_a(double m1, double m2) /adouble r1,r2;double u1,u2;r1=log10(get_r(m1);r2=log10(get_r(m2);u1=log10(get_u(m1);u2=log10(get_u(m2);/*double m;m=log10(r1);*/return(r1-r2)/(u2-u1);double get_c(double m1,double m2,) /cdouble u=pow(get_u(m1),get_a(m1,m2);double no=get_r(m1)*u/get_a(m1,m2);return( no);double get_weidu(double m1, double m2)double wd=asin(get_a( m1, m2);return(wd*180/pi);double get_jiajiao(double m1, double m2,double y)return(get_a(m1,m2)*y);double get_p(double m1,double m2,double x)return(get_c(m1,m2)/pow(get_u(x),get_a(m1,m2);double get_x(double ps,double m1, double m2,double x,double y)return(ps-get_p(m1,m2,x,m)*cos(get_jiajiao(m1,m2,y);double get_y(double m1, double m2,int m,double x,double y)return(get_p(m1,m2,x)*sin(get_jiajiao(m1,m2,y);/*/ /*/ /*/ int main()double m1=27.00; /double m2=45.00; /double ps=6378940.341; /cout<<<<endl;cout<<u1=<<get_u(m1)<<endl;cout<<u2=<<get_u(m2)<<endl;cout<<a<<get_a(m1,m2)<<endl;cout<<c<<get_c(m1,m2)<<endl;cout<<p<< get_pm1,m2,x <<endl;cout<<*<<endl; cout<<36<=x<=43,72<=y<=138<<endl; double x14=2,6,10,14,18,22,26,30,34,38,42,46,50,54;double y12=72,78,84,90,96,102,108,114,120,126,132,138;cout<<<< <<<< <<x<< <<y<<endl;/x,yfor(int i=0;i<14;i+)for(int j=0;j<12;j+)cout<<xi<< <<yj<< <<get_x(ps,m1,m2,xi,yj)<< <<get_y(m1,m2,xi,yj)<<endl;cout&