jjzjj

java - 如何确定 GPS 坐标是否位于矩形区域内？

coder 2024-04-02 原文

我在 Stackoverflow 和其他网站上经历了许多类似的问题，我的解决方案基于这些答案，但我仍然无法让它工作......

我的问题:我想确定某个 GPS 位置 P 是否位于由四个给定 GPS 坐标 A 包围的矩形区域内>, B, C, D.

目前我正在计算三角形ABP、BCP、CDP 和DAP 的面积。如果这些区域中的任何一个大于零(请不要生气，数学家)，该点位于我的矩形之外，如解释的那样 here .

代码:

private static double triangleArea(Location a, Location b, Location c) {
    // (C.x*B.y-B.x*C.y)-(C.x*A.y-A.x*C.y)+(B.x*A.y-A.x*B.y)
    double result = (c.getLongitude()*b.getLatitude()-b.getLongitude()*c.getLatitude())-(c.getLongitude()*a.getLatitude()-a.getLongitude()*c.getLatitude())+(b.getLongitude()*a.getLatitude()-a.getLongitude()*b.getLatitude());
    return result;
}

public static boolean isInsideSquare(Location a, Location b, Location c, Location d, Location p) {
    if (triangleArea(a,b,p)>0 || triangleArea(b,c,p)>0 || triangleArea(c,d,p)>0 || triangleArea(d,a,p)>0) {
        return false;
    }
    return true;
}

但是当我调用这个函数时，它总是返回 false。即使使用以下坐标(纬度、经度):

A: (50.884706, 4.714151)
B: (50.884944, 4.716149)
C: (50.884679, 4.716228)
D: (50.884441, 4.714230)
P: (50.884538, 4.714615)

然而，当我在 map 上绘制这些点时(例如 here)，我可以看到点 P 位于 ABCD 内......(见下文)

我已经阅读了一些简单的解决方案，您只需减去点的 x、y 坐标，但这显然只适用于平行于 x 轴和 y 轴的矩形，而我的矩形可以在任何方向上定向角度。

谁能告诉我我做错了什么或提出更好的解决方案来解决这个问题？

非常感谢!

最佳答案

如果点在矩形内。在飞机上。对于数学家或大地测量学 (GPS) 坐标

  • 让我们延长矩形的边。所以我们有 4 条直线 lAB、lBC、lCD、lDA，或者，为了简短起见, l1, l2, l3, l4.

  • 对每个 li 做一个方程。等式排序:

    fi(P)=0。

P 是一个点。对于属于 li 的点，等式成立。

  • 我们需要等式左边的函数。它们是f1, f2, f3, f4。
  • 注意，对于来自 li 一侧的每个点，函数 fi 都大于 0，对于来自另一侧的点 fi 小于 0。
  • 因此，如果我们要检查 P 是否在矩形中，我们只需要 p 位于所有四条直线的正确边上。因此，我们必须检查四个函数的标志。
  • 但是直线的哪一侧是正确的，矩形属于哪一侧？它是侧面，不属于该线的矩形顶点所在的位置。为了进行检查，我们可以选择两个不属于的顶点中的任何一个。

  • 所以，我们必须检查一下:

    fAB(P) fAB(C) >= 0

    fBC(P) fBC(D) >= 0

    fCD(P) fCD(A) >= 0

    fDA(P) fDA(B) >= 0

不等式并不严格，因为如果一个点在边界上，它也属于矩形。如果你不需要边界上的点，你可以改变严格的不等式。但是，当您从事浮点运算时，选择是无关紧要的。

  • 对于矩形中的一个点，所有四个不等式都为真。请注意，它也适用于每个凸多边形，只是线/方程的数量会有所不同。

  • 唯一剩下的事情就是得到一条穿过两点的直线的方程式。这是一个著名的线性方程。让我们为直线 AB 和点 P 编写它:

    fAB(P) ≡ (xA-xB) (yP-yB) - (yA-yB) (xP-xB)

可以简化检查 - 让我们沿着矩形顺时针 - A、B、C、D、A。然后所有正确的边都会在线的右侧。因此，我们不必与另一个顶点所在的边进行比较。我们需要检查一组更短的不等式:

fAB(P) >= 0

fBC(P) >= 0

fCD(P) >= 0

fDA(P) >= 0

但这对于正常的数学家坐标集来说是正确的，其中 X 位于右侧，Y 位于顶部。对于在 GPS 中使用的大地测量学 坐标，其中 X 位于顶部，Y 位于右侧，我们必须对不等式进行转换:

fAB(P) <=>

fBC(P) <=>

fCD(P) <=>

fDA(P) <=>

如果您不确定轴的方向，请小心进行此简化检查 - 如果您选择了正确的不等式，请手动检查一个点。

顺便说一句，球体上的矩形也有意义。它是两条平行线和两条子午线之间的区域，当然，主轴的任何方向。但是根据评论，您不需要它。

