Beskrivning

Returnerar en geometri som representerar punktuppsättningens skärningspunkt mellan två geometrier. Med andra ord, den del av geometri A och geometri B som delas mellan de två geometrierna.

Om geometrierna inte har några gemensamma punkter (dvs. är disjunkta) returneras en tom atomgeometri av lämplig typ.

Om den valfria parametern gridSize anges (GEOS-3.9.0 eller senare krävs) garanteras att alla resultatnoder hamnar på ett snäppavrundat rutnät med angiven storlek. Observera att åtgärder på ett rutnät kan innehålla små artefakter som uppkommer vid justering mot rutnätet, se ST_ReducePrecision.

ST_Intersection i kombination med ST_Intersects är användbart för att klippa geometrier, t.ex. i frågor om bounding box, buffert eller region där du bara behöver den del av en geometri som ligger inom ett land eller en region av intresse.

Rumsliga predikat som ST_Intersects och överlagringsfunktioner som ST_Intersection beräknas med olika algoritmer. För mycket nära eller nästan sammanfallande koordinater kan flyttalsprecision och geometrisk robusthet göra att ett predikat returnerar true trots att den beräknade skärningen är tom, eller placera en skärningspunkt en bit från indatalinjen. För reproducerbara överlagringsresultat bör indata först snäppas till ett lämpligt rutnät, till exempel med ST_ReducePrecision, eller parametern gridSize användas när alla indatanoder redan ligger på rutnätet.

För geografin är detta ett tunt hölje runt geometriimplementeringen. Den bestämmer först den bästa SRID som passar avgränsningsrutan för de 2 geografiobjekten (om geografiobjekt är inom en halv zon UTM men inte samma UTM väljer en av dem) (gynnar UTM eller Lambert Azimuthal Equal Area (LAEA) nord / sydpol och faller tillbaka på Mercator i värsta fall) och sedan intersektion i den bästa passande plana spatiala ref och retransformerar tillbaka till WGS84-geografi.

	Denna funktion kommer att släppa M-koordinatvärdena om de finns.

	Vid arbete med 3D-geometrier kan du använda den SFCGAL-baserade CG_3DIntersection, som utför en riktig 3D-skärning för 3D-geometrier. Även om den här funktionen fungerar med Z-koordinater beräknar den ett medelvärde av Z-koordinaten.

Utförs av GEOS-modulen

Förbättrad: 3.1.0 accepterar en gridSize-parameter

Kräver GEOS >= 3.9.0 för att använda parametern gridSize

Ändrat: 3.0.0 är inte beroende av SFCGAL.

Tillgänglighet: 1.5 stöd för datatypen Geography infördes.

Denna metod implementerar OGC:s implementeringsspecifikation för enkla funktioner för SQL 1.1. s2.1.1.3

Denna metod implementerar SQL/MM-specifikationen. SQL-MM 3: 5.1.18

Denna funktion stöder 3d och kommer inte att tappa z-index. Resultatet beräknas dock endast med hjälp av XY. De resulterande Z-värdena kopieras, medelvärdesbildas eller interpoleras.

Exempel

SELECT ST_Intersection('POINT(0 0)'::geometry, 'LINESTRING ( 2 0,0 2 )'::geometry);

PUNKT EMPTY

SELECT ST_Intersection('POINT(0 0)'::geometry, 'LINESTRING ( 0 0,0 2 )'::geometry);

PUNKT(0 0)

010100000000000000000000000000000000000000

Figur

Klipp ut alla linjer (spår) efter land. Här antar vi att landets geom är POLYGON eller MULTIPOLYGON. OBS: vi behåller bara korsningar som resulterar i en LINESTRING eller MULTILINESTRING eftersom vi inte bryr oss om spår som bara delar en punkt. Dumpen behövs för att expandera en geometrisamling till enskilda MULT*-delar. Nedanstående är ganska generisk och kommer att fungera för polys, etc. genom att bara ändra var klausulen.

select clipped.gid, clipped.f_name, clipped_geom
from (
         select trails.gid, trails.f_name,
             (ST_Dump(ST_Intersection(country.geom, trails.geom))).geom clipped_geom
         from country
              inner join trails on ST_Intersects(country.geom, trails.geom)
     ) as clipped
where ST_Dimension(clipped.clipped_geom) = 1;

För polygoner, till exempel polygonlandmärken, kan du även använda det ibland snabbare knepet att buffra allt med 0,0. Förutom polygoner resulterar det i en tom geometrisamling. En geometrisamling med polygoner, linjer och punkter som buffras med 0,0 lämnar därför bara polygonerna kvar och löser upp samlingens skal.

select poly.gid,
    ST_Multi(ST_Buffer(ST_Intersection(country.geom, poly.geom),
            0.0
        )
    ) clipped_geom
from country
     inner join poly on ST_Intersects(country.geom, poly.geom)
where not ST_IsEmpty(ST_Buffer(ST_Intersection(country.geom, poly.geom), 0.0));

Observera att detta inte är en verklig 3D-skärning. Den använder 2,5D-geometrier (geometrier med Z-koordinater, men där beräkningar utförs i 2D). Samma indata jämförs med CG_3DIntersection, som utför en verklig 3D-skärning.

WITH data AS (
  SELECT 'LINESTRING Z (2 2 6,1.5 1.5 7,1 1 8,0.5 0.5 8,0 0 10)'::geometry AS input_linestring,
         'POLYGON Z ((0 0 8,0 1 8,1 1 8,1 0 8,0 0 8))'::geometry AS input_polygon
)
SELECT input_linestring AS input_linestring,
       input_polygon AS input_polygon,
       ST_Intersection(input_linestring, input_polygon) AS intersection_2d,
       CG_3DIntersection(input_linestring, input_polygon) AS intersection_3d
FROM data;

LINESTRING Z (2 2 6,1.5 1.5 7,1 1 8,0.5 0.5 8,0 0 10) | POLYGON Z ((0 0 8,0 1 8,1 1 8,1 0 8,0 0 8)) | LINESTRING Z (1 1 8,0.5 0.5 8,0 0 10) | LINESTRING Z (1 1 8,0.5 0.5 8)

010200008005000000000000000000004000000000000000400000000000001840000000000000F83F000000000000F83F0000000000001C40000000000000F03F000000000000F03F0000000000002040000000000000E03F000000000000E03F0000000000002040000000000000000000000000000000000000000000002440 | 010300008001000000050000000000000000000000000000000000000000000000000020400000000000000000000000000000F03F0000000000002040000000000000F03F000000000000F03F0000000000002040000000000000F03F00000000000000000000000000002040000000000000000000000000000000000000000000002040 | 010200008003000000000000000000F03F000000000000F03F0000000000002040000000000000E03F000000000000E03F0000000000002040000000000000000000000000000000000000000000002440 | 010200008002000000000000000000F03F000000000000F03F0000000000002040000000000000E03F000000000000E03F0000000000002040

Figur

Se även

CG_3DIntersection, ST_Difference, ST_Union, ST_ClipByBox2D, ST_Dimension, ST_Dump, ST_Force2D, ST_SymDifference, ST_Intersects, ST_Multi, ST_ReducePrecision