Openlayers使用Proj4js坐标转换cgcs2000分带（zone）投影的问题

问题源于使用4528投影按照官网示例定义view时地图一片空白，于是开始了对Proj4js和投影转换的抽丝剥茧，这里记录下了解的知识和最后的结论。

官网示例和Proj4js的使用

OL示例见wms-image-custom-proj
这里简单介绍下proj4，我认知的proj4是基于一些约定参数组成的用于定义坐标参考的一组参数，和WKT相比更简约，省略了WKT里的完整参数定义。比如4528的wkt

PROJCRS["CGCS2000 / 3-degree Gauss-Kruger zone 40",BASEGEOGCRS["China Geodetic Coordinate System 2000",DATUM["China 2000",ELLIPSOID["CGCS2000",6378137,298.257222101,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],ID["EPSG",4490]],CONVERSION["3-degree Gauss-Kruger zone 40",METHOD["Transverse Mercator",ID["EPSG",9807]],PARAMETER["Latitude of natural origin",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8801]],PARAMETER["Longitude of natural origin",120,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["Scale factor at natural origin",1,SCALEUNIT["unity",1],ID["EPSG",8805]],PARAMETER["False easting",40500000,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["northing (X)",north,ORDER[1],LENGTHUNIT["metre",1]],AXIS["easting (Y)",east,ORDER[2],LENGTHUNIT["metre",1]],USAGE[SCOPE["Cadastre, engineering survey, topographic mapping (large scale)."],AREA["China - onshore between 118°30'E and 121°30'E."],BBOX[24.43,118.5,53.33,121.5]],ID["EPSG",4528]]

它的proj4

+proj=tmerc +lat_0=0 +lon_0=120 +k=1 +x_0=40500000 +y_0=0 +ellps=GRS80 +units=m +no_defs

详细介绍见参考
我对该示例做了简化，去掉了范围，更换图层为OSM，在vue里测试

import Map from "ol/Map";
import View from "ol/View";
import * as Layer from "ol/layer";
import * as Source from "ol/source";
import Projection from "ol/proj/Projection";
import { fromLonLat } from "ol/proj";
import proj4 from "proj4";
import { register } from "ol/proj/proj4";

  private mounted() {// 这里定义编码和proj4参数proj4.defs("EPSG:4549","+proj=tmerc +lat_0=0 +lon_0=120 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs");// 注册之后就可以使用epsg编码调用register(proj4);// 定义了底图 用的OSMconst base = new Layer.Tile({source: new Source.OSM(),});const projection = new Projection({code: "EPSG:4549",units: "m",});// 初始化地图对象const map = new Map({layers: [base],target: "map",view: new View({projection,center: fromLonLat([120.015, 36.408], "EPSG:4549"),zoom: 8,}),});}

获取Proj4的方法

比较方便的方法是在https://epsg.io/里查询，但我一直访问不了，作为替代可以在QGIS里查询
点击右下角的坐标参考按钮，很直观，亲测和网站里获取的是一样的

2000带投影的Proj4的问题

如果使用epsg:4528的proj4和编码，地图会一片空白还不报错。我尝试了其他坐标系，有的行有的不行，对比之后发现凡是以‘CGCS2000 / 3-degree Gauss-Kruger zone’开头的都不行，而同样是两千带投影的‘CGCS2000 / 3-degree Gauss-Kruger CM’系列就没事

这两个坐标系范围啥的都一样，我对比了下他俩的proj4

4528(CGCS2000 / 3-degree Gauss-Kruger zone 40)
+proj=tmerc +lat_0=0 +lon_0=120 +k=1 +x_0=40500000 +y_0=0 +ellps=GRS80 +units=m +no_defs
4549(CGCS2000 / 3-degree Gauss-Kruger CM 120E)
+proj=tmerc +lat_0=0 +lon_0=120 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs

除了+x_0的值，其他都一毛一样，所以问题八成出在这个参数上

如何治标不治本

研究了下，4528(CGCS2000 / 3-degree Gauss-Kruger zone 40)的+x_0的值40500000，代表40带，偏移500000米，如果去掉前边代表带号的40，就和4549的完全一样了。
估计是这种带号+偏移量的赋值方式不被支持，所以在使用4528的proj4时地图空白。
那么一种治标不治本的解决方法就来了，直接使用4549代替4528，或者使用4528，但是移除‘+x_0=40500000’里代表带号的‘40’。

如果有更好的建议欢迎提出，也欢迎指出错误，毕竟当初学的东西基本忘光了，很多都是基于猜测。