在4D numpy网格数据数组中查找不规则区域（lat / lon）

I have a large 4-dimensional dataset of Temperatures [time,pressure,lat,lon]. I need to find all grid points within a region defined by lat/lon indices and calculate an average over the region to leave me with a 2-dimensional array.I have a large 4-dimensional dataset of Tempera