contrib/earthdistance/expected/earthdistance.out

   1 --
   2 --  Test earth distance functions
   3 --
   4 --
   5 -- first, define the datatype.  Turn off echoing so that expected file
   6 -- does not depend on contents of earthdistance.sql or cube.sql.
   7 --
   8 SET client_min_messages = warning;
   9 \set ECHO none
  10 RESET client_min_messages;
  11 --
  12 -- The radius of the Earth we are using.
  13 --
  14 SELECT earth()::numeric(20,5);
  15      earth
  16 ---------------
  17  6378168.00000
  18 (1 row)
  19
  20 --
  21 -- Convert straight line distances to great circle distances.
  22 --
  23 SELECT (pi()*earth())::numeric(20,5);
  24     numeric
  25 ----------------
  26  20037605.73216
  27 (1 row)
  28
  29 SELECT sec_to_gc(0)::numeric(20,5);
  30  sec_to_gc
  31 -----------
  32    0.00000
  33 (1 row)
  34
  35 SELECT sec_to_gc(2*earth())::numeric(20,5);
  36    sec_to_gc
  37 ----------------
  38  20037605.73216
  39 (1 row)
  40
  41 SELECT sec_to_gc(10*earth())::numeric(20,5);
  42    sec_to_gc
  43 ----------------
  44  20037605.73216
  45 (1 row)
  46
  47 SELECT sec_to_gc(-earth())::numeric(20,5);
  48  sec_to_gc
  49 -----------
  50    0.00000
  51 (1 row)
  52
  53 SELECT sec_to_gc(1000)::numeric(20,5);
  54  sec_to_gc
  55 ------------
  56  1000.00000
  57 (1 row)
  58
  59 SELECT sec_to_gc(10000)::numeric(20,5);
  60   sec_to_gc
  61 -------------
  62  10000.00102
  63 (1 row)
  64
  65 SELECT sec_to_gc(100000)::numeric(20,5);
  66   sec_to_gc
  67 --------------
  68  100001.02426
  69 (1 row)
  70
  71 SELECT sec_to_gc(1000000)::numeric(20,5);
  72    sec_to_gc
  73 ---------------
  74  1001027.07131
  75 (1 row)
  76
  77 --
  78 -- Convert great circle distances to straight line distances.
  79 --
  80 SELECT gc_to_sec(0)::numeric(20,5);
  81  gc_to_sec
  82 -----------
  83    0.00000
  84 (1 row)
  85
  86 SELECT gc_to_sec(sec_to_gc(2*earth()))::numeric(20,5);
  87    gc_to_sec
  88 ----------------
  89  12756336.00000
  90 (1 row)
  91
  92 SELECT gc_to_sec(10*earth())::numeric(20,5);
  93    gc_to_sec
  94 ----------------
  95  12756336.00000
  96 (1 row)
  97
  98 SELECT gc_to_sec(pi()*earth())::numeric(20,5);
  99    gc_to_sec
 100 ----------------
 101  12756336.00000
 102 (1 row)
 103
 104 SELECT gc_to_sec(-1000)::numeric(20,5);
 105  gc_to_sec
 106 -----------
 107    0.00000
 108 (1 row)
 109
 110 SELECT gc_to_sec(1000)::numeric(20,5);
 111  gc_to_sec
 112 ------------
 113  1000.00000
 114 (1 row)
 115
 116 SELECT gc_to_sec(10000)::numeric(20,5);
 117  gc_to_sec
 118 ------------
 119  9999.99898
 120 (1 row)
 121
 122 SELECT gc_to_sec(100000)::numeric(20,5);
 123   gc_to_sec
 124 -------------
 125  99998.97577
 126 (1 row)
 127
 128 SELECT gc_to_sec(1000000)::numeric(20,5);
 129   gc_to_sec
 130 --------------
 131  998976.08618
 132 (1 row)
 133
 134 --
 135 -- Set coordinates using latitude and longitude.
 136 -- Extract each coordinate separately so we can round them.
 137 --
 138 SELECT cube_ll_coord(ll_to_earth(0,0),1)::numeric(20,5),
 139  cube_ll_coord(ll_to_earth(0,0),2)::numeric(20,5),
 140  cube_ll_coord(ll_to_earth(0,0),3)::numeric(20,5);
 141  cube_ll_coord | cube_ll_coord | cube_ll_coord
 142 ---------------+---------------+---------------
 143  6378168.00000 |       0.00000 |       0.00000
 144 (1 row)
 145
 146 SELECT cube_ll_coord(ll_to_earth(360,360),1)::numeric(20,5),
 147  cube_ll_coord(ll_to_earth(360,360),2)::numeric(20,5),
 148  cube_ll_coord(ll_to_earth(360,360),3)::numeric(20,5);
 149  cube_ll_coord | cube_ll_coord | cube_ll_coord
 150 ---------------+---------------+---------------
 151  6378168.00000 |       0.00000 |       0.00000
 152 (1 row)
 153
 154 SELECT cube_ll_coord(ll_to_earth(180,180),1)::numeric(20,5),
 155  cube_ll_coord(ll_to_earth(180,180),2)::numeric(20,5),
 156  cube_ll_coord(ll_to_earth(180,180),3)::numeric(20,5);
 157  cube_ll_coord | cube_ll_coord | cube_ll_coord
 158 ---------------+---------------+---------------
 159  6378168.00000 |       0.00000 |       0.00000
 160 (1 row)
 161
 162 SELECT cube_ll_coord(ll_to_earth(180,360),1)::numeric(20,5),
 163  cube_ll_coord(ll_to_earth(180,360),2)::numeric(20,5),
 164  cube_ll_coord(ll_to_earth(180,360),3)::numeric(20,5);
 165  cube_ll_coord  | cube_ll_coord | cube_ll_coord
 166 ----------------+---------------+---------------
 167  -6378168.00000 |       0.00000 |       0.00000
 168 (1 row)
 169
 170 SELECT cube_ll_coord(ll_to_earth(-180,-360),1)::numeric(20,5),
 171  cube_ll_coord(ll_to_earth(-180,-360),2)::numeric(20,5),
 172  cube_ll_coord(ll_to_earth(-180,-360),3)::numeric(20,5);
 173  cube_ll_coord  | cube_ll_coord | cube_ll_coord
 174 ----------------+---------------+---------------
 175  -6378168.00000 |       0.00000 |       0.00000
 176 (1 row)
 177
 178 SELECT cube_ll_coord(ll_to_earth(0,180),1)::numeric(20,5),
 179  cube_ll_coord(ll_to_earth(0,180),2)::numeric(20,5),
 180  cube_ll_coord(ll_to_earth(0,180),3)::numeric(20,5);
 181  cube_ll_coord  | cube_ll_coord | cube_ll_coord
 182 ----------------+---------------+---------------
 183  -6378168.00000 |       0.00000 |       0.00000
 184 (1 row)
 185
 186 SELECT cube_ll_coord(ll_to_earth(0,-180),1)::numeric(20,5),
 187  cube_ll_coord(ll_to_earth(0,-180),2)::numeric(20,5),
 188  cube_ll_coord(ll_to_earth(0,-180),3)::numeric(20,5);
 189  cube_ll_coord  | cube_ll_coord | cube_ll_coord
 190 ----------------+---------------+---------------
 191  -6378168.00000 |       0.00000 |       0.00000
 192 (1 row)
 193
 194 SELECT cube_ll_coord(ll_to_earth(90,0),1)::numeric(20,5),
 195  cube_ll_coord(ll_to_earth(90,0),2)::numeric(20,5),
 196  cube_ll_coord(ll_to_earth(90,0),3)::numeric(20,5);
 197  cube_ll_coord | cube_ll_coord | cube_ll_coord
 198 ---------------+---------------+---------------
 199        0.00000 |       0.00000 | 6378168.00000
 200 (1 row)
 201
 202 SELECT cube_ll_coord(ll_to_earth(90,180),1)::numeric(20,5),
 203  cube_ll_coord(ll_to_earth(90,180),2)::numeric(20,5),
 204  cube_ll_coord(ll_to_earth(90,180),3)::numeric(20,5);
 205  cube_ll_coord | cube_ll_coord | cube_ll_coord
 206 ---------------+---------------+---------------
 207        0.00000 |       0.00000 | 6378168.00000
 208 (1 row)
 209
 210 SELECT cube_ll_coord(ll_to_earth(-90,0),1)::numeric(20,5),
 211  cube_ll_coord(ll_to_earth(-90,0),2)::numeric(20,5),
 212  cube_ll_coord(ll_to_earth(-90,0),3)::numeric(20,5);
 213  cube_ll_coord | cube_ll_coord | cube_ll_coord
 214 ---------------+---------------+----------------
 215        0.00000 |       0.00000 | -6378168.00000
 216 (1 row)
 217
 218 SELECT cube_ll_coord(ll_to_earth(-90,180),1)::numeric(20,5),
 219  cube_ll_coord(ll_to_earth(-90,180),2)::numeric(20,5),
 220  cube_ll_coord(ll_to_earth(-90,180),3)::numeric(20,5);
 221  cube_ll_coord | cube_ll_coord | cube_ll_coord
 222 ---------------+---------------+----------------
 223        0.00000 |       0.00000 | -6378168.00000
 224 (1 row)
 225
 226 --
 227 -- Test getting the latitude of a location.
 228 --
 229 SELECT latitude(ll_to_earth(0,0))::numeric(20,10);
 230    latitude
 231 --------------
 232  0.0000000000
 233 (1 row)
 234
 235 SELECT latitude(ll_to_earth(45,0))::numeric(20,10);
 236    latitude
 237 ---------------
 238  45.0000000000
 239 (1 row)
 240
 241 SELECT latitude(ll_to_earth(90,0))::numeric(20,10);
 242    latitude
 243 ---------------
 244  90.0000000000
 245 (1 row)
 246
 247 SELECT latitude(ll_to_earth(-45,0))::numeric(20,10);
 248     latitude
 249 ----------------
 250  -45.0000000000
 251 (1 row)
 252
 253 SELECT latitude(ll_to_earth(-90,0))::numeric(20,10);
 254     latitude
 255 ----------------
 256  -90.0000000000
 257 (1 row)
 258
 259 SELECT latitude(ll_to_earth(0,90))::numeric(20,10);
 260    latitude
 261 --------------
 262  0.0000000000
 263 (1 row)
 264
 265 SELECT latitude(ll_to_earth(45,90))::numeric(20,10);
 266    latitude
 267 ---------------
 268  45.0000000000
 269 (1 row)
 270
 271 SELECT latitude(ll_to_earth(90,90))::numeric(20,10);
 272    latitude
 273 ---------------
 274  90.0000000000
 275 (1 row)
 276
 277 SELECT latitude(ll_to_earth(-45,90))::numeric(20,10);
 278     latitude
 279 ----------------
 280  -45.0000000000
 281 (1 row)
 282
 283 SELECT latitude(ll_to_earth(-90,90))::numeric(20,10);
 284     latitude
 285 ----------------
 286  -90.0000000000
 287 (1 row)
 288
 289 SELECT latitude(ll_to_earth(0,180))::numeric(20,10);
 290    latitude
 291 --------------
 292  0.0000000000
 293 (1 row)
 294
 295 SELECT latitude(ll_to_earth(45,180))::numeric(20,10);
 296    latitude
 297 ---------------
 298  45.0000000000
 299 (1 row)
 300
 301 SELECT latitude(ll_to_earth(90,180))::numeric(20,10);
 302    latitude
 303 ---------------
 304  90.0000000000
 305 (1 row)
 306
 307 SELECT latitude(ll_to_earth(-45,180))::numeric(20,10);
 308     latitude
 309 ----------------
 310  -45.0000000000
 311 (1 row)
 312
 313 SELECT latitude(ll_to_earth(-90,180))::numeric(20,10);
 314     latitude
 315 ----------------
 316  -90.0000000000
 317 (1 row)
 318
 319 SELECT latitude(ll_to_earth(0,-90))::numeric(20,10);
 320    latitude
 321 --------------
 322  0.0000000000
 323 (1 row)
 324
 325 SELECT latitude(ll_to_earth(45,-90))::numeric(20,10);
 326    latitude
 327 ---------------
 328  45.0000000000
 329 (1 row)
 330
 331 SELECT latitude(ll_to_earth(90,-90))::numeric(20,10);
 332    latitude
 333 ---------------
 334  90.0000000000
 335 (1 row)
 336
 337 SELECT latitude(ll_to_earth(-45,-90))::numeric(20,10);
 338     latitude
 339 ----------------
 340  -45.0000000000
 341 (1 row)
 342
 343 SELECT latitude(ll_to_earth(-90,-90))::numeric(20,10);
 344     latitude
 345 ----------------
 346  -90.0000000000
 347 (1 row)
 348
 349 --
 350 -- Test getting the longitude of a location.
 351 --
 352 SELECT longitude(ll_to_earth(0,0))::numeric(20,10);
 353   longitude
 354 --------------
 355  0.0000000000
 356 (1 row)
 357
 358 SELECT longitude(ll_to_earth(45,0))::numeric(20,10);
 359   longitude
 360 --------------
 361  0.0000000000
 362 (1 row)
 363
 364 SELECT longitude(ll_to_earth(90,0))::numeric(20,10);
 365   longitude
 366 --------------
 367  0.0000000000
 368 (1 row)
 369
 370 SELECT longitude(ll_to_earth(-45,0))::numeric(20,10);
 371   longitude
 372 --------------
 373  0.0000000000
 374 (1 row)
 375
 376 SELECT longitude(ll_to_earth(-90,0))::numeric(20,10);
 377   longitude
 378 --------------
 379  0.0000000000
 380 (1 row)
 381
 382 SELECT longitude(ll_to_earth(0,90))::numeric(20,10);
 383    longitude
 384 ---------------
 385  90.0000000000
 386 (1 row)
 387
 388 SELECT longitude(ll_to_earth(45,90))::numeric(20,10);
 389    longitude
 390 ---------------
 391  90.0000000000
 392 (1 row)
 393
 394 SELECT longitude(ll_to_earth(90,90))::numeric(20,10);
 395    longitude
 396 ---------------
 397  90.0000000000
 398 (1 row)
 399
 400 SELECT longitude(ll_to_earth(-45,90))::numeric(20,10);
 401    longitude
 402 ---------------
 403  90.0000000000
 404 (1 row)
 405
 406 SELECT longitude(ll_to_earth(-90,90))::numeric(20,10);
 407    longitude
 408 ---------------
 409  90.0000000000
 410 (1 row)
 411
 412 SELECT longitude(ll_to_earth(0,180))::numeric(20,10);
 413    longitude
 414 ----------------
 415  180.0000000000
 416 (1 row)
 417
 418 SELECT longitude(ll_to_earth(45,180))::numeric(20,10);
 419    longitude
 420 ----------------
 421  180.0000000000
 422 (1 row)
 423
 424 SELECT longitude(ll_to_earth(90,180))::numeric(20,10);
 425    longitude
 426 ----------------
 427  180.0000000000
 428 (1 row)
 429
 430 SELECT longitude(ll_to_earth(-45,180))::numeric(20,10);
 431    longitude
 432 ----------------
 433  180.0000000000
 434 (1 row)
 435
 436 SELECT longitude(ll_to_earth(-90,180))::numeric(20,10);
 437    longitude
 438 ----------------
 439  180.0000000000
 440 (1 row)
 441
 442 SELECT longitude(ll_to_earth(0,-90))::numeric(20,10);
 443    longitude
 444 ----------------
 445  -90.0000000000
 446 (1 row)
 447
 448 SELECT longitude(ll_to_earth(45,-90))::numeric(20,10);
 449    longitude
 450 ----------------
 451  -90.0000000000
 452 (1 row)
 453
 454 SELECT longitude(ll_to_earth(90,-90))::numeric(20,10);
 455    longitude
 456 ----------------
 457  -90.0000000000
 458 (1 row)
 459
 460 SELECT longitude(ll_to_earth(-45,-90))::numeric(20,10);
 461    longitude
 462 ----------------
 463  -90.0000000000
 464 (1 row)
 465
 466 SELECT longitude(ll_to_earth(-90,-90))::numeric(20,10);
 467    longitude
 468 ----------------
 469  -90.0000000000
 470 (1 row)
 471
 472 --
 473 -- For the distance tests the following is some real life data.
 474 --
 475 -- Chicago has a latitude of 41.8 and a longitude of 87.6.
 476 -- Albuquerque has a latitude of 35.1 and a longitude of 106.7.
 477 -- (Note that latitude and longitude are specified differently
 478 -- in the cube based functions than for the point based functions.)
 479 --
 480 --
 481 -- Test getting the distance between two points using earth_distance.
 482 --
 483 SELECT earth_distance(ll_to_earth(0,0),ll_to_earth(0,0))::numeric(20,5);
 484  earth_distance
 485 ----------------
 486         0.00000
 487 (1 row)
 488
 489 SELECT earth_distance(ll_to_earth(0,0),ll_to_earth(0,180))::numeric(20,5);
 490  earth_distance
 491 ----------------
 492  20037605.73216
 493 (1 row)
 494
 495 SELECT earth_distance(ll_to_earth(0,0),ll_to_earth(90,0))::numeric(20,5);
 496  earth_distance
 497 ----------------
 498  10018802.86608
 499 (1 row)
 500
 501 SELECT earth_distance(ll_to_earth(0,0),ll_to_earth(0,90))::numeric(20,5);
 502  earth_distance
 503 ----------------
 504  10018802.86608
 505 (1 row)
 506
 507 SELECT earth_distance(ll_to_earth(0,0),ll_to_earth(0,1))::numeric(20,5);
 508  earth_distance
 509 ----------------
 510    111320.03185
 511 (1 row)
 512
 513 SELECT earth_distance(ll_to_earth(0,0),ll_to_earth(1,0))::numeric(20,5);
 514  earth_distance
 515 ----------------
 516    111320.03185
 517 (1 row)
 518
 519 SELECT earth_distance(ll_to_earth(30,0),ll_to_earth(30,1))::numeric(20,5);
 520  earth_distance
 521 ----------------
 522     96405.66962
 523 (1 row)
 524
 525 SELECT earth_distance(ll_to_earth(30,0),ll_to_earth(31,0))::numeric(20,5);
 526  earth_distance
 527 ----------------
 528    111320.03185
 529 (1 row)
 530
 531 SELECT earth_distance(ll_to_earth(60,0),ll_to_earth(60,1))::numeric(20,5);
 532  earth_distance
 533 ----------------
 534     55659.48608
 535 (1 row)
 536
 537 SELECT earth_distance(ll_to_earth(60,0),ll_to_earth(61,0))::numeric(20,5);
 538  earth_distance
 539 ----------------
 540    111320.03185
 541 (1 row)
 542
 543 SELECT earth_distance(ll_to_earth(41.8,87.6),ll_to_earth(35.1,106.7))::numeric(20,5);
 544  earth_distance
 545 ----------------
 546   1819303.21265
 547 (1 row)
 548
 549 SELECT (earth_distance(ll_to_earth(41.8,87.6),ll_to_earth(35.1,106.7))*
 550       100./2.54/12./5280.)::numeric(20,5);
 551   numeric
 552 ------------
 553  1130.46261
 554 (1 row)
 555
 556 --
 557 -- Test getting the distance between two points using geo_distance.
 558 --
 559 SELECT geo_distance('(0,0)'::point,'(0,0)'::point)::numeric(20,5);
 560  geo_distance
 561 --------------
 562       0.00000
 563 (1 row)
 564
 565 SELECT geo_distance('(0,0)'::point,'(180,0)'::point)::numeric(20,5);
 566  geo_distance
 567 --------------
 568   12436.77274
 569 (1 row)
 570
 571 SELECT geo_distance('(0,0)'::point,'(0,90)'::point)::numeric(20,5);
 572  geo_distance
 573 --------------
 574    6218.38637
 575 (1 row)
 576
 577 SELECT geo_distance('(0,0)'::point,'(90,0)'::point)::numeric(20,5);
 578  geo_distance
 579 --------------
 580    6218.38637
 581 (1 row)
 582
 583 SELECT geo_distance('(0,0)'::point,'(1,0)'::point)::numeric(20,5);
 584  geo_distance
 585 --------------
 586      69.09318
 587 (1 row)
 588
 589 SELECT geo_distance('(0,0)'::point,'(0,1)'::point)::numeric(20,5);
 590  geo_distance
 591 --------------
 592      69.09318
 593 (1 row)
 594
 595 SELECT geo_distance('(0,30)'::point,'(1,30)'::point)::numeric(20,5);
 596  geo_distance
 597 --------------
 598      59.83626
 599 (1 row)
 600
 601 SELECT geo_distance('(0,30)'::point,'(0,31)'::point)::numeric(20,5);
 602  geo_distance
 603 --------------
 604      69.09318
 605 (1 row)
 606
 607 SELECT geo_distance('(0,60)'::point,'(1,60)'::point)::numeric(20,5);
 608  geo_distance
 609 --------------
 610      34.54626
 611 (1 row)
 612
 613 SELECT geo_distance('(0,60)'::point,'(0,61)'::point)::numeric(20,5);
 614  geo_distance
 615 --------------
 616      69.09318
 617 (1 row)
 618
 619 SELECT geo_distance('(87.6,41.8)'::point,'(106.7,35.1)'::point)::numeric(20,5);
 620  geo_distance
 621 --------------
 622    1129.18983
 623 (1 row)
 624
 625 SELECT (geo_distance('(87.6,41.8)'::point,'(106.7,35.1)'::point)*5280.*12.*2.54/100.)::numeric(20,5);
 626     numeric
 627 ---------------
 628  1817254.87730
 629 (1 row)
 630
 631 --
 632 -- Test getting the distance between two points using the <@> operator.
 633 --
 634 SELECT ('(0,0)'::point <@> '(0,0)'::point)::numeric(20,5);
 635  numeric
 636 ---------
 637  0.00000
 638 (1 row)
 639
 640 SELECT ('(0,0)'::point <@> '(180,0)'::point)::numeric(20,5);
 641    numeric
 642 -------------
 643  12436.77274
 644 (1 row)
 645
 646 SELECT ('(0,0)'::point <@> '(0,90)'::point)::numeric(20,5);
 647   numeric
 648 ------------
 649  6218.38637
 650 (1 row)
 651
 652 SELECT ('(0,0)'::point <@> '(90,0)'::point)::numeric(20,5);
 653   numeric
 654 ------------
 655  6218.38637
 656 (1 row)
 657
 658 SELECT ('(0,0)'::point <@> '(1,0)'::point)::numeric(20,5);
 659  numeric
 660 ----------
 661  69.09318
 662 (1 row)
 663
 664 SELECT ('(0,0)'::point <@> '(0,1)'::point)::numeric(20,5);
 665  numeric
 666 ----------
 667  69.09318
 668 (1 row)
 669
 670 SELECT ('(0,30)'::point <@> '(1,30)'::point)::numeric(20,5);
 671  numeric
 672 ----------
 673  59.83626
 674 (1 row)
 675
 676 SELECT ('(0,30)'::point <@> '(0,31)'::point)::numeric(20,5);
 677  numeric
 678 ----------
 679  69.09318
 680 (1 row)
 681
 682 SELECT ('(0,60)'::point <@> '(1,60)'::point)::numeric(20,5);
 683  numeric
 684 ----------
 685  34.54626
 686 (1 row)
 687
 688 SELECT ('(0,60)'::point <@> '(0,61)'::point)::numeric(20,5);
 689  numeric
 690 ----------
 691  69.09318
 692 (1 row)
 693
 694 SELECT ('(87.6,41.8)'::point <@> '(106.7,35.1)'::point)::numeric(20,5);
 695   numeric
 696 ------------
 697  1129.18983
 698 (1 row)
 699
 700 SELECT (('(87.6,41.8)'::point <@> '(106.7,35.1)'::point)*5280.*12.*2.54/100.)::numeric(20,5);
 701     numeric
 702 ---------------
 703  1817254.87730
 704 (1 row)
 705
 706 --
 707 -- Test getting a bounding box around points.
 708 --
 709 SELECT cube_ll_coord(earth_box(ll_to_earth(0,0),112000),1)::numeric(20,5),
 710        cube_ll_coord(earth_box(ll_to_earth(0,0),112000),2)::numeric(20,5),
 711        cube_ll_coord(earth_box(ll_to_earth(0,0),112000),3)::numeric(20,5),
 712        cube_ur_coord(earth_box(ll_to_earth(0,0),112000),1)::numeric(20,5),
 713        cube_ur_coord(earth_box(ll_to_earth(0,0),112000),2)::numeric(20,5),
 714        cube_ur_coord(earth_box(ll_to_earth(0,0),112000),3)::numeric(20,5);
 715  cube_ll_coord | cube_ll_coord | cube_ll_coord | cube_ur_coord | cube_ur_coord | cube_ur_coord
 716 ---------------+---------------+---------------+---------------+---------------+---------------
 717  6266169.43896 | -111998.56104 | -111998.56104 | 6490166.56104 |  111998.56104 |  111998.56104
 718 (1 row)
 719
 720 SELECT cube_ll_coord(earth_box(ll_to_earth(0,0),pi()*earth()),1)::numeric(20,5),
 721        cube_ll_coord(earth_box(ll_to_earth(0,0),pi()*earth()),2)::numeric(20,5),
 722        cube_ll_coord(earth_box(ll_to_earth(0,0),pi()*earth()),3)::numeric(20,5),
 723        cube_ur_coord(earth_box(ll_to_earth(0,0),pi()*earth()),1)::numeric(20,5),
 724        cube_ur_coord(earth_box(ll_to_earth(0,0),pi()*earth()),2)::numeric(20,5),
 725        cube_ur_coord(earth_box(ll_to_earth(0,0),pi()*earth()),3)::numeric(20,5);
 726  cube_ll_coord  |  cube_ll_coord  |  cube_ll_coord  | cube_ur_coord  | cube_ur_coord  | cube_ur_coord
 727 ----------------+-----------------+-----------------+----------------+----------------+----------------
 728  -6378168.00000 | -12756336.00000 | -12756336.00000 | 19134504.00000 | 12756336.00000 | 12756336.00000
 729 (1 row)
 730
 731 SELECT cube_ll_coord(earth_box(ll_to_earth(0,0),10*earth()),1)::numeric(20,5),
 732        cube_ll_coord(earth_box(ll_to_earth(0,0),10*earth()),2)::numeric(20,5),
 733        cube_ll_coord(earth_box(ll_to_earth(0,0),10*earth()),3)::numeric(20,5),
 734        cube_ur_coord(earth_box(ll_to_earth(0,0),10*earth()),1)::numeric(20,5),
 735        cube_ur_coord(earth_box(ll_to_earth(0,0),10*earth()),2)::numeric(20,5),
 736        cube_ur_coord(earth_box(ll_to_earth(0,0),10*earth()),3)::numeric(20,5);
 737  cube_ll_coord  |  cube_ll_coord  |  cube_ll_coord  | cube_ur_coord  | cube_ur_coord  | cube_ur_coord
 738 ----------------+-----------------+-----------------+----------------+----------------+----------------
 739  -6378168.00000 | -12756336.00000 | -12756336.00000 | 19134504.00000 | 12756336.00000 | 12756336.00000
 740 (1 row)
 741
 742 --
 743 -- Test for points that should be in bounding boxes.
 744 --
 745 SELECT earth_box(ll_to_earth(0,0),
 746        earth_distance(ll_to_earth(0,0),ll_to_earth(0,1))*1.00001) @>
 747        ll_to_earth(0,1);
 748  ?column?
 749 ----------
 750  t
 751 (1 row)
 752
 753 SELECT earth_box(ll_to_earth(0,0),
 754        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.1))*1.00001) @>
 755        ll_to_earth(0,0.1);
 756  ?column?
 757 ----------
 758  t
 759 (1 row)
 760
 761 SELECT earth_box(ll_to_earth(0,0),
 762        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.01))*1.00001) @>
 763        ll_to_earth(0,0.01);
 764  ?column?
 765 ----------
 766  t
 767 (1 row)
 768
 769 SELECT earth_box(ll_to_earth(0,0),
 770        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.001))*1.00001) @>
 771        ll_to_earth(0,0.001);
 772  ?column?
 773 ----------
 774  t
 775 (1 row)
 776
 777 SELECT earth_box(ll_to_earth(0,0),
 778        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.0001))*1.00001) @>
 779        ll_to_earth(0,0.0001);
 780  ?column?
 781 ----------
 782  t
 783 (1 row)
 784
 785 SELECT earth_box(ll_to_earth(0,0),
 786        earth_distance(ll_to_earth(0,0),ll_to_earth(0.0001,0.0001))*1.00001) @>
 787        ll_to_earth(0.0001,0.0001);
 788  ?column?
 789 ----------
 790  t
 791 (1 row)
 792
 793 SELECT earth_box(ll_to_earth(45,45),
 794        earth_distance(ll_to_earth(45,45),ll_to_earth(45.0001,45.0001))*1.00001) @>
 795        ll_to_earth(45.0001,45.0001);
 796  ?column?
 797 ----------
 798  t
 799 (1 row)
 800
 801 SELECT earth_box(ll_to_earth(90,180),
 802        earth_distance(ll_to_earth(90,180),ll_to_earth(90.0001,180.0001))*1.00001) @>
 803        ll_to_earth(90.0001,180.0001);
 804  ?column?
 805 ----------
 806  t
 807 (1 row)
 808
 809 --
 810 -- Test for points that shouldn't be in bounding boxes. Note that we need
 811 -- to make points way outside, since some points close may be in the box
 812 -- but further away than the distance we are testing.
 813 --
 814 SELECT earth_box(ll_to_earth(0,0),
 815        earth_distance(ll_to_earth(0,0),ll_to_earth(0,1))*.57735) @>
 816        ll_to_earth(0,1);
 817  ?column?
 818 ----------
 819  f
 820 (1 row)
 821
 822 SELECT earth_box(ll_to_earth(0,0),
 823        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.1))*.57735) @>
 824        ll_to_earth(0,0.1);
 825  ?column?
 826 ----------
 827  f
 828 (1 row)
 829
 830 SELECT earth_box(ll_to_earth(0,0),
 831        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.01))*.57735) @>
 832        ll_to_earth(0,0.01);
 833  ?column?
 834 ----------
 835  f
 836 (1 row)
 837
 838 SELECT earth_box(ll_to_earth(0,0),
 839        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.001))*.57735) @>
 840        ll_to_earth(0,0.001);
 841  ?column?
 842 ----------
 843  f
 844 (1 row)
 845
 846 SELECT earth_box(ll_to_earth(0,0),
 847        earth_distance(ll_to_earth(0,0),ll_to_earth(0,0.0001))*.57735) @>
 848        ll_to_earth(0,0.0001);
 849  ?column?
 850 ----------
 851  f
 852 (1 row)
 853
 854 SELECT earth_box(ll_to_earth(0,0),
 855        earth_distance(ll_to_earth(0,0),ll_to_earth(0.0001,0.0001))*.57735) @>
 856        ll_to_earth(0.0001,0.0001);
 857  ?column?
 858 ----------
 859  f
 860 (1 row)
 861
 862 SELECT earth_box(ll_to_earth(45,45),
 863        earth_distance(ll_to_earth(45,45),ll_to_earth(45.0001,45.0001))*.57735) @>
 864        ll_to_earth(45.0001,45.0001);
 865  ?column?
 866 ----------
 867  f
 868 (1 row)
 869
 870 SELECT earth_box(ll_to_earth(90,180),
 871        earth_distance(ll_to_earth(90,180),ll_to_earth(90.0001,180.0001))*.57735) @>
 872        ll_to_earth(90.0001,180.0001);
 873  ?column?
 874 ----------
 875  f
 876 (1 row)
 877
 878 --
 879 -- Test the recommended constraints.
 880 --
 881 SELECT is_point(ll_to_earth(0,0));
 882 ERROR:  function is_point(earth) does not exist
 883 LINE 1: SELECT is_point(ll_to_earth(0,0));
 884                ^
 885 HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
 886 SELECT cube_dim(ll_to_earth(0,0)) <= 3;
 887  ?column?
 888 ----------
 889  t
 890 (1 row)
 891
 892 SELECT abs(cube_distance(ll_to_earth(0,0), '(0)'::cube) / earth() - 1) <
 893        '10e-12'::float8;
 894  ?column?
 895 ----------
 896  t
 897 (1 row)
 898
 899 SELECT is_point(ll_to_earth(30,60));
 900 ERROR:  function is_point(earth) does not exist
 901 LINE 1: SELECT is_point(ll_to_earth(30,60));
 902                ^
 903 HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
 904 SELECT cube_dim(ll_to_earth(30,60)) <= 3;
 905  ?column?
 906 ----------
 907  t
 908 (1 row)
 909
 910 SELECT abs(cube_distance(ll_to_earth(30,60), '(0)'::cube) / earth() - 1) <
 911        '10e-12'::float8;
 912  ?column?
 913 ----------
 914  t
 915 (1 row)
 916
 917 SELECT is_point(ll_to_earth(60,90));
 918 ERROR:  function is_point(earth) does not exist
 919 LINE 1: SELECT is_point(ll_to_earth(60,90));
 920                ^
 921 HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
 922 SELECT cube_dim(ll_to_earth(60,90)) <= 3;
 923  ?column?
 924 ----------
 925  t
 926 (1 row)
 927
 928 SELECT abs(cube_distance(ll_to_earth(60,90), '(0)'::cube) / earth() - 1) <
 929        '10e-12'::float8;
 930  ?column?
 931 ----------
 932  t
 933 (1 row)
 934
 935 SELECT is_point(ll_to_earth(-30,-90));
 936 ERROR:  function is_point(earth) does not exist
 937 LINE 1: SELECT is_point(ll_to_earth(-30,-90));
 938                ^
 939 HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
 940 SELECT cube_dim(ll_to_earth(-30,-90)) <= 3;
 941  ?column?
 942 ----------
 943  t
 944 (1 row)
 945
 946 SELECT abs(cube_distance(ll_to_earth(-30,-90), '(0)'::cube) / earth() - 1) <
 947        '10e-12'::float8;
 948  ?column?
 949 ----------
 950  t
 951 (1 row)
 952