from scipy.interpolate import interp2d f = interp2d(longw,latw,Zw)
where longw,latw,Zw are the world coordinates of the NETCDF file and the 2D variable respectively.
Then I applied f(lon,lat)
to all my irregular observations y iteration to interpolate. The lon and lat are the observed coordinates of the local region.