python读取shp地图文件

Jinne

你好，为啥我在绘制浙江省地图的时候没法正确白化，直接全图变白，其他省都ok

jl2587t

Jinne 发表于 2021-10-13 22:11
你好，为啥我在绘制浙江省地图的时候没法正确白化，直接全图变白，其他省都ok

{:eb513:}不知道

yuejun0026

import numpy as np
import xarray as xr
import cartopy.crs as ccrs
import cartopy.feature as cfeat
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
from cartopy.io.shapereader import Reader
import matplotlib.pyplot as plt
import matplotlib.ticker as mticker
import maskout
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter
shp_path='/home/yuejun0026/shp/Sichuan/Sichuan_city.shp'
proj= ccrs.PlateCarree()  # 简写投影
fig = plt.figure(figsize=(4, 4), dpi=400)  # 创建画布
ax = fig.subplots(1, 1, subplot_kw={'projection': proj})  # 创建子图
extent=[96,110,26,35]#限定绘图范围
reader = Reader(shp_path)
sichuancity = cfeat.ShapelyFeature(reader.geometries(), proj, edgecolor='k', facecolor='none')
ax.add_feature(sichuancity, linewidth=0.7)#添加市界细节
ax.set_extent(extent, crs=proj)

大佬  不知道你在使用shp文件过程中有没有遇到这个情况勒？前面的是我的程序
---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
<ipython-input-4-ad6d726e4bd0> in <module>
      1 extent=[96,110,26,35]#限定绘图范围
      2 reader = Reader(shp_path)
----> 3 enshicity = cfeat.ShapelyFeature(reader.geometries(), proj, edgecolor='k', facecolor='none')
      4 ax.add_feature(enshicity, linewidth=0.7)#添加市界细节
      5 ax.set_extent(extent, crs=proj)

~/miniconda3/envs/python/lib/python3.8/site-packages/cartopy/feature/__init__.py in __init__(self, geometries, crs, **kwargs)
    219         """
    220         super(ShapelyFeature, self).__init__(crs, **kwargs)
--> 221         self._geoms = tuple(geometries)
    222
    223     def geometries(self):

~/miniconda3/envs/python/lib/python3.8/site-packages/cartopy/io/shapereader.py in geometries(self)
    167             # Skip the shape that can not be represented as geometry.
    168             if shape.shapeType != shapefile.NULL:
--> 169                 yield sgeom.shape(shape)
    170
    171     def records(self):

~/miniconda3/envs/python/lib/python3.8/site-packages/shapely/geometry/geo.py in shape(context)
     91     True
     92     """
---> 93     if hasattr(context, "__geo_interface__"):
     94         ob = context.__geo_interface__
     95     else:

~/miniconda3/envs/python/lib/python3.8/site-packages/shapefile.py in __geo_interface__(self)
    507
    508         else:
--> 509             raise Exception('Shape type "%s" cannot be represented as GeoJSON.' % SHAPETYPE_LOOKUP[self.shapeType])
    510
    511     @staticmethod

KeyError: -1091265284

zzj

求一份各省shp文件，需要精确到市县的，麻烦楼主了。786535492@qq.com

jl2587t

Jinne 发表于 2021-10-13 22:11
你好，为啥我在绘制浙江省地图的时候没法正确白化，直接全图变白，其他省都ok

这个我不清楚哎，不行你就看看我的代码，套用一下，我画过浙江的

□_□

求吉林省的SHP，能精确到市县的更好，麻烦楼主了，谢谢！。460738506@qq.com

jl2587t

□_□ 发表于 2021-11-8 16:32
求吉林省的SHP，能精确到市县的更好，麻烦楼主了，谢谢！。

给你了

□_□

jl2587t 发表于 2021-11-9 19:37
给你了

谢谢！

那那耶

您好，全国各省市shp文件可以分享一下嘛，谢谢
邮箱：Xxxxuena@163.com

jl2587t

那那耶 发表于 2021-11-18 15:49
您好，全国各省市shp文件可以分享一下嘛，谢谢
邮箱：

发给你啦