底图和轮廓，Python

Question

我有一组卫星数据file here，我为纬度和经度创建了一个网格，为臭氧值创建了一个二维数组。

我知道为了在地图中绘制数据的轮廓，我需要投影坐标，但由于我的网格不是正方形 (144x24)，我无法找到解决方法。我正在覆盖地理区域（0 到 360；-30 到 30），我需要方形像素。

发布的数据很长，但这是我到目前为止的代码，

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import seaborn as sns
from mpl_toolkits.basemap import Basemap, cm
%matplotlib inline

path = '/home/rafaella/month_files_CSV/O3_COLUMNS_MATCHED_fv0005_200306.csv'
df = pd.read_csv(path, skiprows=1)

df = pd.read_csv(path, delim_whitespace=True)
lat = np.array(df['AVG_LAT'])
lon = np.array(df['AVG_LON'])
toc = np.array(df['TROP_COL'])

#new grid for lon[0,360] lat[-30,30]

lomin = 0
lomax = 360
lamin = -30
lamax = 30
stp   = 2.5

loc_lon = np.zeros(int((lomax-lomin)/stp))
loc_lat = np.zeros(int((lamax-lamin)/stp))

for i in range(0,len(loc_lon)):
    loc_lon[i] = i*stp +lomin

for j in range(0,len(loc_lat)):
    loc_lat[j] = j*stp +lamin

mtoc_local = np.zeros((len(loc_lon),len(loc_lat))) 
sdtoc_local = np.zeros((len(loc_lon),len(loc_lat))) 
mtoc_local[:,:] = np.nan      
sdtoc_local[:,:] = np.nan

for i in range (0, len(loc_lon)):
    for j in range (0,len(loc_lat)):
        ix = np.where((lon>=loc_lon[i])& (lat>=loc_lat[j]) & (lon<loc_lon[i]+stp) & (lat<loc_lat[j]+stp))[0] 
        mtoc_local[i,j]=np.nanmean(toc[ix]) 
        sdtoc_local[i,j]=np.nanstd(toc[ix])

fig = plt.figure(figsize=(20, 5))
map = Basemap(llcrnrlon=0,llcrnrlat=-30, urcrnrlon=360.,urcrnrlat=30.,\
        rsphere=(6378137.00,6356752.3142),\
        resolution='l',projection='merc',\
        lat_0=0,lon_0=-30.,lat_ts=30.)

map.drawcoastlines()
# draw parallels
map.drawparallels(np.arange(-30,30,10),labels=[1,1,0,1])
# draw meridians
map.drawmeridians(np.arange(-180,180,20),labels=[1,1,0,1])

map = plt.contourf(loc_lon, loc_lat , mtoc_local.T, vmin=210, vmax=350,  cmap='RdPu')
plt.colorbar(orientation='horizontal', ticks=[200, 220, 240, 260, 280, 300, 320, 340] )
plt.title('Tropical TOC monthly mean 06,2009')

plt.show()

它可以很好地绘制地图或数据，但不能同时绘制两者。这里分别是两者的图像 map

real data

我是python的新手，一个月前开始的，所以对所有的函数和库我还是不熟悉。

Answer 1

您的代码有两个问题。首先，您必须在您的坐标上应用投影，这是使用x,y = map(lon, lat) 完成的。但是，这会在您的情况下引发错误，因为 loc_lon 和 loc_lat 的尺寸不同。您可以传递与 z (mtoc_local.T) 形状相同的数组，而不是将 x 和 y 向量传递给 contourf 函数。您可以使用np.meshgrid 来创建它们。长话短说，用以下三行替换contourf命令的行

X, Y = np.meshgrid(loc_lon, loc_lat)
x,y = map(X,Y)
map = plt.contourf(x, y , mtoc_local.T, vmin=210, vmax=350,  cmap='RdPu')

结果看起来像这样