algorithm/line_vacuate/task_algo_c.cpp

/**
******************************************************************************
* @file           : task_alog_c.cpp
* @author         : wyj
* @brief          : C语言语法测试
* @attention      : None
* @date           : 2025/5/9
******************************************************************************
*/

#ifdef __cplusplus
extern "C"
{
#endif

#include "geography/studio_proj_c.h"
#include "geometry/studio_geo_algo_c.h"

#ifdef __cplusplus
}
#endif

#include <stdio.h>
#include <stdio.h>
#include "geometry/studio_geo_utils.h"

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);
        studio_line_c_add_point(c_line, p);
    }
}

void c_line_to_studio_line(const studio_line_c *c_line, studio_line &line)
{
    for (size_t i = 0; i < c_line->size; ++i)
    {
        const studio_point_c *p = studio_line_c_get_point(c_line, i);
        line.push_back(studio_point(p->x, p->y));
    }
}

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 gauss_line;
        double central = static_cast<int>(coll.m_line[0].x / 3) * 3;
        for (auto &point : coll.m_line)
        {
            double gx = 0.0;
            double gy = 0.0;
            lonlat_to_gauss(central, point.x, point.y, &gx, &gy);
            gauss_line.push_back(studio_point(gx, gy));
        }
        studio_line_c gauss_line_c = studio_line_c_init();
        studio_line_c vac_gauss_line_c = studio_line_c_init();
        studio_line_to_c_line(gauss_line, &gauss_line_c);

        // 简化线段，目标点数为28个
        int max_points = 28;
        double epsilon = 1.0;

        bool res = line_vacuate_c(&gauss_line_c, max_points, epsilon, &vac_gauss_line_c);
        if (!res)
        {
            std::cout << "Failed to simplify line." << std::endl;
            return 1;
        }

        // 高斯投影在转回经纬度
        studio_line simplified_line;
        for (int i = 0; i < vac_gauss_line_c.size; i++)
        {
            double lon = 0.0;
            double lat = 0.0;
            const studio_point_c *p = studio_line_c_get_point(&vac_gauss_line_c, i);
            gauss_to_lonlat(central, p->x, p->y, &lon, &lat);
            simplified_line.push_back(studio_point(lon, lat));
        }

        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(&gauss_line_c);
        studio_line_c_destroy(&vac_gauss_line_c);

        break;
    }
    std::string output_path;
#ifdef IS_WINDOWS
    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/fitting_examples_res_1_c.geojson";
#endif
    silly::geo::utils::write_geo_coll(output_path, res_collections);

    silly::geo::utils::destroy_gdal_env();

    return 0;
}