路径压缩支持大地坐标系

line_vacuate/studio_geo_algo_c.c

@@ -142,7 +142,6 @@ bool line_vacuate_c(const studio_line_c *line, const int max_points, const doubl
 {
     // 初始化状态变量
     bool status = false;
-
     // 如果输入的线段为空，则打印错误信息并返回false
     if (line->size == 0)
     {
@@ -194,6 +193,46 @@ bool line_vacuate_c(const studio_line_c *line, const int max_points, const doubl
     return status;
 }
 
+bool line_vacuate_gcs_c(const studio_line_c *line, const int max_points, const double epsilon,
+                        studio_line_c *vacuate_line)
+{
+    bool status = false;
+    if (!line || !vacuate_line || line->size <= 1) // 添加空指针和长度校验
+    {
+        printf("line or vacuate_line is null or line size <= 1");
+        return status;
+    }
+    // 经纬度转高斯克吕格投影
+    studio_line_c gauss_line = studio_line_c_init();
+    double central = (int) (line->data[0].x) / 3 * 3;
+    for (unsigned int i = 0; i < line->size; i++)
+    {
+        double gx, gy;
+        lonlat_to_gauss(central, line->data[i].x, line->data[i].y, &gx, &gy);
+        studio_point_c point = studio_point_init(gx, gy);
+        studio_line_c_add_point(&gauss_line, point);
+    }
+    studio_line_c vac_gauss_line = studio_line_c_init();
+    bool res = line_vacuate_c(&gauss_line, max_points, epsilon, &vac_gauss_line);
+    if (!res)
+    {
+        printf("Failed to simplify line.");
+        return status;
+    }
+    studio_line_c_destroy(&gauss_line);
+    // 高斯克吕格转经纬度
+    for (unsigned int i = 0; i < vac_gauss_line.size; i++)
+    {
+        double lon, lat;
+        gauss_to_lonlat(central, vac_gauss_line.data[i].x, vac_gauss_line.data[i].y, &lon, &lat);
+        studio_point_c point = studio_point_init(lon, lat);
+        studio_line_c_add_point(vacuate_line, point);
+    }
+    studio_line_c_destroy(&vac_gauss_line);
+    status = true;
+    return status;
+}
+
 
 bool remove_outliers_c(const studio_line_c *line, double distance, studio_line_c *outliers_line)
 {
@@ -216,7 +255,7 @@ bool remove_outliers_c(const studio_line_c *line, double distance, studio_line_c
     }
 
     // 离群点检测 ------------------------------------------------
-    if (gauss_line.size<= 1)
+    if (gauss_line.size <= 1)
     {
         printf("gauss_line size <= 1");
         return status;
@@ -229,10 +268,9 @@ bool remove_outliers_c(const studio_line_c *line, double distance, studio_line_c
         if (d_sq <= dist_square)
         {
             studio_line_c_add_point(&gs_filtered_line, gauss_line.data[i]);
-        }
-        else
+        } else
         {
-            int a=0;
+            int a = 0;
         }
     }
     studio_line_c_destroy(&gauss_line); // 释放原始数据的高斯投影

line_vacuate/studio_geo_algo_c.h

@@ -39,13 +39,25 @@ void douglas_peucker_c(const studio_line_c *points, size_t start, size_t end, do
                        unsigned int *index_count);
 
 /// 矢量线段抽稀算法
-/// \param line 线段
+/// \param line 原始线段
 /// \param max_points 压缩后最大点
 /// \param epsilon 抽稀容差,若某点到当前线段的最大垂直距离超过epsilon，则保留该点,单位与适量坐标系的单位一致
 /// \param vacuate_line 抽稀后的线段
 /// \return
 bool line_vacuate_c(const studio_line_c *line, const int max_points, const double epsilon, studio_line_c *vacuate_line);
 
+/// 矢量线段抽稀算法 大地坐标系(GCS), 坐标点为十进制经纬度
+/// @param line 原始线段
+/// @param max_points 压缩后最大点
+/// @param epsilon 抽稀容差,若某点到当前线段的最大垂直距离超过epsilon，这里单位是米
+/// @param vacuate_line 抽稀后的线段
+/// @return
+bool line_vacuate_gcs_c(const studio_line_c *line, const int max_points, const double epsilon,
+                        studio_line_c *vacuate_line);
+
+
+/////////////////// 移除线段中的离群点 ///////////////////
+
 /// 计算哈希值
 /// \param value
 /// \param size
@@ -58,9 +70,6 @@ unsigned int hash(unsigned int value, unsigned int size);
 void remove_duplicates(unsigned int **indices, unsigned int *size);
 
 
-/////////////////// 移除线段中的离群点 ///////////////////
-
-
 /// 移除线段中的离群点
 /// @param line 原始线段
 /// @param distance 离群点距离

line_vacuate/task_algo_c.cpp

@@ -26,8 +26,6 @@ extern "C"
 
 void studio_line_to_c_line(const studio_line &line, studio_line_c *c_line)
 {
-    //    c_line->size = line.size();
-    c_line->data = (studio_point_c *) malloc(line.size() * sizeof(studio_point_c));
     for (size_t i = 0; i < line.size(); ++i)
     {
         studio_point_c p = studio_point_init(line[i].x, line[i].y);
@@ -44,7 +42,61 @@ void c_line_to_studio_line(const studio_line_c *c_line, studio_line &line)
     }
 }
 
-///////////////////// 移除离群点 ////////////////////////////////
+// ///////////////////// 移除离群点 ////////////////////////////////
+// int main()
+// {
+//     printf("\n\n===================== %s =====================\n\n", __FILE__);
+//     silly::geo::utils::init_gdal_env();
+//     std::string path;
+// #ifdef IS_WINDOWS
+//     path = "D:/5_file/2_readfile/geojson/multi_point/fitting_examples.geojson";
+// #else
+//     path = "/home/wyj/myself/2_data/2_geojson/multi_point/fitting_examples.geojson";
+// #endif
+//
+//     std::cout << "path: " << path << std::endl;
+//     std::vector<studio_geo_coll> res_collections;
+//     std::vector<studio_geo_coll> collections;
+//     silly::geo::utils::read_geo_coll(path, collections);
+//
+//     for (auto &coll: collections)
+//     {
+//         studio_line_c line_c = studio_line_c_init();
+//         studio_line_to_c_line(coll.m_line, &line_c);
+//         std::cout << "原始数据量：" << coll.m_line.size() << std::endl;
+//
+//         studio_line_c outliers_line = studio_line_c_init();
+//         remove_outliers_c(&line_c, 1, &outliers_line);
+//
+//         studio_line simplified_line;
+//         c_line_to_studio_line(&outliers_line, simplified_line);
+//         std::cout << "简化后数据量：" << simplified_line.size() << std::endl;
+//
+//         studio_geo_coll temp;
+//         temp.m_type = enum_geometry_type::egtLineString;
+//         temp.m_line = simplified_line;
+//         res_collections.push_back(temp);
+//
+//         studio_line_c_destroy(&line_c);
+//         studio_line_c_destroy(&outliers_line);
+//
+//         break;
+//     }
+//
+//
+//     std::string output_path;
+// #ifdef IS_WINDOWS
+//     output_path = "D:/5_file/2_readfile/geojson/multi_point/remove_outliers.geojson";
+// #else
+//     output_path = "/home/wyj/myself/2_data/2_geojson/multi_point/remove_outliers.geojson";
+// #endif
+//     silly::geo::utils::write_geo_coll(output_path, res_collections);
+//      silly::geo::utils::destroy_gdal_env();
+//     return 0;
+// }
+
+
+// //////////////////// 路径压缩算法测试 //////////////////////
 int main()
 {
     printf("\n\n===================== %s =====================\n\n", __FILE__);
@@ -60,19 +112,23 @@ int main()
     std::vector<studio_geo_coll> res_collections;
     std::vector<studio_geo_coll> collections;
     silly::geo::utils::read_geo_coll(path, collections);
-
     for (auto &coll: collections)
     {
         studio_line_c line_c = studio_line_c_init();
+        studio_line_c vac_line_c = studio_line_c_init();
         studio_line_to_c_line(coll.m_line, &line_c);
-        std::cout << "原始数据量：" << coll.m_line.size() << std::endl;
-
-        studio_line_c outliers_line = studio_line_c_init();
-        remove_outliers_c(&line_c, 1, &outliers_line);
 
+        // 简化线段，目标点数为28个
+        int max_points = 28;
+        double epsilon = 1.0;
+        bool res = line_vacuate_gcs_c(&line_c, max_points, epsilon, &vac_line_c);
+        if (!res)
+        {
+            std::cout << "Failed to simplify line." << std::endl;
+            return 1;
+        }
         studio_line simplified_line;
-        c_line_to_studio_line(&outliers_line, simplified_line);
-        std::cout << "简化后数据量：" << simplified_line.size() << std::endl;
+        c_line_to_studio_line(&vac_line_c, simplified_line);
 
         studio_geo_coll temp;
         temp.m_type = enum_geometry_type::egtLineString;
@@ -80,25 +136,23 @@ int main()
         res_collections.push_back(temp);
 
         studio_line_c_destroy(&line_c);
-        studio_line_c_destroy(&outliers_line);
+        studio_line_c_destroy(&vac_line_c);
 
         break;
     }
-
-
     std::string output_path;
 #ifdef IS_WINDOWS
-    output_path = "D:/5_file/2_readfile/geojson/multi_point/remove_outliers.geojson";
+    output_path = "D:/5_file/2_readfile/geojson/multi_point/fitting_examples_res_1.geojson";
 #else
-    output_path = "/home/wyj/myself/2_data/2_geojson/multi_point/remove_outliers.geojson";
+    output_path = "/home/wyj/myself/2_data/2_geojson/multi_point/fitting_examples_res_2_c.geojson";
 #endif
     silly::geo::utils::write_geo_coll(output_path, res_collections);
+
     silly::geo::utils::destroy_gdal_env();
+
     return 0;
 }
 
-
-// //////////////////// 路径压缩算法测试 //////////////////////
 // int main()
 // {
 //     printf("\n\n===================== %s =====================\n\n", __FILE__);