Remove empty geometries from geodataframe
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ margin-bottom:0;
}
3
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
asked 58 mins ago
ImanolUrImanolUr
381212
add a comment |
3
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
asked 58 mins ago
ImanolUrImanolUr
381212
add a comment |
3
3
3
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
asked 58 mins ago
ImanolUrImanolUr
381212
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
python geopandas
asked 58 mins ago
ImanolUrImanolUr
381212
asked 58 mins ago
ImanolUrImanolUr
381212
asked 58 mins ago
ImanolUrImanolUr
381212
asked 58 mins ago
ImanolUrImanolUr
381212
asked 58 mins ago
ImanolUrImanolUr
381212
381212
add a comment |
add a comment |
1 Answer
1
active
oldest
votes
3
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
answered 45 mins ago
jorisjoris
1,175613
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
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
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
3
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
answered 45 mins ago
jorisjoris
1,175613
add a comment |
3
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
answered 45 mins ago
jorisjoris
1,175613
add a comment |
3
3
3
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
answered 45 mins ago
jorisjoris
1,175613
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
answered 45 mins ago
jorisjoris
1,175613
answered 45 mins ago
jorisjoris
1,175613
answered 45 mins ago
jorisjoris
1,175613
answered 45 mins ago
jorisjoris
1,175613
1,175613
add a comment |
add a comment |
draft saved
draft discarded
Thanks for contributing an answer to Geographic Information Systems Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
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
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
Required, but never shown
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
Required, but never shown
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
Required, but never shown
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
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
