- 积分
- 3442
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-7-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
有没有大神知道哇~~~
按照萝卜大佬的步骤做了白化,但是结果出来之后只有地图没有等值线图了,哭泣
附上代码
fig = plt.figure()
proj = ccrs.PlateCarree(central_longitude=105)
ax = fig.add_subplot(111,projection=proj)
leftlon, rightlon, lowerlat, upperlat = (73,135,10,70)
ax.set_extent([leftlon, rightlon, lowerlat, upperlat, crs=ccrs.PlateCarree())
area = shpreader.Reader('liaoning.dbf').geometries()
ax.add_geometries(area, ccrs.PlateCarree(),facecolor='none', edgecolor='black',zorder = 1)
ax.set_xticks(np.arange(leftlon,rightlon,0.1), crs=ccrs.PlateCarree())
ax.set_yticks(np.arange(lowerlat,upperlat,0.1), crs=ccrs.PlateCarree())
lon_formatter = cticker.LongitudeFormatter()
lat_formatter = cticker.LatitudeFormatter()
ax.xaxis.set_major_formatter(lon_formatter)
ax.yaxis.set_major_formatter(lat_formatter)
m = ax.contourf(lon,lat, I, levels=np.arange(0,3.2,0.5), zorder=0, extend = 'both',transform=ccrs.PlateCarree(), cmap=cmaps.MPL_YlGn)
clip = maskout.shp2clip(m,ax,'liaoning.shp',210000)
plt.colorbar(m)
plt.show()
|
|