切图算法简述2——WGS 84坐标系下XYZ（WMTS）、TMS切片文件坐标生成

• 简述

  基于WGS 84 地理坐标系就WMTS（XYZ）服务、TMS服务切图文件坐标生成进行简述，并且简单以C# 代码进行实现。下面将所有代码贴出：

• 具体代码实现（C#） 

/// <summary>
        /// 通过文件后缀名判断该文件是否为图片格式
        /// </summary>
        /// <param name="fileExtension">文件后缀名（.tif、.bil、.jpg、.raster、.bt、.png）</param>
        /// <returns></returns>
        public bool IsImage(string fileExtension)
        {
            return !string.IsNullOrEmpty(fileExtension) && (fileExtension.Equals(".tif") || fileExtension.Equals(".bil") || (fileExtension.Equals(".jpg") || fileExtension.Equals(".raster")) || (fileExtension.Equals(".bt") || fileExtension.Equals(".png")));
        }

        /// <summary>
        /// 通过栅格数据后缀名获取坐标文件后缀名
        /// </summary>
        /// <param name="fileExtension">栅格文件后缀名（.tif、.bil、.jpg、.raster、.bt、.png）</param>
        /// <returns></returns>
        public string GetCoordinateFileExtension(string fileExtension)
        {
            if (string.IsNullOrEmpty(fileExtension))
                return null;
            if (fileExtension.Equals(".tif"))
                return ".tfw";
            else if (fileExtension.Equals(".bil"))
                return ".blw";
            else if (fileExtension.Equals(".jpg"))
                return ".jgw";
            else if (fileExtension.Equals(".raster"))
                return ".rasterw";
            else if (fileExtension.Equals(".bt"))
                return ".btw";
            else if (fileExtension.Equals(".png"))
                return ".pgw";
            return null;
        }


        /// <summary>
        /// 创建WGS84地理坐标系文件（.prj文件），这儿需要引入ProjNet.dll类库（可用NuGET获取）
        /// </summary>
        /// <param name="fileName">文件名称，我这儿不标准（仅仅适用于作者本人），可根据自己情况重写。
        /// （不带后缀名的文件路径：例如：D:\\sample\\image）</param>
        public void CreateCoordinateSystemWktFile(string fileName)
        {
            if (string.IsNullOrEmpty(fileName))
                return;
            IGeographicCoordinateSystem wgS84 = GeographicCoordinateSystem.WGS84;
            try
            {
                StreamWriter streamWriter = new StreamWriter(fileName + ".prj", false);
                streamWriter.WriteLine(wgS84.WKT);
                streamWriter.Close();
                Console.WriteLine(fileName + ".prj:创建成功！");
            }
            catch
            {
                Console.WriteLine(fileName + ".prj:创建失败！");
            }
        }

        /// <summary>
        /// 基于TMS为栅格数据创建坐标文件
        /// </summary>
        /// <param name="row">行序号</param>
        /// <param name="column">列序号</param>
        /// <param name="zoom">地图层级</param>
        /// <param name="fullFileName">坐标文件名称（整个文件名称，例如：D:\\smaple\\image.pgw）</param>
        public void CreateTMSCoordinateFile(int row,int column,int zoom,string fullFileName)
        {
            double unitsPerPixel=UnitsPerPixel(256, zoom);
            //栅格数据左上角经度
            double leftTopLng = TMSLongitude(column, zoom);
            //栅格数据左上角纬度
            double leftTopLat = TMSLatitude(row + 1, zoom);
            double rotate_1 = 0.0;
            double rotate_2 = 0.0;
            string messageCoor = string.Format("{0}\r\n{1}\r\n{2}\r\n{3}\r\n{4}\r\n{5}\r\n", unitsPerPixel, rotate_1.ToString("0.00000000000000"), rotate_2.ToString("0.00000000000000"), -unitsPerPixel, leftTopLng, leftTopLat);

            try
            {
                StreamWriter streamWriter = new StreamWriter(fullFileName, false);
                streamWriter.WriteLine(messageCoor);
                streamWriter.Close();
                Console.WriteLine(fullFileName + ":创建成功！");
            }
            catch
            {
                Console.WriteLine(fullFileName + ":创建失败！");
            }
        }

        /// <summary>
        /// 基于XYZ为栅格数据创建坐标文件
        /// </summary>
        /// <param name="row">行序号</param>
        /// <param name="column">列序号</param>
        /// <param name="zoom">地图层级</param>
        /// <param name="fullFileName">坐标文件名称（整个文件名称，例如：D:\\smaple\\image.pgw）</param>
        public void CreateXYZCoordinateFile(int row, int column, int zoom, string fullFileName)
        {
            double unitsPerPixel = UnitsPerPixel(256, zoom);
            //栅格数据左上角经度
            double leftTopLng = XYZLongitude(column, zoom);
            //栅格数据左上角纬度
            double leftTopLat = XYZLatitude(row, zoom);
            double rotate_1 = 0.0;
            double rotate_2 = 0.0;
            string messageCoor = string.Format("{0}\r\n{1}\r\n{2}\r\n{3}\r\n{4}\r\n{5}\r\n", unitsPerPixel, rotate_1.ToString("0.00000000000000"), rotate_2.ToString("0.00000000000000"), -unitsPerPixel, leftTopLng, leftTopLat);

            try
            {
                StreamWriter streamWriter = new StreamWriter(fullFileName, false);
                streamWriter.WriteLine(messageCoor);
                streamWriter.Close();
                Console.WriteLine(fullFileName + ":创建成功！");
            }
            catch
            {
                Console.WriteLine(fullFileName + ":创建失败！");
            }
        }