高斯投影正算
高斯投影正、反算代码
//高斯投影正、反算
//////6度带宽 54年北京坐标系
//高斯投影由经纬度(Unit:DD)反算大地坐标(含带号，Unit:Metres)
void GaussProjCal(double longitude, double latitude, double *X, double *Y)
{
int ProjNo=0; int ZoneWide; ////带宽
double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval;
double a,f, e2,ee, NN, T,C,A, M, iPI;
iPI = 0.0174532925199433; ////3.1415926535898/180.0;
ZoneWide = 6; ////6度带宽
a=6378245.0; f=1.0/298.3; //54年北京坐标系参数
////a=6378140.0; f=1/298.257; //80年西安坐标系参数
ProjNo = (int)(longitude / ZoneWide) ;
longitude0 = ProjNo * ZoneWide + ZoneWide / 2;
longitude0 = longitude0 * iPI ;
latitude0=0;
longitude1 = longitude * iPI ; //经度转换为弧度
latitude1 = latitude * iPI ; //纬度转换为弧度
e2=2*f-f*f;
ee=e2*(1.0-e2);
NN=a/sqrt(1.0-e2*sin(latitude1)*sin(latitude1));
T=tan(latitude1)*tan(latitude1);
C=ee*cos(latitude1)*cos(latitude1);
A=(longitude1-longitude0)*cos(latitude1);
M=a*((1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2 /32+45*e2*e2
*e2/1024)*sin(2*latitude1)
+(15*e2*e2/256+45*e2*e2*e2/1024)*sin(4*latitude1)-(35*e2*e2*e2/30 72)*sin(6*l
atitude1));
xval = NN*(A+(1-T+C)*A*A*A/6+(5-18*T+T*T+72*C-58*ee)*A*A*A*A*A/120);
yval = M+NN*tan(latitude1)*(A*A/2+(5-T+9*C+4*C*C)*A*A*A*A/24
+(61-58*T+T*T+600*C-330*ee)*A*A*A*A*A*A/720);
X0 = 1000000L*(ProjNo+1)+500000L;
Y0 = 0;
xval = xval+X0; yval = yval+Y0;
*X = xval;
*Y = yval;
}
//高斯投影由大地坐标(Unit:Metres)反算经纬度(Unit:DD)
void GaussProjInvCal(double X, double Y, double *longitude, double *latitude)
{
int ProjNo; int ZoneWide; ////带宽
double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval;
double e1,e2,f,a, ee, NN, T,C, M, D,R,u,fai, iPI;
iPI = 0.0174532925199433; ////3.1415926535898/180.0;
a = 6378245.0; f = 1.0/298.3; //54年北京坐标系参数
////a=6378140.0; f=1/298.257; //80年西安坐标系参数
ZoneWide = 6; ////6度带宽
ProjNo = (int)(X/1000000L) ; //查找带号
longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2;
longitude0 = longitude0 * iPI ; //中央经线
X0 = ProjNo*1000000L+500000L;
Y0 = 0;
xval = X-X0; yval = Y-Y0; //带内大地坐标
e2 = 2*f-f*f;
e1 = (1.0-sqrt(1-e2))/(1.0+sqrt(1-e2));
ee = e2/(1-e2);
M = yval;
u = M/(a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256));
fai = u+(3*e1/2-27*e1*e1*e1/32)*sin(2*u)+(21*e1*e1/16-55*e1*e1*e1*e1/32)*si n(
4*u)
+(151*e1*e1*e1/96)*sin(6*u)+(1097*e1*e1*e1*e1/512)*sin(8*u);
C = ee*cos(fai)*cos(fai);
T = tan(fai)*tan(fai);
NN = a/sqrt(1.0-e2*sin(fai)*sin(fai));
R = a*(1-e2)/sqrt((1-e2*sin(fai)*sin(fai))*(1-e2*sin(fai)*sin(fai))*(1-e2 *sin
(fai)*sin(fai)));
D = xval/NN;
//计算经度(Longitude) 纬度(Latitude)
longitude1 = longitude0+(D-(1+2*T+C)*D*D*D/6+(5-2*C+28*T-3*C*C+8*ee+24*T*T)*D
*D*D*D*D/120)/cos(fai);
latitude1 = fai -(NN*tan(fai)/R)*(D*D/2-(5+3*T+10*C-4*C*C-9*ee)*D*D*D*D/24
+(61+90*T+298*C+45*T*T-256*ee-3*C*C)*D*D*D*D*D*D/720);
//转换为度 DD
*longitude = longitude1 / iPI;
*latitude = latitude1 / iPI;
}。
高斯投影坐标正算名词解释
高斯投影坐标正算名词解释
高斯投影坐标正算，也称为高斯投影正算，是指已知椭球面上某点的大地坐标（B，L），通过一定的数学公式和投影变换，求该点在高斯投影平面上的直角坐标（X，Y）。
这是大地测量和地图制作中常用的坐标转换方法。
高斯投影正算的实现通常依赖于高斯-克吕格投影，其转换公式涉及到大地纬度B、经度L、地球椭球参数、投影常数等多个因素。
通过这一转换，可以将地理坐标（经纬度）转换为更易于进行地图制作和平面测量处理的投影坐标。
在实际应用中，需要充分考虑地球椭球模型、投影方式和精度要求等因素，以选择适当的数学模型和计算方法进行高斯投影坐标正算。
高斯投影正反算python代码
高斯投影正反算1. 什么是高斯投影高斯投影是一种常用的地图投影方法，它将地球表面的经纬度坐标转换为平面坐标，常用于地理信息系统（GIS）和测绘工程中。
高斯投影分为正算和反算两个过程。
•正算：将经纬度坐标转换为平面坐标。
•反算：将平面坐标转换为经纬度坐标。
2. 高斯投影正算2.1 原理高斯投影正算的原理是根据椭球体上某一点处的曲率半径、子午线弧长和东西方向上的距离，计算该点在平面上的x、y坐标。
2.2 具体步骤高斯投影正算的具体步骤如下：1.根据给定的椭球体参数（长半轴a、短半轴b），计算椭球体第一偏心率e。
2.根据给定的中央子午线经度λ0，计算λ - λ0 的差值Δλ。
3.计算曲率半径N和子午线弧长A0。
4.根据给定的纬度φ和经度λ，计算Δφ和Δλ。
5.计算子午线弧长A1、A2、A3和A4。
6.计算平面坐标x和y。
2.3 Python代码实现下面是使用Python实现高斯投影正算的示例代码：import math# 输入参数a = 6378137.0 # 长半轴b = 6356752.314245 # 短半轴e = math.sqrt(1 - (b/a)**2) # 第一偏心率λ0 = math.radians(120) # 中央子午线经度，单位为弧度# 输入经纬度坐标φ = math.radians(30) # 纬度，单位为弧度λ = math.radians(121) # 经度，单位为弧度# 计算ΔλΔλ = λ - λ0# 计算曲率半径N和子午线弧长A0N = a / math.sqrt(1 - e**2 * math.sin(φ)**2)A0 = a * (1 - e**2) / (1 - e**2 * math.sin(φ)**2)**1.5# 计算Δφ和ΔλΔφ = φ - φ0# 计算子午线弧长A1、A2、A3和A4A1 = A0 + N * math.tan(φ) / 2 * Δλ**2 * math.cos(φ)A2 = A0 + N * math.tan(φ) / 24 * (5 - math.tan(φ)**2 + 9 * e2 * math.cos(φ)**2 + 4 * e2**2 * math.cos(φ)**4) * Δλ**4 * math.cos(φ)A3 = A0 + N * math.tan(φ) / 720 * (61 - 58 * math.tan(φ)**2 + math.tan(φ)** 4) * Δλ**6 * math.cos(φ)A4 = A0 + N * math.tan(φ) / 40320 * (1385 - 3111*math.tan(φ)**2 + 543*math.tan(φ)**4 - math.tan(φ)**6) \* Δλ**8 \* math.cos(phi)# 计算平面坐标x和yx = A1 + A2 + A3 + A4y = N / math.cos(phi) \(Δλ - Δλ**3/6*(1+math.tan(phi))**2/N/A0^2 \+ Δλ^5/120*(5+28*math.tan(phi)^2+24*math.tan(phi)^4)*N/A0^4/N/A0^3)# 输出结果print("平面坐标(x, y):", x, y)3. 高斯投影反算3.1 原理高斯投影反算的原理是根据平面坐标和中央子午线经度，计算对应的经纬度坐标。
高斯投影坐标正算公式
高斯投影坐标正算公式高斯投影坐标正反算公式2.2.2. 1高斯投影坐标正算公式： B, x,y高斯投影必须满足以下三个条件：⑴中央子午线投影后为直线；⑵中央子午线投影后长度不变；⑶投影具有正形性质，即正形投影条件。
由第一条件知中央子午线东西两侧的投影必然对称于中央子午线，即式中，x为的偶函数，y为的奇函数；，即，如展开为的级数，收敛。
（2-10）式中是待定系数，它们都是纬度B的函数。
由第三个条件知：分别对和q求偏导数并代入上式(2-11)上两式两边相等，其必要充分条件是同次幂前的系数应相等，即(2-12)(2-12)是一种递推公式，只要确定了就可依次确定其余各系数。
由第二条件知:位于中央子午线上的点，投影后的纵坐标x应等于投影前从赤道量至该点的子午线弧长X，即(2-10)式第一式中，当时有：(2-13)顾及(对于中央子午线)得：(2-14,15)(2-16)依次求得并代入(2-10)式，得到高斯投影正算公式(2-17)2.2.2. 2高斯投影坐标反算公式x,y B,投影方程：(2-18)高斯投影坐标反算公式推导要复杂些。
⑴由x求底点纬度(垂足纬度),对应的有底点处的等量纬度，求x,y与的关系式，仿照式有，由于y和椭球半径相比较小(1/16.37)，可将展开为y的幂级数；又由于是对称投影，q必是y的偶函数，必是y的奇函数。
(2-19)是待定系数，它们都是x的函数.由第三条件知：，， (2-20)(2-19)式分别对x和y求偏导数并代入上式上式相等必要充分条件，是同次幂y前的系数相等，第二条件，当y=0时，点在中央子午线上，即x=X，对应的点称为底点，其纬度为底点纬度，也就是x=X时的子午线弧长所对应的纬度，设所对应的等量纬度为。
也就是在底点展开为y的幂级数。
由(2-19)1式依次求得其它各系数(2-21)(2-21)1…………将代入(2-19)1式得(2-22)1(2-22)将代入(2-19)2式得(2-23)2式。
高斯平面直角坐标与大地坐标的相互转换——高斯投影的正算与反算.
昆明冶金高等专科学校测绘学院 （4）计算公式
3 2 2 2 4 ( 5 3 t 9 t ) y f f f f 2M f N f 2 4M f N 3 f tf 2 4 6 (6 1 9 0t f 4 5t f ) y 7 2 0M f N 5 f 1 1 2 2 3 l y (1 2t f f ) y 3 N f co s B f 6 N f co s B f 1 2 5 (5 2 8t 2 t4 2 2 f 24 f 6 f 8 f t f )y 5 1 2 0N f co s B f B Bf tf y2 tf
式中：
2 e 2 cos2 B
t 2 tan2 B l （L L0） X为B对应子午线弧长 N为卯酉圈曲率半径 20626 5
昆明冶金高等专科学校测绘学院
2
高斯投影坐标反算公式
（1）高斯投影反算：
已知某点 x, y ，求该点 L, B ，即 x, y ( L, B) 的坐标变换。 （2）投影变换必须满足的条件
昆明冶金高等专科学校测绘学院
二、高斯投影坐标正反算得实用公式及算例
1 高斯投影坐标正算公式 （1）高斯投影正算： 已知某点的 L, B ，求该点的 x, y ，即 （2）投影变换必须满足的条件： 中央子午线投影后为直线； 中央子午线投影后长度不变； 投影具有正形性质，即正形投影条件。 （3）投影过程 在椭球面上有对称于中央子午线的两点 P1 和 P2 ,它们的大地坐标 分别为 ( L1 , B1 )或（l1 , B1）及 （L2 , B2）或(l2 , B2 ) 式中 l 为椭球面上点的经 度与中央子午线 ( L0 ) 的经度差：l L L0 ,点在中央子午线之东, l 为正，在西则为负，则投影后的平面坐标一定为P1 ( x1 , y1 ) 和 P2 ( x 2 , y 2 ) 。
高斯投影正反算
高斯投影正、反算及换带程序执行条件※数组投影选择T、换算点个数“Z=0 F≠0”、=0正算0、≠0反算※坐标系选择“54 ≠54”、=54换算为1954年北京坐标系输入54、≠54换算为1988年西安坐标系M、中央子午线经度（°′″）输入※大地坐标I、序列号B、L：大地纬度和经度（地理坐标）（°′″）※高斯平面坐标轴子午线I、序列号X、Y：高斯平面坐标（m） Z、轴子午线（°）输出※大地坐标子午收敛角N、序列号B、L：大地纬度和经度（地理坐标）（°′″） R、子午收敛角（°′″）※高斯平面坐标子午收敛角N、序列号X、Y：高斯平面坐标（m） R、子午收敛角（°′″）注：1、程序执行前必须进行数组定位。
如：Defm 10 T×2=5×2=102、Y坐标值要去掉带号及避免出现负值的500公里；4、本程序运算时，各已知数据、观测变量不会随之变化，可非常方便地进行各数据的核对;5、本程序在进行换带计算时采用的是间接换带计算法。
Prog GSXYDefm 10:TA“Z=0 F≠0”G“54 ≠54”Z:Fixm:I=0:「b」0:I=I+1◢J=2I-1:M=Z[J:L=Z[J+1:A=0=>Prog“3”:B=M:M=L+Z:Prog“3”:L=M:{BL}:M=B:Prog“2”: B=M:M=L:Prog“2”:L=M-Z:≠>X=M:Y=L:{XY}:B=X:L=Y⊿Z[J]=B:Z[J+1]=L:I<T=>Goto 0⊿G=54=>C=6399698.90178271:E=.006738525414684:≠>C=6399596.65198801:E=.006 739501819473⊿I=0:「b」0:I“N”=I+1◢J=2I-1:B=Z[J:L=Z[J+1:A≠0=>X=B:Y=L:Goto 2⊿S=sin B:G=54=>F=111134.8611B-(32 005.7799S+133.9238S∧3+.6973S∧5+.0039S∧7)cos B:≠>F=111133.0047B-(32009.857 S+133.9602S∧3+.6976S∧5+.0039S∧7)cos B⊿U=√Ecos B:V=√(1+U2:N=C÷V:W=tan B: M=cos B(Lπ÷180:X=F+NW(.5M2+1┛24(5-W2+9U2+4U∧4)M∧4+1┛720(61-58W2+W∧4)M∧6◢Y=N(M+1┛6(1-W 2+U 2)M ∧3+1┛120(5-18W 2+W ∧4+14U 2-58U 2W 2)M ∧5◢M=W ┛π(180M+60(1+3U 2+2U ∧4)M ∧3+12(2-W 2)M ∧5:Goto 3:「b 」2:W=E ﹣6X-3:G=54=>F=27.11115372595+9.024********W-.00579740442W 2-4.3532572E ﹣4W ∧3+4.857285E ﹣5W ∧4+2.15727E ﹣6W ∧5-1.9399E ﹣7W ∧6:≠>F=27.11162289465+9.024********W-.00579850656W2-4.3540029E ﹣4W ∧3+4.858357E ﹣5W ∧4+2.15769E ﹣6W ∧5-1.9404E ﹣7W ∧6⊿U=√Ecos F:V=√(1+U 2:Q=YV ÷C:W=tan F:M=F-(1+U 2)W ┛π(90Q 2-7.5(5+3W 2+U 2-9U 2W 2)Q ∧4+.25(61+90W 2+45W ∧4)Q ∧6:Prog “3”:B=M ◢M=Z+1┛(πcos F)(180Q-30(1+2W 2+U 2)Q ∧3+1.5(5+28W 2+24W ∧4)Q ∧5:Prog “3”:L=M ◢M=W ┛π(180Q-60(1+W 2-U 2)Q ∧3+12(2+5W 2+3W ∧4)Q ∧5:「b 」3:Prog “3”:R=M ◢ I<T=>Goto 1⊿“END ”概要说明：我国的经度范围西边自73°起，东边至135°，可分成6°带共11带或3°共22带。
高斯投影正反算原理
高斯投影正反算原理高斯投影是一种常用于地图制图的投影方式，也被广泛应用于其他领域的空间数据处理。
高斯投影正反算是对于已知的地球坐标系上的位置（经纬度），通过计算得到该点的平面坐标（东、北坐标），或者对于已知的平面坐标（东、北坐标），通过计算得到该点的地球坐标系上的位置（经纬度）的过程。
本文将详细介绍高斯投影正反算的原理。
一、高斯投影简介高斯投影是一种圆锥投影，其投影面在地球表面的某个经线上，也就是说，投影面是以该经线为轴的圆锥面。
经过对圆锥体的调整后，使其切于地球椭球面，在该经线上进行投影，同时保持沿该经线方向的比例尺一致，从而达到地图上各点在包括该经线的垂直面上映射的目的。
这种投影方式在某一特定区域内得到高精度的结果，因此广泛应用于地图制图。
二、高斯投影数学模型对于高斯投影正反算，需要先建立高斯投影坐标系与地球坐标系的转换模型。
1.高斯投影坐标系的建立高斯投影坐标系的建立需要确定圆锥面的基本参数，首先需要确定其所处的中央子午线，再确定该子午线上的经度为零点，并利用该经线上某一点的经度和该点的高度来确定该点所在的圆锥体。
圆锥体的底面包括所有与地球椭球面相切的圆面，通过对这些圆面进行调整，使得圆锥体转动后能够在中央子午线上进行投影。
在此基础上，可建立高斯投影坐标系，其中投影面为圆锥面，且中央子午线与投影面的交点称为该投影坐标系的中心，投影面的上端点和下端点分别对应正北方向和正南方向。
2.地球坐标系的建立地球坐标系是以地球椭球体为基础建立的，其坐标系原点确定为地球椭球体上的一个特定点。
在已知该点经纬度和高度的前提下，可确定以该点为中心的地球椭球体，并可根据它与地球坐标系之间的转换关系得到平面坐标系。
3.高斯投影坐标系与地球坐标系之间的转换关系由于高斯投影坐标系与地球坐标系存在不同的坐标体系和基准面，因此需要通过数学关系式来建立它们之间的转换关系。
（1）高斯投影坐标系转地球坐标系：已知高斯投影坐标系中任意一点的东北坐标（N，E），以及所属的中央子午线经度λ0、椭球参数a和e，则可通过以下公式求出该点的地球坐标系经纬度（φ，λ）和高度H：A0为以地球椭球体中心为原点，高斯投影坐标系中心投影坐标为（0，0）的点到椭球面的距离。
高斯投影坐标计算
l 4
4
)
2、高斯投影坐标反算公式
已知高斯平面坐标(x,y)，求椭球面上的大地坐标(B，L）的 问题称高斯投影坐标反算。 B 1 ( x, y ) 函数式： l 2 ( x, y ) 同正算一样，对投影函数提出三个条件 (1) x 坐标轴投影成中央子午线，是投影的对称轴； (2) x 轴上长度投影保持不变； (3) 正形投影条件。
B ) cos
B
N b2 b3 b 4 b5
N
f
y cos B f cos
2
( 0 . 5 0 . 00336975
B f ) sin B f cos B f 0 . 001123 cos
2 2 2 2
0 . 333333 ( 0 . 1666667
高斯投影的特点分析
(1)当l等于常数时，随着B的增加x值增大，y值 减小；无论B值为正或负，y值不变。这就是说， 椭球面上除中央子午线外，其他子午线投影后， 均向中央子午线弯曲，并向两极收敛，同时还 对称于中央子午线和赤道。 (2)当B等于常数时，随着l的增加，x值和y值都 增大。所以在椭球面上对称于赤道的纬圈，投 影后仍成为对称的曲线，同时与子午线的投影 曲线互相垂直凹向两极。
(1)中央子午线投影后为直线； (2)中央子午线投影后长度不变； (3)投影具有正形性质，即正形投影 条件。
高斯投影坐标正算
l =3／ρ=0.052
1) 由第一个条件（中央子午线投影后为直线） 可知，由于地球椭球体是一个旋转椭球体，即 中央子午线东西两侧的投影必然对称于中央子 午线。 x 为 l 的偶函数，而y 则为 l 的奇函数。
B f ) cos
2
2
高斯投影坐标计算
B
d B dq
2
dX dq dq
c
(
cos B dV V dB
2
dB dq
sin B dB V dq
2
)
2
d B dq
2
cos B c ( tan B V
2 2
3
V
sin B cos B
)
N sin B cos B
同理得
d X dq
3
N cos B ( 1
3
3
2
0
l
L
L
0
高斯投影坐标正算的函数式：
x y
l 是以弧度为单位的经度差。
F B , l F B , l
1 2
一 高斯投影坐标正算公式计算
如图，椭球面上一点投影 到平面后为d点，椭球面上 该点的平行圈（B或q为一 常数）与中央子午线的交 点为e点，若将上式中的展 开点z0设为e处，则很据高 斯投影条件，中央子午线 的长度比m=1，且纵坐标x 等于从赤道起到该平行圈 间的子午线弧长X。此时 可以写出下列方程：
4 2
二、高斯投影坐标反算公式
最后得到坐标反算的公式为：
B B
f
2M
f
t
f
y N
f
2
t 24 M
2 f
f
f
f
N
4 f
3 f
5 3 t
6
2 f
2 f
9 f t
2
2 f
y
4
t
第四章 7高斯投影坐标正反算
——正形投影的一般条件 ——高斯投影坐标正算 ——高斯投影坐标反算 ——高斯投影几何解释
提前在黑板上写出四个m2
上一讲应掌握的内容
1、地图(数学)投影：将椭球面上元素(包括坐标，方位和 距离)按一定的数学法则投影到可展平面上。 x F1 ( L, B) 坐标投影公式： y F2 ( L, B) 2、地图投影变形几个概念： 长度比，主方向，变形椭圆 3、四种投影变形： 长度变形，方向变形，角度变形，面积变形
x m0 m 2 l 2 m 4 l 4 y m1l m3 l 3 m5 l 5
分别对l 和q 求偏导数
2) 由第三个条件正形投影条件
y x x y 和 l q l q
dm0 dm2 2 dm4 4 2 4 m1 3m3 l 5m5 l dq dq l dq l 2m l 4m l 3 dm1 l dm3 l 3 ）
将各系数代入，略去高次项，得高斯投影坐标正算公式 精度为0.001m
xX N N sin B cos Bl 2 + sin B cos 3 B(5 - t 9 2 4 4 )l 4 + 2 24
N sin B cos 5 B(61 - 58t 2 t 4 )l 6 720
dl tan Adq
2 2 2 2 E ( dq ) 2 F tan A ( dq ) G tan A ( dq ) m2 2 2 2 r2 ( dq ) tan A ( dq )
E 2 F tan A G tan 2 A = r 2 sec 2 A E cos 2 A 2 F sin A cos A G sin 2 A = r2
高斯投影高斯投影正算公式
① 公式推导 a)级数展开
展开条件：经差 l 不大，在0～3.5°（0.061rad）
以内，展开后的形式（ l的幂级数 ） ：
xm0m1lm2l2m3l3m4l4..... yn0n1ln2l2n3l3n4l4......
由于展开点为P0（B， 0），则式中待定系数是等量 纬度q（或大地纬度B）的函数
纬线投影为平 行直线，经线 投影为与纬线 垂直而且间隔 相等的平行直 线，两经线间 的距离与相应 的经差成正比。
x f
y c
a dx Rd
b cd rd
dx c Rd r
dx c d cos
x
c
ln
tan
45
2
y c m c r
c=R切圆柱投影
如为椭 球呢？
② 高斯投影的几何概念高－斯等投角影横平切面椭圆柱投影
－
两带坐标重叠40km相
当于的经差22′。
3、高斯平面直角坐标系
Gauss Plane Rectangular Coordinate System
① 定义
x
纵坐标轴x：中央子午线的投
影线；
横坐标轴y：赤道的投影线；
o
y
原点o：中央子午线与赤道的
交点的投影点；
② 自然坐标与通用坐标
500km 中央子午线投影
④ UTM
长度比0.9996的选 择可以使6°带的 中央经线与边缘经 线的长度变形的绝 对值大致相等；
两条无长度变形的 割线的位置距中央 经线以东以西各 180km，相当于经 差约± 1°40′。
2、高斯投影的分带
Zone-dividing of Gauss Projection