标签：GIS

“我算的面积怎么差了十倍？”——Shapely 的投影陷阱你是不是也遇到过这种情况：用 Shapely 算出来的多边形面积，跟 ArcGIS 或 QGIS 里显示的结果完全对不上？甚至差出几个数量级？别慌，这不是代码写错了，而是你忘了最关键一步——投影转换。我在参与某省耕地保护项目时，就曾因为没转投影，导致计算的地块面…

别再手动查经纬度了！用Geopy让Python自动搞定地理编码你是不是还在Excel里一列列复制粘贴地址，然后打开百度地图一个个查经纬度？或者写爬虫时被反爬机制卡住，项目进度一拖再拖？我在参与某省智慧交通项目时，就曾因手动处理5万条POI数据差点错过交付节点——直到我遇见了Geopy。地理编码不是魔法，是“地址⇄坐标”…

“Fiona打不开GDB？图层名都不知道怎么填！”——你的报错我太懂了上周一个在规划院实习的研究生小张私信我：“Dr. Gis，我用Fiona.open()打开本地.gdb文件，报错说‘Layer not found’，但我压根不知道里面有哪些图层名啊！总不能一个个猜吧？”——这简直是GIS自动化脚本新手的“经典噩梦三…

你写的Reproject为啥总报错？可能是坐标系理解错了上周一位读者在后台留言：‘Dr. Gis，我用Rasterio做reproject，结果图像全黑，控制台还跳出CRS mismatch警告……’——这太典型了。我在参与全国耕地遥感监测项目时，也栽过这个坑：明明代码照着文档敲，输出却像被泼了墨。问题根源往往不在语法…

你不是一个人：GDAL安装报错的“地狱级”体验“pip install gdal” —— 看似简单的一行命令，却让无数GIS新手在深夜对着满屏红色报错欲哭无泪。别慌，Dr. Gis当年第一次部署遥感处理环境时，也在cmd窗口前枯坐三小时，最后靠一杯美式续命才搞定。今天这篇，就是为你写的‘GDAL安装生存指南’。”“GD…

为什么你的空间连接总是报错？90%的人没搞懂sjoin的“潜规则”上周一个研究生私信我：“Dr. Gis，我用GeoPandas做空间连接，结果返回空DataFrame，代码也没报错，是不是库坏了？”——这太典型了！不是库坏了，是你没搞懂sjoin函数底层在“想什么”。我在国土空间规划项目中，曾因忽略坐标系一致性，导致…

NDVI一片黑？别慌，可能是你波段顺序搞反了上周一个研究生私信我：‘老师，我用Python算出来的NDVI全黑，拉伸也没用，是不是代码写错了？’——结果一查，红光和近红外波段位置颠倒了。这问题太常见了！今天我就手把手带你从原理到实战，彻底搞懂NDVI计算中那些‘数组操作’的坑。NDVI本质：一场植被与光的‘恋爱指数’N…

为什么你用Rasterio读影像总报错？可能是没搞懂这三层结构“明明路径没错，代码也没拼写错误，怎么一运行就报‘IndexError: band index out of range’？”——这是我后台收到最多的提问之一。别慌，这不是你的错，而是你还没理解遥感影像在Python眼中的“三维本质”。我在参与某省耕地遥感监…

你是不是也遇到过：地图打开全是“锟斤拷”？别慌，这不是你的电脑中病毒了，也不是数据坏了——这是典型的编码不匹配。我在参与某省国土空间规划项目时，第一次用GeoPandas读县级行政区划.shp文件，屏幕上蹦出一串“涓浗鍖哄煙”，差点以为硬盘烧了。后来才明白，这其实是UTF-8和GBK在打架。乱码的本质，是你用错了解码…

为什么你加的坐标系总报错？——自定义投影不是填空游戏上周一位在国土局工作的读者私信我：‘Dr. Gis，我在QGIS里加载地方测绘院给的Shapefile，属性表坐标动辄上亿，地图全跑到太平洋去了！’——这不是数据坏了，而是缺了最关键的一环：自定义投影。很多人以为‘选个EPSG代码就完事’，结果一遇到地方独立坐标系、工…