关于java - 如何确定 GPS 坐标是否位于矩形区域内？，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/21787743/

位于javasubcodestrongalgorithmmathgpscoordinates

有关java - 如何确定 GPS 坐标是否位于矩形区域内？的更多相关文章

  1. ruby - 如何使用 Nokogiri 的 xpath 和 at_xpath 方法 - 2

    我正在学习如何使用Nokogiri，根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div

  2. ruby - 如何从 ruby​​ 中的字符串运行任意对象方法？ - 2

    总的来说，我对ruby​​还比较陌生，我正在为我正在创建的对象编写一些rspec测试用例。许多测试用例都非常基础，我只是想确保正确填充和返回值。我想知道是否有办法使用循环结构来执行此操作。不必为我要测试的每个方法都设置一个assertEquals。例如:describeitem,"TestingtheItem"doit"willhaveanullvaluetostart"doitem=Item.new#HereIcoulddotheitem.name.shouldbe_nil#thenIcoulddoitem.category.shouldbe_nilendend但我想要一些方法来使用

  3. python - 如何使用 Ruby 或 Python 创建一系列高音调和低音调的蜂鸣声？ - 2

    关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点？我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等)，但我确实想创建一个输出文件。

  4. ruby-on-rails - 如何验证 update_all 是否实际在 Rails 中更新 - 2

    给定这段代码defcreate@upgrades=User.update_all(["role=?","upgraded"],:id=>params[:upgrade])redirect_toadmin_upgrades_path,:notice=>"Successfullyupgradeduser."end我如何在该操作中实际验证它们是否已保存或未重定向到适当的页面和消息？ 最佳答案 在Rails3中，update_all不返回任何有意义的信息，除了已更新的记录数(这可能取决于您的DBMS是否返回该信息)。http://ar.ru

  5. ruby-on-rails - 'compass watch' 是如何工作的/它是如何与 rails 一起使用的 - 2

    我在我的项目目录中完成了compasscreate.和compassinitrails。几个问题:我已将我的.sass文件放在public/stylesheets中。这是放置它们的正确位置吗？当我运行compasswatch时，它不会自动编译这些.sass文件。我必须手动指定文件:compasswatchpublic/stylesheets/myfile.sass等。如何让它自动运行？文件ie.css、print.css和screen.css已放在stylesheets/compiled。如何在编译后不让它们重新出现的情况下删除它们？我自己编译的.sass文件编译成compiled/t

  6. ruby - 如何将脚本文件的末尾读取为数据文件(Perl 或任何其他语言) - 2

    我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚

  7. ruby - 如何指定 Rack 处理程序 - 2

    Rackup通过Rack的默认处理程序成功运行任何Rack应用程序。例如:classRackAppdefcall(environment)['200',{'Content-Type'=>'text/html'},["Helloworld"]]endendrunRackApp.new但是当最后一行更改为使用Rack的内置CGI处理程序时，rackup给出“NoMethodErrorat/undefinedmethod`call'fornil:NilClass”:Rack::Handler::CGI.runRackApp.newRack的其他内置处理程序也提出了同样的反对意见。例如Rack

  8. ruby - 如何每月在 Heroku 运行一次 Scheduler 插件？ - 2

    在选择我想要运行操作的频率时，唯一的选项是“每天”、“每小时”和“每10分钟”。谢谢!我想为我的Rails3.1应用程序运行调度程序。 最佳答案 这不是一个优雅的解决方案，但您可以安排它每天运行，并在实际开始工作之前检查日期是否为当月的第一天。 关于ruby-如何每月在Heroku运行一次Scheduler插件？，我们在StackOverflow上找到一个类似的问题： https://stackoverflow.com/questions/8692687/

  9. ruby-on-rails - 如何从 format.xml 中删除 <hash></hash> - 2

    我有一个对象has_many应呈现为xml的子对象。这不是问题。我的问题是我创建了一个Hash包含此数据，就像解析器需要它一样。但是rails自动将整个文件包含在.........我需要摆脱type="array"和我该如何处理？我没有在文档中找到任何内容。 最佳答案 我遇到了同样的问题；这是我的XML:我在用这个:entries.to_xml将散列数据转换为XML，但这会将条目的数据包装到中所以我修改了:entries.to_xml(root:"Contacts")但这仍然将转换后的XML包装在“联系人”中，将我的XML代码修改为

  10. ruby - 如何使用文字标量样式在 YAML 中转储字符串？ - 2

    我有一大串格式化数据(例如JSON)，我想使用Psychinruby​​同时保留格式转储到YAML。基本上，我希望JSON使用literalstyle出现在YAML中:---json:|{"page":1,"results":["item","another"],"total_pages":0}但是，当我使用YAML.dump时，它不使用文字样式。我得到这样的东西:---json:!"{\n\"page\":1,\n\"results\":[\n\"item\",\"another\"\n],\n\"total_pages\":0\n}\n"我如何告诉Psych以想要的样式转储标量？解

随机推荐