寻找一种方法，可以根据给定的纬度和经度（如果必要还有国际代码），获取电话区号列表。请注意，我所说的不是国际拨号前缀，而是其中的区号。
例如，科罗拉多州丹佛市由303和720区号覆盖。它位于39.739 -104.985，属于NANP 1。因此，给定39.739，-104.985，1，我想得到[303,720]。
图书馆、网络服务、数据库或需要解析为数据库的原始数据，例如形状点的网页，都可以接受，覆盖面越广越好，但仅限于NANP 1也会很有帮助。
请注意，我已经使用了MaxMind，并可以将纬度和经度转换为虚假IP，并将其用作查找键，但MaxMind声称仅支持美国地区码（无论它们是否真正意味着美国，或者实际上是NANP，我尚未测试），似乎每个位置仅限一个区号（例如丹佛市仅限303）。因此这是一种可能性，但并不是太好的选择。
更新: 我找到了一些更相关的信息，但没有明确的解决方案，所以我在这里列出来，而不是在答案中：
例如，科罗拉多州丹佛市由303和720区号覆盖。它位于39.739 -104.985，属于NANP 1。因此，给定39.739，-104.985，1，我想得到[303,720]。
图书馆、网络服务、数据库或需要解析为数据库的原始数据，例如形状点的网页，都可以接受，覆盖面越广越好，但仅限于NANP 1也会很有帮助。
请注意，我已经使用了MaxMind，并可以将纬度和经度转换为虚假IP，并将其用作查找键，但MaxMind声称仅支持美国地区码（无论它们是否真正意味着美国，或者实际上是NANP，我尚未测试），似乎每个位置仅限一个区号（例如丹佛市仅限303）。因此这是一种可能性，但并不是太好的选择。
更新: 我找到了一些更相关的信息，但没有明确的解决方案，所以我在这里列出来，而不是在答案中：
我能够找到两个美国数据库http://www.area-codes.com/area-code-database.asp 和 http://www.nationalnanpa.com/area_codes/index.html （页面下方50%，MS Access文件）。前者包括$450的纬度/经度信息，而后者需要最近邻匹配，就像KeithS所说的那样（它可能是NANPA City Query发现的相同DB）。
此外，我发现Teleatlas拥有区号边界地图，ESRI随ArcGIS副本一起提供区号形状文件。 Maponics似乎有可用的数据：Maponics' 的数据在Google Maps上实现，网址为http://www.usnaviguide.com/areacode.htm。