cartopy白化分享_2

坎坷

是笔者上个帖子的脚本资源，大家支持一下哈！

map_f.py (1.85 KB, 下载次数: 23, 售价: 1 贡献)

2020-8-18 22:25 上传

点击文件名下载附件
maskout_new.py (1.93 KB, 下载次数: 34, 售价: 1 贡献)

2020-8-18 22:25 上传

点击文件名下载附件

songzhh16

请问楼主有没有散点图白化的方法呀

LIUCAS

请问楼主可以免费发一份吗？我的贡献值不够，谢谢啊，我的邮箱是2514330828@qq.com

Divine

'utf-8' codec can't decode byte 0xb1 in position 0: invalid start byte

楼主大大，这个问题怎么解决啊

edwardli

Divine 发表于 2021-8-25 09:30
'utf-8' codec can't decode byte 0xb1 in position 0: invalid start byte

楼主大大，这个问题怎么解决 ...

在mask的reader中传递shp文件的编码。可以通过notepad++打开shp，右下角有编码方式。
比如encoding=‘ANSI’  或者‘GBK’

Divine

edwardli 发表于 2021-8-25 09:51
在mask的reader中传递shp文件的编码。可以通过notepad++打开shp，右下角有编码方式。
比如encoding=‘AN ...

ANSI和GBK都试了还是不行，报出了同样的错误，还有什么编码吗？

edwardli

Divine 发表于 2021-8-25 14:47
ANSI和GBK都试了还是不行，报出了同样的错误，还有什么编码吗？

"可以通过notepad++打开shp，右下角有编码方式。"

Divine

edwardli 发表于 2021-8-25 15:05
"可以通过notepad++打开shp，右下角有编码方式。"

右下角是ANSI，但是内容乱码，写进去也还是同样的错误，愁死我了

edwardli

Divine 发表于 2021-8-25 15:53
右下角是ANSI，但是内容乱码，写进去也还是同样的错误，愁死我了

1. 
   

复制代码

我打开了Province.shp，内容和你的一样。


def shp2clip(originfig,ax,shpfile,region, clabel = None, vcplot = None):
    sf = shapefile.Reader(shpfile,encoding='ANSI')
    vertices = []
    codes = []
    for shape_rec in sf.shapeRecords():
        ####这里需要找到和region匹配的唯一标识符，record[]中必有一项是对应的。
        if shape_rec.record[1] == region:   #####在country1.shp上，对中国以外的其他国家或地区进行maskout
        # if shape_rec.record[7] in region:  #####在bou2_4p.shp上，对中国的某几个省份或地区之外的部分进行maskout
            pts = shape_rec.shape.points
            prt = list(shape_rec.shape.parts) + [len(pts)]
            for i in range(len(prt) - 1):
                for j in range(prt, prt[i+1]):
                    vertices.append((pts[j][0], pts[j][1]))
                codes += [Path.MOVETO]
                codes += [Path.LINETO] * (prt[i+1] - prt -2)
                codes += [Path.CLOSEPOLY]
            clip = Path(vertices, codes)
            plate_carree_transform = ccrs.PlateCarree()._as_mpl_transform(ax)
            clip = PathPatch(clip, transform=plate_carree_transform)
    if vcplot:
        if isinstance(originfig,Iterable):
            for ivec in originfig:
                ivec.set_clip_path(clip)
        else:
            originfig.set_clip_path(clip)
    else:
        for contour in originfig.collections:
        # for contour in originfig:
            contour.set_clip_path(clip)
            print(0)
    if  clabel:
        clip_map_shapely = ShapelyPolygon(vertices)
        for text_object in clabel:
            if not clip_map_shapely.contains(ShapelyPoint(text_object.get_position())):
                text_object.set_visible(False)   
    return clip



clip = shp2clip(cs, ax, 'D:/pycode/shp/Province.shp','河北省')

皮叶克尼斯

感谢楼主
！
maskout文件很有用