Python画色斑图和等值线的问题

yx1230

研究了两天，终于用Python结合matplotlib画出了色斑图和等值线，使用shape进行白化，但是遇到个问题，等值线裁剪后，部分等值线的标注不能显示裁剪后的等值线上，大家都是怎么处理的？

yx1230

补下裁剪前后裁剪后的效果

yx1230

当未来

你要先画等值线，再标注，再裁剪才行吧。。不知道你裁剪怎么裁的

yx1230

当未来 发表于 2016-10-10 15:04
你要先画等值线，再标注，再裁剪才行吧。。不知道你裁剪怎么裁的

你有试过吗，画等值，裁剪吗。我的就是先裁剪后，再标注的

阿阿飞飞

有个解决办法，就是生成一个比现在的shape文件范围稍微大一点的shape文件，然后判断mask掉这个shape文件以外的数据，然后再绘制，再用老的shape文件白化

heidicat

from shapely.geometry import Polygon as ShapelyPolygon
from shapely.geometry import Point as ShapelyPoint
import shapefile
cn =m.contour(olon,olat,data_new,levels=levels,colors="k",linewidths=3)
cn_label=plt.clabel(cn,cn.levels,inline=True,fmt="%.1f",fontsize=20)
        sf = shapefile.Reader(r'files\CHN_adm1.shp') # note, no suffix, all 3 files are used
        recs = sf.records()
        shapes = sf.shapes()
        for i in range(len(recs)):
                if recs[4]=='Fujian':
                        fj_area=ShapelyPolygon(shapes.points)
                        break
        for text_object in cn_label:
                #print (text_object.get_position())
                if not fj_area.contains(ShapelyPoint(text_object.get_position())):
                        text_object.set_visible(False)

############################
把剪裁的地图区域读到polygon 然后判断label的位置是否在地图区域内
如果不是就设置为不可见

TangYh

可以的zsbd