def st_polars():
gdf = pst.read_file("/Users/TU99MK0/Downloads/ellipse_r2000.shp")
gdf = gdf.with_columns(pst.geom("geometry").st.to_wkt().alias("wkt"))
gdf = gdf.select(pl.exclude(['geometry']))
print(gdf.head(1))
t1 = conn.sql("""
select *,st_geomfromtext(wkt) as geometry
from gdf
""")
conn.sql("""
select * from t1
""").show()
conn.sql("""
COPY t1
TO '1.geojson'
WITH(FORMAT GDAL, DRIVER 'GeoJSON', SRS 'EPSG:25833')""")
conn.sql("""
COPY (select color,geometry from t1)
TO '2.geojson'
WITH(FORMAT GDAL, DRIVER 'GeoJSON', SRS 'EPSG:25833')""")
df = df.with_columns(pl.struct(pl.all()).map_elements(get_pms_task, return_dtype=pl.String).alias("历史任务信息"))
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。