Python与GIS-北京尧图网络

Python矢量与栅格数据处理原创孟祥帅祥帅的小屋2022-10-21 18:55发表于山东持续更新ing，欢迎各位点赞关注哦！python地理信息基础知识，基本图形的创建https://github.com/CNFeffery/DataScienceStudyNotesgeopandas建立缓冲区与地图处理、合并、交叉、缓冲https://blog.csdn.net/fengdu78/article/details/1070542491. Python处理栅格影像1.1 读取tifimport rasteriofrom rasterio.plot import showfrom matplotlib import colors, cmrs = rasterio.open(r'C:\Users\lenovo\Desktop\pzh_map_dispose\sf1.tif','r')result1=rs.read()1.2 替换Nodata数据#满足条件，则替换，否则保持原样result=np.where(result1==result1.min(),np.nan,result1)result# np.unique(rss[0])out:array([[[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],...,[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan],[nan, nan, nan, ..., nan, nan, nan]]])1.2 显示TIFimport geopandas as gpdshp = gpd.read_file(r"C:\Users\lenovo\Desktop\pzh_map_dispose\pzh_city.shp")fig, ax = plt.subplots(figsize=(5,9))shp.plot(ax=ax,color='none')show(result, transform=rs_mask.transform,ax=ax, cmap='gist_earth')fig.colorbar(cm.ScalarMappable(norm=colors.Normalize(vmin=np.nanmin(result), vmax=np.nanmax(result)), cmap='gist_earth'), ax=ax,extend='both',fraction=0.05)2.从头开始的一个例子——完整版2.1创建shp面，并写入文件import osimport geopandasfrom shapely import geometryimport matplotlib.pyplot as pltx1,y1=30,30x2,y2=50,50# 对应shapely.geometry中的Polygon，用于表示面，下面我们创建一个由若干Polygon对象组成cq = geopandas.GeoSeries([geometry.Polygon([(x1,y1), (x2,y1), (x2,y2), (x1,y2)]),geometry.Polygon([(x2,y1),(55,40), (x2,y2)])],index=['1', '2'], # 构建一个索引字段crs='EPSG:4326', # 坐标系是：WGS 1984)cq.to_file(r'simple_poly.shp',driver='ESRI Shapefile',encoding='utf-8')cq2.2 使用geopandas读取gdf=geopandas.read_file(r'simple_poly.shp')gdfgdf.plot(column='index')2.3 面转为栅格——根据字段转栅格# 矢量数据转栅格:不需要模板def shpToRaster2(shp_path,raster_path,cellsize,field):'''矢量数据转栅格:不需要模板:param shp_path:需要转换的矢量数据:param raster_path:输出后栅格保存路径:param cellsize:栅格大小

Python与GIS

相关推荐

Oracle 10g 服务端/客户端安装

首字响应280ms+动态慢思考：M3.1-Flash-Preview编程模型实测

基于PCA9422与PIC18F46K40的低功耗电源管理方案设计与实践

更多资讯

C# StringBuilder 详解：字符串拼接别只会用 `+`，从原理到实战一次讲透

L3开源静态评测｜Facebook Memlab：浏览器内存泄漏工具源码尽调（证据可复现）

【老计带你懂概率统计】03：条件概率与贝叶斯定理，从线索反推真相的神器

Vue 3 页面性能排查与优化

不用复杂命令，Windows 快速跑通 OpenClaw 本地 AI

Kiro 本身是亚马逊云科技于 2025 年 7 月推出的 AI 原生集成开发环境（Agentic IDE）

今日更新

MATLAB缩进快捷键与批量格式化实战指南

用Poco实现Qt风格信号与槽的跨平台C++方案

ESXi裸机虚拟化安装指南：从硬件兼容到虚拟机部署

本周精选

金融AI落地指南：从调研报告到技术选型与工程实践

大模型上下文成本与注意力优化技术深度解析——为何O(n²)平方复杂度无法被常规优化消除｜TaoToken 统一 Key 通道实测

基于GMM的语音识别系统原理与MATLAB实现详解

本月热门

我发现了一个新思路：用 Remotion + Claude Code 像写代码一样自动化生成短视频

Windows下 Codex 中 Chrome 和 Computer Use 插件不可用问题排查及解决参考方式：TaoToken 统一 Key 配置与验证

2026 大模型集体涨价：用 Python 做企业 Token 成本测算与选型避坑（附配置）

读完文章，想聊聊您的网站？