7

假设我的数据库中有许多几何图形。我想在 PostgreSQL 中创建一个函数，借此我可以传入一个纬度/经度，并让它返回给定半径的距离（线性和道路）最近的几何图形。

我是这方面的新手，所以任何建议都值得赞赏。

我正在运行以下版本：

  • PostgreSQL：9.2
  • 地理信息系统：2.0

这是架构：

-- ----------------------------
--  Table structure for "cities-usa"
-- ----------------------------
DROP TABLE IF EXISTS "cities-usa";
CREATE TABLE "cities-usa" (
"gid" int4 NOT NULL DEFAULT nextval('"cities-usa_gid_seq"'::regclass),
"st_fips" varchar(4),
"sfips" varchar(2),
"county_fip" varchar(4),
"cfips" varchar(4),
"pl_fips" varchar(7),
"id" varchar(20),
"name" varchar(39),
"elevation" varchar(60),
"pop_1990" numeric,
"population" varchar(30),
"st" varchar(6),
"state" varchar(16),
"warngenlev" varchar(16),
"warngentyp" varchar(16),
"watch_warn" varchar(3),
"zwatch_war" float8,
"prog_disc" int4,
"zprog_disc" float8,
"comboflag" float8,
"land_water" varchar(16),
"recnum" float8,
"lon" float8,
"lat" float8,
"geom" "geometry"
)
WITH (OIDS=FALSE);
ALTER TABLE "cities-usa" OWNER TO "postgres";

-- ----------------------------
--  Primary key structure for table "cities-usa"
-- ----------------------------
ALTER TABLE "cities-usa" ADD CONSTRAINT "cities-usa_pkey" PRIMARY KEY ("gid") NOT DEFERRABLE INITIALLY IMMEDIATE;
4

1 回答 1

4

逻辑可以包装到一个函数中，但我建议使用查询来测试逻辑。您将遇到的困难是线性距离（英里）与角坐标（纬度和经度）相结合。此外，上面指定的几何类型没有 SRID 或几何类型。我猜应该geometry(Point,4326)是

您可以尝试一些策略，例如使用geography类型，它会自动使用计量器来执行ST_Distance, 或ST_DWithin. 下面的示例仅使用ST_Distance_Sphere从英里到米的转换来帮助您前进。或者，如果您需要性能，您可以尝试索引最近邻搜索<->的功能。

您可以尝试以下方法：

CREATE OR REPLACE FUNCTION WhatAmINear(lat float8, lon float8,
                                       radius_mi float8, num int DEFAULT 10)
    RETURNS SETOF "cities-usa" AS
$body$
SELECT *
FROM "cities-usa"
WHERE ST_Distance_Sphere(geom, ST_MakePoint(lon, lat)) <= radius_mi * 1609.34
LIMIT num;
$body$
LANGUAGE sql VOLATILE;

接着：

SELECT WhatAmINear(44.9, -93.1, 100);
于 2013-01-02T22:34:28.087 回答