OK list, here is one more philosophical GIS question, with more than one answer possible.<div>Task: create one layer with street names, house numbers and city name</div><div>Set up:</div><div>- vector layers are all shapefiles</div>
<div>- one layer with points, attribute filed is "house_number" (presumably, someone will input this data point by point, by digitizing or GPS surveying...)</div><div>- the second layer is "street buffer" - created from "streets", contains "street_name" field, among others</div>
<div><br></div><div>If I perform "intersect" I get a new layer with point within polygons, that now contain both "street_name" and "house_number" among others, _except_ for the houses/points that are on the crossroads, because they intersect with more than one polygon (street buffer). Let's say there are maybe 5%, or less, points that are on the crossroads, so my solution for now is to add them manually later after intersection.</div>
<div><br></div><div>Of course, the resulting shapefile contains more fields (columns), that are not needed, but cannot be erased in QGIS, or at least I haven't figured out how...</div><div><br></div><div>The second approach that crosses my mind is to import all of the shapefiles into PostGIS, but I do not have very much experience with "spatial queries" per se - is it possible to refine this process in PostGIS?</div>
<div><br></div><div>If anyone has any ideas/pointers/suggestions/experience with similar tasks I would be very grateful :-).</div><div><br></div><div><br></div><div>BR,</div><div><br clear="all">--- <br>Zoran Jankovic<br>
ZIS-Izrada softvera i savjetovanje / ZIS - Software Development and Consulting<br><br><a href="http://www.zisis.hr">http://www.zisis.hr</a><br><br>M: 00 385 98 682 902<br>T:  00 385 44 683 374<br><br>A. Senoe 4<br>44320 Kutina<br>
Croatia<br>
</div>