问题

我想在Matlab中检测两条WSG84线的交点，已知以下信息。

红线：latStart0, lonStart0, latEnd0, lonEnd0

蓝线：latStart1, lonStart1, latEnd1, lonEnd1

英文:

I would like to detect intersection of two WSG84 lines in Matlab given the information below known.

RedLine: latStart0, lonStart0, latEnd0, lonEnd0

BlueLine: latStart1, lonStart1, latEnd1, lonEnd1

How can i achieve this ?

答案1

得分: 0

以下是解决方案：

function intersectionPoint = findIntersection(line1_start, line1_end, line2_start, line2_end)
    line1_start_rad = deg2rad(line1_start);
    line1_end_rad = deg2rad(line1_end);
    line2_start_rad = deg2rad(line2_start);
    line2_end_rad = deg2rad(line2_end);
    
    % 计算线段的分母
    denom = (line1_end_rad(1) - line1_start_rad(1)) * (line2_end_rad(2) - line2_start_rad(2)) - ...
        (line1_end_rad(2) - line1_start_rad(2)) * (line2_end_rad(1) - line2_start_rad(1));
    if abs(denom) < eps
        intersectionPoint = struct('latitude', NaN, 'longitude', NaN);
        return;
    end
    numer = (line2_end_rad(1) - line2_start_rad(1)) * (line1_start_rad(2) - line2_start_rad(2)) - ...
        (line2_end_rad(2) - line2_start_rad(2)) * (line1_start_rad(1) - line2_start_rad(1));
    t1 = numer / denom;
    
    if t1 < 0
        intersectionPoint = struct('latitude', NaN, 'longitude', NaN);
        return;
    end
    
    % 计算交点
    lonIntersection_rad = line1_start_rad(2) + (line1_end_rad(2) - line1_start_rad(2)) * t1;
    latIntersection_rad = line1_start_rad(1) + (line1_end_rad(1) - line1_start_rad(1)) * t1;
    latIntersection_deg = rad2deg(latIntersection_rad);
    lonIntersection_deg = rad2deg(lonIntersection_rad);
    intersectionPoint = struct('latitude', latIntersection_deg, 'longitude', lonIntersection_deg);
end

英文:

Here it is the solution:

function intersectionPoint = findIntersection(line1_start, line1_end, line2_start, line2_end)
    line1_start_rad = deg2rad(line1_start);
    line1_end_rad = deg2rad(line1_end);
    line2_start_rad = deg2rad(line2_start);
    line2_end_rad = deg2rad(line2_end);
    
    % Calculate the denominator for the lines
    denom = (line1_end_rad(1) - line1_start_rad(1)) * (line2_end_rad(2) - line2_start_rad(2)) - ...
        (line1_end_rad(2) - line1_start_rad(2)) * (line2_end_rad(1) - line2_start_rad(1));
    if abs(denom) &lt; eps
        intersectionPoint = struct(&#39;latitude&#39;, NaN, &#39;longitude&#39;, NaN);
        return;
    end
    numer = (line2_end_rad(1) - line2_start_rad(1)) * (line1_start_rad(2) - line2_start_rad(2)) - ...
        (line2_end_rad(2) - line2_start_rad(2)) * (line1_start_rad(1) - line2_start_rad(1));
    t1 = numer / denom;
    
    if t1 &lt; 0
        intersectionPoint = struct(&#39;latitude&#39;, NaN, &#39;longitude&#39;, NaN);
        return;
    end
    
    % Calculate the intersection point
    lonIntersection_rad = line1_start_rad(2) + (line1_end_rad(2) - line1_start_rad(2)) * t1;
    latIntersection_rad = line1_start_rad(1) + (line1_end_rad(1) - line1_start_rad(1)) * t1;
    latIntersection_deg = rad2deg(latIntersection_rad);
    lonIntersection_deg = rad2deg(lonIntersection_rad);
    intersectionPoint = struct(&#39;latitude&#39;, latIntersection_deg, &#39;longitude&#39;, lonIntersection_deg);
end

通过集体智慧和协作来改善编程学习和解决问题的方式。致力于成为全球开发者共同参与的知识库，让每个人都能够通过互相帮助和分享经验来进步。

在Matlab中检测两条WGS84线的交点。

问题

答案1

Reading a Fortran Data File in Python

在Python中无法读取.mat文件。

解决二次阻力耦合微分方程。

在MATLAB中，如何在x轴上绘制离散点数组 [x1, … xn]？

如何在Playwright视觉比较中屏蔽多个定位器？

在C++中，可以使用可变模板参数来检索类型的内部类型。

selenium.common.exceptions.StaleElementReferenceException: Message: stale element reference: stale element not found

Creating and opening a URL to log in to Website via Basic Auth with Robot Framework/Selenium (Python)

AG Grid 在上下文菜单中以大文本形式打开

What's the correct way to type hint an empty list as a literal in python?

如何在Highcharts Gantt中更改本地化的星期名称

如何在同一个流中使用多个过滤器和映射函数？

如何使用Map/Set来将代码优化到O(n)？

.NET MAUI Android在GitHub Actions上构建失败，错误代码为1。