学步园

计算的方法如下图：

相应的代码如下（100W，50W，25W，10W，5W几种比例尺，其它比例尺同理）：

/// <summary> 从新图幅号计算角点坐标 </summary><br /> /// <param name="newNumber">输入的新图幅号(请先验证图幅号格式合法，再)</param><br /> /// <returns>角点列表</returns><br /> public static List<ESRI.ArcGIS.Geometry.IPoint> CalPointListFromNewNo(string newNumber)<br /> {<br /> List<ESRI.ArcGIS.Geometry.IPoint> resultLst = new List<ESRI.ArcGIS.Geometry.IPoint>();<br /> double delX = 0, delY = 0; //该比例尺下的图幅经度差和纬度差<br /> int H = 0, L = 0; //H为1:100万图幅纬度带字符对应的数字（A对应1，以此类推），<br /> //L为1:100万图幅经度带的数字码，例如J50中的50<br /> int h = 0, l = 0; //h为所求比例尺地形图所在的行号，例如I50B002001中的002<br /> //l为所求比例尺地形图所在的列号，例如I50B002001中的001</p> <p> #region 读取比例尺及行列号，并为经差及纬差赋值<br /> if (newNumber.Length == 3) //100w<br /> {<br /> delX = 6;<br /> delY = 4;<br /> //100w比例尺的情况下，没有后边六位的行列号，直接认为h和l是0即可<br /> }<br /> else<br /> {<br /> //其它比例尺下，提取行列号<br /> if (!int.TryParse(newNumber.Substring(4, 3), out h)<br /> || !int.TryParse(newNumber.Substring(7, 3), out l))<br /> return null;<br /> }<br /> if (!int.TryParse(newNumber.Substring(1, 2), out L)) return null;<br /> H = (int)newNumber[0] - (int)'A' + 1;<br /> switch (newNumber[3])<br /> {<br /> case 'B': //50w<br /> delX = 3;<br /> delY = 2;<br /> break;<br /> case 'C': //25w<br /> delX = 1.5;<br /> delY = 1;<br /> break;<br /> case 'D': //10w<br /> delX = 0.5;<br /> delY = 1.0 / 3;<br /> break;<br /> case 'E': //5w<br /> delX = 0.25;<br /> delY = 1.0 / 6;<br /> break;<br /> default:<br /> return null; //其它比例尺，暂没有写<br /> }<br /> #endregion<br /> //左下角点<br /> double x = (L - 31) * 6 + (l - 1) * delX;<br /> double y = (H - 1) * 4 + (4 / delY - h) * delY;<br /> ESRI.ArcGIS.Geometry.IPoint mPnt = new ESRI.ArcGIS.Geometry.PointClass();<br /> mPnt.PutCoords(x, y);<br /> mPnt.Z = 1; //这里赋Z值并没有实际作用，只是为了防止出现“the geometry has null z value”错误<br /> resultLst.Add(mPnt);<br /> mPnt = new ESRI.ArcGIS.Geometry.PointClass();<br /> mPnt.PutCoords(x+delX, y);<br /> mPnt.Z = 1;<br /> resultLst.Add(mPnt);<br /> mPnt = new ESRI.ArcGIS.Geometry.PointClass();<br /> mPnt.PutCoords(x+delX, y+delY);<br /> mPnt.Z = 1;<br /> resultLst.Add(mPnt);<br /> mPnt = new ESRI.ArcGIS.Geometry.PointClass();<br /> mPnt.PutCoords(x, y+delY);<br /> mPnt.Z = 1;<br /> resultLst.Add(mPnt);<br /> return resultLst;<br /> }