1-- 2-- Regression tests that were filed as cases in bug tickets, 3-- referenced by bug number for historical interest. 4-- 5 6SET postgis.backend = 'sfcgal'; 7 8-- NOTE: some tests _require_ spatial_ref_sys entries. 9-- In particular, the GML output ones want auth_name and auth_srid too, 10-- so we provide one for EPSG:4326 11DELETE FROM spatial_ref_sys; 12INSERT INTO spatial_ref_sys ( srid, proj4text ) VALUES ( 32611, '+proj=utm +zone=11 +ellps=WGS84 +datum=WGS84 +units=m +no_defs' ); 13INSERT INTO spatial_ref_sys ( auth_name, auth_srid, srid, proj4text ) VALUES ( 'EPSG', 4326, 4326, '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs' ); 14INSERT INTO spatial_ref_sys ( srid, proj4text ) VALUES ( 32602, '+proj=utm +zone=2 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ' ); 15INSERT INTO spatial_ref_sys ( srid, proj4text ) VALUES ( 32702, '+proj=utm +zone=2 +south +ellps=WGS84 +datum=WGS84 +units=m +no_defs ' ); 16INSERT INTO spatial_ref_sys ( srid, proj4text ) VALUES ( 3395, '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ' ); 17 18INSERT INTO spatial_ref_sys (srid,proj4text) VALUES (32707,'+proj=utm +zone=7 +south +ellps=WGS84 +datum=WGS84 +units=m +no_defs '); 19 20-- #2 -- 21SELECT '#2', ST_AsText(ST_Union(g)) FROM 22( VALUES 23('SRID=4326;MULTIPOLYGON(((1 1,1 2,2 2,2 1,1 1)))'), 24('SRID=4326;MULTIPOLYGON(((2 1,3 1,3 2,2 2,2 1)))') 25) AS v(g); 26 27-- #11 -- 28SELECT '#11', round(ST_Distance (a.g, ST_Intersection(b.g, a.g))) AS distance 29FROM (SELECT '01020000000200000050E8303FC2E85141B017CFC05A825541000000E0C0E85141000000205C825541'::geometry AS g) a, 30 (SELECT 'LINESTRING(4694792.35840419 5638508.89950758,4694793.20840419 5638506.34950758)'::geometry AS g) b; 31 32-- #21 -- 33SELECT '#21', ST_AsEWKT(ST_Locate_Along_Measure(g, 4566)) FROM 34( VALUES 35(ST_GeomFromEWKT('SRID=31293;LINESTRINGM( 6193.76 5337404.95 4519, 6220.13 5337367.145 4566, 6240.889 5337337.383 4603 )')) 36) AS v(g); 37 38-- #22 -- 39SELECT ST_Within(g, 'POLYGON((0 0,10 0,20 10,10 20,0 20,-10 10,0 0))') FROM 40(VALUES 41('POLYGON((4 9,6 9,6 11,4 11,4 9))') 42) AS v(g); 43 44-- #33 -- 45CREATE TABLE road_pg (ID INTEGER, NAME VARCHAR(32)); 46SELECT '#33', AddGeometryColumn( '', 'public', 'road_pg','roads_geom', 330000, 'POINT', 2 ); 47DROP TABLE road_pg; 48 49-- #44 -- 50SELECT '#44', ST_Relate(g1, g2, 'T12101212'), ST_Relate(g1, g2, 't12101212') FROM 51(VALUES 52('POLYGON((0 0, 2 0, 2 2, 0 2, 0 0))', 'POLYGON((1 1, 3 1, 3 3, 1 3, 1 1))') 53) AS v(g1, g2); 54 55-- #58 -- 56SELECT '#58', round(ST_xmin(g)),round(ST_ymin(g)),round(ST_xmax(g)),round(ST_ymax(g)) FROM (SELECT ST_Envelope('CIRCULARSTRING(220268.439465645 150415.359530563,220227.333322076 150505.561285879,220227.353105332 150406.434743975)'::geometry) as g) AS foo; 57 58-- #65 -- 59SELECT '#65', ST_AsGML(ST_GeometryFromText('CURVEPOLYGON(CIRCULARSTRING(4 2,3 -1.0,1 -1,-1.0 4,4 2))')); 60 61-- #66 -- 62SELECT '#66', ST_AsText((ST_Dump(ST_GeomFromEWKT('CIRCULARSTRING(0 0,1 1,2 2)'))).geom); 63 64-- #68 -- 65SELECT '#68a', ST_AsText(ST_ShiftLongitude(ST_GeomFromText('MULTIPOINT(1 3, 4 5)'))); 66SELECT '#68b', ST_AsText(ST_ShiftLongitude(ST_GeomFromText('CIRCULARSTRING(1 3, 4 5, 6 7)'))); 67 68-- #69 -- 69SELECT '#69', ST_AsText(ST_Translate(ST_GeomFromText('CIRCULARSTRING(220268 150415,220227 150505,220227 150406)'),1,2)); 70 71-- #70 -- 72SELECT '#70', ST_NPoints(ST_LinetoCurve(ST_Buffer('POINT(1 2)',3))); 73 74-- #73 -- 75SELECT '#73', ST_AsText(ST_ForceCollection(ST_GeomFromEWKT('CIRCULARSTRING(1 1, 2 3, 4 5, 6 7, 5 6)'))); 76 77-- #80 -- 78SELECT '#80', ST_AsText(ST_Multi('MULTILINESTRING((0 0,1 1))')); 79 80-- #83 -- 81SELECT '#83', ST_AsText(ST_Multi(ST_GeomFromText('CIRCULARSTRING(220268 150415,220227 150505,220227 150406)'))); 82 83-- #85 -- 84SELECT '#85', ST_Distance(ST_GeomFromText('CIRCULARSTRING(220268 150415,220227 150505,220227 150406)'), ST_Point(220268, 150415)); 85 86-- #112 -- 87SELECT '#112', ST_AsText(ST_CurveToLine('GEOMETRYCOLLECTION(POINT(-10 50))'::geometry)); 88 89-- #113 -- 90SELECT '#113', ST_Locate_Along_Measure('LINESTRING(0 0 0, 1 1 1)', 0.5); 91 92-- #116 -- 93SELECT '#116', ST_AsText('010300000000000000'); 94 95-- #122 -- 96SELECT '#122', ST_AsText(ST_SnapToGrid(ST_GeomFromText('CIRCULARSTRING(220268 150415,220227 150505,220227 150406)'), 0.1)); 97 98-- #124 -- 99SELECT '#124a', ST_AsText(ST_GeomFromEWKT('COMPOUNDCURVE(CIRCULARSTRING(0 0,1 1,1 0),(1 0,30 5),CIRCULARSTRING(30 5,34 56,67 89))')); 100SELECT '#124b', ST_AsText(ST_GeomFromEWKT('COMPOUNDCURVE(CIRCULARSTRING(0 0,1 1,1 0),(1 0,30 6),CIRCULARSTRING(30 5,34 56,67 89))')); 101 102-- #145 -- 103SELECT '#145a', ST_AsText(ST_Buffer(ST_GeomFromText('LINESTRING(-116.93414544665981 34.16033385105459,-116.87777514700957 34.10831080544884,-116.86972224705954 34.086748622072776,-116.9327074288116 34.08458099517253,-117.00216369088065 34.130329331330216,-117.00216369088065 34.130329331330216)', 4326), 0)); 104SELECT '#145b', ST_Area(ST_Buffer(ST_GeomFromText('LINESTRING(-116.93414544665981 34.16033385105459,-116.87777514700957 34.10831080544884,-116.86972224705954 34.086748622072776,-116.9327074288116 34.08458099517253,-117.00216369088065 34.130329331330216,-117.00216369088065 34.130329331330216)', 4326), 0)); 105 106-- #146 -- 107SELECT '#146', ST_Distance(g1,g2), ST_Dwithin(g1,g2,0.01), ST_AsEWKT(g2) FROM (SELECT ST_geomFromEWKT('LINESTRING(1 2, 2 4)') As g1, ST_Collect(ST_GeomFromEWKT('LINESTRING(0 0, -1 -1)'), ST_GeomFromEWKT('MULTIPOINT(1 2,2 3)')) As g2) As foo; 108 109-- #156 -- 110SELECT '#156', ST_AsEWKT('0106000000010000000103000000010000000700000024213D12AA7BFD40945FF42576511941676A32F9017BFD40B1D67BEA7E511941C3E3C640DB7DFD4026CE38F4EE531941C91289C5A7EFD40017B8518E3531941646F1599AB7DFD409627F1F0AE521941355EBA49547CFD407B14AEC74652194123213D12AA7BFD40945FF42576511941'); 111 112-- #157 -- 113SELECT 114 '#157', 115 ST_GeometryType(g) As newname, 116 GeometryType(g) as oldname 117FROM ( VALUES 118 (ST_GeomFromText('POLYGON((-0.25 -1.25,-0.25 1.25,2.5 1.25,2.5 -1.25,-0.25 -1.25), (2.25 0,1.25 1,1.25 -1,2.25 0),(1 -1,1 1,0 0,1 -1))') ), 119 ( ST_Point(1,2) ), 120 ( ST_Buffer(ST_Point(1,2), 3) ), 121 ( ST_LineToCurve(ST_Buffer(ST_Point(1,2), 3)) ) , 122 ( ST_LineToCurve(ST_Boundary(ST_Buffer(ST_Point(1,2), 3))) ) 123 ) AS v(g); 124 125-- #168 -- 126SELECT '#168', ST_NPoints(g), ST_AsText(g), ST_isValidReason(g) 127FROM ( VALUES 128('01060000C00100000001030000C00100000003000000E3D9107E234F5041A3DB66BC97A30F4122ACEF440DAF9440FFFFFFFFFFFFEFFFE3D9107E234F5041A3DB66BC97A30F4122ACEF440DAF9440FFFFFFFFFFFFEFFFE3D9107E234F5041A3DB66BC97A30F4122ACEF440DAF9440FFFFFFFFFFFFEFFF'::geometry) 129) AS v(g); 130 131-- #175 -- 132SELECT '#175', ST_AsEWKT(ST_GeomFromEWKT('SRID=26915;POINT(482020 4984378.)')); 133 134-- #178 -- 135SELECT '#178a', ST_XMin(ST_MakeBox2D(ST_Point(5, 5), ST_Point(0, 0))); 136SELECT '#178b', ST_XMax(ST_MakeBox2D(ST_Point(5, 5), ST_Point(0, 0))); 137 138-- #179 -- 139SELECT '#179a', ST_MakeLine(ARRAY[NULL,NULL,NULL,NULL]); 140SELECT '#179b', ST_MakeLine(ARRAY[NULL,NULL,NULL,NULL]); 141 142-- #183 -- 143SELECT '#183', ST_AsText(ST_SnapToGrid(ST_LineToCurve(ST_LineMerge(ST_Collect(ST_CurveToLine(ST_GeomFromEWKT('CIRCULARSTRING(0 0, 1 1, 1 0)')),ST_GeomFromEWKT('LINESTRING(1 0, 0 1)') ))), 1E-10)); 144 145-- #210 -- 146SELECT '#210a', ST_Union(ARRAY[NULL,NULL,NULL,NULL]) ; 147SELECT '#210b', ST_MakeLine(ARRAY[NULL,NULL,NULL,NULL]) ; 148 149-- #213 -- 150SELECT '#213', round(ST_Perimeter(ST_CurveToLine(ST_GeomFromEWKT('CURVEPOLYGON(COMPOUNDCURVE(CIRCULARSTRING(0 0,2 0, 2 1, 2 3, 4 3),(4 3, 4 5, 1 4, 0 0)))')))); 151 152-- #234 -- 153SELECT '#234', ST_AsText(ST_GeomFromText('COMPOUNDCURVE( (0 0,1 1) )')); 154 155-- #239 -- 156--SELECT '#239', ST_AsSVG('010700002031BF0D0000000000'); 157 158-- #241 -- 159CREATE TABLE c ( the_geom GEOMETRY); 160INSERT INTO c SELECT ST_MakeLine(ST_Point(-10,40),ST_Point(40,-10)) As the_geom; 161INSERT INTO c SELECT ST_MakeLine(ST_Point(-10,40),ST_Point(40,-10)) As the_geom; 162SELECT '#241', sum(ST_LineCrossingDirection(the_geom, ST_GeomFromText('LINESTRING(1 2,3 4)'))) FROM c; 163DROP TABLE c; 164 165-- #254 -- 166SELECT '#254', ST_AsText(ST_Segmentize(ST_GeomFromText('GEOMETRYCOLLECTION EMPTY'), 0.5)); 167 168-- #259 -- 169SELECT '#259', ST_Distance(ST_GeographyFromText('SRID=4326;POLYGON EMPTY'), ST_GeographyFromText('SRID=4326;POINT(1 2)')); 170 171-- #260 -- 172SELECT '#260', round(ST_Distance(ST_GeographyFromText('SRID=4326;POINT(-10 40)'), ST_GeographyFromText('SRID=4326;POINT(-10 55)'))); 173 174-- #261 -- 175SELECT '#261', ST_Distance(ST_GeographyFromText('SRID=4326;POINT(-71.0325022849392 42.3793285830812)'),ST_GeographyFromText('SRID=4326;POLYGON((-71.0325022849392 42.3793285830812,-71.0325745928559 42.3793012556699,-71.0326708728343 42.3794450989722,-71.0326045866257 42.3794706688942,-71.0325022849392 42.3793285830812))')); 176 177-- #262 -- 178SELECT '#262', ST_AsText(pt.the_geog) As wkt_pt, ST_Covers(poly.the_geog, pt.the_geog) As geog, 179 ST_Covers( 180 ST_Transform(CAST(poly.the_geog As geometry),32611), 181 ST_Transform(CAST(pt.the_geog As geometry),32611)) As utm, 182 ST_Covers( 183 CAST(poly.the_geog As geometry), 184 CAST(pt.the_geog As geometry) 185 ) As pca 186FROM (SELECT ST_GeographyFromText('SRID=4326;POLYGON((-119.5434 34.9438,-119.5437 34.9445,-119.5452 34.9442,-119.5434 34.9438))') As the_geog) 187 As poly 188 CROSS JOIN 189 (VALUES 190 ( ST_GeographyFromText('SRID=4326;POINT(-119.5434 34.9438)') ) , 191 ( ST_GeographyFromText('SRID=4326;POINT(-119.5452 34.9442)') ) , 192 ( ST_GeographyFromText('SRID=4326;POINT(-119.5434 34.9438)') ), 193 ( ST_GeographyFromText('SRID=4326;POINT(-119.5438 34.9443)') ) 194 ) 195 As pt(the_geog); 196 197-- #263 -- 198SELECT '#263', ST_AsEWKT(geometry(geography(pt.the_geom))) As wkt, 199 ST_Covers( 200 ST_Transform(poly.the_geom,32611), 201 ST_Transform(pt.the_geom,32611)) As utm, 202 ST_Covers( 203 poly.the_geom, 204 pt.the_geom) 205 As pca, 206 ST_Covers(geometry(geography(poly.the_geom)), 207 geometry(geography(pt.the_geom))) As gm_to_gg_gm_pca 208 209FROM (SELECT ST_GeomFromEWKT('SRID=4326;POLYGON((-119.5434 34.9438,-119.5437 34.9445,-119.5452 34.9442,-119.5434 34.9438))') As the_geom) 210 As poly 211 CROSS JOIN 212 (VALUES 213 ( ST_GeomFromEWKT('SRID=4326;POINT(-119.5434 34.9438)') ) , 214 ( ST_GeomFromEWKT('SRID=4326;POINT(-119.5452 34.9442)') ) , 215 ( ST_GeomFromEWKT('SRID=4326;POINT(-119.5434 34.9438)') ), 216 ( ST_GeomFromEWKT('SRID=4326;POINT(-119.5438 34.9443)') ) 217 ) 218 As pt(the_geom); 219 220-- #271 -- 221SELECT '#271', ST_Covers( 222'POLYGON((-9.123456789 50,51 -11.123456789,-10.123456789 50,-9.123456789 50))'::geography, 223'POINT(-10.123456789 50)'::geography 224); 225 226-- #272 -- 227SELECT '#272', ST_LineCrossingDirection(foo.line1, foo.line2) As l1_cross_l2 , 228 ST_LineCrossingDirection(foo.line2, foo.line1) As l2_cross_l1 229FROM (SELECT 230 ST_GeomFromText('LINESTRING(25 169,89 114,40 70,86 43)') As line1, ST_GeomFromText('LINESTRING(2.99 90.16,71 74,20 140,171 154)') As line2 ) As foo; 231 232-- #277 -- 233SELECT '#277', ST_AsGML(2, ST_GeomFromText('POINT(1 1e308)')); 234 235-- #299 -- 236SELECT '#299', round(ST_Y(geometry(ST_Intersection(ST_GeographyFromText('POINT(1.2456 2)'), ST_GeographyFromText('POINT(1.2456 2)'))))); 237 238-- #304 -- 239 240SELECT '#304'; 241 242CREATE OR REPLACE FUNCTION utmzone(geometry) 243 RETURNS integer AS 244$BODY$ 245DECLARE 246 geomgeog geometry; 247 zone int; 248 pref int; 249 250BEGIN 251 geomgeog:= ST_Transform($1,4326); 252 253 IF (ST_Y(geomgeog))>0 THEN 254 pref:=32600; 255 ELSE 256 pref:=32700; 257 END IF; 258 259 zone:=floor((ST_X(geomgeog)+180)/6)+1; 260 IF ( zone > 60 ) THEN zone := 60; END IF; 261 262 RETURN zone+pref; 263END; 264$BODY$ LANGUAGE 'plpgsql' IMMUTABLE 265 COST 100; 266 267CREATE TABLE utm_dots ( the_geog geography, utm_srid integer); 268INSERT INTO utm_dots SELECT geography(ST_SetSRID(ST_Point(i*10,j*10),4326)) As the_geog, utmzone(ST_SetSRID(ST_Point(i*10,j*10),4326)) As utm_srid FROM generate_series(-17,17) As i CROSS JOIN generate_series(-8,8) As j; 269 270SELECT ST_AsText(the_geog) as the_pt, 271 utm_srid, 272 ST_Area(ST_Buffer(the_geog,10)) As the_area, 273 ST_Area(geography(ST_Transform(ST_Buffer(ST_Transform(geometry(the_geog),utm_srid),10),4326))) As geog_utm_area 274FROM utm_dots 275WHERE ST_Area(ST_Buffer(the_geog,10)) NOT between 307 and 315 276LIMIT 10; 277 278SELECT '#304.a', Count(*) FROM utm_dots WHERE ST_DWithin(the_geog, 'POINT(0 0)'::geography, 3000000); 279 280CREATE INDEX utm_dots_gix ON utm_dots USING GIST (the_geog); 281SELECT '#304.b', Count(*) FROM utm_dots WHERE ST_DWithin(the_geog, 'POINT(0 0)'::geography, 300000); 282 283DROP FUNCTION utmzone(geometry); 284DROP TABLE utm_dots; 285 286-- #408 -- 287SELECT '#408', substring(st_isvalidreason('0105000020E0670000010000000102000020E06700000100000016DA52BA62A04141FFF3AD290B735241') from E'.*\\['); 288SELECT '#408.1', st_isvalid('01050000000100000001020000000100000000000000000000000000000000000000'); 289SELECT '#408.2', st_isvalidreason('01020000000100000000000000000000000000000000000000'); 290SELECT '#408.3', st_isvalid('0106000020BB0B000001000000010300000005000000D6000000000000C0F1A138410AD7A3103190524114AE4721F7A138410000000030905241713D0A57FAA1384185EB51982C9052417B14AE87FAA13841000000402A905241AE47E13AFBA1384114AE474128905241EC51B81EFDA138413D0AD7632690524152B81E85FFA13841D7A3707D2590524152B81E4500A23841713D0A072590524100000000FFA13841713D0AF724905241CDCCCCCCFFA13841CDCCCC0C249052419A99991901A238419A999919249052411F85EBD101A23841D7A3704D23905241E17A14AE02A23841A4703D3A2290524185EB51F808A23841000000D021905241E17A14EE0BA238413333335321905241666666A615A23841B81E85AB1F905241713D0AD710A23841B81E858B1E905241666666260CA2384114AE4731209052410AD7A37009A23841CDCCCCEC20905241AE47E13A06A23841EC51B87E219052419A9999D903A23841EC51B80E21905241EC51B8DE08A238418FC2F5681F905241666666660EA23841A4703D9A1D9052413D0AD7A30FA2384114AE47A11C9052413333333311A23841000000101C90524148E17A1414A23841333333931A905241EC51B81E1AA23841EC51B89E1890524185EB51381CA23841D7A370BD17905241295C8FC220A2384185EB51C816905241C3F528DC2BA238418FC2F5E814905241E17A14EE2EA23841EC51B8FE13905241AE47E1FA2BA23841EC51B8EE119052415C8FC2F52BA2384185EB51C810905241333333332EA2384114AE4701109052417B14AEC73BA238410AD7A3300D905241CDCCCCCC43A238417B14AE870B905241CDCCCC4C45A23841A4703D3A0B9052413D0AD7A34DA238411F85EB810A905241295C8F8250A23841CDCCCCBC09905241000000C053A238410AD7A36009905241A4703D0A5CA2384114AE4751099052418FC2F5A866A23841A4703D4A0A9052410000008069A2384114AE47B10A9052415C8FC2B573A23841B81E851B0D905241AE47E1FA76A23841666666C60D9052417B14AE077BA238411F85EB510E905241AE47E1FA7FA23841666666E60E9052415C8FC2358DA23841000000F010905241EC51B85E95A2384148E17A14119052417B14AE4795A23841CDCCCCCC1090524114AE47E18CA2384185EB51F80F90524148E17AD480A2384114AE47A10E9052413D0AD76383A23841713D0AF70C9052417B14AEC784A23841A4703D8A0C9052410000000086A23841A4703D0A0B9052411F85EB9187A23841B81E851B09905241F6285CCF86A2384100000010099052417B14AE4787A2384185EB51180890524114AE472186A23841CDCCCCFC079052410AD7A3B086A23841C3F5280C07905241333333B385A23841CDCCCC0C07905241000000407EA2384152B81E15079052419A9999197BA23841A4703D2A079052416666666673A2384114AE47710790524114AE47A16EA238413D0AD7F30690524114AE472170A238418FC2F5D8059052413D0AD72372A2384185EB511806905241A4703DCA71A23841D7A3707D05905241EC51B8DE7EA238411F85EBD104905241E17A14AE76A23841333333E30390524185EB51B875A23841713D0A47029052415C8FC2F576A23841D7A3703D029052419A99999976A2384114AE475101905241333333B378A23841EC51B83E019052410AD7A3307CA238410AD7A320019052418FC2F52885A23841295C8FD20090524148E17A9483A238416666664600905241295C8FC276A238411F85EB9100905241CDCCCC4C76A2384100000070FF8F5241AE47E13A76A2384133333343FF8F52417B14AE0783A23841D7A370FDFE8F5241E17A14AE81A23841C3F5283CFE8F52415C8FC2357DA2384114AE47A1FE8F5241CDCCCC8C7BA238415C8FC275FD8F524148E17AD47FA238417B14AE17FD8F5241000000807EA238413D0AD7A3FC8F5241713D0AD776A2384133333363FD8F5241CDCCCC8C73A2384152B81EF5FD8F5241E17A14AE70A2384114AE4711FE8F524185EB51B870A23841B81E852BFF8F524152B81E056CA2384114AE4731FF8F5241B81E856B6BA2384152B81E25FE8F5241AE47E13A69A23841E17A142EFE8F52419A99995962A238415C8FC275FE8F52415C8FC27562A238411F85EB710090524114AE47A15BA23841EC51B86E00905241EC51B85E5BA238419A9999E902905241295C8F025BA238417B14AE8704905241A4703D8A54A23841713D0A7704905241C3F5285C51A2384148E17A44049052418FC2F5E84CA23841A4703D4A049052410AD7A3304CA2384148E17A94049052419A9999994CA238413D0AD7F305905241AE47E1FA49A2384114AE471106905241EC51B85E49A23841AE47E1FA07905241AE47E13A43A238415C8FC255099052416666662644A23841D7A3707D0A9052415C8FC2F53FA23841CDCCCC2C0B90524185EB517839A2384185EB51380C90524152B81E4532A23841D7A370DD099052419A9999592FA2384152B81E4509905241D7A3707D2AA23841B81E857B089052415C8FC27526A238416666667608905241EC51B81E26A238417B14AE7707905241295C8F0226A23841F6285C1F07905241666666A61CA23841CDCCCC6C07905241713D0A1716A23841B81E856B07905241B81E85EB15A2384100000040079052411F85EB1108A2384114AE4741079052415C8FC23505A2384148E17AC40790524152B81E0506A23841EC51B84E09905241713D0A5707A23841AE47E10A0A905241EC51B81E12A2384152B81EA50890524152B81E4521A23841F6285CDF0890524148E17A1427A23841EC51B80E099052415C8FC2B52AA2384185EB5168099052413D0AD7E32FA23841CDCCCC3C0A9052413D0AD7A331A23841333333930A90524185EB51F833A238410AD7A3B00B9052411F85EB5135A23841EC51B8AE0B90524185EB513836A23841333333530D9052410000004030A238415C8FC2B50E905241000000802FA238415C8FC2750E9052410000008028A238419A9999B90E905241F6285C0F28A23841333333630E90524114AE47E123A23841333333730E905241D7A3707D23A2384152B81E750E9052410000008023A2384185EB51E80E905241D7A370BD1EA23841713D0A170F9052419A9999991EA23841333333830F9052419A9999991CA238411F85EB810F90524152B81E450EA2384152B81E750F9052417B14AE470AA23841EC51B89E0F905241B81E85EB05A238410AD7A3F00F905241D7A3707DFFA138417B14AEF71090524185EB51F8FDA138413D0AD72311905241295C8F02FAA13841EC51B83E11905241713D0A17EFA13841E17A145E129052417B14AE87CBA1384152B81E05179052418FC2F528C0A13841D7A3709D17905241D7A370BDC0A13841F6285C5F18905241295C8F82CBA138419A9999F917905241AE47E1BAD7A138415C8FC255169052411F85EB91E0A13841000000401590524114AE47A1E2A1384185EB5108159052415C8FC235ECA13841C3F5286C13905241C3F528DCF4A13841000000701290524152B81E45F5A138415C8FC2A51290524100000040F6A13841CDCCCC7C129052417B14AE47FBA138410000000012905241EC51B85E00A2384185EB51A8119052413D0AD7A305A238419A99990911905241EC51B81E0AA238411F85EB611090524148E17A940EA23841A4703D0A1090524152B81E0514A2384152B81EE50F905241F6285CCF15A238410AD7A3E00F9052417B14AEC718A23841CDCCCCDC0F9052419A9999191AA238415C8FC2E50F905241EC51B81E20A238415C8FC2E50F90524148E17A1420A23841333333C30F9052417B14AE8729A23841713D0AA70F905241000000002AA238413D0AD733119052411F85EB9129A23841C3F528EC12905241E17A142E28A238419A999919149052418FC2F5A825A23841CDCCCC1C15905241C3F5285C21A2384114AE473116905241F6285C4F1BA2384152B81E4517905241D7A3707D1AA2384148E17A1417905241C3F5285C19A23841D7A3702D179052413333333318A23841B81E85CB17905241000000C017A2384185EB5108189052413D0AD72310A238419A9999F91790524148E17A5410A2384148E17AD41690524148E17AD40DA2384152B81ED516905241B81E85AB09A23841F6285C1F179052418FC2F5A809A2384185EB51C8179052418FC2F56808A238417B14AEC7179052410000008008A2384185EB515818905241AE47E1BA16A23841F6285C6F189052410AD7A33017A23841C3F5289C18905241000000C012A23841AE47E10A1A9052411F85EB5110A238413D0AD7131B905241000000800DA238415C8FC2451C90524148E17A140DA23841A4703D4A1C9052413D0AD7230BA2384148E17A141E905241295C8F4205A23841295C8F421E905241F6285C4F00A238419A9999691E905241C3F5289CFDA13841AE47E16A1E905241D7A370FDFDA1384114AE47711F90524185EB51F801A23841B81E855B1F905241A4703D0A05A23841CDCCCCAC1F905241E17A14AE02A23841EC51B89E20905241295C8F82FFA138419A9999E9219052419A999919FEA1384185EB5188239052410AD7A330FAA13841E17A14DE25905241C3F5281CF9A13841F6285CFF26905241AE47E13AF8A138415C8FC2E527905241CDCCCC8CF7A1384185EB51B82890524185EB51F8F7A138413333334329905241C3F5281CF6A138418FC2F5182C905241A4703D4AF5A13841EC51B8BE2C905241CDCCCC0CF4A13841713D0AF72E90524185EB5138F3A13841000000502F905241000000C0F1A138410AD7A3103190524111000000CDCCCC4C6AA23841713D0A170A905241CDCCCC4C6FA238411F85EB710790524148E17A947AA2384114AE47410890524185EB51387BA23841D7A370ED07905241666666267BA23841713D0A7707905241EC51B85E7EA2384114AE476107905241E17A14AE7EA2384114AE473108905241713D0A5785A23841CDCCCCBC089052410AD7A33083A23841F6285C9F0B9052413333337382A23841666666460C90524148E17AD481A23841333333D30C90524114AE472181A23841333333730D905241AE47E1BA7FA23841CDCCCC3C0E9052418FC2F56876A238415C8FC2050D90524185EB513872A238418FC2F5D80B90524185EB513870A23841EC51B82E0B905241CDCCCC4C6AA23841713D0A170A90524113000000A4703DCA6BA238415C8FC295FF8F5241B81E856B74A238410AD7A380FF8F5241EC51B89E74A23841E17A149E00905241AE47E1BA74A23841F6285C4F019052419A99995974A23841295C8F520290524148E17AD473A2384114AE47B10390524152B81EC572A2384114AE4701059052417B14AE476DA23841AE47E1DA04905241F6285C8F6DA23841333333A303905241000000006CA2384152B81E9503905241EC51B81E69A23841A4703D7A039052410AD7A3B068A23841B81E85AB049052415C8FC2F55CA23841CDCCCC8C0490524152B81E455DA238411F85EB0103905241333333B35DA23841AE47E1DA00905241EC51B85E67A23841C3F528FC009052410AD7A3B067A23841CDCCCCCCFF8F524185EB51B86BA23841333333C3FF8F5241A4703DCA6BA238415C8FC295FF8F5241020000003D0AD7A346A238413D0AD7E3099052413D0AD7A346A238413D0AD7E309905241100000007B14AEC74DA238410AD7A3B008905241EC51B85E4EA23841B81E850B079052413D0AD7234FA23841C3F528DC0490524185EB51385CA23841EC51B8EE04905241B81E852B60A23841A4703DDA049052410000008066A23841A4703DFA049052411F85EB116AA23841CDCCCC0C05905241CDCCCC8C6EA23841A4703D5A05905241A4703D0A6CA238417B14AE37079052411F85EBD16AA23841A4703DCA07905241AE47E1FA68A23841E17A14FE08905241295C8F0268A238410AD7A320099052419A99999963A23841AE47E10A099052415C8FC2B55AA23841333333C30890524114AE476156A23841CDCCCCBC089052417B14AEC74DA238410AD7A3B008905241'); 291SELECT '#408.4', st_isvalidreason('0106000020BB0B000001000000010300000005000000D6000000000000C0F1A138410AD7A3103190524114AE4721F7A138410000000030905241713D0A57FAA1384185EB51982C9052417B14AE87FAA13841000000402A905241AE47E13AFBA1384114AE474128905241EC51B81EFDA138413D0AD7632690524152B81E85FFA13841D7A3707D2590524152B81E4500A23841713D0A072590524100000000FFA13841713D0AF724905241CDCCCCCCFFA13841CDCCCC0C249052419A99991901A238419A999919249052411F85EBD101A23841D7A3704D23905241E17A14AE02A23841A4703D3A2290524185EB51F808A23841000000D021905241E17A14EE0BA238413333335321905241666666A615A23841B81E85AB1F905241713D0AD710A23841B81E858B1E905241666666260CA2384114AE4731209052410AD7A37009A23841CDCCCCEC20905241AE47E13A06A23841EC51B87E219052419A9999D903A23841EC51B80E21905241EC51B8DE08A238418FC2F5681F905241666666660EA23841A4703D9A1D9052413D0AD7A30FA2384114AE47A11C9052413333333311A23841000000101C90524148E17A1414A23841333333931A905241EC51B81E1AA23841EC51B89E1890524185EB51381CA23841D7A370BD17905241295C8FC220A2384185EB51C816905241C3F528DC2BA238418FC2F5E814905241E17A14EE2EA23841EC51B8FE13905241AE47E1FA2BA23841EC51B8EE119052415C8FC2F52BA2384185EB51C810905241333333332EA2384114AE4701109052417B14AEC73BA238410AD7A3300D905241CDCCCCCC43A238417B14AE870B905241CDCCCC4C45A23841A4703D3A0B9052413D0AD7A34DA238411F85EB810A905241295C8F8250A23841CDCCCCBC09905241000000C053A238410AD7A36009905241A4703D0A5CA2384114AE4751099052418FC2F5A866A23841A4703D4A0A9052410000008069A2384114AE47B10A9052415C8FC2B573A23841B81E851B0D905241AE47E1FA76A23841666666C60D9052417B14AE077BA238411F85EB510E905241AE47E1FA7FA23841666666E60E9052415C8FC2358DA23841000000F010905241EC51B85E95A2384148E17A14119052417B14AE4795A23841CDCCCCCC1090524114AE47E18CA2384185EB51F80F90524148E17AD480A2384114AE47A10E9052413D0AD76383A23841713D0AF70C9052417B14AEC784A23841A4703D8A0C9052410000000086A23841A4703D0A0B9052411F85EB9187A23841B81E851B09905241F6285CCF86A2384100000010099052417B14AE4787A2384185EB51180890524114AE472186A23841CDCCCCFC079052410AD7A3B086A23841C3F5280C07905241333333B385A23841CDCCCC0C07905241000000407EA2384152B81E15079052419A9999197BA23841A4703D2A079052416666666673A2384114AE47710790524114AE47A16EA238413D0AD7F30690524114AE472170A238418FC2F5D8059052413D0AD72372A2384185EB511806905241A4703DCA71A23841D7A3707D05905241EC51B8DE7EA238411F85EBD104905241E17A14AE76A23841333333E30390524185EB51B875A23841713D0A47029052415C8FC2F576A23841D7A3703D029052419A99999976A2384114AE475101905241333333B378A23841EC51B83E019052410AD7A3307CA238410AD7A320019052418FC2F52885A23841295C8FD20090524148E17A9483A238416666664600905241295C8FC276A238411F85EB9100905241CDCCCC4C76A2384100000070FF8F5241AE47E13A76A2384133333343FF8F52417B14AE0783A23841D7A370FDFE8F5241E17A14AE81A23841C3F5283CFE8F52415C8FC2357DA2384114AE47A1FE8F5241CDCCCC8C7BA238415C8FC275FD8F524148E17AD47FA238417B14AE17FD8F5241000000807EA238413D0AD7A3FC8F5241713D0AD776A2384133333363FD8F5241CDCCCC8C73A2384152B81EF5FD8F5241E17A14AE70A2384114AE4711FE8F524185EB51B870A23841B81E852BFF8F524152B81E056CA2384114AE4731FF8F5241B81E856B6BA2384152B81E25FE8F5241AE47E13A69A23841E17A142EFE8F52419A99995962A238415C8FC275FE8F52415C8FC27562A238411F85EB710090524114AE47A15BA23841EC51B86E00905241EC51B85E5BA238419A9999E902905241295C8F025BA238417B14AE8704905241A4703D8A54A23841713D0A7704905241C3F5285C51A2384148E17A44049052418FC2F5E84CA23841A4703D4A049052410AD7A3304CA2384148E17A94049052419A9999994CA238413D0AD7F305905241AE47E1FA49A2384114AE471106905241EC51B85E49A23841AE47E1FA07905241AE47E13A43A238415C8FC255099052416666662644A23841D7A3707D0A9052415C8FC2F53FA23841CDCCCC2C0B90524185EB517839A2384185EB51380C90524152B81E4532A23841D7A370DD099052419A9999592FA2384152B81E4509905241D7A3707D2AA23841B81E857B089052415C8FC27526A238416666667608905241EC51B81E26A238417B14AE7707905241295C8F0226A23841F6285C1F07905241666666A61CA23841CDCCCC6C07905241713D0A1716A23841B81E856B07905241B81E85EB15A2384100000040079052411F85EB1108A2384114AE4741079052415C8FC23505A2384148E17AC40790524152B81E0506A23841EC51B84E09905241713D0A5707A23841AE47E10A0A905241EC51B81E12A2384152B81EA50890524152B81E4521A23841F6285CDF0890524148E17A1427A23841EC51B80E099052415C8FC2B52AA2384185EB5168099052413D0AD7E32FA23841CDCCCC3C0A9052413D0AD7A331A23841333333930A90524185EB51F833A238410AD7A3B00B9052411F85EB5135A23841EC51B8AE0B90524185EB513836A23841333333530D9052410000004030A238415C8FC2B50E905241000000802FA238415C8FC2750E9052410000008028A238419A9999B90E905241F6285C0F28A23841333333630E90524114AE47E123A23841333333730E905241D7A3707D23A2384152B81E750E9052410000008023A2384185EB51E80E905241D7A370BD1EA23841713D0A170F9052419A9999991EA23841333333830F9052419A9999991CA238411F85EB810F90524152B81E450EA2384152B81E750F9052417B14AE470AA23841EC51B89E0F905241B81E85EB05A238410AD7A3F00F905241D7A3707DFFA138417B14AEF71090524185EB51F8FDA138413D0AD72311905241295C8F02FAA13841EC51B83E11905241713D0A17EFA13841E17A145E129052417B14AE87CBA1384152B81E05179052418FC2F528C0A13841D7A3709D17905241D7A370BDC0A13841F6285C5F18905241295C8F82CBA138419A9999F917905241AE47E1BAD7A138415C8FC255169052411F85EB91E0A13841000000401590524114AE47A1E2A1384185EB5108159052415C8FC235ECA13841C3F5286C13905241C3F528DCF4A13841000000701290524152B81E45F5A138415C8FC2A51290524100000040F6A13841CDCCCC7C129052417B14AE47FBA138410000000012905241EC51B85E00A2384185EB51A8119052413D0AD7A305A238419A99990911905241EC51B81E0AA238411F85EB611090524148E17A940EA23841A4703D0A1090524152B81E0514A2384152B81EE50F905241F6285CCF15A238410AD7A3E00F9052417B14AEC718A23841CDCCCCDC0F9052419A9999191AA238415C8FC2E50F905241EC51B81E20A238415C8FC2E50F90524148E17A1420A23841333333C30F9052417B14AE8729A23841713D0AA70F905241000000002AA238413D0AD733119052411F85EB9129A23841C3F528EC12905241E17A142E28A238419A999919149052418FC2F5A825A23841CDCCCC1C15905241C3F5285C21A2384114AE473116905241F6285C4F1BA2384152B81E4517905241D7A3707D1AA2384148E17A1417905241C3F5285C19A23841D7A3702D179052413333333318A23841B81E85CB17905241000000C017A2384185EB5108189052413D0AD72310A238419A9999F91790524148E17A5410A2384148E17AD41690524148E17AD40DA2384152B81ED516905241B81E85AB09A23841F6285C1F179052418FC2F5A809A2384185EB51C8179052418FC2F56808A238417B14AEC7179052410000008008A2384185EB515818905241AE47E1BA16A23841F6285C6F189052410AD7A33017A23841C3F5289C18905241000000C012A23841AE47E10A1A9052411F85EB5110A238413D0AD7131B905241000000800DA238415C8FC2451C90524148E17A140DA23841A4703D4A1C9052413D0AD7230BA2384148E17A141E905241295C8F4205A23841295C8F421E905241F6285C4F00A238419A9999691E905241C3F5289CFDA13841AE47E16A1E905241D7A370FDFDA1384114AE47711F90524185EB51F801A23841B81E855B1F905241A4703D0A05A23841CDCCCCAC1F905241E17A14AE02A23841EC51B89E20905241295C8F82FFA138419A9999E9219052419A999919FEA1384185EB5188239052410AD7A330FAA13841E17A14DE25905241C3F5281CF9A13841F6285CFF26905241AE47E13AF8A138415C8FC2E527905241CDCCCC8CF7A1384185EB51B82890524185EB51F8F7A138413333334329905241C3F5281CF6A138418FC2F5182C905241A4703D4AF5A13841EC51B8BE2C905241CDCCCC0CF4A13841713D0AF72E90524185EB5138F3A13841000000502F905241000000C0F1A138410AD7A3103190524111000000CDCCCC4C6AA23841713D0A170A905241CDCCCC4C6FA238411F85EB710790524148E17A947AA2384114AE47410890524185EB51387BA23841D7A370ED07905241666666267BA23841713D0A7707905241EC51B85E7EA2384114AE476107905241E17A14AE7EA2384114AE473108905241713D0A5785A23841CDCCCCBC089052410AD7A33083A23841F6285C9F0B9052413333337382A23841666666460C90524148E17AD481A23841333333D30C90524114AE472181A23841333333730D905241AE47E1BA7FA23841CDCCCC3C0E9052418FC2F56876A238415C8FC2050D90524185EB513872A238418FC2F5D80B90524185EB513870A23841EC51B82E0B905241CDCCCC4C6AA23841713D0A170A90524113000000A4703DCA6BA238415C8FC295FF8F5241B81E856B74A238410AD7A380FF8F5241EC51B89E74A23841E17A149E00905241AE47E1BA74A23841F6285C4F019052419A99995974A23841295C8F520290524148E17AD473A2384114AE47B10390524152B81EC572A2384114AE4701059052417B14AE476DA23841AE47E1DA04905241F6285C8F6DA23841333333A303905241000000006CA2384152B81E9503905241EC51B81E69A23841A4703D7A039052410AD7A3B068A23841B81E85AB049052415C8FC2F55CA23841CDCCCC8C0490524152B81E455DA238411F85EB0103905241333333B35DA23841AE47E1DA00905241EC51B85E67A23841C3F528FC009052410AD7A3B067A23841CDCCCCCCFF8F524185EB51B86BA23841333333C3FF8F5241A4703DCA6BA238415C8FC295FF8F5241020000003D0AD7A346A238413D0AD7E3099052413D0AD7A346A238413D0AD7E309905241100000007B14AEC74DA238410AD7A3B008905241EC51B85E4EA23841B81E850B079052413D0AD7234FA23841C3F528DC0490524185EB51385CA23841EC51B8EE04905241B81E852B60A23841A4703DDA049052410000008066A23841A4703DFA049052411F85EB116AA23841CDCCCC0C05905241CDCCCC8C6EA23841A4703D5A05905241A4703D0A6CA238417B14AE37079052411F85EBD16AA23841A4703DCA07905241AE47E1FA68A23841E17A14FE08905241295C8F0268A238410AD7A320099052419A99999963A23841AE47E10A099052415C8FC2B55AA23841333333C30890524114AE476156A23841CDCCCCBC089052417B14AEC74DA238410AD7A3B008905241'); 292 293-- #457 -- 294SELECT '#457.1', st_astext(st_collectionExtract('POINT(0 0)', 1)); 295SELECT '#457.2', st_astext(st_collectionExtract('POINT(0 0)', 2)); 296SELECT '#457.3', st_astext(st_collectionExtract('POINT(0 0)', 3)); 297SELECT '#457.4', st_astext(st_collectionExtract('LINESTRING(0 0, 1 1)', 1)); 298SELECT '#457.5', st_astext(st_collectionExtract('LINESTRING(0 0, 1 1)', 2)); 299SELECT '#457.6', st_astext(st_collectionExtract('LINESTRING(0 0, 1 1)', 3)); 300SELECT '#457.7', st_astext(st_collectionExtract('POLYGON((0 0, 1 0, 1 1, 0 1, 0 0))', 1)); 301SELECT '#457.8', st_astext(st_collectionExtract('POLYGON((0 0, 1 0, 1 1, 0 1, 0 0))', 2)); 302SELECT '#457.9', st_astext(st_collectionExtract('POLYGON((0 0, 1 0, 1 1, 0 1, 0 0))', 3)); 303 304-- #835 -- 305SELECT '#835.1', st_astext(st_collectionExtract('POLYGON EMPTY', 1)); 306SELECT '#835.2', st_astext(st_collectionExtract('POLYGON EMPTY', 2)); 307SELECT '#835.3', st_astext(st_collectionExtract('POLYGON EMPTY', 3)); 308SELECT '#835.4', st_astext(st_collectionExtract('LINESTRING EMPTY', 1)); 309SELECT '#835.5', st_astext(st_collectionExtract('LINESTRING EMPTY', 2)); 310SELECT '#835.6', st_astext(st_collectionExtract('LINESTRING EMPTY', 3)); 311SELECT '#835.7', st_astext(st_collectionExtract('POINT EMPTY', 1)); 312SELECT '#835.8', st_astext(st_collectionExtract('POINT EMPTY', 2)); 313SELECT '#835.9', st_astext(st_collectionExtract('POINT EMPTY', 3)); 314SELECT '#835.10', st_astext(st_collectionExtract('GEOMETRYCOLLECTION EMPTY', 1)); 315SELECT '#835.11', st_astext(st_collectionExtract('GEOMETRYCOLLECTION EMPTY', 2)); 316SELECT '#835.12', st_astext(st_collectionExtract('GEOMETRYCOLLECTION EMPTY', 3)); 317 318-- #650 -- 319SELECT '#650', ST_AsText(ST_Collect(ARRAY[ST_MakePoint(0,0), ST_MakePoint(1,1), null, ST_MakePoint(2,2)])); 320 321-- #662 -- 322--SELECT '#662', ST_MakePolygon(ST_AddPoint(ST_AddPoint(ST_MakeLine(ST_SetSRID(ST_MakePointM(i+m,j,m),4326),ST_SetSRID(ST_MakePointM(j+m,i-m,m),4326)),ST_SetSRID(ST_MakePointM(i,j,m),4326)),ST_SetSRID(ST_MakePointM(i+m,j,m),4326))) As the_geom FROM generate_series(-10,50,20) As i CROSS JOIN generate_series(50,70, 20) As j CROSS JOIN generate_series(1,2) As m ORDER BY i, j, m, i*j*m LIMIT 1; 323 324-- #667 -- 325WITH vs AS (SELECT c[1]::int AS major, c[2]::int AS minor, c[3]::int AS patch 326 FROM (SELECT regexp_split_to_array(postgis_sfcgal_version(),'\.') AS c) AS r ) 327SELECT '#667', 328CASE WHEN vs.major >= 1 AND vs.minor >= 0 AND vs.patch >=5 THEN 329ST_AsEWKT(ST_LineToCurve(ST_Buffer(ST_SetSRID(ST_Point(i,j),4326), j))) 330ELSE 331'SRID=4326;CURVEPOLYGON(CIRCULARSTRING(30 40,-50 39.9999999999999,30 40))' 332END As the_geom 333FROM vs, generate_series(-10,50,10) As i CROSS JOIN generate_series(40,70, 20) As j ORDER BY i, j, i*j LIMIT 1; 334 335-- #677 -- 336SELECT '#677',round(ST_DistanceSpheroid(ST_GeomFromEWKT('MULTIPOLYGON(((-10 40,-10 55,-10 70,5 40,-10 40)))'), ST_GeomFromEWKT('MULTIPOINT(20 40,20 55,20 70,35 40,35 55,35 70,50 40,50 55,50 70)'), 'SPHEROID["GRS_1980",6378137,298.257222101]')) As result; 337 338-- #680 -- 339SELECT '#680', encode(ST_AsBinary(geography(foo1.the_geom), 'NDR'),'hex') As result FROM ((SELECT ST_SetSRID(ST_MakePointM(i,j,m),4326) As the_geom FROM generate_series(-10,50,10) As i CROSS JOIN generate_series(50,70, 20) AS j CROSS JOIN generate_series(1,2) As m ORDER BY i, j, i*j*m)) As foo1 LIMIT 1; 340 341-- #681 -- 342SELECT '#681a', ST_AsGML(ST_GeomFromText('POINT EMPTY', 4326)); 343SELECT '#681b', ST_AsGML(ST_GeomFromText('POLYGON EMPTY', 4326)); 344SELECT '#681c', ST_AsGML(ST_GeomFromText('LINESTRING EMPTY', 4326)); 345SELECT '#681d', ST_AsGML(ST_GeomFromText('MULTIPOINT EMPTY', 4326)); 346SELECT '#681e', ST_AsGML(ST_GeomFromText('MULTILINESTRING EMPTY', 4326)); 347SELECT '#681f', ST_AsGML(ST_GeomFromText('MULTIPOLYGON EMPTY', 4326)); 348SELECT '#681g', ST_AsGML(ST_GeomFromText('GEOMETRYCOLLECTION EMPTY', 4326)); 349 350-- #682 -- 351SELECT '#682', ST_AsText(ST_Buffer(ST_GeomFromText('POLYGON EMPTY',4326) , 0.5)); 352 353-- #683 -- 354SELECT '#683', ST_AsText(ST_BuildArea(ST_GeomFromText('POINT EMPTY',4326))); 355 356-- #684,#2109 -- 357SELECT '#684,#2109', ST_AsEWKT(ST_Centroid(ST_GeomFromText('POLYGON EMPTY',4326))); 358SELECT '#2109', ST_AsEWKT(ST_Centroid(ST_GeomFromText('MULTILINESTRING ZM EMPTY',3395))); 359 360-- #685 -- 361SELECT '#685', ST_AsText(ST_ConvexHull(ST_GeomFromText('POLYGON EMPTY',4326))); 362 363-- #686 -- 364SELECT '#686', St_AsText(ST_COLLECT(ST_GeomFromText('POLYGON EMPTY',4326),ST_GeomFromText('TRIANGLE EMPTY',4326))); 365 366-- #687 -- 367SELECT '#687', ST_DFullyWithin(ST_GeomFromText('LINESTRING(-10 50,50 -10)',4326), ST_GeomFromText('POLYGON EMPTY',4326),5); 368 369-- #689 -- 370SELECT '#689', ST_CoveredBy(ST_GeomFromText('POLYGON EMPTY'), ST_GeomFromText('LINESTRING(-10 50,50 -10)')); 371 372-- #690 -- 373SELECT '#690'; 374SELECT encode(ST_AsEWKB(ST_MakeLine(ST_GeomFromText('POINT(-11.1111111 40)'), ST_GeomFromText('LINESTRING(-11.1111111 70,70 -11.1111111)')),'ndr'),'hex') As result; 375 376-- #693 -- 377SELECT '#693a', ST_AsText(ST_GeomFromEWKT('SRID=4326;POLYGONM((-71.1319 42.2503 1,-71.132 42.2502 3,-71.1323 42.2504 -2,-71.1322 42.2505 1,-71.1319 42.2503 0))')); 378SELECT '#693b', ST_AsText(ST_GeomFromEWKT('SRID=4326;POLYGONM((-71.1319 42.2512 0,-71.1318 42.2511 20,-71.1317 42.2511 -20,-71.1317 42.251 5,-71.1317 42.2509 4,-71.132 42.2511 6,-71.1319 42.2512 30))')); 379 380-- #694 -- 381SELECT '#694'; 382SELECT ST_MakePolygon('POINT(1 2)'::geometry); 383 384-- #695 -- 385SELECT '#695'; 386SELECT ST_RemovePoint('POINT(-11.1111111 40)'::geometry, 1); 387 388-- #696 -- 389SELECT '#696',St_AsText(ST_Segmentize(ST_GeomFromEWKT('PolyhedralSurface( ((0 0 0, 0 0 1, 0 1 1, 0 1 0, 0 0 0)), ((0 0 0, 0 1 0, 1 1 0, 1 0 0, 0 0 0)), ((0 0 0, 1 0 0, 1 0 1, 0 0 1, 0 0 0)), ((1 1 0, 1 1 1, 1 0 1, 1 0 0, 1 1 0)), ((0 1 0, 0 1 1, 1 1 1, 1 1 0, 0 1 0)), ((0 0 1, 1 0 1, 1 1 1, 0 1 1, 0 0 1)) )'), 0.5)); 390 391-- #720 -- 392SELECT '#720', ST_AsText(ST_SnapTogrid(ST_Transform(ST_GeomFromText('MULTIPOINT(-10 40,-10 55,-10 70,5 40,5 55,5 70,20 40,20 55,20 70,35 40,35 55,35 70,50 40,50 55,50 70)',4326), 3395), 0.01)); 393 394-- #723 -- 395SELECT '#723', 396 ST_AsText(ST_SnapToGrid( ST_Intersection(a.geog, b.geog)::geometry, 0.00001)) 397FROM (VALUES (ST_GeogFromText('SRID=4326;POINT(-11.1111111 40)') ), (ST_GeogFromText('SRID=4326;POINT(-11.1111111 55)') ) ) As a(geog) CROSS JOIN ( VALUES (ST_GeogFromText('SRID=4326;POINT(-11.1111111 40)') ), (ST_GeogFromText('SRID=4326;POINT(-11.1111111 55)') )) As b(geog); 398 399-- #729 -- 400--SELECT '#729',ST_MakeLine(foo1.the_geom) As result FROM ((SELECT ST_GeomFromText('POINT EMPTY',4326) As the_geom UNION ALL SELECT ST_GeomFromText('MULTIPOINT EMPTY',4326) As the_geom UNION ALL SELECT ST_GeomFromText('MULTIPOLYGON EMPTY',4326) As the_geom UNION ALL SELECT ST_GeomFromText('LINESTRING EMPTY',4326) As the_geom UNION ALL SELECT ST_GeomFromText('MULTILINESTRING EMPTY',4326) As the_geom ) ) As foo1; 401 402-- #804 403SELECT '#804', ST_AsGML(3, 'SRID=4326;POINT(0 0)'::geometry, 0, 1); 404 405-- #845 406SELECT '#845', ST_Intersects('POINT(169.69960846592 -46.5061209281002)'::geometry, 'POLYGON((169.699607857174 -46.5061218662,169.699607857174 -46.5061195965597,169.699608806526 -46.5061195965597,169.699608806526 -46.5061218662,169.699607857174 -46.5061218662))'::geometry); 407 408-- #834 409SELECT '#834', ST_AsEWKT(ST_Intersection('LINESTRING(0 0,0 10,10 10,10 0)', 'LINESTRING(10 10 4,10 0 5,0 0 5)')); 410 411-- #884 -- 412CREATE TABLE foo (id integer, the_geom geometry); 413INSERT INTO foo VALUES (1, st_geomfromtext('MULTIPOLYGON(((-113.6 35.4,-113.6 35.8,-113.2 35.8,-113.2 35.4,-113.6 35.4),(-113.5 35.5,-113.3 35.5,-113.3 35.7,-113.5 35.7,-113.5 35.5)))')); 414INSERT INTO foo VALUES (2, st_geomfromtext('MULTIPOLYGON(((-113.7 35.3,-113.7 35.9,-113.1 35.9,-113.1 35.3,-113.7 35.3),(-113.6 35.4,-113.2 35.4,-113.2 35.8,-113.6 35.8,-113.6 35.4)),((-113.5 35.5,-113.5 35.7,-113.3 35.7,-113.3 35.5,-113.5 35.5)))')); 415 416select '#884', id, ST_Within( 417ST_GeomFromText('POINT (-113.4 35.6)'), the_geom 418) from foo; 419 420select '#938', 'POLYGON EMPTY'::geometry::box2d; 421 422DROP TABLE foo; 423 424-- #668 -- 425select '#668',box2d('CIRCULARSTRING(10 2,12 2,14 2)'::geometry) as b; 426 427-- #711 -- 428select '#711', ST_GeoHash(ST_GeomFromText('POLYGON EMPTY',4326)); 429 430-- #712 -- 431SELECT '#712',ST_IsValid(ST_GeomFromText('POLYGON EMPTY',4326)); 432 433-- #756 434WITH inp AS ( SELECT 'LINESTRING(0 0, 1 1)'::geometry as s, 435 'LINESTRING EMPTY'::geometry as e ) 436 SELECT '#756.1', ST_Equals(s, st_multi(s)), 437 ST_Equals(s, st_collect(s, e)) 438 FROM inp; 439 440-- #1023 -- 441select '#1023', 'POINT(10 4)'::geometry = 'POINT(10 4)'::geometry; 442select '#1023.a', 'POINT(10 4)'::geometry = 'POINT(10 5)'::geometry; 443select '#1023.b', postgis_addbbox('POINT(10 4)'::geometry) = 'POINT(10 4)'::geometry; 444 445-- #1069 -- 446select '#1060', ST_Relate(ST_GeomFromText('POINT EMPTY',4326), ST_GeomFromText('POINT EMPTY',4326)) As result; 447 448-- #1273 -- 449WITH p AS ( SELECT 'POINT(832694.188 816254.625)'::geometry as g ) 450SELECT '#1273', st_equals(p.g, postgis_addbbox(p.g)) from p; 451 452-- Another for #1273 -- 453WITH p AS ( SELECT 'MULTIPOINT((832694.188 816254.625))'::geometry as g ) 454SELECT '#1273.1', st_equals(p.g, postgis_dropbbox(p.g)) from p; 455 456-- #877, #818 457create table t(g geometry); 458select '#877.1', ST_EstimatedExtent('t','g'); 459analyze t; 460select '#877.2', ST_EstimatedExtent('public', 't','g'); 461select '#877.2.deprecated', ST_Estimated_Extent('public', 't','g'); 462insert into t(g) values ('LINESTRING(-10 -50, 20 30)'); 463 464-- #877.3 465select '#877.3', 466 round(st_xmin(e.e)::numeric, 5), round(st_xmax(e.e)::numeric, 5), 467 round(st_ymin(e.e)::numeric, 5), round(st_ymax(e.e)::numeric, 5) 468from ( select ST_EstimatedExtent('t','g') as e limit 1 ) e; 469 470-- #877.4 471analyze t; 472select '#877.4', 473 round(st_xmin(e.e)::numeric, 5), round(st_xmax(e.e)::numeric, 5), 474 round(st_ymin(e.e)::numeric, 5), round(st_ymax(e.e)::numeric, 5) 475from ( select ST_EstimatedExtent('t','g') as e limit 1 ) e; 476 477-- #877.5 478truncate t; 479select '#877.5', 480 round(st_xmin(e.e)::numeric, 5), round(st_xmax(e.e)::numeric, 5), 481 round(st_ymin(e.e)::numeric, 5), round(st_ymax(e.e)::numeric, 5) 482from ( select ST_EstimatedExtent('t','g') as e limit 1 ) e; 483drop table t; 484 485-- #1292 486SELECT '#1292', ST_AsText(ST_SnapToGrid(ST_GeomFromText( 487 'GEOMETRYCOLLECTION(POINT(180 90),POLYGON((140 50,150 50,180 50,140 50),(140 60,150 60,180 60,140 60)))' 488 , 4326), 0.00001)::geography); 489 490-- #1292.1 491SELECT '#1292.1', ST_AsText(ST_GeomFromText('POINT(180.00000000001 95)')::geography), 492 ST_AsText(ST_GeomFromText('POINT(185 90.00000000001)')::geography); 493 494-- #1320 495SELECT '<#1320>'; 496CREATE TABLE A ( geom geometry(MultiPolygon, 4326), 497 geog geography(MultiPolygon, 4326) ); 498-- Valid inserts 499INSERT INTO a(geog) VALUES('SRID=4326;MULTIPOLYGON (((0 0, 10 0, 10 10, 0 0)))'::geography); 500INSERT INTO a(geom) VALUES('SRID=4326;MULTIPOLYGON (((0 0, 10 0, 10 10, 0 0)))'::geometry); 501SELECT '#1320.geog.1', geometrytype(geog::geometry), st_srid(geog::geometry) FROM a where geog is not null; 502SELECT '#1320.geom.1', geometrytype(geom), st_srid(geom) FROM a where geom is not null; 503-- Type mismatches is not allowed 504INSERT INTO a(geog) VALUES('SRID=4326;POLYGON ((0 0, 10 0, 10 10, 0 0))'::geography); 505INSERT INTO a(geom) VALUES('SRID=4326;POLYGON ((0 0, 10 0, 10 10, 0 0))'::geometry); 506SELECT '#1320.geog.2', geometrytype(geog::geometry), st_srid(geog::geometry) FROM a where geog is not null; 507SELECT '#1320.geom.2', geometrytype(geom), st_srid(geom) FROM a where geom is not null; 508-- Even if it's a trigger changing the type 509CREATE OR REPLACE FUNCTION triga() RETURNS trigger AS 510$$ BEGIN 511 NEW.geom = ST_GeometryN(New.geom,1); 512 NEW.geog = ST_GeometryN(New.geog::geometry,1)::geography; 513 RETURN NEW; 514END; $$ language plpgsql VOLATILE; 515CREATE TRIGGER triga_before 516 BEFORE INSERT ON a FOR EACH ROW 517 EXECUTE PROCEDURE triga(); 518INSERT INTO a(geog) VALUES('SRID=4326;MULTIPOLYGON (((0 0, 10 0, 10 10, 0 0)))'::geography); 519INSERT INTO a(geom) VALUES('SRID=4326;MULTIPOLYGON (((0 0, 10 0, 10 10, 0 0)))'::geometry); 520SELECT '#1320.geog.3', geometrytype(geog::geometry), st_srid(geog::geometry) FROM a where geog is not null; 521SELECT '#1320.geom.3', geometrytype(geom), st_srid(geom) FROM a where geom is not null; 522DROP TABLE A; 523DROP FUNCTION triga(); 524SELECT '</#1320>'; 525 526-- st_AsText POLYGON((0 0,10 0,10 10,0 0)) 527 528-- #1344 529select '#1344', octet_length(ST_AsEWKB(st_makeline(g))) FROM ( values ('POINT(0 0)'::geometry ) ) as foo(g); 530 531-- #1385 532SELECT '#1385', ST_Extent(g) FROM ( select null::geometry as g ) as foo; 533 534-- #657 535SELECT '#657.1',Round(ST_X(ST_Project('POINT(175 10)'::geography, 2000000, 3.1415/2)::GEOMETRY)::numeric,2); 536SELECT '#657.2',Round(ST_Distance(ST_Project('POINT(10 10)'::geography, 10, 0), 'POINT(10 10)'::geography)::numeric,2); 537SELECT '#657.3',ST_DWithin(ST_Project('POINT(10 10)'::geography, 2000, pi()/2), 'POINT(10 10)'::geography, 2000); 538 539-- #1305 540SELECT '#1305.1',ST_AsText(ST_Project('POINT(10 10)'::geography, 0, 0)); 541WITH pts AS ( SELECT 'POINT(0 45)'::geography AS s, 'POINT(45 45)'::geography AS e ) 542SELECT '#1305.2',abs(ST_Distance(e, ST_Project(s, ST_Distance(s, e), ST_Azimuth(s, e)))) < 0.001 FROM pts; 543SELECT '#1305.3',ST_Azimuth('POINT(0 45)'::geography, 'POINT(0 45)'::geography) IS NULL; 544 545-- #1445 546SELECT '01060000400200000001040000400100000001010000400000000000000000000000000000000000000000000000000101000040000000000000F03F000000000000F03F000000000000F03F'::geometry; 547SELECT '01050000400200000001040000400100000001010000400000000000000000000000000000000000000000000000000101000040000000000000F03F000000000000F03F000000000000F03F'::geometry; 548SELECT '01040000400200000001040000400100000001010000400000000000000000000000000000000000000000000000000101000040000000000000F03F000000000000F03F000000000000F03F'::geometry; 549SELECT '01090000400200000001040000400100000001010000400000000000000000000000000000000000000000000000000101000040000000000000F03F000000000000F03F000000000000F03F'::geometry; 550SELECT '010B0000400200000001040000400100000001010000400000000000000000000000000000000000000000000000000101000040000000000000F03F000000000000F03F000000000000F03F'::geometry; 551SELECT '010C0000400200000001040000400100000001010000400000000000000000000000000000000000000000000000000101000040000000000000F03F000000000000F03F000000000000F03F'::geometry; 552 553-- #1453 554SELECT '#1453.1', ST_OrderingEquals('POINT EMPTY', 'POINT EMPTY'); 555SELECT '#1453.2', ST_OrderingEquals('POINT EMPTY', 'POINT Z EMPTY'); 556 557-- #1454 558with inp as ( select 'MULTILINESTRING((0 0, 2 0))'::geometry as g ) 559SELECT '#1454', st_orderingequals(g,g) from inp; 560 561-- #1414 562SELECT '#1414', st_astext(st_Force3DZ('CURVEPOLYGON EMPTY')); 563 564-- #1478 565SELECT '#1478', ST_AsText('SRID=1;POINT EMPTY'::geometry::text::geometry); 566 567-- #745 568SELECT '#745', ST_AsEWKT(ST_Split('POLYGON((-72 42 1,-70 43 1,-71 41 1,-72 42 1))', 569 'LINESTRING(-10 40 1,-9 41 1)')); 570 571-- #1450 572SELECT '#1450', GeometryType('POINT(0 0)'::geography), GeometryType('POLYGON EMPTY'::geography); 573 574-- #1482 575select '#1482', ST_Srid('POINT(0 0)'::geography(point, 0)::geometry); 576 577-- #852 578CREATE TABLE cacheable (id int, g geometry); 579COPY cacheable FROM STDIN; 5801 POINT(0.5 0.5000000000001) 5812 POINT(0.5 0.5000000000001) 582\. 583select '#852.1', id, -- first run is not cached, consequent are cached 584 st_intersects(g, 'POLYGON((0 0, 10 10, 1 0, 0 0))'::geometry), 585 st_intersects(g, 'POLYGON((0 0, 1 1, 1 0, 0 0))'::geometry) from cacheable; 586UPDATE cacheable SET g = 'POINT(0.5 0.5)'; 587-- New select, new cache 588select '#852.2', id, -- first run is not cached, consequent are cached 589 st_intersects(g, 'POLYGON((0 0, 10 10, 1 0, 0 0))'::geometry), 590 st_intersects(g, 'POLYGON((0 0, 1 1, 1 0, 0 0))'::geometry) from cacheable; 591DROP TABLE cacheable; 592 593-- #1489 594with inp AS ( SELECT 595 st_multi('POINT EMPTY'::geometry) as mp, 596 st_multi('LINESTRING EMPTY'::geometry) as ml, 597 st_multi('POLYGON EMPTY'::geometry) as ma, 598 st_multi('GEOMETRYCOLLECTION EMPTY'::geometry) as mm 599) select '#1489', 600 st_astext(mp), st_numgeometries(mp), 601 st_astext(ml), st_numgeometries(ml), 602 st_astext(ma), st_numgeometries(ma), 603 st_astext(mm), st_numgeometries(mm) 604FROM inp; 605 606-- #1150 607insert into spatial_ref_sys (srid, proj4text) values (500001,NULL); 608insert into spatial_ref_sys (srid, proj4text) values (500002, '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'); 609select '#1150', st_astext(st_transform('SRID=500002;POINT(0 0)'::geometry,500001)); 610 611-- #1038 612select '#1038', ST_AsSVG('POLYGON EMPTY'::geometry); 613 614-- #1042 615select '#1042',round((st_ymax(st_minimumboundingcircle('LINESTRING(-1 -1, 1 1)')) * st_xmax(st_minimumboundingcircle('LINESTRING(-1 -1, 1 1)')))::numeric,0); 616 617-- #1170 -- 618SELECT '#1170', ST_Y(ST_Intersection( ST_GeogFromText( 'POINT(0 90)'), ST_GeogFromText( 'POINT(0 90)' ))::geometry); 619 620-- #1264 -- 621SELECT '#1264', ST_DWithin('POLYGON((-10 -10, -10 10, 10 10, 10 -10, -10 -10))'::geography, 'POINT(0 0)'::geography, 0); 622 623-- #1398 624select '#1398a', st_astext(st_snaptogrid(st_project('POINT(-120 45)'::geography, 100000, radians(45))::geometry,0.000001)); 625select '#1398b', st_astext(st_snaptogrid(st_project('POINT(20 85)'::geography, 2000000, radians(0.1))::geometry,0.000001)); 626 627-- #1543 628with inp as ( select 629'0105000000020000000102000000040000000000000000000000000000000000000000000000000024400000000000000000000000000000244000000000000024400000000000000000000000000000000001020000000100000000000000000000000000000000000000' 630::geometry as g ) 631select '#1543', st_astext(g), st_astext(st_buildarea(g)) from inp; 632 633-- #1580 634select '#1580.1', ST_Summary(ST_Transform('SRID=4326;POINT(0 0)'::geometry, 3395)); 635select '#1580.2', ST_Transform('SRID=4326;POINT(180 90)'::geometry, 3395); -- fails 636select '#1580.3', ST_Summary(ST_Transform('SRID=4326;POINT(0 0)'::geometry, 3395)); 637 638-- #1596 -- 639CREATE TABLE road_pg (ID INTEGER, NAME VARCHAR(32)); 640SELECT '#1596.1', AddGeometryColumn( 'road_pg','roads_geom', 3395, 'POINT', 2 ); 641SELECT '#1596.2', UpdateGeometrySRID( 'road_pg','roads_geom', 330000); 642SELECT '#1596.3', srid FROM geometry_columns 643 WHERE f_table_name = 'road_pg' AND f_geometry_column = 'roads_geom'; 644SELECT '#1596.4', UpdateGeometrySRID( 'road_pg','roads_geom', 999000); 645SELECT '#1596.5', srid FROM geometry_columns 646 WHERE f_table_name = 'road_pg' AND f_geometry_column = 'roads_geom'; 647SELECT '#1596.6', UpdateGeometrySRID( 'road_pg','roads_geom', -1); 648SELECT '#1596.7', srid FROM geometry_columns 649 WHERE f_table_name = 'road_pg' AND f_geometry_column = 'roads_geom'; 650DROP TABLE road_pg; 651 652-- #1596 653WITH inp AS ( SELECT 654 'POLYGON((-176 -22,-176 -21,-175 -21,-175 -22,-176 -22))'::geography as a, 655 'POINT(-176 -22)'::geography as p 656) SELECT '#1596', ST_Summary(ST_Intersection(a,p)) FROM inp; 657 658-- #1695 659SELECT '#1695', ST_AsEWKT(ST_SnapToGrid('MULTIPOLYGON(((0 0, 10 0, 10 10, 0 10, 0 0)))'::geometry, 20)); 660 661-- #1697 -- 662CREATE TABLE eg(g geography, gm geometry); 663CREATE INDEX egi on eg using gist (g); 664CREATE INDEX egind on eg using gist (gm gist_geometry_ops_nd); 665INSERT INTO eg (g, gm) 666 select 'POINT EMPTY'::geography, 667 'POINT EMPTY'::geometry 668 from generate_series(1,1024); 669SELECT '#1697.1', count(*) FROM eg WHERE g && 'POINT(0 0)'::geography; 670SELECT '#1697.2', count(*) FROM eg WHERE gm && 'POINT(0 0)'::geometry; 671SELECT '#1697.3', count(*) FROM eg WHERE gm ~= 'POINT EMPTY'::geometry; 672DROP TABLE eg; 673 674-- #1734 -- 675create table eg (g geography); 676create index egi on eg using gist (g); 677INSERT INTO eg(g) VALUES (NULL); 678INSERT INTO eg (g) VALUES ('POINT(0 0)'::geography); 679INSERT INTO eg (g) select 'POINT(0 0)'::geography 680 FROM generate_series(1,1024); 681SELECT '#1734.1', count(*) FROM eg; 682DROP table eg; 683 684-- #1755 -- 685select '#1755', ST_AsText(st_geographyFromText('SRID=4326;Point(85 35 0)')); 686 687-- #1776 -- 688with inp as ( SELECT 689 'POLYGON EMPTY'::geometry as A, 690 'POLYGON((0 0, 10 0, 10 10, 0 0))'::geometry as B ) 691SELECT '#1776', 692 ST_AsText(ST_SymDifference(A,B)), ST_AsText(ST_SymDifference(B, A)) 693FROM inp; 694 695-- #1780 -- 696SELECT '#1780',ST_GeoHash('POINT(4 4)'::geometry) = ST_GeoHash('POINT(4 4)'::geography); 697 698-- #1791 -- 699with inp as ( SELECT 700 '010100000000000000004065C0041AD965BE5554C0'::geometry as a, 701 '010100000001000000004065C0041AD965BE5554C0'::geometry as b 702) SELECT '#1791', round(ST_Azimuth(a,b)*10)/10 from inp; 703 704-- #1799 -- 705SELECT '#1799', ST_Segmentize('LINESTRING(0 0, 10 0)'::geometry, 0); 706 707-- #1936 -- 708select st_astext(st_geomfromgml( 709 '<gml:Polygon xmlns:gml="http://www.opengis.net/gml/3.2" 710 gml:id="HPA.15449990010" srsName="urn:ogc:def:crs:EPSG::4326" 711 srsDimension="2"> 712 <gml:exterior> 713 <gml:Ring> 714 <gml:curveMember> 715 <gml:LineString gml:id="HPA.15449990010.1"> 716 <gml:posList>711540.35 1070163.61 711523.82 1070166.54 711521.30 1070164.14 711519.52 1070162.44 711518.57 1070164.62 712154.47 1070824.94</gml:posList> 717 </gml:LineString> 718 </gml:curveMember> 719 <gml:curveMember> 720 <gml:Curve gml:id="HPA.15449990010.2"> 721 <gml:segments><gml:ArcString> 722 <gml:posList>712154.47 1070824.94 712154.98 1070826.04 712154.41 1070827.22</gml:posList> 723 </gml:ArcString> 724 </gml:segments> 725 </gml:Curve> 726 </gml:curveMember> 727 <gml:curveMember> 728 <gml:LineString gml:id="HPA.15449990010.3"> 729 <gml:posList>712154.41 1070827.22 712160.31 1070837.07 712160.92 1070835.36 712207.89 1071007.95</gml:posList> 730 </gml:LineString> 731 </gml:curveMember> 732 <gml:curveMember> 733 <gml:Curve gml:id="HPA.15449990010.4"><gml:segments><gml:ArcString><gml:posList>712207.89 1071007.95 712207.48 1071005.59 712208.38 1071001.28</gml:posList></gml:ArcString></gml:segments></gml:Curve></gml:curveMember><gml:curveMember><gml:LineString gml:id="HPA.15449990010.5"><gml:posList>712208.38 1071001.28 712228.74 1070949.67 712233.98 1070936.15 712124.93 1070788.72</gml:posList></gml:LineString></gml:curveMember><gml:curveMember><gml:Curve gml:id="HPA.15449990010.6"><gml:segments><gml:ArcString><gml:posList>712124.93 1070788.72 712124.28 1070785.87 712124.63 1070783.38</gml:posList></gml:ArcString></gml:segments></gml:Curve></gml:curveMember><gml:curveMember><gml:LineString gml:id="HPA.15449990010.7"><gml:posList>712124.63 1070783.38 712141.04 1070764.12 712146.60 1070757.01 711540.35 1070163.61</gml:posList></gml:LineString></gml:curveMember></gml:Ring></gml:exterior> 734 <gml:interior> 735 <gml:LinearRing> 736 <gml:posList>713061.62 1070354.46 713053.59 1070335.12 713049.58 1070315.92 713049.65 1070298.33 713061.62 1070354.46</gml:posList> 737 </gml:LinearRing> 738 </gml:interior> 739 </gml:Polygon>')); 740 741-- #1978 -- 742SELECT '#1978', round(ST_Length(ST_GeomFromText('CIRCULARSTRING(0 0,1 0,0 0)',0))::numeric,4); 743 744-- #1996 -- 745SELECT '#1996', ST_AsGeoJSON(ST_GeomFromText('POINT EMPTY')); 746 747-- #2001 -- 748SELECT '#2001', ST_AsText(ST_CurveToLine(ST_GeomFromText('CURVEPOLYGON((0 0, 0 1, 1 1, 0 0))'), 2)); 749 750-- #2028 -- 751SELECT '#2028', ST_AsText(ST_Multi('TRIANGLE((0 0, 0 1, 1 1, 0 0))')); 752 753-- #2035 START ------------------------------------------------------------ 754 755-- Simple geographic table, with single point. 756CREATE TABLE "city" ( 757 "id" integer, 758 "name" varchar(30) NOT NULL, 759 "point" geometry(POINT,4326) NOT NULL 760); 761CREATE INDEX "city_point_id" ON "city" USING GIST ( "point" ); 762 763-- Initial data, with points around the world. 764INSERT INTO "city" (id, name, point) VALUES (1, 'Houston', 'SRID=4326;POINT(-95.363151 29.763374)'); 765INSERT INTO "city" (id, name, point) VALUES (2, 'Dallas', 'SRID=4326;POINT(-95.363151 29.763374)'); 766INSERT INTO "city" (id, name, point) VALUES (3, 'Oklahoma City', 'SRID=4326;POINT(-97.521157 34.464642)'); 767INSERT INTO "city" (id, name, point) VALUES (4, 'Wellington', 'SRID=4326;POINT(174.783117 -41.315268)'); 768INSERT INTO "city" (id, name, point) VALUES (5, 'Pueblo', 'SRID=4326;POINT(-104.609252 38.255001)'); 769INSERT INTO "city" (id, name, point) VALUES (6, 'Lawrence', 'SRID=4326;POINT(-95.23506 38.971823)'); 770INSERT INTO "city" (id, name, point) VALUES (7, 'Chicago', 'SRID=4326;POINT(-87.650175 41.850385)'); 771INSERT INTO "city" (id, name, point) VALUES (8, 'Victoria', 'SRID=4326;POINT(-123.305196 48.462611)'); 772 773-- This query, or COUNT(*), does not return anything; should return 6 cities, 774-- excluding Pueblo and Victoria. The Polygon is a simple approximation of 775-- Colorado. 776SELECT '#2035a', Count(*) FROM "city" 777 WHERE "city"."point" >> ST_GeomFromEWKT('SRID=4326;POLYGON ((-109.060253 36.992426, -109.060253 41.003444, -102.041524 41.003444, -102.041524 36.992426, -109.060253 36.992426))'); 778 779-- However, when a LIMIT is placed on statement, the query suddenly works. 780SELECT '#2035b', Count(*) FROM "city" 781 WHERE "city"."point" >> ST_GeomFromEWKT('SRID=4326;POLYGON ((-109.060253 36.992426, -109.060253 41.003444, -102.041524 41.003444, -102.041524 36.992426, -109.060253 36.992426))') LIMIT 6; 782 783DROP TABLE "city"; 784-- #2035 END -------------------------------------------------------------- 785 786-- #2084 -- 787SELECT '#2048', num, ST_Within('POINT(-54.394 56.522)', "the_geom"), ST_CoveredBy('POINT(-54.394 56.522)', "the_geom") 788FROM ( VALUES 789(1, '0103000000010000000E00000051C6F7C5A5324BC02EB69F8CF13F4C40F12EA4C343364BC0326AA2CF47434C402BC1A8A44E364BC02A50E10852434C407F2990D959364BC0A0D1730B5D434C404102452C62364BC0ECF335CB65434C400903232F6B364BC0F635E84B6F434C40BD0CC51D6F364BC0D2805EB873434C40B9E6E26F7B364BC0F20B93A982434C40D9FAAF73D3344BC0FE84D04197444C40BD5C8AABCA344BC0CED05CA791444C4023F2237EC5344BC02A84F23E8E444C40BDCDD8077B324BC0C60FB90F01434C409FD1702E65324BC04EF1915C17404C4051C6F7C5A5324BC02EB69F8CF13F4C40'::geometry), 790(2, '0103000000010000001C00000003F25650F73B4BC098477F523E3E4C40C9A6A344CE3C4BC0C69698653E3E4C40BDD0E979373E4BC0081FA0FB723E4C400FD252793B3E4BC01A137F14753E4C40537170E998414BC070D3BCE314414C4023FC51D499474BC0D4D100DE024F4C40638C47A984454BC024130D52F0504C40B9442DCDAD404BC03A29E96168554C40C7108DEE20404BC07C7C26FBE7554C40195D6BEF533F4BC0E20391459A564C40239FE40E9B344BC08C1ADB6B41514C40132D3F7095314BC0BA2ADF33124F4C409DB91457952D4BC02C7B681F2B4C4C4089DC60A8C32C4BC07C5C3810924B4C40D7ED409DF22A4BC0F64389963C4A4C405D1EF818AC2A4BC00EC84274084A4C401B48A46DFC294BC0B271A8DF85494C40E78AA6B393294BC01ED0EFFB37494C4081C64B3789294BC0DC5BE7DF2E494C409B23329287294BC0F0D6974E2D494C40CD22D5D687294BC0844316D72C494C40F5229D4FE2294BC002F19825AB484C40A3D0BD5AE9294BC06C0776A9A2484C409FD1702E65324BC04EF1915C17404C409F860AA7BD324BC0162CA390E33F4C40539A5C1C23334BC0FE86B04EB03F4C4081511DFF90334BC088FF36D4873F4C4003F25650F73B4BC098477F523E3E4C40'::geometry), 791(3, '010300000001000000100000008D57CD101A214BC0AECDD34E072C4C400DBB72E6EC274BC0A8088D60E32C4C40CF8FD7E6734E4BC0B22695BE4A324C40BFA74213934F4BC020BE505D4C354C4057CD4BEE454E4BC0BA6CF3940F3D4C40E7BDC5FD263E4BC09A4B297D5B484C4073A46A86701C4BC0B287F08D93364C4045501F86701C4BC05EBDB78D93364C40A37DB6586D1C4BC0841E7D2891364C409FBF445F6D1C4BC01E225C5690364C40D1BA97726D1C4BC06E2AF7EA8D364C4019B60C9B751C4BC0D2FD702575364C40FDE4394B5E1F4BC08C40F231CC2F4C402343DF40F51F4BC022008E3D7B2E4C400BB57B45F9204BC0908CE2EA3A2C4C408D57CD101A214BC0AECDD34E072C4C40'::geometry) 792) AS f(num, the_geom); 793 794-- #2112 -- Start 795SELECT '#2112a', ST_3DDistance(a,b), ST_ASEWKT(ST_3DShortestLine(a,b)) 796FROM (SELECT 'POLYGON((1 1 1, 5 1 1,5 5 1, 1 5 1,1 1 1))'::geometry as a, 'LINESTRING(0 0 2, 0 0 0,5 5 2)'::geometry as b 797 ) as foo; 798 799SELECT '#2112b', ST_3DDistance(a,b), ST_ASEWKT(ST_3DShortestLine(a,b)) 800FROM (SELECT 'POLYGON((1 1 1, 5 1 1,5 5 1, 1 5 1,1 1 1))'::geometry as a, 'LINESTRING(1 0 2, 1 0 0,5 5 -1)'::geometry as b 801 ) as foo; 802-- 2112 -- End 803 804SELECT '#2108', ST_AsEWKT(ST_LineInterpolatePoint('SRID=3395;LINESTRING M EMPTY'::geometry, 0.5)); 805SELECT '#2117', ST_AsEWKT(ST_PointOnSurface('SRID=3395;MULTIPOLYGON M EMPTY'::geometry)); 806 807SELECT '#2110.1', 'POINT(0 0)'::geometry = 'POINT EMPTY'::geometry; 808SELECT '#2110.2', 'POINT EMPTY'::geometry = 'POINT EMPTY'::geometry; 809SELECT '#2110.3', 'POINT(0 0)'::geometry = 'POINT(0 0)'::geometry; 810 811SELECT '#2145', 812round(ST_Length(St_Segmentize(ST_GeographyFromText('LINESTRING(-89.3000030518 28.2000007629,-89.1999969482 89.1999969482,-89.1999969482 89.1999969482)'), 10000))::numeric,0); 813 814-- Clean up 815DELETE FROM spatial_ref_sys; 816