UTM坐标的遥感图，Python终于把它扶正啦(手动狗头)

Value  Description
cea  Cylindrical Equal Area
mbtfpq  McBryde-Thomas Flat-Polar Quartic
aeqd  Azimuthal Equidistant
sinu  Sinusoidal
poly  Polyconic
omerc  Oblique Mercator
gnom  Gnomonic
moll  Mollweide
lcc  Lambert Conformal
tmerc  Transverse Mercator
nplaea  North-Polar Lambert Azimuthal
gall  Gall Stereographic Cylindrical
npaeqd  North-Polar Azimuthal Equidistant
mill  Miller Cylindrical
merc  Mercator
stere  Stereographic
eqdc  Equidistant Conic
rotpole  Rotated Pole
cyl  Cylindrical Equidistant
npstere  North-Polar Stereographic
spstere  South-Polar Stereographic
hammer  Hammer
geos  Geostationary
nsper  Near-Sided Perspective
eck4  Eckert IV
aea  Albers Equal Area
kav7  Kavrayskiy VII
spaeqd  South-Polar Azimuthal Equidistant
ortho  Orthographic
cass  Cassini-Soldner
vandg  van der Grinten
laea  Lambert Azimuthal Equal Area
splaea  South-Polar Lambert Azimuthal
robin  Robinson

rsphere radius of the sphere used to define map projection (default6370997 meters, close to the arithmetic mean radius of the earth). If given as a sequence, the first two elements are interpreted as the radii of the major and minor axes of an ellipsoid. Note: sometimes an ellipsoid is specified by the major axis and an inverse flattening parameter (if). The minor axis (b) can be computed from the major axis (a) and the inverse flattening parameter using the formula if = a/(a-b).

no_rot only used by oblique mercator. If settoTrue, the map projection coordinates will not be rotated totrue North. DefaultisFalse (projection coordinates are automatically rotated).

patches = [PolygonPatch(feature, edgecolor="red", facecolor="none", linewidth=2) for feature in features]
ax.add_collection(mpl.collections.PatchCollection(patches, match_original=True))

p1 = pyproj.Proj(init='epsg:32650')
for i in range(len(projlatlon)):
    coor[i,0],coor[i,1]= p1(projlatlon[i,0],projlatlon[i,1],inverse=True)