M1 Mac下Django GEOS创建Polygon等几何时持续触发GEOS C错误

该问题由Apple M1系列芯片下Homebrew安装的高版本GEOS库与Django 3.2内置的GEOS C接口绑定存在ctypes参数类型不匹配引发，特征完全匹配描述的现象：

• 仅创建需要传入子几何/坐标环数组的复合几何类型（Polygon/MultiPoint/MultiLineString/MultiPolygon）时抛出该异常，Point/LineString这类不需要传数组参数的简单几何类型创建正常
• 通过WKT字符串调用GEOSGeometry创建对象时不触发异常，因为该路径走GEOS内置的WKT解析逻辑，不经过Python层ctypes构造几何对象的数组传参流程
• Intel芯片Mac无该问题，原因是M1架构的指针长度、内存对齐规则和x86架构存在差异，Django 3.2早期版本的ctypes参数声明没有适配该差异
• 切换Python版本、重装Homebrew依赖无效，因为问题出在Django框架层的GEOS绑定逻辑，并非Python环境或依赖包安装损坏

按以下步骤可100%确认问题：

1. 执行命令查看当前GEOS版本：

   geos-config --version
   

   若返回版本号大于等于3.10.0，基本可判定命中该兼容问题。
2. 进入Django shell执行以下测试代码，验证接口调用结果：

   from ctypes import c_void_p
   from django.contrib.gis.geos.libgeos import GEOSFuncFactory, geos_version_info
   from django.contrib.gis.geos.coordseq import CoordSeq
   
   # 打印当前加载的GEOS版本信息
   print(geos_version_info())
   
   # 测试线性环创建接口
   create_linear_ring = GEOSFuncFactory("GEOSGeom_createLinearRing_r")
   create_linear_ring.restype = c_void_p
   cs = CoordSeq([(4,4),(4,4),(4,4),(4,4)])
   ring_ptr = create_linear_ring(cs.ptr)
   print("线性环对象指针：", ring_ptr)
   
   # 测试多边形创建接口
   create_polygon = GEOSFuncFactory("GEOSGeom_createPolygon_r")
   create_polygon.restype = c_void_p
   poly_ptr = create_polygon(ring_ptr, None, 0)
   print("多边形对象指针：", poly_ptr)
   

   若最终多边形指针打印结果为None，即可确认是ctypes参数类型不匹配导致GEOS接口调用失败返回空指针，触发框架层异常。

按优先级从高到低选择即可：

• 升级Django版本（推荐）
  优先升级到Django 3.2最新补丁版本（3.2.15及以上），或直接升级到Django 4.0+正式版。上述版本官方已经合入M1架构下GEOS ctypes绑定的适配修复，升级后原有代码无需修改即可正常运行。
• 临时补丁（无法升级Django时使用）
  若项目暂时无法升级Django版本，在项目settings.py的最顶部（必须提前于所有django.contrib.gis相关模块的导入）添加如下补丁代码，手动修正ctypes函数的参数类型声明：

  import ctypes
  from django.contrib.gis.geos.libgeos import GEOSFuncFactory
  
  # 修正多边形创建接口参数类型
  geos_create_poly = GEOSFuncFactory("GEOSGeom_createPolygon_r")
  geos_create_poly.argtypes = [
      ctypes.c_void_p,
      ctypes.c_void_p,
      ctypes.POINTER(ctypes.c_void_p),
      ctypes.c_uint
  ]
  geos_create_poly.restype = ctypes.c_void_p
  
  # 修正所有多几何类型、几何集合创建接口参数类型
  for func_name in [
      "GEOSGeom_createMultiPoint_r",
      "GEOSGeom_createMultiLineString_r",
      "GEOSGeom_createMultiPolygon_r",
      "GEOSGeom_createGeometryCollection_r"
  ]:
      geos_create_multi = GEOSFuncFactory(func_name)
      geos_create_multi.argtypes = [
          ctypes.c_void_p,
          ctypes.c_uint,
          ctypes.POINTER(ctypes.c_void_p)
      ]
      geos_create_multi.restype = ctypes.c_void_p
  

• 降级GEOS版本（不推荐）
  若以上方案均无法实施，可通过Homebrew安装3.9.x版本的GEOS库替换当前高版本，该版本不存在与Django 3.2早期版本的绑定兼容问题，但会限制其他依赖高版本GEOS的空间分析功能使用。

内容的提问来源于stack exchange，提问作者arielkaluzhny