projection oracle,Oracle Spatial构建自定义投影坐标系

--PROJCS[--"Xian 1980 / 3-degree Gauss-Kruger CM 118.83E",--GEOGCS [--"Xian 1980",--DATUM [--"Xian 1980 (EPSG ID 6610)",--SPHEROID ["Xian 1980 (EPSG ID 7049)", 6378140, 298.257]],--PRIMEM [ "Greenwich", 0.000000 ],--UNIT ["Decimal Degree", 0.01745329251994328]],--PROJECTION ["Transverse Mercator"],--PARAMETER ["Latitude_Of_Origin", 0],--PARAMETER ["Central_Meridian", 118.8333],--PARAMETER ["Scale_Factor", 1],--PARAMETER ["False_Easting", 500000],--PARAMETER ["False_Northing", 0],--UNIT ["Meter", 1]]--首先，基础地理坐标系是Xian 1980， 它的EPSG ID为6610--然后找到与计算PROJECTION["Transverse Mercator"]ESPG等式

selectcoord_op_method_id,

legacy_namefromsdo_coord_op_methodswhere

not legacy_name is null

order bycoord_op_method_id;--结果显示--COORD_OP_METHOD_ID LEGACY_NAME---------------------------------------------9802 Lambert Conformal Conic--9803 Lambert Conformal Conic (Belgium 1972)--9805 Mercator--9806 Cassini--9807 Transverse Mercator--9829 Polar Stereographic---------------根据PROJECTION ["Transverse Mercator"]，选择9807--建立坐标转换规则，编号96170

insert intoMDSYS.SDO_COORD_OPS(

COORD_OP_ID,

COORD_OP_NAME,

COORD_OP_TYPE,

SOURCE_SRID,

TARGET_SRID,

COORD_TFM_VERSION,

COORD_OP_VARIANT,

COORD_OP_METHOD_ID,

UOM_ID_SOURCE_OFFSETS,

UOM_ID_TARGET_OFFSETS,

INFORMATION_SOURCE,

DATA_SOURCE,

SHOW_OPERATION,

IS_LEGACY,

LEGACY_CODE,

REVERSE_OP,

IS_IMPLEMENTED_FORWARD,

IS_IMPLEMENTED_REVERSE)VALUES(96170,'3-degree Gauss-Kruger CM 118.83E (EPSG OP 96170)','CONVERSION',NULL,NULL,NULL,NULL,9807,NULL,NULL,'EPSG','EPSG',1,'FLASE',NULL,1,1,1);--看看需要添加哪些参数

select

use.parameter_id ||':'||

use.legacy_param_namefromsdo_coord_op_param_useuse

where use.coord_op_method_id=9807

--结果：--8801:Latitude_Of_Origin--8802:Central_Meridian--8805:Scale_Factor--8806:False_Easting--8807:False_Northing--检查一下完全吻合啊!!!!!--现在查查经常使用的单位：

selectUOM_ID||':'||UNIT_OF_MEAS_NAMEfromsdo_units_of_measurewhereuom_idin (9001,9101,9102,9201)order byuom_id;--结果--9001:metre--9101:radian--9102:degree--9201:unity

--现在，插入参数--8801:Latitude_Of_Origin

insert intoMDSYS.Sdo_Coord_Op_Param_Vals(

COORD_OP_ID,

COORD_OP_METHOD_ID,

PARAMETER_ID,

PARAMETER_VALUE,

PARAM_VALUE_FILE_REF,

UOM_ID

)VALUES(96170,9807,8801,0,NULL,9102);--8802:Central_Meridian

insert intoMDSYS.SDO_COORD_OP_PARAM_VALS(

COORD_OP_ID,

COORD_OP_METHOD_ID,

PARAMETER_ID,

PARAMETER_VALUE,

PARAM_VALUE_FILE_REF,

UOM_ID)VALUES(96170,9807,8802,118.833333333333,NULL,9102);--8805:Scale_Factor

insert intoMDSYS.SDO_COORD_OP_PARAM_VALS(

COORD_OP_ID,

COORD_OP_METHOD_ID,

PARAMETER_ID,

PARAMETER_VALUE,

PARAM_VALUE_FILE_REF,

UOM_ID)VALUES(96170,9807,8805,1,NULL,9201);--8806:False_Easting

insert intoMDSYS.SDO_COORD_OP_PARAM_VALS(

COORD_OP_ID,

COORD_OP_METHOD_ID,

PARAMETER_ID,

PARAMETER_VALUE,

PARAM_VALUE_FILE_REF,

UOM_ID)VALUES(96170,9807,8806,500000,NULL,9001);--8807:False_Northing

insert intoMDSYS.SDO_COORD_OP_PARAM_VALS(

COORD_OP_ID,

COORD_OP_METHOD_ID,

PARAMETER_ID,

PARAMETER_VALUE,

PARAM_VALUE_FILE_REF,

UOM_ID)VALUES(96170,9807,8807,0,NULL,9001);--现在建立投影坐标系。首先要确定GEOG_CRS_DATUM_ID和COORD_SYS_ID。--GEOG_CRS_DATUM_ID是地理坐标系对应的datum_id--select datum_id from sdo_coord_ref_sys where srid=4610--GEOG_CRS_DATUM_ID为6610--CORD_SYS_ID是4530，CORD_SYS_NAME为Cartesian 2D CS. Axes: northing, easting (X,Y). Orientations: north, east. UoM: m.--现在开始建立投影坐标系

insert intoMDSYS.SDO_COORD_REF_SYSTEM(

SRID,

COORD_REF_SYS_NAME,

COORD_REF_SYS_KIND,

COORD_SYS_ID,

DATUM_ID,

GEOG_CRS_DATUM_ID,

SOURCE_GEOG_SRID,

PROJECTION_CONV_ID,

CMPD_HORIZ_SRID,

CMPD_VERT_SRID,

INFORMATION_SOURCE,

DATA_SOURCE,

IS_LEGACY,

LEGACY_CODE,

LEGACY_WKTEXT,

LEGACY_CS_BOUNDS,

IS_VALID,

SUPPORTS_SDO_GEOMETRY

)VALUES(9992010,'Xian 1980 / 3-degree Gauss-Kruger CM 118.83E','PROJECTED',4530,NULL,6610,4610,96170,NULL,NULL,'EPSG','EPSG','FALSE',NULL,NULL,NULL,'TRUE','TRUE');