小白画浓度分布图求助

Ronin.

本人python小白，最近老师让尝试画一个污染物的浓度分布图，在查找资料和视频后，瞎鼓捣成下面这样。不太懂为什么背景图是紫色的，然后颜色填充的好像也不对，希望有大佬可以帮忙解答一下。
import xarray as xr
import matplotlib.pyplot as plt
import numpy as np
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTER
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus']=False

with open(r'D:\Users\hp\Desktop\2\CN-border-La.dat') as src:
    context = src.read()
    blocks = [cnt for cnt in context.split('>') if len(cnt) > 0]
    borders = [np.fromstring(block, dtype=float, sep=' ') for block in blocks]
proj = ccrs.PlateCarree()
fig = plt.figure(figsize=(10,8),dpi=200)
ax = fig.add_axes([0.5,0.5,0.9,0.9], projection=proj)
ax.set_extent([97,132,17,47])
for line in borders:
    ax.plot(line[0::2], line[1::2], '-', lw=1.5, color='k',
            transform=ccrs.Geodetic())
ds = xr.open_dataset(r'D:\Users\hp\Desktop\1\ECHAP_PM10_V3.nc')
lon = ds.Lon
lat = ds.Lat
pm10 = ds.PM10
fig.canvas.draw()
cf = ax.contourf(lon,lat,pm10,transform=proj,camp='Spectral_r',aloha=0.8)
ax.set_xticks(np.arange(97,132,5))
ax.set_yticks(np.arange(17,47,5),crs=proj)
ax.xaxis.set_major_formatter(LONGITUDE_FORMATTER)
ax.yaxis.set_major_formatter(LATITUDE_FORMATTER)
fig.canvas.draw()
#设置刻度值，画经纬网格
xticks = [97,102,107,112,117,122,127]
yticks = [17,22,27,32,37,42,47]
ax.gridlines(xlocs=xticks, ylocs=yticks,linestyle='--',lw=1,color='dimgrey')
ax.coastlines()
ax.add_feature(cfeature.OCEAN.with_scale('50m'))
ax.add_feature(cfeature.LAND.with_scale('50m'))
plt.title('pm10浓度分布图')
plt.show()

gjlxr

真厉害，佩服佩服

扇贝单词

边缘没有白化呀，可以参考一下这个帖子，4.cartopy的白化
http://bbs.06climate.com/forum.p ... &fromuid=130924
(出处: 气象家园)
并且，可以设置缩小一些地图范围

Ronin.

扇贝单词 发表于 2022-4-11 23:05
边缘没有白化呀，可以参考一下这个帖子，4.cartopy的白化
http://bbs.06climate.com/forum.php?mod=viewth ...

就是感觉那个中间的紫色方块把陆地和海洋的颜色都遮住了，这是没有白化的原因么，因为我看白化好像只显示选中的区域

松鼠算天气

紫色那部分是你的数据范围吧 可以把紫色换成透明色

来自苹果APP客户端

Ronin.

松鼠算天气 发表于 2022-4-11 23:59
紫色那部分是你的数据范围吧 可以把紫色换成透明色

可以问下怎么换么，因为是小白还不太懂

扇贝单词

Ronin. 发表于 2022-4-11 23:11
就是感觉那个中间的紫色方块把陆地和海洋的颜色都遮住了，这是没有白化的原因么，因为我看白化好像只显示 ...

对，白化是框定范围的，你这是相当于加了一层把之前的地图覆盖了，白化可以保留一部分地图不被覆盖，那些紫色部分我感觉都是一种缺值数据，感觉白化应该不要紧，不会影响你的结果

Ronin.

扇贝单词 发表于 2022-4-12 13:46
对，白化是框定范围的，你这是相当于加了一层把之前的地图覆盖了，白化可以保留一部分地图不被覆盖，那些 ...

多谢大佬，紫色部分的确是缺测值，全部替换为nan后紫色就消失了

Ronin.

紫色部分是缺测值，在大佬提醒和询问师兄后加了两句就除掉了紫色，接下来就是加色条和让图变得更好看，学习还是任重而道远
pm10 = ds.PM10.values
pm10=np.where(pm10==-999,np.nan,pm10)

chaoschaos

ddddddddd字数补丁