using ExportMap.db;
|
using ExportMap.Models;
|
using System;
|
using System.Collections.Generic;
|
using System.IO;
|
using System.Linq;
|
using System.Web;
|
|
namespace ExportMap.cs
|
{
|
public class PyLasUtils
|
{
|
/// <summary>
|
/// 获取路径
|
/// </summary>
|
public static string GetPath(int id)
|
{
|
return Path.Combine(SGUtils.LFData, "3d\\3dtiles\\las", id.ToString());
|
}
|
|
/// <summary>
|
/// 获取发布地址
|
/// </summary>
|
public static string GetReleaseUrl(string path)
|
{
|
return "{host}/LFData/" + path.Replace("\\", "/");
|
}
|
|
/// <summary>
|
/// 生成
|
/// </summary>
|
/// <param name="args">XYZ参数</param>
|
/// <param name="err">错误信息</param>
|
/// <returns>数据发布ID集合</returns>
|
public static List<int> Generate(XYZArgs args, ref string err)
|
{
|
try
|
{
|
List<SysMeta> list = XYZUtils.SelectMetas(args.ids, "and type in ('las', 'laz')");
|
if (null == list || list.Count == 0) return null;
|
|
string pyPath = Tools.GetSetting("pythonFolder");
|
string uploadFolder = Tools.GetSetting("uploadFolder");
|
|
List<int> ids = new List<int>();
|
foreach (SysMeta meta in list)
|
{
|
string lasPath = Path.Combine(uploadFolder, meta.path);
|
if (!File.Exists(lasPath)) continue;
|
|
string outPath = GetPath(meta.id);
|
if (Directory.Exists(outPath)) Tools.DelPath(outPath);
|
if (!Directory.Exists(outPath)) Directory.CreateDirectory(outPath);
|
|
//int idx = args.ids.IndexOf(meta.id); // args.srids[idx], args.zs[idx])
|
List<string> cmds = new List<string>();
|
cmds.Add(string.Format("set PYTHONHOME={0}", pyPath));
|
cmds.Add(string.Format("{0}\\Scripts\\py3dtiles.exe convert --overwrite {1} --out {2}", pyPath, lasPath, outPath));
|
|
SysTask task = TaskDBHelper.CreateTask(args, meta, "LAS", "点云数据(LAS)");
|
err = Tools.ExecCmd(task, cmds, false);
|
err = null;
|
|
string jsonFile = Path.Combine(outPath, "tileset.json");
|
if (File.Exists(jsonFile))
|
{
|
string path = jsonFile.Replace(Tools.GetSetting("lfData") + "\\", "");
|
int pubid = InsertToDB(meta, args, path);
|
if (pubid > 0) ids.Add(pubid);
|
}
|
}
|
|
return ids;
|
}
|
catch (Exception ex)
|
{
|
LogOut.Error(ex.Message + "\r\n" + ex.StackTrace);
|
err = ex.Message;
|
return null;
|
}
|
}
|
|
/// <summary>
|
/// 插入数据库
|
/// </summary>
|
private static int InsertToDB(SysMeta meta, XYZArgs args, string path)
|
{
|
//if (PubDBHelper.IsPublish(meta.id)) return 1;
|
int pubid = PubDBHelper.GetPushlishId(meta.id);
|
if (pubid > 0)
|
{
|
PubDBHelper.UpdatePublish(pubid, meta.name, args.userId, null);
|
return pubid;
|
}
|
|
SysPublish sys = Tools.NewPublish(meta, args, GetReleaseUrl(path), path);
|
sys.json = GetLasGeomJSON(Path.Combine(Tools.GetSetting("uploadFolder"), meta.path), meta.id);
|
|
pubid = PubDBHelper.InsertPublish(sys);
|
if (pubid > 0)
|
{
|
sys.id = pubid;
|
PubDBHelper.InsertLayer(sys, meta);
|
PubDBHelper.InsertMetaPub(meta.id, pubid, args.userId);
|
}
|
|
return pubid;
|
}
|
|
/// <summary>
|
/// 坐标转换 *
|
/// </summary>
|
public static void CsTransform(string epsg, double x, double y, double z = 0)
|
{
|
List<string> cmds = new List<string>();
|
|
string gdalPath = Tools.GetSetting("gdalPath");
|
cmds.Add(string.Format("\"{0}\\gdaltransform.exe\" -s_srs EPSG:{1} -t_srs EPSG:4326", gdalPath, epsg.Replace("EPSG:", "")));
|
cmds.Add(string.Format("echo {0} {1} {2}", x, y, z));
|
|
string rs = null;
|
string str = Tools.ExecCmd(cmds, ref rs);
|
|
//return rs;
|
}
|
|
/// <summary>
|
/// 获取las坐标JSON
|
/// </summary>
|
public static string GetLasGeomJSON(string lasFile, int id)
|
{
|
try
|
{
|
string cmd = string.Format("\"{0}\\lasinfo.exe\" \"{1}\"", Tools.GetSetting("tilerPath"), lasFile);
|
string rs = Tools.ExecCmd(new List<string> { cmd });
|
if (string.IsNullOrEmpty(rs)) return null;
|
|
string[] strs = rs.Split(new string[] { "\r\n" }, StringSplitOptions.RemoveEmptyEntries);
|
string epsg = GetEpsgFromStrs(strs);
|
string minStr = GetCoordFromStrs(strs, "min x y z:");
|
string maxStr = GetCoordFromStrs(strs, "max x y z:");
|
string scaleStr = GetCoordFromStrs(strs, "scale factor x y z:");
|
if (string.IsNullOrEmpty(epsg) || string.IsNullOrEmpty(minStr) || string.IsNullOrEmpty(maxStr) || string.IsNullOrEmpty(scaleStr)) return null;
|
|
string[] minXYZ = minStr.Split(new string[] { " " }, StringSplitOptions.RemoveEmptyEntries);
|
string[] maxXYZ = maxStr.Split(new string[] { " " }, StringSplitOptions.RemoveEmptyEntries);
|
string[] scales = scaleStr.Split(new string[] { " " }, StringSplitOptions.RemoveEmptyEntries);
|
double minX = double.Parse(minXYZ[0]), minY = double.Parse(minXYZ[1]), minZ = double.Parse(minXYZ[2]);
|
double maxX = double.Parse(maxXYZ[0]), maxY = double.Parse(maxXYZ[1]), maxZ = double.Parse(maxXYZ[2]);
|
double x = (minX + maxX) / 2, y = (minY + maxY) / 2, z = maxZ, scale = Math.Round(1 / double.Parse(scales[0]), 2) / 10;
|
|
string url = ExportUtil.LFServer + "/publish/selectCsTransform?token=" + Tools.Token + "&x=" + x + "&y=" + y + "&epsg=" + epsg;
|
rs = ExportUtil.GetData(url);
|
if (string.IsNullOrEmpty(rs)) return null;
|
|
strs = rs.Replace("\"", "").Split(new string[] { "," }, StringSplitOptions.None);
|
x = double.Parse(strs[0]);
|
y = double.Parse(strs[1]);
|
|
return string.Format("{{\"lon\":{0},\"lat\":{1},\"height\":{2},\"yaw\":0,\"alpha\":1,\"modelid\":\"DbId\",\"scale\":{3},\"lasId\":{4}}}", x, y, z, scale, id);
|
}
|
catch (Exception ex)
|
{
|
LogOut.Error(ex.StackTrace);
|
return null;
|
}
|
}
|
|
/// <summary>
|
/// 从字符串数组中获取EPSG
|
/// </summary>
|
private static string GetEpsgFromStrs(string[] strs)
|
{
|
foreach (string str in strs)
|
{
|
if (str.Contains("key 3072 "))
|
{
|
int start = str.IndexOf("value_offset") + 12;
|
int end = str.IndexOf("-", start);
|
|
return str.Substring(start, end - start).Replace(" ", "");
|
}
|
}
|
|
return null;
|
}
|
|
/// <summary>
|
/// 从字符串数组中获取坐标
|
/// </summary>
|
private static string GetCoordFromStrs(string[] strs, string key)
|
{
|
foreach (string str in strs)
|
{
|
if (str.Contains(key))
|
{
|
return str.Replace(key, "");
|
}
|
}
|
|
return null;
|
}
|
}
|
}
|
