1

13693261870

2023-09-23 ca69628172dbbc8e576ae55d588e99262bb8ac19

 ExportMap/cs/TerrainUtils.cs

@@ -13,6 +13,8 @@
    {
        private static int terrainMaxLevel = 0;

        private const string EPSG4326 = "EPSG:4326";

        /// <summary>
        /// 地形最大级别
        /// </summary>
@@ -62,10 +64,9 @@
        /// <returns>数据发布ID集合</returns>
        public static List<int> Generate(XYZArgs args, ref string err)
        {
            string tifFile = null;
            try
            {
                List<SysMeta> metas = XYZUtils.SelectMetas(args.ids, "and type in ('tif', 'tiff')");
                List<SysMeta> metas = XYZUtils.SelectMetas(args.ids, "and type in ('tif', 'tiff', 'dem')");
                if (null == metas || metas.Count == 0) return PrintInfo("找不到元数据");

                string dirPath = GetTerrainPath(args.dircode);
@@ -77,37 +78,34 @@
                int pubid = PubDBHelper.GetPushlishId(args.dircode, "DEM");
                List<int> metsIds = PubDBHelper.GetPublishMetaId(pubid);

                List<string> files = new List<string>();
                string uploadFolder = Tools.GetSetting("uploadFolder");
                foreach (SysMeta meta in metas)
                {
                    if (metsIds.Contains(meta.id)) return PrintInfo("元数据[" + meta.id + "] 已发布。");

                    if (metsIds.Contains(meta.id))
                    {
                        PrintInfo("元数据[" + meta.id + "] 已发布。");
                        continue;
                    }
                    string sourceFile = Path.Combine(uploadFolder, meta.path);
                    if (!File.Exists(sourceFile)) return PrintInfo("元数据[" + meta.id + "] 不存在。");
                    if (!File.Exists(sourceFile))
                    {
                        PrintInfo("元数据[" + meta.id + "] 不存在。");
                        continue;
                    }

                    string targetFile = Path.Combine(subPath, meta.name);
                    Reproject(sourceFile, targetFile, "", "EPSG:4326");
                    string targetFile = Path.Combine(subPath, meta.id + ".tif"); // meta.path.Split(new string[] { "\\", "//" }, StringSplitOptions.None)[1]
                    ConvertRaster(sourceFile, targetFile);
                    if (File.Exists(targetFile))
                    {
                        files.Add(targetFile);
                        CreateTerrain(args, targetFile, ref err);
                    }
                }



                /*tifFile = Merge(metas, args, ref err);
                if (!File.Exists(tifFile))
                {
                    LogOut.Info("TerraUtils：找不到tifFile数据。");
                    return null;
                }

                CreateLayerJson(args, files, ref err);
                string json = Path.Combine(dirPath, "layer.json");
                if (File.Exists(json)) File.Delete(json);
                //Generate(args, tifFile, ref err);
                CreateTerrain(args, tifFile, ref err);

                if (!File.Exists(json))
                {
                    LogOut.Info("TerraUtils：找不到layer.json文件。");
                    return null;
                }*/
                if (!File.Exists(json)) return PrintInfo("找不到layer.json文件");

                Complement(args);
                pubid = InsertToDB(metas, args);
@@ -131,25 +129,71 @@
            return null;
        }

        /// <summary>
        /// 转换栅格文件
        /// </summary>
        private static void ConvertRaster(string sourceFile, string targetFile)
        {
            if (File.Exists(targetFile)) return;

            string epsg = Tools.GetEPSG(sourceFile);
            if (epsg == EPSG4326)
            {
                File.Copy(sourceFile, targetFile);
                return;
            }

            Reproject(sourceFile, targetFile, epsg, EPSG4326);
        }

        /// <summary>
        /// 创建高程切片
        /// </summary>
        private static void CreateTerrain(XYZArgs args, string tifFile, ref string err)
        {
            string ctbPath = Tools.GetSetting("ctbPath");
            string ctbPath = Tools.GetSetting("ctbJL");
            string dirPath = GetTerrainPath(args.dircode);
            int maxLevel = GetTerrainMaxLevel(args, tifFile);
            //int maxLevel =  GetTerrainMaxLevel(args, tifFile);

            string gdal_data = string.Format("set GDAL_DATA={0}\\data", ctbPath);
            // -N 顶点法线, -C 强制创建缺失根瓦片, -R 不覆盖现有文件
            string createMesh = string.Format("{0}\\ctb-tile.exe -R -C -s {4} -o \"{1}\" -f Mesh \"{2}\\{3}.tif\"", ctbPath, dirPath, dirPath, args.dircode, maxLevel);
            //string createLayer = string.Format("{0}\\ctb-tile.exe -l -s {4} -o \"{1}\" -f Mesh \"{2}\\{3}.tif\"", ctbPath, dirPath, dirPath, args.dircode, maxLevel);
            // -N 顶点法线, -C 强制创建缺失根瓦片, -R 不覆盖现有文件，-f Mesh
            string createMesh = string.Format("{0}\\ctb-tile.exe -s {1} -o \"{2}\" \"{3}\"", ctbPath, TERRAIN_MAX_LEVEL, dirPath, tifFile);

            List<string> list = new List<string>() { gdal_data, createMesh };
            SysTask task = TaskDBHelper.CreateTask(args, "DEM", "高程数据(DEM)");
            err = Tools.ExecCmd(task, list);
        }

        /// <summary>
        /// 创建layer.json
        /// </summary>
        private static string CreateLayerJson(XYZArgs args, List<string> files, ref string err)
        {
            string dirPath = GetTerrainPath(args.dircode);
            string vrt = Path.Combine(dirPath, "subs", "dem.vrt");
            CreateVrt(files, vrt);
            if (!File.Exists(vrt)) return null;

            string ctbPath = Tools.GetSetting("ctbJL");
            string gdal_data = string.Format("set GDAL_DATA={0}\\data", ctbPath);
            string createLayer = string.Format("{0}\\ctb-tile.exe -l -s {1} -o \"{2}\" \"{3}\"", ctbPath, TERRAIN_MAX_LEVEL, dirPath, vrt);
            Tools.ExecCmd(new List<string>() { gdal_data, createLayer });

            return vrt;
        }

        /// <summary>
        /// 创建VRT
        /// </summary>
        private static void CreateVrt(List<string> files, string vrt)
        {
            string str = string.Join("\r\n", files);
            string filelist = vrt.Replace("dem.vrt", "list.txt");
            File.WriteAllText(filelist, str);

            string gdalPath = Tools.GetSetting("gdalPath");
            string cmd = string.Format("\"{0}\\gdalbuildvrt.exe\" -input_file_list \"{1}\" \"{2}\" -overwrite", gdalPath, filelist, vrt);
            Tools.ExecCmd(new List<string> { cmd });
        }

        /// <summary>
@@ -160,15 +204,16 @@
            FileInfo fi = new FileInfo(tifFile);
            if (fi.Length > DEFAULT_MAX_SIZE) return TERRAIN_MAX_LEVEL;

            string ctbPath = Tools.GetSetting("ctbPath");
            string ctbPath = Tools.GetSetting("ctbJL");
            string dirPath = GetTerrainPath(args.dircode);
            string subPath = Path.Combine(dirPath, "subs");

            string gdal_data = string.Format("set GDAL_DATA={0}\\data", ctbPath);
            string createLayer = string.Format("{0}\\ctb-tile.exe -l -o \"{1}\" -f Mesh \"{2}\\{3}.tif\"", ctbPath, dirPath, dirPath, args.dircode);
            string createLayer = string.Format("{0}\\ctb-tile.exe -l -o \"{1}\" \"{2}\"", ctbPath, subPath, tifFile);

            Tools.ExecCmd(new List<string>() { gdal_data, createLayer });

            string layerJson = Path.Combine(dirPath, "layer.json");
            string layerJson = Path.Combine(subPath, "layer.json");
            if (!File.Exists(layerJson)) return TERRAIN_MAX_LEVEL;

            string[] lines = File.ReadAllLines(layerJson, Encoding.UTF8);
@@ -186,7 +231,7 @@
        /// <summary>
        /// 补充文件
        /// </summary>
        private static void Complement(XYZArgs args)
        public static void Complement(XYZArgs args)
        {
            string dirPath = GetTerrainPath(args.dircode);
            string p_0_0 = Path.Combine(dirPath, "0", "0");
@@ -195,13 +240,12 @@
            string p_0_1 = Path.Combine(dirPath, "0", "1");
            if (!Directory.Exists(p_0_1)) Directory.CreateDirectory(p_0_1);

            string s_0_0_0 = Path.Combine(SGUtils.LFData, "dem", "0", "0", "0.terrain");
            string s_0 = Tools.BaseDir + "\\Sources\\0.terrain";
            string d_0_0_0 = Path.Combine(dirPath, "0", "0", "0.terrain");
            if (!File.Exists(d_0_0_0)) File.Copy(s_0_0_0, d_0_0_0, true);
            if (!File.Exists(d_0_0_0)) File.Copy(s_0, d_0_0_0, true);

            string s_0_1_0 = Path.Combine(SGUtils.LFData, "dem", "0", "1", "0.terrain");
            string d_0_1_0 = Path.Combine(dirPath, "0", "1", "0.terrain");
            if (!File.Exists(d_0_1_0)) File.Copy(s_0_1_0, d_0_1_0, true);
            if (!File.Exists(d_0_1_0)) File.Copy(s_0, d_0_1_0, true);

            string layerJson = Path.Combine(dirPath, "layer.json");
            string[] lines = File.ReadAllLines(layerJson, Encoding.UTF8);
@@ -233,7 +277,7 @@
                    if (!ids.Contains(m.id)) PubDBHelper.InsertMetaPub(m.id, pubid, args.userId);
                }
                string geom = GetPointZ(args);
                PubDBHelper.UpdatePublish(pubid, args.name, args.userId, geom);
                int rows = PubDBHelper.UpdatePublish(pubid, args.name, args.userId, geom);

                return pubid;
            }
@@ -268,11 +312,47 @@
        /// <summary>
        /// 获取中心点
        /// </summary>
        private static string GetPointZ(XYZArgs args)
        public static string GetPointZ(XYZArgs args)
        {
            double x = 0, y = 0;
            string dirPath = GetTerrainPath(args.dircode);
            string vrt = Path.Combine(dirPath, "subs", "dem.vrt");
            if (!File.Exists(vrt)) return null;

            return string.Format("ST_GeomFromText('POINT Z ({0} {1} {2})')", x, y, 12);
            string gdalPath = Tools.GetSetting("gdalPath");
            string cmd = string.Format("\"{0}\\gdalinfo.exe\" \"{1}\"", gdalPath, vrt);

            string rs = null;
            Tools.ExecCmd(new List<string> { cmd }, ref rs);
            string center = GetCenter(rs);
            if (string.IsNullOrEmpty(rs)) return null;

            string[] strs = center.Split(new string[] { "," }, StringSplitOptions.None);
            double x = double.Parse(strs[0]);
            double y = double.Parse(strs[1]);

            return string.Format("ST_GeomFromText('POINT Z ({0} {1} {2})')", y, x, 12);
        }

        /// <summary>
        /// 获取中心坐标
        /// </summary>
        private static string GetCenter(string rs)
        {
            if (string.IsNullOrEmpty(rs)) return null;

            string[] strs = rs.Split(new string[] { "\r\n" }, StringSplitOptions.RemoveEmptyEntries);
            foreach (string str in strs)
            {
                if (str.StartsWith("Center "))
                {
                    int start = str.IndexOf("(");
                    int end = str.IndexOf(")");

                    return str.Substring(start + 1, end - start - 1).Replace(" ", "");
                }
            }

            return null;
        }

        /// <summary>
@@ -281,13 +361,13 @@
        public static void Reproject(string sourceFile, string targetFile, string sourceSrs, string targetSrs)
        {
            string gdalPath = Tools.GetSetting("gdalPath");
            string cmd = string.Format("{0}\\gdalwarp.exe -s_srs {1} -t_srs {2} -r near -of GTiff {3} {4}", gdalPath, sourceSrs, targetSrs, sourceFile, targetFile);
            string cmd = string.Format("\"{0}\\gdalwarp.exe\" -s_srs {1} -t_srs {2} -r near -of GTiff \"{3}\" \"{4}\"", gdalPath, sourceSrs, targetSrs, sourceFile, targetFile);

            Tools.ExecCmd(new List<string> { cmd });
        }

        /// <summary>
        /// 指定投影
        /// 指定投影 *
        /// </summary>
        public static void Project(string sourceFile, string targetSrs)
        {