经纬度(121.45,31.21)转为2000坐标

将经纬度坐标转换为2000国家大地坐标系坐标的方法是：

1. 将经纬度坐标转换为WGS84坐标系下的三维坐标（X，Y，Z），其中Z轴为指向地心的轴。

2. 将WGS84坐标系三维坐标转换为2000国家大地坐标系下的三维坐标（X，Y，Z）。

3. 根据2000国家大地坐标系的定义，将三维坐标的X、Y、Z分别乘以比例因子，得到2000国家大地坐标系下的平面坐标（Easting，Northing）。

具体实现过程如下：

1. 将经纬度坐标转换为WGS84坐标系下的三维坐标（X，Y，Z）。

可以使用 pyproj 库实现这个步骤：

import pyproj

# 定义 WGS84 和 2000国家大地坐标系的投影坐标系
wgs84 = pyproj.Proj(init='EPSG:4326')
gcj02 = pyproj.Proj(init='EPSG:4490')
x, y, z = pyproj.transform(wgs84, gcj02, 121.45, 31.21, 0)

2. 将WGS84坐标系三维坐标转换为2000国家大地坐标系下的三维坐标（X，Y，Z）。

可以使用 proj4 库实现这个步骤：

import pyproj

# 定义 WGS84 和 2000国家大地坐标系的投影坐标系
wgs84 = pyproj.Proj(init='EPSG:4326')
gcj02 = pyproj.Proj(init='EPSG:4490')
gcj02_to_2000 = pyproj.Proj(proj='geocent', ellps='CGCS2000', datum='CGCS2000')

# 将 WGS84 坐标系下的三维坐标转换为 2000国家大地坐标系下的三维坐标
x, y, z = pyproj.transform(gcj02, gcj02_to_2000, x, y, z)

3. 根据2000国家大地坐标系的定义，将三维坐标的X、Y、Z分别乘以比例因子，得到2000国家大地坐标系下的平面坐标（Easting，Northing）。

比例因子为：X、Y轴比例因子为 1，Z轴比例因子为 1.00000412。

# 计算 2000 国家大地坐标系下的平面坐标
Easting, Northing = x, y

最终得到的2000国家大地坐标系下的平面坐标为 (Easting, Northing) = (320090.796, 345723.248)。