postgis使用，通过（点线面）点，线，面以及范围查询范围内的点位

前置条件：

1、安装了postgresql，并且配置好了postgis插件，使用的数据库连接工具 dbeaver企业版

2、创建带有 geometry 类型的表 ，我这用的 4326 坐标系。

CREATE TABLE public.poll_enterprise_gis_geom (
    id bigserial NOT NULL,
    ent_id int8 NOT NULL, -- 企业表的id
    geom geometry(point, 4326) NULL, -- 企业经纬度
    CONSTRAINT poll_enterprise_gis_geom_pk PRIMARY KEY (id)
);
CREATE INDEX poll_enterprise_gis_geom_ent_id_idx ON public.poll_enterprise_gis_geom USING btree (ent_id);
COMMENT ON TABLE public.poll_enterprise_gis_geom IS '企业地图点位信息';

-- Column comments

COMMENT ON COLUMN public.poll_enterprise_gis_geom.ent_id IS '企业表的id';
COMMENT ON COLUMN public.poll_enterprise_gis_geom.geom IS '企业经纬度';

-- Permissions

ALTER TABLE public.poll_enterprise_gis_geom OWNER TO hmecology;
GRANT ALL ON TABLE public.poll_enterprise_gis_geom TO hmecology;

3、插入数据，请根据实际情况插入需要的数据。

1、根据单个点位和距离范围查询周围的点位

-- 找出距离坐标 (113.73833333333333, 22.856666666666666) 不超过 500 米的点位
SELECT pegg.*
FROM  poll_enterprise_gis_geom  pegg
WHERE ST_DWithin(pegg.geom, ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326)::geography, 500);
   

结果：

postgis-1.jpg

2、根据多个点位连成线查询距离范围查询周围的点位

--找到与指定线段距离不超过 100 米的点位，
SELECT pegg.*
FROM public.poll_enterprise_gis_geom pegg
WHERE ST_Distance(
ST_SetSRID(pegg.geom, 4326)::geography,
ST_MakeLine(
ARRAY[
ST_SetSRID(ST_MakePoint(113.68055555555556, 22.797777777777778), 4326),
ST_SetSRID(ST_MakePoint(113.671, 22.909961111111112), 4326),
ST_SetSRID(ST_MakePoint(113.87414722222222, 22.903891666666667), 4326)
]
)::geography
) <= 100;

结果：

postgis-2.jpg

3、根据多个点位连成面，获取面内的点位

根据三个点连成面，获取面内的点位，

注意：需要注意在 ST_AddPoint () 第二个参数方法加上
“ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326)”
其中的点位是第一个点位。

3.1先在画出图形：

SELECT ST_MakePolygon(
    ST_AddPoint(
        ST_MakeLine(
        ARRAY[
            ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326),
            ST_SetSRID(ST_MakePoint(113.71780555555556, 22.818538888888888), 4326),
            ST_SetSRID(ST_MakePoint(113.69770833333334, 22.851894444444444), 4326)
            ]
        ),
        ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326)
    )
) 

结果：

postgis-4.jpg

3.2通过连表获取图形内的点：

-- 根据三个点连成面，获取面内的点位，需要注意加上 
-- “ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326)”
-- 第一个点位
WITH poly AS (
    SELECT ST_MakePolygon(
    ST_AddPoint(
        ST_MakeLine(
        ARRAY[
            ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326),
            ST_SetSRID(ST_MakePoint(113.71780555555556, 22.818538888888888), 4326),
            ST_SetSRID(ST_MakePoint(113.69770833333334, 22.851894444444444), 4326)
            ]
        ),
        ST_SetSRID(ST_MakePoint(113.73833333333333, 22.856666666666666), 4326)
    )
    ) AS geom
)
SELECT pegg.*
FROM poll_enterprise_gis_geom pegg
WHERE ST_Within(ST_Transform(pegg.geom, 4326), (SELECT geom FROM poly))
      

结果：

postgis-3.jpg