将形状文件转换为地理坐标

Question

对geospatial数据非常陌生-我目前正在努力获得一组在其正确地理区域内标记的邮政编码，如NUTS1所定义。我想最好的办法是：

1. 得到NUTS1区域的坐标
2. 把它们转换成多边形
3. 获取邮政编码的坐标并查看它位于哪个多边形内

然而，我在第（1）步中遇到了一个障碍——从我收集的数据来看，国家统计局关于坚果区域的数据不是经纬度格式的（见下面的来源和错误）。根据类似的问题，我发现ogr2ogr是转换的方法-所以我下载了OSGeo4W，它有{}，但是我并不完全是我应该在我的OSGeo4Wpowershell中运行以进行适当的转换。在

我也参考了下面的备忘单，但还是没有运气。在

我想geopandas对我的问题有一个非常简单的答案-但是，在尝试安装这个软件包两天后，我仍然没有任何运气，所以我暂时退出了这条路（另一条注意事项是，如果有人有幸用python3.7安装geopandas，一定要给我指出正确的方向。我尝试了conda安装和下载所有的whl文件单独，似乎都没有工作）。在

我想知道是否有人能为一个新手说清楚我需要运行的命令来解决我的问题。在

提前谢谢！在

数据源：http://geoportal.statistics.gov.uk/datasets/nuts-level-1-january-2018-full-extent-boundaries-in-the-united-kingdom/data

Error: ValueError: shapefile must have lat/lon vertices

看起来这个有顶点 地图投影坐标。可以将形状文件转换为地理位置 使用shapelib工具中的shpproj实用程序进行坐标测量 （http://shapelib.maptools.org/shapelib-tools.html）

备忘单：http://www.bostongis.com/PrinterFriendly.aspx?content_name=ogr_cheatsheet

Answer 1

第一步是使用ogr2ogrcli程序的-t_srs选项将原始数据转换为无处不在的WGS84坐标系。阅读有关ogr2ogrhere的所有当前选项的详细信息。在

ogr2ogr boundaries.shp NUTS_Level_1_January_2018_Full_Extent_Boundaries_in_the_United_Kingdom.shp -t_srs epsg:4326

数据集已经是“ST_Polygon”形式，因此第2步没有实际意义。在

我将把它作为练习留给您自己来解决步骤3；提示：您可以使用-sql -dialect sqlite与ogr2ogr进行空间连接。在