如何计算区域平均SIC

sic_son = xr.open_dataset('G:/NSIDC-SIC-SH-V4/00-10/00-11.nc')['nsidc_nt_seaice_conc'].sel(time=slice('2010-09','2011-04'))

1. 首先，观察SIC数据，发现没有海冰的地方SIC为0，大陆为nan值。
   
2. 例如，要计算东罗斯海区域的平均海冰情况

#计算区域平均的SIC
time_len = len(sic_son)
ross_sie_nt = np.zeros((time_len),'float')

for i in range(time_len):
#    print(i) 
    mask_ross_nt = (lons>-180.)& (lons<-130.) & (lats<-70.) & (sic_son[i]>=0)
    sic_ross_nt = np.where(mask_ross_nt,sic_son[i],np.nan)
    ross_sie_nt[i] = np.nanmean(sic_ross_nt)

先筛选出东罗斯海的区域，然后其他地方为np.nan。然后在东罗斯海区域内计算区域平均SIC，包括SIC为0的地方。