1.绘制兰勃脱投影的中国区域(包含南海子图):
import numpy as np
import xarray as xr
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from copy import copy
from cartopy.mpl.gridliner import LATITUDE_FORMATTER, LONGITUDE_FORMATTER
import shapely.geometry as sgeom
def find_side(ls, side):
"""
Given a shapely LineString which is assumed to be rectangular, return the
line corresponding to a given side of the rectangle.
"""
minx, miny, maxx, maxy = ls.bounds
points = {left: [(minx, miny), (minx, maxy)],
ight: [(maxx, miny), (maxx, maxy)],
ottom: [(minx, miny), (maxx, miny)],
op: [(minx, max