मेरे पास जीपीएस डेटा है जो मैंने बस मार्गों से लिया है, और अब मैं इसे अपने सड़क नेटवर्क पर स्नैप करना चाहूंगा। PostGIS DB में दोनों लेयर लाइन लेयर हैं। मैं या तो QGIS या PostGIS का उपयोग करना चाहता हूं, लेकिन अगर मुझे GRASS या ArcMap के साथ जाना है, तो यह ठीक है, भी। धन्यवाद!
स्पष्ट करने के लिए, मैं लाइनों को इंगित करने के लिए लाइनों को इंगित करने की कोशिश कर रहा हूं, लाइनों की ओर नहीं।
जवाबों:
मैं इस समारोह मैं इस्तेमाल किया है। सावधान रहें कि यह मौजूदा बिंदु तालिका में ज्यामिति को बदलता है। मैंने लंबे समय तक इसका उपयोग नहीं किया लेकिन ऐसा लगता है कि इसे काम करना चाहिए। जहाँ तक मुझे याद है यह अच्छी तरह से काम करता है यदि आपके पास दोनों तालिका में स्थानिक सूचकांक हैं। आदेश में यह फोन करने के लिए
चयन करें Snap_point_to_line ('points_table', 'line_table', 500)। यह आपके प्रक्षेपण प्रणाली की इकाई होने के नाते 500, 500 की सहिष्णुता के साथ स्नैप करेगा। मैं लैम्बर्ट के साथ काम करता था।
CREATE OR REPLACE FUNCTION snap_point_to_line(points_table character varying, line_table character varying, tolerance double precision)
RETURNS boolean AS
$BODY$
DECLARE
srid integer;
i integer;
row record;
row_1 record;
closest_distance double precision;
query varchar;
snapped_point geometry;
BEGIN
--Get the srid of the points table
FOR row IN EXECUTE 'select getsrid(the_geom) as srid from '||points_table||' where gid = (select min(gid) from '||points_table||')' LOOP
END LOOP;
srid := row.srid;
-- Add a column in which it will store the closest nodes from the line
FOR row IN EXECUTE 'SELECT the_geom FROM '||points_table LOOP
query := 'SELECT ST_Transform(the_geom,'||srid||') as the_geom, ST_Distance(GeometryFromText('''||ST_AsText(row.the_geom)||''','||srid||'), ST_Transform(the_geom,'||srid||')) as distance FROM ' ||line_table||' ORDER BY ST_Distance(GeometryFromText('''||ST_AsText(row.the_geom)||''','||srid||'), ST_Transform(the_geom,'||srid||')) LIMIT 1';
RAISE NOTICE '%',query;
FOR row_1 IN EXECUTE query LOOP
closest_distance := row_1.distance;
--If below the distance threeshold, then snap the point
IF closest_distance < tolerance THEN
snapped_point := ST_line_interpolate_point(ST_LineMerge(row_1.the_geom),ST_line_locate_point(ST_LineMerge(row_1.the_geom), row.the_geom));
--UPDATE the_geometry
EXECUTE 'UPDATE '||points_table||' SET the_geom = GeometryFromText('''||ST_AsText(snapped_point)||''','||srid||') WHERE ST_AsText(the_geom) = '''||ST_AsText(row.the_geom)||'''';
END IF;
END LOOP;
END LOOP;
RETURN true;
END;
$BODY$
LANGUAGE 'plpgsql' VOLATILE STRICT
COST 100;
ALTER FUNCTION snap_point_to_line(character varying, character varying, double precision) OWNER TO yourowner;