postgis修复几何，修复自相交

Posted 2022-02-02 JerFer

1、报错问题

一般进行空间查询分析的时候，数据库报错如下：

NOTICE:  Ring Self-intersection at or near point 98.56186929203426 25.024940000380752
NOTICE:  Ring Self-intersection at or near point 98.840470640108038 24.256657960096845
NOTICE:  Ring Self-intersection at or near point 101.0043257882046 24.190313156127331

表示有图形自相交了。
图形如图所示：

2、处理方式如下

-- 把自相交的数据放在一个临时表里面
CREATE TABLE temp1 as select * from zd where ST_IsValid(shape) = false;
-- 原表中删除自相交的数据
delete from  zd where ST_IsValid(shape) = false;
select st_geometryType(st_multi(ST_Buffer(shape, 0.0))),st_geometryType(shape) from temp1;
-- 处理自相交数据
update temp1 set shape =st_multi(ST_Buffer(shape, 0.0))
-- 写入
insert into zd select * from temp1