使用geopandas命令\u crs epsg“无列显示”时出错

2024-04-27 11:03:39 发布

您现在位置:Python中文网/ 问答频道 /正文

我有这个代码,我运行它,一切都很好,但当我使用.head()列中没有显示raw,其余代码没有运行,因为raw是emputy。。你知道吗

这是密码

import geopandas as gpd
import json
import numpy as np
from shapely.geometry import LineString, Point, box
import ast
from pyproj import Proj
paths = road_features.SHAPE.map(lambda x: np.array(ast.literal_eval(x)["paths"][0]))
pathLineStrings = paths.map(LineString)
gdf = gpd.GeoDataFrame(road_features,geometry=pathLineStrings)
#gdf.crs = {'init': 'epsg:3857'}
gdf.crs = {'init': 'epsg:4326'}
gdf = gdf.to_crs({'init': 'epsg:4326'})
x0 = -112.8
x1 = -111.0
y0 = 39.6
y1 = 41.6
SLC = box(x0,y0,x1,y1)

slc_df = gdf[gdf.intersects(SLC)]

我跑了

slc_df.head()

显示这个

segment_id
OBJECTID_1
Join_Count
TARGET_FID
Join_Count_1
TARGET_FID_1
Join_Count_12
TARGET_FID_12
OBJECTID
pre_dir
...
road_orient_approx
proximity_to_signal
proximity_to_billboard
proximity_to_nearest_intersection
proximity_to_major_road
population_density
station_id
SHAPE
accident_counts
geometry

0 rows × 32 columns

我错过什么了吗? 谢谢你的新年时光


Tags: to代码importtargetinitcountepsgpaths