Hi All I have a need to identify roads (in say streets_ny) neighboring the state of new jersey which I have in streets_nj. I used to do this with a function that looked for the new york streets near/touching the new jersey roads, then looped over these roads finding roads near them. I have an updated data set and this runs much slower than it did before and I was concerned about getting duplicates with this method so I'm looking for alternatives. I tried the following but it seems to actually be worse
SELECT * FROM streets_ny ns, streets_ny s WHERE (( st_dWithin(ns.the_geom, s.the_geom, 0.08) AND s.func_class = '1' ) OR ( st_dWithin(ns.the_geom, s.the_geom, 0.032) AND s.func_class = '2' ) OR ( st_dWithin(ns.the_geom, s.the_geom, 0.008) AND ( s.func_class = '3' OR s.func_class = '4') )) AND ns.gid in ( SELECT nns.gid from streets_nj s, streets_ny nns WHERE s.the_geom && nns.the_geom AND st_dWithin(s.the_geom, nns.the_geom, 0.001) ) Basically I want longer stretches of the major roads and smaller stretches of secondary/side roads in the neighboring states. Is there a faster way to accomplish this? I think it would be simpler if I had a polygon defining the state of nj but not sure it should be necessary to do this faster. thanks, charles _______________________________________________ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users