SHP转STL【原创】| ①用QGIS裁剪Shape file

用Python实现快速的SHP转STL三维建模，所需的库如下；脚本原型由Gemini提供：

import geopandas
import trimesh
from shapely.geometry import Polygon, MultiPolygon 
import numpy as np
import pandas as pd 
from matplotlib import pyplot as plt

因原SHP文件太大，需要先裁剪出区域的SHP文件再进行STL转换处理，所以下载了QGIS来进行处理，详细步骤及截图如下所示。

① QGIS是从官网下载的、适用于win64安装的【QGISQT6-OSGeo4W-3.40.3-1.msi】，直接双击、管理员权限安装即可。

②打开QGIS，在左侧导航栏打开原shp文件，加载图层；在上方工具栏选择【Layer】→【Create Layer】→【New Temporary Scratch Layer】创建新图层，创建前选【Geometry type】为【Polygon】，坐标系统【Project CRS】选择了与原shp文件一样的CRS，点击创建。

③临时图层创建好后，在上方工具栏选择【Edit】→【Add polygon feature】，然后点击鼠标左键开始在原shp文件图层上选取想要的区域，点击鼠标右键后完成选区，这时所选区为填色区域。

④在上方工具栏选择【Vector】→【Clip】，然后【Input layer】选原shp文件，【Output layer】选创建的临时图层，然后请【Run】，若成功则完成！若不成功，可参考我：在同一个选项卡里打开【Algorithm Settings】，在【Invalid Feature Filtering】中选择【Skip (ignore) Feature with Invalid Geometries】，再【Run】即完成裁剪。

【PS. 我的原shp文件的Geometry应该是没什么问题的，也不知道为什么直接clip会报Geometry not valid？然后选了【skip】裁剪后出来的文件geometry好像也没发现什么问题，有恰当的经纬度坐标，所以不知道啦属于傻瓜流程分享!】 

⑤【Clip】顺利完成后会多出一个图层【Clipped】，选中该图层后鼠标右键有【Export】→【Save Feature As..】选项，可导出为Shape file，我依然沿用原文件的CRS；选中图层后鼠标右键也有【Properties】可以查看裁剪出来的图层信息完不完整。