python basemap readshapefile二三事
|
今天要用到basemap读取shp文件报错，查了很多资料，都没有解决。 先是： fig,ax = plt.subplots(figsize=(15,10)) from mpl_toolkits.basemap import Basemap m = Basemap(llcrnrlon=-128.,llcrnrlat=30.,urcrnrlon=-110.,urcrnrlat=43.,projection=‘merc‘,lat_1=20.,lat_2=40.,lon_0=-60.,resolution =‘l‘,area_thresh=1000.,ax = ax) fm = np.vectorize(m,otypes=[np.int32,np.int32]) x,y = fm(housing.longitude,housing.latitude) m.scatter(x,y,s=housing[‘population‘]/50,alpha=0.4,label=‘population‘,c=housing[‘median_house_value‘],cmap=plt.get_cmap(‘jet‘)) m.readshapefile(‘World_GIS_data/country‘,‘country‘)
成功了。 可再用另一中国地图来加载时就出错了： 报错： ValueError: shapefile must have lat/lon vertices - it looks like this one has vertices in map projection coordinates. You can convert the shapefile to geographic coordinates using the shpproj utility from the shapelib tools (http://shapelib.maptools.org/shapelib-tools.html)我 我按着提示去找了shapelib，源码下载到了，需要编译。 我也没有vc++环境，太麻烦了，也不知道能成不。 想了好多办法都没成功，包括用pyshpelib，这个也还是要用到上面的shapelib。 然后找到了： https://github.com/matplotlib/basemap/issues/285 说可以用geopandas来处理。 然后就是安装，这又是一个坑。 安装数次报错，找不着门道。 最后看到一篇文章。该作者可能也被坑惨了，总结写的很好（墙裂推荐）： https://geoffboeing.com/2014/09/using-geopandas-windows/ 跟着上面这篇文章的终于将geopandas安装好了 然后改shapefile的crs为空字典。本来以为就好了，结果还是不行！！！ 还是不行！ ? 然后又看到这篇下面的解答，提到了ogr2ogr https://stackoverflow.com/questions/54508676/convert-a-shapefile-to-geographic-coordinates 弄半天才明白的是安装的GDAL中的可执行文件。 后面又找到一个库，调用的ogr2ogr: 地址： https://github.com/charlie1347/basemaps 最后：
大功告成！！！ 最后在解决这个问题过程中发现一个资源： http://geoportal.statistics.gov.uk/datasets/nuts-level-1-january-2018-full-extent-boundaries-in-the-united-kingdom/data （编辑：李大同） 【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容! |
