中国领土的四端为：
最东端：
在黑龙江和乌苏里江的主航道中心线的相交处（135°2′30′′E）；
最西端：
在帕米尔高原附近（73°29′59.79′′E）；
最南端：
在立地暗沙（3°31′00′′N，112°17′09′E），（英语：Lidi Ansha或Lydi Shoal）为中国南海南沙群岛区域的一座暗沙，是实际上的中国领土的最南端（非位于其东北约15海里的曾母暗沙）。按中华人民共和国行政区划，立地暗沙属于海南省三沙市管辖；
最北端：
在漠河以北黑龙江主航道的中心线上（53°33′N，124°20′E）。
中国东西跨越经度60多度，最东端的乌苏里江畔和最西端的帕米尔高原高原相差5个时区。
中国领土辽阔广大，陆地总面积约960万平方千米，仅次于俄罗斯、加拿大，居世界第3位，第四位为美国。中国南北相距5500千米，东西相距5200千米。
EPSG是什么？
European Peroleum Survey Group，欧洲石油勘探组织。对，简单的来说就是它搞的一套标准，大家以他为权威，他收录了很多坐标系，并给这些坐标系分配到详细的“EPSG代码”，比如：EPSG:3857、EPSG:4326、EPSG:4527 等。。。
CGCS2000、北京54、西安80、WGS84 是什么？有什么关系？
详细介绍可以参考这篇文章：
(36条消息) 北京54、西安80、CGCS2000、WGS84坐标系及区别54 80 2000坐标系区别大作家佚名的博客-CSDN博客
现在我国已经停用北京54和西安80坐标系，改用CGCS2000坐标系了。北京54和西安80是参心坐标系，CGCS2000是地心坐标系与WGS84一样，只是椭球扁率微小差异，实际坐标结果几乎一致，约有毫米级别的差异 。
EPSG:3857是什么？
最早是google在web上使用的伪墨卡托投影坐标系（Web Mercator），一开始EPSG组织由于它“伪”拒绝收录，后来google自己定制了一个ID为EPSG:900913（非官方认证的ID），900913-google长得很像。。。，再然后，EPSG官方在2008年收录了并指定代号为EPSG:3785，再然后代号被弃用重新命名为EPSG:3857，也就是现在我们常用的Web 地图的坐标系，并且给定官方命名 “WGS 84 / Pseudo-Mercator“？？？
EPSG:4326 和 EPSG:3857 ？
这俩都是投影，世界级别（WGS 84），只不过EPSG:4326标注的单位是度，EPSG:3857标注的单位是米。
所以 EPSG:4326 叫：“WGS 84”；EPSG:3857 叫：“WGS 84 / Pseudo-Mercator”
举例对比：
EPSG对照表（参照中国领土的四端）：
坐标系对应EPSG代号、经度范围、中央经线 - 豆丁网 (docin.com)
坐标在线转换：
https://epsg.io/transform
转换举例对比：
image.png
CGCS2000 / 3-degree Gauss-Kruger zone 39，解释为：【CGCS2000坐标系（国家2000）、三度带（经度跨度单元）、高斯克吕格（投影算法）、时区 39】对应 EPSG的代号为 4527；
参心坐标系由于参考中心不同、距离参考中心远近不同导致有明显的误差的；
加代号和不加代号少了两位数而已。
扩展研究（选看）：
image.png
image.png
在postgis数据库中可以查询到国内常用的几个坐标系的空间参考系信息：
-- 六度带：zone 属于 [13,23]，CM 属于 [75E, 135E]，数量：11
-- 三度带：zone 属于 [25,45]，CM 属于 [75E, 135E]，数量：21
-- 1. 2000坐标系：
-- CGCS2000 6度带 高斯投影 加代号 srid 属于 [4491,4501]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["CGCS2000 / Gauss-Kruger zone%' order by srid asc;
-- CGCS2000 6度带 高斯投影 不加代号 srid 属于 [4502,4512]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["CGCS2000 / Gauss-Kruger CM%' order by srid asc;
-- CGCS2000 3度带 高斯投影 加代号 srid 属于 [4513,4533]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["CGCS2000 / 3-degree Gauss-Kruger zone%' order by srid asc;
-- CGCS2000 3度带 高斯投影 不加代号 srid 属于 [4534,4554]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["CGCS2000 / 3-degree Gauss-Kruger CM%' order by srid asc;
-- 2. 西安80坐标系：
-- Xian 1980 6度带 高斯投影 加代号 srid 属于 [2327,2337]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Xian 1980 / Gauss-Kruger zone%' order by srid asc;
-- Xian 1980 6度带 高斯投影 不加代号 srid 属于 [2338,2348]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Xian 1980 / Gauss-Kruger CM%' order by srid asc;
-- Xian 1980 3度带 高斯投影 加代号 srid 属于 [2349,2369]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Xian 1980 / 3-degree Gauss-Kruger zone%' order by srid asc;
-- Xian 1980 3度带 高斯投影 不加代号 srid 属于 [2370,2390]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Xian 1980 / 3-degree Gauss-Kruger CM%' order by srid asc;
-- 3. 北京54坐标系：
-- Beijing 1954 6度带 高斯投影 加代号 srid 属于 [21413,21423]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Beijing 1954 / Gauss-Kruger zone%' order by srid asc;
-- Beijing 1954 6度带 高斯投影 不加代号 srid 属于 [21453,21463]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Beijing 1954 / Gauss-Kruger CM%' order by srid asc;
-- Beijing 1954 3度带 高斯投影 加代号 srid 属于 [2401,2421]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Beijing 1954 / 3-degree Gauss-Kruger zone%' order by srid asc;
-- Beijing 1954 3度带 高斯投影 不加代号 srid 属于 [2422,2442]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["Beijing 1954 / 3-degree Gauss-Kruger CM%' order by srid asc;
-- 1. 新北京54坐标系：
-- （新北京54）New Beijing 6度带 高斯投影 加代号 srid 属于 [4568,4578]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["New Beijing / Gauss-Kruger zone%' order by srid asc;
-- （新北京54）New Beijing 6度带 高斯投影 不加代号 srid 属于 [4579,4589]，数量：11
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["New Beijing / Gauss-Kruger CM%' order by srid asc;
-- （新北京54）New Beijing 3度带 高斯投影 加代号 srid 属于 [4652,4656]并[4766,4781]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["New Beijing / 3-degree Gauss-Kruger zone%' order by srid asc;
-- （新北京54）New Beijing 3度带 高斯投影 不加代号 srid 属于 [4782,4800]并[4812]并[4822]，数量：21
SELECT * FROM "spatial_ref_sys" where srtext like 'PROJCS["New Beijing / 3-degree Gauss-Kruger CM%' order by srid asc;
附：python根据给定的投影坐标x,y值，自动识别epsg代码：
适用于 CGCS2000，Xian80，Beijing54，NewBeijing54 四大crs
from math import inf
from pprint import pprint
import pyproj
# 中国领土的四端为：
# 最东端：
# 在黑龙江和乌苏里江的主航道中心线的相交处（135°2'30"E）；
# 最西端：
# 在帕米尔高原附近（73°29'59.79"E）；
# 最南端：
# 在立地暗沙（3°31'00"N，112°17'09"E），（英语：Lidi Ansha或Lydi Shoal）为中国南海南沙群岛区域的一座暗沙，是实际上的中国领土的最南端（非位于其东北约15海里的曾母暗沙）。按中华人民共和国行政区划，立地暗沙属于海南省三沙市管辖；
# 最北端：
# 在漠河以北黑龙江主航道的中心线上（53°33'N，124°20'E）。
# minx, miny, maxx, maxy = (73°29'59.79"E, 3°31'00"N, 135°2'30"E, 53°33'N)
# # 3857
# minx, miny, maxx, maxy = (8181976.079668645, 391719.56668402057, 15032769.569208318, 7085388.165495079)
# # 4326
ChinaBbox4326 = [73.4999417, 3.5166667, 135.0416667, 53.55]
# # cgcs2000-4491
# minx, miny, maxx, maxy = (13333307.577437889, 388991.6394066638, 17137769.453445204, 7744052.21995985)
# # cgcs2000-4501
# minx, miny, maxx, maxy = (14778679.655726733, 820784.1522103351, 23502761.70636996, 5935831.474453285)
# # cgcs2000-3degree-4513
# minx, miny, maxx, maxy = (25333307.57743789, 388991.6394066638, 29137769.453445204, 7744052.21995985)
# # cgcs2000-3degree-4533
# minx, miny, maxx, maxy = (36778679.65572673, 820784.1522103351, 45502761.70636996, 5935831.474453285)
# class SpatialRefSys:
# srid: int = 4326
# auth_name: str = 'EPSG'
# auth_srid: int = 4326
# srtext: str = 'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]'
# proj4text: str = '+proj=longlat +datum=WGS84 +no_defs '
# 已知 web墨卡托
WebMercator3857 = 'WebMercator3857'
# 已知 wgs84 4326
WGS84 = 'WGS84'
# 已知 中国四大坐标系得 zone 和 cm 范围：
CGCS2000 = 'CGCS2000'
Xian80 = 'Xian80'
Beijing54 = 'Beijing54'
NewBeijing54 = 'NewBeijing54'
zoneRangesDegree6, zoneRangesDegree3 = (range(13, 23 + 1, 1), range(25, 45 + 1, 1))
zoneRangesDegree6Dict, zoneRangesDegree3Dict = ({item: index for index, item in enumerate(zoneRangesDegree6)},
{item: index for index, item in enumerate(zoneRangesDegree3)})
cmRangesDegree6, cmRangesDegree3 = (range(75, 135 + 6, 6), range(75, 135 + 3, 3))
chinaSrsDegree6List1323: list = []
chinaSrsDegree3List2545: list = []
class ChinaSrs:
crs: str = 'CGCS2000' # 坐标系
isZone: bool = False # 是否加代号
is3DegreeBelt: bool = False # 是否3度带
epsg: str = 'EPSG:4491' # epsg代码
zone: int = 13 # 时区/代号
cm: str = '75E' # 中央经线
convertCoordinate: list[float] = [0.0, 0.0] # 转换结果
def __init__(self, crs, isZone, is3DegreeBelt, epsg, zone, cm):
self.crs = crs
self.isZone = isZone
self.is3DegreeBelt = is3DegreeBelt
self.epsg = epsg
self.zone = zone
self.cm = cm
def __repr__(self) -> str:
return f"(crs={self.crs!r}, isZone={self.isZone!r}, is3DegreeBelt={self.is3DegreeBelt!r}, epsg={self.epsg!r}, zone={self.zone!r}, cm={self.cm!r}, convertCoordinate={self.convertCoordinate!r})"
def initCrs():
if len(chinaSrsDegree6List1323) > 0:
print('===》chinaSrs 缓存...')
return
print('===》chinaSrs 初始化...')
is3DegreeBelt = False
for index, zone in enumerate(zoneRangesDegree6):
centralMeridian = cmRangesDegree6.start + cmRangesDegree6.step * index
chinaSrsDegree6List1323.append((
# CGCS2000 6度带 高斯投影 加代号 srid 属于 [4491,4501]，数量：11
# CGCS2000 6度带 高斯投影 不加代号 srid 属于 [4502,4512]，数量：11
ChinaSrs(crs=CGCS2000, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{4491 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=CGCS2000, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{4502 + index}', zone=zone,
cm=f'{centralMeridian}E'),
# Xian 1980 6度带 高斯投影 加代号 srid 属于 [2327,2337]，数量：11
# Xian 1980 6度带 高斯投影 不加代号 srid 属于 [2338,2348]，数量：11
ChinaSrs(crs=Xian80, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{2327 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=Xian80, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{2338 + index}', zone=zone,
cm=f'{centralMeridian}E'),
# Beijing 1954 6度带 高斯投影 加代号 srid 属于 [21413,21423]，数量：11
# Beijing 1954 6度带 高斯投影 不加代号 srid 属于 [21453,21463]，数量：11
ChinaSrs(crs=Beijing54, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{21413 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=Beijing54, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{21453 + index}', zone=zone,
cm=f'{centralMeridian}E'),
# （新北京54）New Beijing 6度带 高斯投影 加代号 srid 属于 [4568,4578]，数量：11
# （新北京54）New Beijing 6度带 高斯投影 不加代号 srid 属于 [4579,4589]，数量：11
ChinaSrs(crs=NewBeijing54, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{4568 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=NewBeijing54, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{4579 + index}',
zone=zone,
cm=f'{centralMeridian}E')
))
is3DegreeBelt = True
for index, zone in enumerate(zoneRangesDegree3):
centralMeridian = cmRangesDegree3.start + cmRangesDegree3.step * index
# 新北京54坐标系-3度带-加代号特殊处理
epsgNewBeijing54Zone = 4652 + index
if epsgNewBeijing54Zone >= 4657:
epsgNewBeijing54Zone = epsgNewBeijing54Zone - 4657 + 4766
# 新北京54坐标系-3度带-不加代号特殊处理
epsgNewBeijing54 = 4782 + index
if epsgNewBeijing54 == 4801:
epsgNewBeijing54 = 4812
elif epsgNewBeijing54 == 4802:
epsgNewBeijing54 = 4822
chinaSrsDegree3List2545.append((
# CGCS2000 3度带 高斯投影 加代号 srid 属于 [4513,4533]，数量：21
# CGCS2000 3度带 高斯投影 不加代号 srid 属于 [4534,4554]，数量：21
ChinaSrs(crs=CGCS2000, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{4513 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=CGCS2000, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{4534 + index}', zone=zone,
cm=f'{centralMeridian}E'),
# Xian 1980 3度带 高斯投影 加代号 srid 属于 [2349,2369]，数量：21
# Xian 1980 3度带 高斯投影 不加代号 srid 属于 [2370,2390]，数量：21
ChinaSrs(crs=Xian80, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{2349 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=Xian80, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{2370 + index}', zone=zone,
cm=f'{centralMeridian}E'),
# Beijing 1954 3度带 高斯投影 加代号 srid 属于 [2401,2421]，数量：21
# Beijing 1954 3度带 高斯投影 不加代号 srid 属于 [2422,2442]，数量：21
ChinaSrs(crs=Beijing54, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{2401 + index}', zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=Beijing54, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{2422 + index}', zone=zone,
cm=f'{centralMeridian}E'),
# （新北京54）New Beijing 3度带 高斯投影 加代号 srid 属于 [4652,4656]并[4766,4781]，数量：21
# （新北京54）New Beijing 3度带 高斯投影 不加代号 srid 属于 [4782,4800]并[4812]并[4822]，数量：21
ChinaSrs(crs=NewBeijing54, isZone=True, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{epsgNewBeijing54Zone}',
zone=zone,
cm=f'{centralMeridian}E'),
ChinaSrs(crs=NewBeijing54, isZone=False, is3DegreeBelt=is3DegreeBelt, epsg=f'EPSG:{epsgNewBeijing54}',
zone=zone,
cm=f'{centralMeridian}E')
))
# print(chinaSrsDegree6List1323)
print('------------------------------')
# print(chinaSrsDegree3List2545)
results: list[ChinaSrs] = []
cgcs2000Results = []
xian80Results = []
beijing54Results = []
newbeijing54Results = []
def identify(_x: float, _y: float, zoneNum, isZone):
if zoneRangesDegree6.start <= zoneNum <= zoneRangesDegree6.stop - 1:
# print('==>六度带')
tmpZoneRangesDict = zoneRangesDegree6Dict
tmpChinaSrsList = chinaSrsDegree6List1323
elif zoneRangesDegree3.start <= zoneNum <= zoneRangesDegree3.stop - 1:
# print('==>三度带')
tmpZoneRangesDict = zoneRangesDegree3Dict
tmpChinaSrsList = chinaSrsDegree3List2545
else:
# print('==>无解，包括x == 24')
return
if zoneNum in tmpZoneRangesDict:
index = tmpZoneRangesDict.get(zoneNum)
chinaSrsTuple: tuple[ChinaSrs] = tmpChinaSrsList[index]
for item in chinaSrsTuple:
if not item.isZone == isZone:
# 如果要求识别 isZone ，则跳过不加代号
# 如果不要求识别 isZone ，则跳过加代号
continue
try:
transformer = pyproj.Transformer.from_crs(item.epsg, "EPSG:4326", always_xy=True)
other_x, other_y = transformer.transform(xx=_x, yy=_y)
if xyTransferValid(xt=other_x, yt=other_y) is False:
continue
item.convertCoordinate = [other_x, other_y]
if item.crs == CGCS2000:
cgcs2000Results.append(item)
elif item.crs == Xian80:
xian80Results.append(item)
elif item.crs == Beijing54:
beijing54Results.append(item)
elif item.crs == NewBeijing54:
newbeijing54Results.append(item)
else:
continue
except Exception as e:
print('e=%s' % e)
continue
# 判断转换后坐标的有效性
def xyTransferValid(xt, yt):
if xt == inf or yt == inf:
return False
elif ChinaBbox4326[0] <= xt <= ChinaBbox4326[2] and ChinaBbox4326[1] <= yt <= ChinaBbox4326[3]:
return True
else:
return False
# 根据坐标计算epsg代号
def identify_crs(_x: float, _y: float, crsName: str = 'unlimited', degreeBelt: str = 'unlimited'):
"""
Args:
_x:
_y:
crsName: 坐标系名称，可选为 'unlimited'、'WebMercator3857'、'CGCS2000'、'Xian80'、'Beijing54'、'NewBeijing54'，默认 'unlimited' 不限
degreeBelt: 坐标系度带，可选为 '6'、'3'、'unlimited'，默认 'unlimited' 不限
Returns: []
"""
if crsName not in (WebMercator3857, CGCS2000, Xian80, Beijing54, NewBeijing54):
crsName = 'unlimited'
if degreeBelt not in ('6', '3'):
degreeBelt = 'unlimited'
print(f'(坐标X值：{_x}，坐标Y值：{_y}，坐标系：{crsName}，度带：{degreeBelt})')
_vo = []
try:
# xInt, yInt = int(_x), int(_y)
xLen, yLen = len(str(int(_x))), len(str(int(_y)))
# 1. 3857坐标系识别
epsg3857Result = None
transformer = pyproj.Transformer.from_crs('EPSG:3857', "EPSG:4326", always_xy=True)
other_x, other_y = transformer.transform(xx=_x, yy=_y)
if xyTransferValid(xt=other_x, yt=other_y) is False:
pass
else:
epsg3857Result = {
'crs': WebMercator3857,
'epsg': 'EPSG:3857',
'convertCoordinate': [other_x, other_y]
}
# 1-1. 3857 坐标特殊处理返回
if crsName == WebMercator3857:
_vo.append(epsg3857Result)
return _vo
# 2-0. 初始化 chinaSrs
initCrs()
# 2-1. 中国四大坐标系识别
if xLen == 8 and yLen == 7:
# 加代号
# print('==>加代号')
zoneNum = int(_x / 1000000)
identify(_x, _y, zoneNum, True)
else:
# 不加代号
# print('==>不加代号')
for zoneNum in zoneRangesDegree6:
identify(_x, _y, zoneNum, False)
for zoneNum in zoneRangesDegree3:
identify(_x, _y, zoneNum, False)
# 整合结果集合，按 crs 分组
results.extend(cgcs2000Results)
results.extend(xian80Results)
results.extend(beijing54Results)
results.extend(newbeijing54Results)
# 2-2. 中国四大坐标统一处理返回
if crsName != 'unlimited' and degreeBelt != 'unlimited':
# 国家2000 6度带
for item in results:
if item.crs == crsName and item.is3DegreeBelt == (False if degreeBelt == '6' else True):
_vo.append(item)
elif crsName != 'unlimited' and degreeBelt == 'unlimited':
# 国家2000 不限
for item in results:
if item.crs == crsName:
_vo.append(item)
elif crsName == 'unlimited' and degreeBelt != 'unlimited':
# 不限 6度带
for item in results:
if item.is3DegreeBelt == (False if degreeBelt == '6' else True):
_vo.append(item)
else:
# 不限 不限
_vo = results
if epsg3857Result is not None:
_vo.append(epsg3857Result)
return _vo
except Exception as e:
print('e=%s' % e)
return _vo
def getOne(_vo):
_voDict = {}
if len(_vo) > 0:
for item in _vo:
if type(item) == ChinaSrs:
_voDict[item.crs] = item.epsg
else:
_voDict[item['crs']] = item['epsg']
else:
return None
if WebMercator3857 in _voDict:
return _voDict.get(WebMercator3857)[5:]
elif CGCS2000 in _voDict:
return _voDict.get(CGCS2000)[5:]
elif Beijing54 in _voDict:
return _voDict.get(Beijing54)[5:]
elif Xian80 in _voDict:
return _voDict.get(Xian80)[5:]
elif NewBeijing54 in _voDict:
return _voDict.get(NewBeijing54)[5:]
else:
return None
if __name__ == '__main__':
# 坐标系自动识别逻辑：
# 已知中心点：
x, y = (39500000.1234, 3500000.5678)
# x, y = (38471609.56449664,3918692.0474597644)
# x, y = (500000, 3500000)
# x, y = (13033287.131065886,4080601.624207251)
vo: list = []
vo = identify_crs(_x=x, _y=y, crsName='unlimited', degreeBelt='unlimited')
print('1. 获取可能的空间参考系epsg代码列表为：')
pprint(vo)
print('---------------------------------------------------------')
print('2. 获取一个可能的空间参考系epsg代码为：')
pprint(getOne(vo))
vo.clear()
#
# x, y = (39440055.461313486,3855857.7584496653)
# vo = identify_crs(_x=x, _y=y, crsName='unlimited', degreeBelt='unlimited')
# pprint(vo)
# pprint(getOne(vo))
# vo.clear()
#
# x, y = (38471609.56449664, 3918692.0474597644)
# vo = identify_crs(_x=x, _y=y, crsName='unlimited', degreeBelt='unlimited')
# pprint(vo)
# pprint(getOne(vo))
# vo.clear()
# vo = identify_crs(_x=x, _y=y, crsName='unlimited', degreeBelt='unlimited')
# vo = identify_crs(_x=x, _y=y, crsName='WebMercator3857', degreeBelt='unlimited')
# vo = identify_crs(_x=x, _y=y, crsName='CGCS2000', degreeBelt='unlimited')
# vo = identify_crs(_x=x, _y=y, crsName='CGCS2000', degreeBelt='6')
# vo = identify_crs(_x=x, _y=y, crsName='CGCS2000', degreeBelt='3')
# vo = identify_crs(_x=x, _y=y, crsName='Xian80', degreeBelt='unlimited')
# vo = identify_crs(_x=x, _y=y, crsName='Xian80', degreeBelt='6')
# vo = identify_crs(_x=x, _y=y, crsName='Xian80', degreeBelt='3')
# vo = identify_crs(_x=x, _y=y, crsName='Beijing54', degreeBelt='unlimited')
# vo = identify_crs(_x=x, _y=y, crsName='Beijing54', degreeBelt='6')
# vo = identify_crs(_x=x, _y=y, crsName='Beijing54', degreeBelt='3')
# vo = identify_crs(_x=x, _y=y, crsName='NewBeijing54', degreeBelt='unlimited')
# vo = identify_crs(_x=x, _y=y, crsName='NewBeijing54', degreeBelt='6')
# vo = identify_crs(_x=x, _y=y, crsName='NewBeijing54', degreeBelt='3')
# vo = identify_crs(_x=x, _y=y, crsName='unlimited', degreeBelt='6')
# vo = identify_crs(_x=x, _y=y, crsName='unlimited', degreeBelt='3')
# pprint(vo)
# pprint(getOne(vo))
# 4326 -180, -90, 180, 90
# 3857 -20037508.342789244, -242528680.94374272, 20037508.342789244, 242528680.94374272
# 国家2000六度带加代号 4491 -13500000, -10001965.729230462, 13500000, 10001965.729230462
# 国家2000六度带加代号 4501 -23500000, -10001965.729230462, 23500000, 10001965.729230462
# 国家2000三度带加代号 4513 -25500000, -10001965.729230462, 25500000, 10001965.729230462
# 国家2000三度带加代号 4533 -45500000, -10001965.729230462, 45500000, 10001965.729230462
pass
测试：
image.png
结果：
image.png