ST_Transform() 関数

以下の変換が有効です。

• 2 つの UNKNOWN (不明) 座標系間 (つまり、UNKNOWN を含む両 SRID の spatial_references 表内の srtext 列)
• 基準となる地理座標系が同一の、投影された座標系と非投影座標系との間
• 基準となる地理座標系が同一の、2 つの投影された座標系間
• 仮原点または座標系単位が異なる、地理座標系が同一の 2 つの座標系間

SE_CreateSrid() 関数で、オーストラリアに適した spatial_references 表エントリを作成していると仮定します。

EXECUTE FUNCTION SE_CreateSrid (110, -45, 156, -10,
                                "Australia: lat/lon coords");

(expression) 

        1002

ここで、オーストラリアに関するすべてのデータをロードします。

CREATE TABLE aus_locns (name varchar(128), locn ST_Point);

INSERT INTO aus_locns VALUES ("Adelaide",    '1002 point(139.14 -34.87)');
INSERT INTO aus_locns VALUES ("Brisbane",    '1002 point(153.36 -27.86)');
INSERT INTO aus_locns VALUES ("Canberra",    '1002 point(148.84 -35.56)');
INSERT INTO aus_locns VALUES ("Melbourne",   '1002 point(145.01 -37.94)');
INSERT INTO aus_locns VALUES ("Perth",       '1002 point(116.04 -32.12)');
INSERT INTO aus_locns VALUES ("Sydney",      '1002 point(151.37 -33.77)');

オーストラリア大陸のデータをすべてロードした後、ノーフォーク島やココス諸島などのいくつかの周辺諸島のデータを入れる必要があることがわかります。ただし、SRID 1002 用に選択した仮座標原点およびスケール ファクタは、この諸島には機能しません。

INSERT INTO aus_locns VALUES ("Norfolk Is.", '1002 point(167.83 -29.24)');
(USE19) - Coordinates out of bounds in ST_PointIn.

INSERT INTO aus_locns VALUES ("Cocos Is.",   '1002 point( 96.52 -12.08)');
(USE19) - Coordinates out of bounds in ST_PointIn.

解決策は、旧のデータと新規データの両方に対応する仮座標原点とスケール ファクタで、spatial_references 表エントリを作成して、旧のデータを更新することです。

EXECUTE FUNCTION SE_CreateSrid (95, -55, 170, -10,
                                "Australia + outer islands: lat/lon coords");

(expression) 

        1003


INSERT INTO aus_locns VALUES ("Norfolk Is.", '1003 point(167.83 -29.24)');
INSERT INTO aus_locns VALUES ("Cocos Is.",   '1003 point( 96.52 -12.08)');

UPDATE aus_locns
   SET locn = ST_Transform(locn, 1003)::ST_Point
   WHERE ST_Srid(locn) = 1002;

まず、非投影データに適した spatial_references 表エントリを作成します。この例では、1983 北米測地基準を使用します。この測地基準は、既知の標準測地基準のため、SE_CreateSrtext() 関数を使用して srtext フィールドを作成できます。

INSERT INTO spatial_references 
  (srid, description, falsex, falsey, xyunits,
   falsez, zunits, falsem, munits, srtext)
VALUES (1004, "Unprojected lat/lon, NAD 83 datum",
        -180, -90, 5000000, 0, 1000, 0, 1000, 
        SE_CreateSrtext(4269));

次に、表を作成してデータをロードします。

CREATE TABLE airports (id    char(4),
                       name  varchar(128),
                       locn  ST_Point);

INSERT INTO airports VALUES(
   'BTM', 'Bert Mooney',          '1004 point(-112.4975 45.9548)');
INSERT INTO airports VALUES(
   'BZN', 'Gallatin Field',       '1004 point(-111.1530 45.7769)');
INSERT INTO airports VALUES(
   'COD', 'Yellowstone Regional', '1004 point(-109.0238 44.5202)');
INSERT INTO airports VALUES(
   'JAC', 'Jackson Hole',         '1004 point(-110.7377 43.6073)');
INSERT INTO airports VALUES(
   'IDA', 'Fanning Field',        '1004 point(-112.0702 43.5146)');

必要な投影法に対して、1 つ以上の spatial_references 表エントリを作成します。基準となる地理座標系 (この場合は、NAD 83) が同一であることを確認します。

INSERT INTO spatial_references 
  (srid, description, falsex, falsey, xyunits,
   falsez, zunits, falsem, munits, srtext)
VALUES (1005, "UTM zone 12N, NAD 83 datum",
        336000, 4760000, 1000, 0, 1000, 0, 1000, 
        SE_CreateSrtext(26912));

必要に応じて、投影座標系にデータを変換します。

SELECT id, ST_Transform(locn, 1005) as utm FROM airports;


id   BTM
utm  1005 POINT (383951.152 5090115.666) 

id   BZN
utm  1005 POINT (488105.331 5069271.419) 

id   COD
utm  1005 POINT (657049.762 4931552.365) 

id   JAC
utm  1005 POINT (521167.881 4828291.447) 

id   IDA
utm  1005 POINT (413500.979 4818519.081) 

SRID が異なる 2 つの形状を比較する場合、ST_Transform() 関数を使用して形状を変換できます。

SELECT * FROM tab1 a, tab2 b WHERE
   ST_Intersects(a.shape, ST_Transform(b.shape, ST_SRID(a.shape)));