标签：Shapely

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

你写的“点在面内”判断，可能悄悄漏掉了百万数据上周一位读者私信我：“Dr. Gis，我用Shapely判断10万个点是否在某个行政区内，结果跑了半小时还报错——明明代码逻辑没问题啊？” 我一看就笑了：你不是第一个踩坑的。这问题背后藏着计算几何的“温柔陷阱”，今天我就带你彻底拆解。“点在面内”不是直觉那么简单——先搞懂“…

缓冲区生成总报错？你可能连“缓冲”的本质都没搞懂上周一位研究生私信我：“Dr. Gis，我用Shapely给点做500米缓冲区，结果出来是个椭圆，代码没报错但结果完全不对！”——这其实是90%初学者踩的第一个大坑：坐标系没转换。缓冲区不是画个圈那么简单，它背后藏着空间参考的“潜规则”。缓冲区的本质：不是画圈，是“空间影…