WKT空間参照の順列定義は本来この順列に従う必要があるらしいとか。。。
ISO19162 の7.5
When is present in the WKT string the descriptions shall be ordered according to the axis order sequence.
このあたりかな??
GEODCRS["JGD2000", DATUM["Japanese Geodetic Datum 2000", ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1.0]]], CS[ellipsoidal,2], AXIS["latitude",north,ORDER[1]], AXIS["longitude",east,ORDER[2]], ANGLEUNIT["degree",0.01745329252], ID["EPSG",4612]]
PROJCRS["JGD2000 / Japan Plane Rectangular CS I", BASEGEODCRS["JGD2000", DATUM["Japanese Geodetic Datum 2000", ELLIPSOID["GRS 1980",6378137,298.257222101,LENGTHUNIT["metre",1.0]]]], CONVERSION["Japan Plane Rectangular CS zone I", METHOD["Transverse Mercator",ID["EPSG",9807]], PARAMETER["Latitude of natural origin",33,ANGLEUNIT["degree",0.01745329252]], PARAMETER["Longitude of natural origin",129.5,ANGLEUNIT["degree",0.01745329252]], PARAMETER["Scale factor at natural origin",0.9999,SCALEUNIT["unity",1.0]], PARAMETER["False easting",0,LENGTHUNIT["metre",1.0]], PARAMETER["False northing",0,LENGTHUNIT["metre",1.0]]], CS[cartesian,2], AXIS["northing (X)",north,ORDER[1]], AXIS["easting (Y)",east,ORDER[2]], LENGTHUNIT["metre",1.0], ID["EPSG",2443]]
平面直角座標系 1 系
X座標 1m
Y座標 2m
緯度 33°00′00.03246″
経度 129°30′00.07705″
33.000009016666666666666666666667
129.50002140277777777777777777778
SELECT ST_AsText( ST_Transform (ST_GeomFromText('POINT(129.5000214027777777777777777777 33.00000901666666666666666666666)', 4612), 2443) ) ;
‘POINT(1.99995816522274 0.999888726019038)’
‘2.1 USE_GEOS=1 USE_PROJ=1 USE_STATS=1’
SELECT ST_AsText( ST_GeographyFromText('SRID=4612;POINT(130 30)') );
‘POINT(30 50)’
現状は、こうだがそのうちXY逆転が必要なのかもな。。。めんど。
まあ空間参照テーブルのWKT定義側のORDER指定で望む側に直せば済む気がするがそういう問題でもないのか?。