location RECORD;
relation_members TEXT[];
+ centroid GEOMETRY;
+
addr_street TEXT;
addr_place TEXT;
-- Speed up searches - just use the centroid of the feature
-- cheaper but less acurate
NEW.centroid := ST_PointOnSurface(NEW.geometry);
- -- For searching near features rather use the centroid
- NEW.postcode := null;
--DEBUG: RAISE WARNING 'Computing preliminary centroid at %',ST_AsText(NEW.centroid);
+ NEW.postcode := null;
+
-- recalculate country and partition
IF NEW.rank_search = 4 AND NEW.address is not NULL AND NEW.address ? 'country' THEN
-- for countries, believe the mapped country code,
and ( relation_members[i+1] != 'side_stream' or NEW.name->'name' = name->'name')
LOOP
UPDATE placex SET linked_place_id = NEW.place_id WHERE place_id = linked_node_id;
+ DELETE FROM search_name WHERE place_id = linked_node_id;
END LOOP;
END IF;
END LOOP;
IF location.place_id is not null THEN
--DEBUG: RAISE WARNING 'Linked %', location;
- -- Use this as the centre point of the geometry
- NEW.centroid := coalesce(location.centroid,
- ST_Centroid(location.geometry));
+ -- Use the linked point as the centre point of the geometry,
+ -- but only if it is within the area of the boundary.
+ centroid := coalesce(location.centroid, ST_Centroid(location.geometry));
+ IF centroid is not NULL AND ST_Within(centroid, NEW.geometry) THEN
+ NEW.centroid := centroid;
+ END IF;
-- Use the address rank of the linked place, if it has one
IF location.rank_address between 5 and 25 THEN
-- mark the linked place (excludes from search results)
UPDATE placex set linked_place_id = NEW.place_id
WHERE place_id = location.place_id;
+ -- ensure that those places are not found anymore
+ DELETE FROM search_name WHERE place_id = location.place_id;
SELECT wikipedia, importance
FROM compute_importance(location.extratags, NEW.country_code,
| ID | osm_type | osm_id |
| 0 | N | 2 |
+ # github #1352
+ Scenario: Do not use linked centroid when it is outside the area
+ Given the named places
+ | osm | class | type | admin | geometry |
+ | R13 | boundary | administrative | 4 | poly-area:0.01 |
+ And the named places
+ | osm | class | type | geometry |
+ | N2 | place | city | 0.1 0.1 |
+ And the relations
+ | id | members | tags+type |
+ | 13 | N2:label | boundary |
+ When importing
+ Then placex contains
+ | object | linked_place_id |
+ | N2 | R13 |
+ And placex contains
+ | object | centroid |
+ | R13 | in geometry |
| object | nameaddress_vector |
| W1 | bonn, new york, smalltown |
+ Scenario: a linked place does not show up in search name
+ Given the named places
+ | osm | class | type | admin | geometry |
+ | R13 | boundary | administrative | 9 | poly-area:0.01 |
+ And the named places
+ | osm | class | type | geometry |
+ | N2 | place | city | 0.1 0.1 |
+ And the relations
+ | id | members | tags+type |
+ | 13 | N2:label | boundary |
+ When importing
+ Then placex contains
+ | object | linked_place_id |
+ | N2 | R13 |
+ And search_name has no entry for N2
+
+ Scenario: a linked waterway does not show up in search name
+ Given the scene split-road
+ And the places
+ | osm | class | type | name | geometry |
+ | W1 | waterway | river | Rhein | :w-2 |
+ | W2 | waterway | river | Rhein | :w-3 |
+ | R13 | waterway | river | Rhein | :w-1 + :w-2 + :w-3 |
+ And the relations
+ | id | members | tags+type |
+ | 13 | W1,W2:main_stream | waterway |
+ When importing
+ Then placex contains
+ | object | linked_place_id |
+ | W1 | R13 |
+ | W2 | R13 |
+ And search_name has no entry for W1
+ And search_name has no entry for W2
return
if column.startswith('centroid'):
- fac = float(column[9:]) if column.startswith('centroid*') else 1.0
- x, y = value.split(' ')
- assert_almost_equal(float(x) * fac, row['cx'], "Bad x coordinate")
- assert_almost_equal(float(y) * fac, row['cy'], "Bad y coordinate")
+ if value == 'in geometry':
+ query = """SELECT ST_Within(ST_SetSRID(ST_Point({}, {}), 4326),
+ ST_SetSRID('{}'::geometry, 4326))""".format(
+ row['cx'], row['cy'], row['geomtxt'])
+ cur = context.db.cursor()
+ cur.execute(query)
+ eq_(cur.fetchone()[0], True, "(Row %s failed: %s)" % (column, query))
+ else:
+ fac = float(column[9:]) if column.startswith('centroid*') else 1.0
+ x, y = value.split(' ')
+ assert_almost_equal(float(x) * fac, row['cx'], msg="Bad x coordinate")
+ assert_almost_equal(float(y) * fac, row['cy'], msg="Bad y coordinate")
elif column == 'geometry':
geom = context.osm.parse_geometry(value, context.scene)
cur = context.db.cursor()