Otdfbt

Add elements inside Array conditionally in JavaScript

Examples where existence is harder than evaluation

Are there vaccine ingredients which may not be disclosed ("hidden", "trade secret", or similar)?

My Sixteen Friendly Students

Justification of physical currency in an interstellar civilization?

Why doesn't increasing the temperature of something like wood or paper set them on fire?

Whose birthyears are canonically established in the MCU?

Employee is self-centered and affects the team negatively

Why doesn't a particle exert force on itself?

Capturing the entire webpage with WebExecute's CaptureImage

99 coins into the sacks

How do I minimise waste on a flight?

Is it safe to keep the GPU on 100% utilization for a very long time?

What happens when the drag force exceeds the weight of an object falling into earth?

Why is the episode called "The Last of the Starks"?

Is it possible to do moon sighting in advance for 5 years with 100% accuracy?

Exactly which act of bravery are Luke and Han awarded a medal for?

GLM: Modelling proportional data - account for variation in total sample size

Is there an application which does HTTP PUT?

Why are thrust reversers not used down to taxi speeds?

Light Switch Neutrals: Bundle all together?

Cyclic queue using an array in C#

Expl3 and recent xparse on overleaf: No expl3 loader detected

When an electron around an atom drops to a lower state, is 100% of the energy converted to a photon?

Remove empty geometries from geodataframe

Which data providers support loading attribute tables without geometry?How to merge polygons using GeoDataFrame?GeoPandas GeoDataFrame plot statistics - how?Create GPS (gpx) routes from GeoPandas geodataframeGeopanda write GeoDataFrame into shapefile or spatialiteIntersect does not remove features, leaves them empty. PythonOverlay two linestring objects in geopandas, accounting for the attributesOverlay Union Geopandas improve performanceGeodataframe not plotting properly in geopandasMultiple shapefiles to one GeoPandas GeoDataFrame

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

6

Since geopandas can't perform overlay/intersection with lines, I found a workaround:

import geopandas as gpd
import fiona
from shapely.geometry import shape, MultiPolygon, LineString, MultiLineString


lines = gpd.read_file('lines.shp')
polyons= gpd.read_file('polygons.shp')
pol_f = fiona.open('polygons.shp')

pol_shp = MultiPolygon([shape(x["geometry"]) for x in pol_f])

lines2 = lines.geometry.apply(lambda x: x.intersection(pol_shp))
lines2.to_file('lines_int.shp')

But it does not let me because there are the empty geometries are GeometryCollections.

I tried this to filter, but it won't work.

lines3 = lines2[lines2.geometry != np.nan]
lines3 = lines2[lines2.geometry == LineString]
lines3 = lines2[lines2.geometry == LineString()]
lines3 = lines2[lines2.geometry == MultiLineString()]
lines3 = lines2[lines2.geometry == MultiLineString]
lines3 = lines2[lines2.geometry != GeometryCollection]

Any ideas?

python geopandas

share|improve this question

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

add a comment | 

6

Since geopandas can't perform overlay/intersection with lines, I found a workaround:

import geopandas as gpd
import fiona
from shapely.geometry import shape, MultiPolygon, LineString, MultiLineString


lines = gpd.read_file('lines.shp')
polyons= gpd.read_file('polygons.shp')
pol_f = fiona.open('polygons.shp')

pol_shp = MultiPolygon([shape(x["geometry"]) for x in pol_f])

lines2 = lines.geometry.apply(lambda x: x.intersection(pol_shp))
lines2.to_file('lines_int.shp')

But it does not let me because there are the empty geometries are GeometryCollections.

I tried this to filter, but it won't work.

lines3 = lines2[lines2.geometry != np.nan]
lines3 = lines2[lines2.geometry == LineString]
lines3 = lines2[lines2.geometry == LineString()]
lines3 = lines2[lines2.geometry == MultiLineString()]
lines3 = lines2[lines2.geometry == MultiLineString]
lines3 = lines2[lines2.geometry != GeometryCollection]

Any ideas?

python geopandas

share|improve this question

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

add a comment | 

Since geopandas can't perform overlay/intersection with lines, I found a workaround:

import geopandas as gpd
import fiona
from shapely.geometry import shape, MultiPolygon, LineString, MultiLineString


lines = gpd.read_file('lines.shp')
polyons= gpd.read_file('polygons.shp')
pol_f = fiona.open('polygons.shp')

pol_shp = MultiPolygon([shape(x["geometry"]) for x in pol_f])

lines2 = lines.geometry.apply(lambda x: x.intersection(pol_shp))
lines2.to_file('lines_int.shp')

But it does not let me because there are the empty geometries are GeometryCollections.

I tried this to filter, but it won't work.

lines3 = lines2[lines2.geometry != np.nan]
lines3 = lines2[lines2.geometry == LineString]
lines3 = lines2[lines2.geometry == LineString()]
lines3 = lines2[lines2.geometry == MultiLineString()]
lines3 = lines2[lines2.geometry == MultiLineString]
lines3 = lines2[lines2.geometry != GeometryCollection]

Any ideas?

python geopandas

share|improve this question

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

import geopandas as gpd
import fiona
from shapely.geometry import shape, MultiPolygon, LineString, MultiLineString


lines = gpd.read_file('lines.shp')
polyons= gpd.read_file('polygons.shp')
pol_f = fiona.open('polygons.shp')

pol_shp = MultiPolygon([shape(x["geometry"]) for x in pol_f])

