#python #arcpy #geopandas #epsg
#python #arcpy #geopandas #epsg
Вопрос:
У меня есть этот код, и я запускаю его, и все хорошо, но когда я использую.head(), в столбцах не отображается raw, а остальные коды не выполняются, потому что raw-файлы пусты..
это код
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
Я что-то пропустил?
спасибо за год
Комментарии:
1. это важная часть моего машинного обучения, если вы не знаете его GEAI…
2. Теги имеют отношение к содержимому вопроса, а не к контексту; возможно, вы очень хорошо разрабатывали это для космического приложения, но это не сделало бы его подходящим для
space-engineering
тега…3. Если результат пустой (без строк), это означает, что ни одна из строк не пересекалась с
SLC
полем. Вы уверены, что они пересекаются? Например, вы можете сначала отобразить их, чтобы проверить это.