Changes between Version 31 and Version 32 of PostGIS_Raster_SoC_Idea_2012/Distance_Analysis_Tools/test


Ignore:
Timestamp:
Aug 29, 2012, 1:32:29 PM (12 years ago)
Author:
qliu
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • PostGIS_Raster_SoC_Idea_2012/Distance_Analysis_Tools/test

    v31 v32  
    3737}}}
    3838
    39  0-a: Test NULL Raster. Should be true.
     39 * 0-a: Test NULL Raster. Should be true.
    4040{{{
    4141#!sql
    … …  
    4343}}}
    4444
    45  0-b: Test Empty Raster. Should be true.
     45 * 0-b: Test Empty Raster. Should be true.
    4646{{{
    4747#!sql
    … …  
    5151== Test 1 ==
    5252=== Small number of sources vs low resolution raster: (Source: 10 points, Raster: 10x10) ===
    53  1-a: Source geometry within the extent of the desired raster:[[BR]]
     53 * 1-a: Source geometry within the extent of the desired raster:[[BR]]
    5454{{{
    5555#!sql
    … …  
    5858CREATE TABLE test_ref_raster_1 AS SELECT 1 AS rid,rast FROM (SELECT ST_MakeEmptyRaster(10,10,0,0,1,1,0,0,4326) AS rast)foo;
    5959}}}
    60   1-a-i: Generate Euclidean Distance from a one-band reference raster
    61    1-a-i-1: Snap is True; Max distance NOT specified:
     60  * 1-a-i: Generate Euclidean Distance from a one-band reference raster
     61   * 1-a-i-1: Snap is True; Max distance NOT specified:
    6262{{{
    6363#!sql
    … …  
    6666    result:[[BR]][[BR]]
    6767    [[Image(test_eudist_1_a_i_1.gif)]][[BR]]
    68    1-a-i-2: Snap is True; Max distance is specified:
     68   * 1-a-i-2: Snap is True; Max distance is specified:
    6969{{{
    7070#!sql
    … …  
    7373    result:[[BR]][[BR]]
    7474    [[Image(test_eudist_1_a_i_2.gif)]][[BR]]
    75    1-a-i-3: Snap is False; Max distance NOT specified:
     75   * 1-a-i-3: Snap is False; Max distance NOT specified:
    7676{{{
    7777#!sql
    … …  
    8080    result:[[BR]][[BR]]
    8181    [[Image(test_eudist_1_a_i_3.gif)]][[BR]]
    82    1-a-i-4: Snap is False; Max distance is specified:
     82   * 1-a-i-4: Snap is False; Max distance is specified:
    8383{{{
    8484#!sql
    … …  
    8787    result:[[BR]][[BR]]
    8888    [[Image(test_eudist_1_a_i_4.gif)]][[BR]]
    89   1-a-ii: Generate Euclidean Distance from a set of raster specifications:
    90    1-a-ii-1: Snap is True; Max distance NOT specified:
     89  * 1-a-ii: Generate Euclidean Distance from a set of raster specifications:
     90   * 1-a-ii-1: Snap is True; Max distance NOT specified:
    9191{{{
    9292#!sql
    … …  
    9494}}}
    9595    result: same as test_eudist_1_a_i_1.gif[[BR]][[BR]]
    96    1-a-ii-2: Snap is True; Max distance is specified:
     96   * 1-a-ii-2: Snap is True; Max distance is specified:
    9797{{{
    9898#!sql
    … …  
    100100}}}
    101101    result: same as test_eudist_1_a_i_2.gif[[BR]][[BR]]
    102    1-a-ii-3: Snap is False; Max distance NOT specified:
     102   * 1-a-ii-3: Snap is False; Max distance NOT specified:
    103103{{{
    104104#!sql
    … …  
    106106}}}
    107107    result: same as test_eudist_1_a_i_3.gif[[BR]][[BR]]
    108    1-a-ii-4: Snap is False; Max distance is specified:
     108   * 1-a-ii-4: Snap is False; Max distance is specified:
    109109{{{
    110110#!sql
    … …  
    113113    result: same as test_eudist_1_a_i_4.gif[[BR]][[BR]]
    114114[[BR]]
    115  1-b: Source geometry completely outside of the extent of the desired raster:
     115 * 1-b: Source geometry completely outside of the extent of the desired raster:
    116116{{{
    117117#!sql
    … …  
    119119CREATE TABLE test_source_geometry_1_b AS SELECT generate_series(1,10) AS id,ST_RandomPoints(the_geom,10) AS the_geom FROM (SELECT ST_SetSRID(ST_Extent(rast::geometry),4326) AS the_geom FROM test_source_raster_1_b)foo;
    120120}}}
    121   1-b-i: Generate Euclidean Distance from a one-band reference raster:
    122    1-b-i-1: Snap is True; Max distance NOT specified:
     121  * 1-b-i: Generate Euclidean Distance from a one-band reference raster:
     122   * 1-b-i-1: Snap is True; Max distance NOT specified:
    123123{{{
    124124#!sql
    … …  
    127127    result:[[BR]][[BR]]
    128128    [[Image(test_eudist_1_b_i_1.gif)]][[BR]]
    129    1-b-i-2: Snap is True; Max distance is specified:
     129   * 1-b-i-2: Snap is True; Max distance is specified:
    130130{{{
    131131#!sql
    … …  
    134134    result:[[BR]][[BR]]
    135135    [[Image(test_eudist_1_b_i_2.gif)]][[BR]]
    136    1-b-i-3: Snap is False; Max distance NOT specified:
     136   * 1-b-i-3: Snap is False; Max distance NOT specified:
    137137{{{
    138138#!sql
    … …  
    141141    result:[[BR]][[BR]]
    142142    [[Image(test_eudist_1_b_i_3.gif)]][[BR]]
    143    1-b-i-4: Snap is False; Max distance is specified:
     143   * 1-b-i-4: Snap is False; Max distance is specified:
    144144{{{
    145145#!sql
    … …  
    148148    result:[[BR]][[BR]]
    149149    [[Image(test_eudist_1_b_i_4.gif)]][[BR]]
    150   1-b-ii: Generate Euclidean Distance from a set of raster specifications:
    151    1-b-ii-1: Snap is True; Max distance NOT specified:
     150  * 1-b-ii: Generate Euclidean Distance from a set of raster specifications:
     151   * 1-b-ii-1: Snap is True; Max distance NOT specified:
    152152{{{
    153153#!sql
    … …  
    155155}}}
    156156    result: same as test_eudist_1_b_i_1.gif[BR]][[BR]]
    157    1-b-ii-2: Snap is True; Max distance is specified:
     157   * 1-b-ii-2: Snap is True; Max distance is specified:
    158158{{{
    159159#!sql
    … …  
    161161}}}
    162162    result: same as test_eudist_1_b_i_2.gif[BR]][[BR]]
    163    1-b-ii-3: Snap is False; Max distance NOT specified:
     163   * 1-b-ii-3: Snap is False; Max distance NOT specified:
    164164{{{
    165165#!sql
    … …  
    167167}}}
    168168    result: same as test_eudist_1_b_i_3.gif[BR]][[BR]]
    169    1-b-ii-4: Snap is False; Max distance is specified:
     169   * 1-b-ii-4: Snap is False; Max distance is specified:
    170170{{{
    171171#!sql
    … …  
    174174    result: same as test_eudist_1_b_i_4.gif[BR]][[BR]]
    175175[[BR]]
    176   1-c: Source geometry partially outside of the extent of the desired raster:
     176  * 1-c: Source geometry partially outside of the extent of the desired raster:
    177177{{{
    178178#!sql
    … …  
    180180CREATE TABLE test_source_geometry_1_c AS SELECT generate_series(1,10) AS id,ST_RandomPoints(the_geom,10) AS the_geom FROM (SELECT ST_SetSRID(ST_Extent(rast::geometry),4326) AS the_geom FROM test_source_raster_1_c)foo;
    181181}}}
    182    1-c-i: Generate Euclidean Distance from a one-band reference raster:
     182   * 1-c-i: Generate Euclidean Distance from a one-band reference raster:
    183183{{{
    184184#!sql
    … …  
    187187    result:[[BR]][[BR]]
    188188    [[Image(test_eudist_1_c_i_1.gif)]][[BR]]
    189    1-c-i-2: Snap is True; Max distance is specified:
     189   * 1-c-i-2: Snap is True; Max distance is specified:
    190190{{{
    191191#!sql
    … …  
    194194    result:[[BR]][[BR]]
    195195    [[Image(test_eudist_1_c_i_2.gif)]][[BR]]
    196    1-c-i-3: Snap is False; Max distance NOT specified:
     196   * 1-c-i-3: Snap is False; Max distance NOT specified:
    197197{{{
    198198#!sql
    … …  
    201201    result:[[BR]][[BR]]
    202202    [[Image(test_eudist_1_c_i_3.gif)]][[BR]]
    203    1-c-i-4: Snap is False; Max distance is specified:
     203   * 1-c-i-4: Snap is False; Max distance is specified:
    204204{{{
    205205#!sql
    … …  
    209209    [[Image(test_eudist_1_c_i_4.gif)]][[BR]]
    210210
    211   1-c-ii: Generate Euclidean Distance from a set of raster specifications:
    212    1-c-ii-1: Snap is True; Max distance NOT specified:
     211  * 1-c-ii: Generate Euclidean Distance from a set of raster specifications:
     212   * 1-c-ii-1: Snap is True; Max distance NOT specified:
    213213{{{
    214214#!sql
    … …  
    216216}}}
    217217    result: same as test_eudist_1_c_i_1.gif[BR]][[BR]]
    218    1-c-ii-2: Snap is True; Max distance is specified:
     218   * 1-c-ii-2: Snap is True; Max distance is specified:
    219219{{{
    220220#!sql
    … …  
    222222}}}
    223223    result: same as test_eudist_1_c_i_2.gif[BR]][[BR]]
    224    1-c-ii-3: Snap is False; Max distance NOT specified:
     224   * 1-c-ii-3: Snap is False; Max distance NOT specified:
    225225{{{
    226226#!sql
    … …  
    228228}}}
    229229    result: same as test_eudist_1_c_i_3.gif[BR]][[BR]]
    230    1-c-ii-4: Snap is False; Max distance is specified:
     230   * 1-c-ii-4: Snap is False; Max distance is specified:
    231231{{{
    232232#!sql
    … …  
    248248CREATE INDEX test_source_geometry_2_gist_idx ON test_source_geometry_2 USING GIST (the_geom);
    249249}}}
    250  2-a: Generate Euclidean Distance from a one-band reference raster
    251   2-a-1: Snap is True; Max distance NOT specified:
     250 * 2-a: Generate Euclidean Distance from a one-band reference raster
     251  * 2-a-1: Snap is True; Max distance NOT specified:
    252252{{{
    253253#!sql
    … …  
    256256    result:[[BR]][[BR]]
    257257    [[Image(test_eudist_2_a.gif)]][[BR]]
    258   2-a-2: Snap is True; Max distance is specified:
     258  * 2-a-2: Snap is True; Max distance is specified:
    259259{{{
    260260#!sql
    … …  
    263263    result:[[BR]][[BR]]
    264264    [[Image(test_eudist_2_a_2.gif)]][[BR]]
    265   2-a-3: Snap is False; Max distance NOT specified:
     265  * 2-a-3: Snap is False; Max distance NOT specified:
    266266{{{
    267267#!sql
    … …  
    270270    result:[[BR]][[BR]]
    271271    [[Image(test_eudist_2_a_3.gif)]][[BR]]
    272   2-a-4: Snap is False; Max distance is specified:
     272  * 2-a-4: Snap is False; Max distance is specified:
    273273{{{
    274274#!sql
    … …  
    277277    result:[[BR]][[BR]]
    278278    [[Image(test_eudist_2_a_4.gif)]][[BR]]
    279  2-b: Generate Euclidean Distance from a set of raster specifications:
    280   2-b-1: Snap is True; Max distance NOT specified:
     279 * 2-b: Generate Euclidean Distance from a set of raster specifications:
     280  * 2-b-1: Snap is True; Max distance NOT specified:
    281281{{{
    282282#!sql
    … …  
    284284}}}
    285285    result: same as test_eudist_2_a.gif[[BR]]
    286   2-b-2: Snap is True; Max distance is specified:
     286  * 2-b-2: Snap is True; Max distance is specified:
    287287{{{
    288288#!sql
    … …  
    290290}}}
    291291    result: same as test_eudist_2_a_2.gif[[BR]]
    292   2-b-3: Snap is False; Max distance NOT specified:
     292  * 2-b-3: Snap is False; Max distance NOT specified:
    293293{{{
    294294#!sql
    … …  
    296296}}}
    297297    result: same as test_eudist_2_a_3.gif[[BR]]
    298   2-b-4: Snap is False; Max distance is specified:
     298  * 2-b-4: Snap is False; Max distance is specified:
    299299{{{
    300300#!sql
    … …  
    311311CREATE INDEX test_source_geometry_3_gist_idx ON test_source_geometry_3 USING GIST (the_geom);
    312312}}}
    313  3-a: Generate Euclidean Distance from a one-band reference raster:
    314   3-a-1: Snap is True; Max distance NOT specified:
     313 * 3-a: Generate Euclidean Distance from a one-band reference raster:
     314  * 3-a-1: Snap is True; Max distance NOT specified:
    315315{{{
    316316#!sql
    … …  
    321321    result zoom in:[[BR]][[BR]]
    322322    [[Image(test_eudist_3_a_1_zoomin.gif)]][[BR]]
    323   3-a-2: Snap is True; Max distance is specified:
     323  * 3-a-2: Snap is True; Max distance is specified:
    324324{{{
    325325#!sql
    326326CREATE TABLE test_eudist_3_a_2 AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_3','the_geom',True,12) AS rast FROM test_ref_raster_3);
    327327}}}
    328   3-a-3: Snap is False; Max distance NOT specified:
     328  * 3-a-3: Snap is False; Max distance NOT specified:
    329329{{{
    330330#!sql
    … …  
    335335    result zoom in:[[BR]][[BR]]
    336336    [[Image(test_eudist_3_a_3_zoomin.gif)]][[BR]]
    337   3-a-4: Snap is False; Max distance is specified:
     337  * 3-a-4: Snap is False; Max distance is specified:
    338338{{{
    339339#!sql
    340340CREATE TABLE test_eudist_3_a_4 AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_3','the_geom',False,12) AS rast FROM test_ref_raster_3);
    341341}}}
    342  3-b: Generate Euclidean Distance from a set of raster specifications:
    343   3-b-1: Snap is True; Max distance NOT specified:
     342 * 3-b: Generate Euclidean Distance from a set of raster specifications:
     343  * 3-b-1: Snap is True; Max distance NOT specified:
    344344{{{
    345345#!sql
    346346CREATE TABLE test_eudist_3_b_1 AS (SELECT 1 AS rid,ST_EuclideanDistance(300,300,0,0,1,1,0,0,4326,NULL,'32BF','public','test_source_geometry_3','the_geom',True) AS rast FROM test_ref_raster_3);
    347347}}}
    348   3-b-2: Snap is True; Max distance is specified:
     348  * 3-b-2: Snap is True; Max distance is specified:
    349349{{{
    350350#!sql
    … …  
    352352..........SELECT 1 AS rid,ST_EuclideanDistance(300,300,0,0,1,1,0,0,4326,NULL,'32BF','public','test_source_geometry_3','the_geom',True,12) AS rast FROM test_ref_raster_3);
    353353}}}
    354   3-b-3: Snap is False; Max distance NOT specified:
     354  * 3-b-3: Snap is False; Max distance NOT specified:
    355355{{{
    356356#!sql
    357357CREATE TABLE test_eudist_3_b_3 AS (SELECT 1 AS rid,ST_EuclideanDistance(300,300,0,0,1,1,0,0,4326,NULL,'32BF','public','test_source_geometry_3','the_geom',False) AS rast FROM test_ref_raster_3);
    358358}}}
    359   3-b-4: Snap is False; Max distance is specified:
     359  * 3-b-4: Snap is False; Max distance is specified:
    360360{{{
    361361#!sql
    … …  
    369369CREATE TABLE test_ref_raster_4 AS SELECT 1 AS rid,rast FROM (SELECT ST_MakeEmptyRaster(300,300,0,0,1,1,0,0,4326) AS rast)foo;
    370370}}}
    371  4-a: Source: 100 points; Raster: 300x300
     371 * 4-a: Source: 100 points; Raster: 300x300
    372372{{{
    373373#!sql
    … …  
    376376CREATE TABLE test_eudist_4_a AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_a','the_geom',True) AS rast FROM test_ref_raster_4);
    377377}}}
    378   * 122325 ms execution time.
     378  122325 ms execution time.
    379379    result:[[BR]][[BR]]
    380380    [[Image(test_eudist_4_a.gif)]][[BR]]
    381  4-b: Source: 1,000 points; Raster: 300x300
     381 * 4-b: Source: 1,000 points; Raster: 300x300
    382382{{{
    383383#!sql
    … …  
    386386CREATE TABLE test_eudist_4_b AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_b','the_geom',True) AS rast FROM test_ref_raster_4);
    387387}}}
    388   * 132139 ms execution time.
     388  132139 ms execution time.
    389389    result:[[BR]][[BR]]
    390390    [[Image(test_eudist_4_b.gif)]][[BR]][[BR]]
    391  4-c: Source: 5,000 points; Raster: 300x300
     391 * 4-c: Source: 5,000 points; Raster: 300x300
    392392{{{
    393393#!sql
    … …  
    396396CREATE TABLE test_eudist_4_c AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_c','the_geom',True) AS rast FROM test_ref_raster_4);
    397397}}}
    398   * 126652 ms execution time.
     398  126652 ms execution time.
    399399    result:[[BR]][[BR]]
    400400    [[Image(test_eudist_4_c.gif)]][[BR]][[BR]]
    401  4-d: Source: 10,000 points; Raster: 300x300
     401 * 4-d: Source: 10,000 points; Raster: 300x300
    402402{{{
    403403#!sql
    … …  
    406406CREATE TABLE test_eudist_4_d AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_d','the_geom',True) AS rast FROM test_ref_raster_4);
    407407}}}
    408   * 143764 ms execution time.
     408  143764 ms execution time.
    409409    result:[[BR]][[BR]]
    410410    [[Image(test_eudist_4_d.gif)]][[BR]][[BR]]
    411  4-e: Source: 20,000 points; Raster: 300x300
     411 * 4-e: Source: 20,000 points; Raster: 300x300
    412412{{{
    413413#!sql
    … …  
    416416CREATE TABLE test_eudist_4_e AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_e','the_geom',True) AS rast FROM test_ref_raster_4);
    417417}}}
    418   * 120016 ms execution time.
     418  120016 ms execution time.
    419419    result:[[BR]][[BR]]
    420420    [[Image(test_eudist_4_e.gif)]][[BR]][[BR]]
    421  4-f: Source: 40,000 points; Raster: 300x300
     421 * 4-f: Source: 40,000 points; Raster: 300x300
    422422{{{
    423423#!sql
    … …  
    426426CREATE TABLE test_eudist_4_f AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_f','the_geom',True) AS rast FROM test_ref_raster_4);
    427427}}}
    428   * 169386 ms execution time.
     428  169386 ms execution time.
    429429    result:[[BR]][[BR]]
    430430    [[Image(test_eudist_4_f.gif)]][[BR]][[BR]]
    431  4-g: Source: 80,000 points; Raster: 300x300
     431 * 4-g: Source: 80,000 points; Raster: 300x300
    432432{{{
    433433#!sql
    … …  
    436436CREATE TABLE test_eudist_4_g AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_g','the_geom',True) AS rast FROM test_ref_raster_4);
    437437}}}
    438   * 164385 ms execution time.
     438  164385 ms execution time.
    439439    result:[[BR]][[BR]]
    440440    [[Image(test_eudist_4_g.gif)]][[BR]][[BR]]
    441  4-h: Source: 160,000 points; Raster: 300x300
     441 * 4-h: Source: 160,000 points; Raster: 300x300
    442442{{{
    443443#!sql
    … …  
    446446CREATE TABLE test_eudist_4_h AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_h','the_geom',True) AS rast FROM test_ref_raster_4);
    447447}}}
    448   * 160800 ms execution time.
     448  160800 ms execution time.
    449449    result:[[BR]][[BR]]
    450450    [[Image(test_eudist_4_h.gif)]][[BR]][[BR]]
    451  4-i: Source: 320,000 points; Raster: 300x300
     451 * 4-i: Source: 320,000 points; Raster: 300x300
    452452{{{
    453453#!sql
    … …  
    456456CREATE TABLE test_eudist_4_i AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_i','the_geom',True) AS rast FROM test_ref_raster_4);
    457457}}}
    458   * 204460 ms execution time.
     458  204460 ms execution time.
    459459    result: too many source point, picture doesn't make too much sense.[[BR]][[BR]]
    460  4-j: Source: 640,000 points; Raster: 300x300
     460 * 4-j: Source: 640,000 points; Raster: 300x300
    461461{{{
    462462#!sql
    … …  
    465465CREATE TABLE test_eudist_4_j AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_4_j','the_geom',True) AS rast FROM test_ref_raster_4);
    466466}}}
    467   * 168992 ms execution time.
     467  168992 ms execution time.
    468468    result: too many source point, picture doesn't make too much sense.[[BR]][[BR]]
    469469[[BR]]
    … …  
    471471=== Scalability test ===
    472472Very large size source points vs very high resolution raster:
    473  5-a: Source: 10,000 points; Raster: 1000x1000
     473 * 5-a: Source: 10,000 points; Raster: 1000x1000
    474474{{{
    475475#!sql
    … …  
    479479CREATE TABLE test_eudist_5_a AS (SELECT 1 AS rid,ST_EuclideanDistance(rast,'32BF','public','test_source_geometry_5_a','the_geom',True) AS rast FROM test_ref_raster_5);
    480480}}}
    481   * 1576434 ms execution time.
     481  1576434 ms execution time.
    482482    result:[[BR]][[BR]]
    483483    [[Image(test_eudist_5_a.gif)]][[BR]][[BR]]