lines2 = lines.geometry.apply(lambda x: x.intersection(pol_shp))
lines2.to_file('lines_int.shp')

lines3 = lines2[lines2.geometry != np.nan]
lines3 = lines2[lines2.geometry == LineString]
lines3 = lines2[lines2.geometry == LineString()]
lines3 = lines2[lines2.geometry == MultiLineString()]
lines3 = lines2[lines2.geometry == MultiLineString]
lines3 = lines2[lines2.geometry != GeometryCollection]

share|improve this question

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

share|improve this question

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

asked Apr 29 at 10:13

ImanolUrImanolUr

396212

1 Answer
1

active

oldest

votes

5

The GeoDataFrame and GeoSeries have an is_empty attribute, similar to shapely geometry objects. You can use this to filter out the empty ones:

In [7]: lines = geopandas.GeoSeries([
 ...: LineString([(0, 0), (1, 1)]), LineString([]), LineString([(1, 1), (2, 2)])])

In [8]: lines 
Out[8]: 
0 LINESTRING (0 0, 1 1)
1 GEOMETRYCOLLECTION EMPTY
2 LINESTRING (1 1, 2 2)
dtype: object

In [9]: lines.is_empty 
Out[9]: 
0 False
1 True
2 False
dtype: bool

In [10]: lines[~lines.is_empty] 
Out[10]: 
0 LINESTRING (0 0, 1 1)
2 LINESTRING (1 1, 2 2)
dtype: object

share|improve this answer

answered Apr 29 at 10:25

jorisjoris

1,205713

add a comment | 

Your Answer

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);

draft saved

draft discarded

Sign up or log in

StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);

Sign up using Google

Sign up using Facebook

Sign up using Email and Password

Post as a guest

Name

Email

Required, but never shown

StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f321115%2fremove-empty-geometries-from-geodataframe%23new-answer', 'question_page');

);

Post as a guest

Name

Email

Required, but never shown

5

The GeoDataFrame and GeoSeries have an is_empty attribute, similar to shapely geometry objects. You can use this to filter out the empty ones:

In [7]: lines = geopandas.GeoSeries([
 ...: LineString([(0, 0), (1, 1)]), LineString([]), LineString([(1, 1), (2, 2)])])

In [8]: lines 
Out[8]: 
0 LINESTRING (0 0, 1 1)
1 GEOMETRYCOLLECTION EMPTY
2 LINESTRING (1 1, 2 2)
dtype: object

In [9]: lines.is_empty 
Out[9]: 
0 False
1 True
2 False
dtype: bool

In [10]: lines[~lines.is_empty] 
Out[10]: 
0 LINESTRING (0 0, 1 1)
2 LINESTRING (1 1, 2 2)
dtype: object

share|improve this answer

answered Apr 29 at 10:25

jorisjoris

1,205713

add a comment | 

5

The GeoDataFrame and GeoSeries have an is_empty attribute, similar to shapely geometry objects. You can use this to filter out the empty ones:

In [7]: lines = geopandas.GeoSeries([
 ...: LineString([(0, 0), (1, 1)]), LineString([]), LineString([(1, 1), (2, 2)])])

In [8]: lines 
Out[8]: 
0 LINESTRING (0 0, 1 1)
1 GEOMETRYCOLLECTION EMPTY
2 LINESTRING (1 1, 2 2)
dtype: object

In [9]: lines.is_empty 
Out[9]: 
0 False
1 True
2 False
dtype: bool

In [10]: lines[~lines.is_empty] 
Out[10]: 
0 LINESTRING (0 0, 1 1)
2 LINESTRING (1 1, 2 2)
dtype: object

share|improve this answer

answered Apr 29 at 10:25

jorisjoris

1,205713

add a comment | 

The GeoDataFrame and GeoSeries have an is_empty attribute, similar to shapely geometry objects. You can use this to filter out the empty ones:

In [7]: lines = geopandas.GeoSeries([
 ...: LineString([(0, 0), (1, 1)]), LineString([]), LineString([(1, 1), (2, 2)])])

In [8]: lines 
Out[8]: 
0 LINESTRING (0 0, 1 1)
1 GEOMETRYCOLLECTION EMPTY
2 LINESTRING (1 1, 2 2)
dtype: object

In [9]: lines.is_empty 
Out[9]: 
0 False
1 True
2 False
dtype: bool

In [10]: lines[~lines.is_empty] 
Out[10]: 
0 LINESTRING (0 0, 1 1)
2 LINESTRING (1 1, 2 2)
dtype: object

share|improve this answer

answered Apr 29 at 10:25

jorisjoris

1,205713

In [7]: lines = geopandas.GeoSeries([
 ...: LineString([(0, 0), (1, 1)]), LineString([]), LineString([(1, 1), (2, 2)])])

In [8]: lines 
Out[8]: 
0 LINESTRING (0 0, 1 1)
1 GEOMETRYCOLLECTION EMPTY
2 LINESTRING (1 1, 2 2)
dtype: object

In [9]: lines.is_empty 
Out[9]: 
0 False
1 True
2 False
dtype: bool

In [10]: lines[~lines.is_empty] 
Out[10]: 
0 LINESTRING (0 0, 1 1)
2 LINESTRING (1 1, 2 2)
dtype: object

share|improve this answer

answered Apr 29 at 10:25

jorisjoris

1,205713

answered Apr 29 at 10:25

jorisjoris

1,205713

answered Apr 29 at 10:25

jorisjoris

1,205713