Kimi | QGIS制作分幅图指南

假设你是QGIS软件专家，我现在有一幅“G:\生态调查评价\Geodata\云南省LUCC_CAS_L2.tif”土地覆被栅格数据，现在我在QGIS中使用Polygonize (raster to vector)工具，把栅格转为了矢量，但是出现了一个问题，由于栅格数据是一个个像元，转矢量之后会出现一个个的锯齿，我现在想要生成平滑的矢量，如何改进操作

思考已完成

用户的问题是关于如何在 QGIS 中对栅格转矢量后的锯齿进行平滑处理。这个问题涉及到技术操作，同时也隐含了用户对结果的期望——生成更美观、更符合实际需求的矢量数据。我需要先明确几个关键点：

首先，我要分析用户的需求和潜在意图。用户已经完成了栅格转矢量的操作，但结果不符合预期，主要是因为锯齿状的边界问题。这说明用户希望得到更平滑的边界，可能是为了更好地展示数据，或者是为后续分析（例如缓冲区分析、叠加分析等）提供更良好的输入数据。用户的意图是明确的：减少锯齿，得到更平滑的矢量。

其次，我要考虑问题的技术背景。栅格转矢量的本质是将像元的边界转化为矢量多边形的边界，因此像元的方格状特性自然会导致锯齿状的边界。这是栅格数据的固有特性，无法完全避免。解决这个问题的关键是通过后续处理来平滑边界，而不是直接减少锯齿的生成。

接着，我需要思考有哪些方法可以实现平滑边界的目标。QGIS 中的平滑工具是一个直接的选择，它可以通过插值和调整多边形的顶点来减少锯齿。但我也意识到，单纯平滑可能会导致一些微妙的几何变形，因此我需要提醒用户注意平滑的程度。除此之外，我想到可以结合缓冲区操作：先给多边形添加一个小的缓冲区，然后再移除缓冲区，这种方法通常能有效减少锯齿，同时保持多边形的总体形状。

在思考如何组织回答时，我决定按步骤清晰地描述操作流程。首先介绍平滑工具的使用方法，其次补充缓冲区操作的技巧，最后提到数据预处理（例如重分类或栅格重采样）作为可能的辅助手段。这样可以给用户提供一个全面的解决方案，同时避免让回答显得过于复杂。

我还需要注意回答的语气和风格。用户的问题是技术性的，因此我的回答要尽量简洁明了，避免过多的术语。同时，我也要确保回答具有指导性，方便用户直接操作。

在 QGIS 中，如果你通过 Polygonize (raster to vector) 工具将栅格数据转为矢量后出现锯齿状边界，可以通过以下方法来平滑矢量边界：

使用平滑工具（Smooth）
- 在 QGIS 菜单栏中，选择 “矢量（Vector）”>“几何工具（Geometry Tools）”>“平滑（Smooth）”。
- 在弹出的 “平滑” 对话框中，选择你刚才通过 Polygonize 工具生成的锯齿状矢量图层作为输入图层（Input layer）。
- 设置迭代次数（Iterations）参数。迭代次数越多，平滑效果越明显，但可能会使边界偏离原始形状较多。你可以先尝试设置为 3 - 5 次，根据实际效果进行调整。
- 设置最大角度（Maximum angle）参数。这个参数用于控制平滑过程中允许的最大角度变化。较小的值会使边界更平滑，但可能会丢失一些重要的边界细节。你可以设置为 10 - 20 度左右，根据实际情况进行微调。
- 设置容差（Tolerance）参数。容差值决定了平滑过程中允许的边界偏移量。较大的容差值会使边界更加平滑，但可能会导致边界与原始位置偏离较大。你可以先设置为 1 - 2 个地图单位（如米或千米，具体取决于你的数据投影），然后根据效果进行调整。
- 运行平滑工具后，就会生成一个新的平滑后的矢量图层。
利用缓冲区（Buffer）操作进行平滑（间接方法）
- 首先，选择 “矢量（Vector）”>“几何工具（Geometry Tools）”>“缓冲区（Buffer）”。
- 在 “缓冲区” 对话框中，将锯齿状矢量图层作为输入图层。
- 设置缓冲区距离为较小的正值（如 1 - 5 个地图单位）。这个距离根据你的数据精度和需求来确定，目的是通过缓冲区操作来填充锯齿间的凹陷部分。
- 运行缓冲区工具，得到缓冲后的矢量图层。
- 然后再对缓冲后的矢量图层进行平滑操作（按照前面提到的平滑工具步骤），这样可以进一步改善边界形状。
数据预处理（在栅格转矢量之前）
- 重分类（Reclassify） ：在进行 Polygonize 操作之前，可以对栅格数据进行重分类。选择 “栅格（Raster）”>“转换（Conversion）”>“重分类值（Reclassify values）”。通过重分类，可以将相近的土地覆被类别合并，减少类别边界的复杂性，从而在转矢量后可能会减少锯齿的产生。例如，如果你的土地覆被类型有多个细分的亚类，而这些亚类在矢量化后边界锯齿较为明显，可以将它们合并为一个主类。
- 栅格重采样（Resample） ：使用 “栅格（Raster）”>“投影（Projections）”>“重采样（Resample）” 工具。通过改变栅格的像元大小，对栅格数据进行重采样。适当增大像元大小可以减少栅格的细节，使后续矢量化后的边界相对平滑。不过要注意，增大像元大小会降低数据的空间分辨率，可能会丢失一些细节信息。

QGIS有没有类似ArcGIS炸开的操作，我现在有一个面要素，有多个部分，我想类似ArcGIS那样把不同位置的部分给炸开，分成独立的

思考已完成

我先仔细琢磨了下用户的需求，就是想在 QGIS 里面把一个多部分的面要素给分开，就像 ArcGIS 里的炸开操作一样，搞成独立的部分。

我立马就想到了 QGIS 有几个现成的工具应该能搞定这事儿。“Multipart to Singleparts”这个工具名字就很直观，感觉就是用来干这个的，能直接把多部分的面要素变成一个个独立的。不过，我也虑到用户可能还想手动挑着分，于是又想到了“按字段分离”和“按表达式分离”，这两个方法可以满足更细粒度的分离需求。

然后我就是想着得把这些方法都给用户讲清楚，让他们能根据自己的情况去选。我就先从最简单的“Multipart to Singleparts”讲起，把操作步骤都列出来，这样用户一目了然。接着又把“按字段分离”和“按表达式分离”的方法加上，把它们的适用场景和操作流程都详细说一说。

我还琢磨着用户可能不太熟悉这些功能在哪找，所以就都把工具的找到路径给写上了，这样用户操作起来就不用瞎找了。最后，我觉得这样应该就能把问题给解决好啦，那就可以开始把这思路整理成回复了。

在 QGIS 中可以通过以下工具实现类似 ArcGIS 炸开的操作，将多部分面要素拆分为独立的单一面要素：

方法一：“Multipart to Singleparts”工具

操作步骤
- 在 QGIS 菜单栏中选择 “矢量（Vector）”>“几何工具（Geometry Tools）”>“Multipart to Singleparts”，会弹出相应的对话框。
- 选择你想要拆分的多部分面要素图层，点击 “确定（OK）”，工具会将原本属于同一多部分要素的各个面部分拆分成独立的面要素，生成新的图层。

方法二：“按字段分离（Split Vector Layer）”工具

操作步骤
- 选择 “矢量（Vector）”>“数据管理工具（Data Management Tools）”>“Split Vector Layer”。
- 在弹出的对话框中，选择要拆分的面要素图层。
- 选择一个字段作为拆分依据，这个字段可以是包含每个面部分唯一标识的字段。例如，如果你的面要素中有 “ID” 字段，每个小面部分都有唯一编码，就可以选择 “ID” 字段。
- 指定输出文件夹，运行工具后，每个面部分会根据所选字段的值生成独立的形状文件。

方法三：“按表达式分离（Split with Expressions）”工具

操作步骤
- 如果没有合适的字段，可以使用表达式来定义拆分规则。选择 “矢量（Vector）”>“数据管理工具（Data Management Tools）”>“Split with Expressions”。
- 在对话框中选择你的面要素图层。
- 输入一个表达式来确定如何将面要素拆分为不同的组。例如，可以根据面要素的某些属性（如名称、类型等）来构建表达式。比如，如果你想根据面要素的 “类型” 字段将不同类型的部分分开，可以使用表达式 “类型 = '某种类型'”，运行工具后，每个符合表达式的面要素部分会生成独立的图层。

现在有一个名为“田坝保护动植物”的点要素图层，里面有“Name”字段，为保护动植物的名称，有“Type”字段，写了是“保护动物”还是“保护植物”，我想把保护植物用绿色方块符号，保护动物用黄色圆点符号，每个不同名称的动植物叠加上不同的序号，请问如何设置符号

思考已完成

用户的问题是如何在 QGIS 中为一个点要素图层设置不同类型的符号，并且在符号上叠加序号标注。我先分析一下用户的需求：用户提到的图层有两个重要字段——“Name” 和 “Type”，并希望根据 “Type” 字段区分符号（保护植物为绿色方块，保护动物为黄色圆点），同时为每个不同名称的动植物叠加序号。这意味着用户的重点是符号的分类和标注的叠加。

我先考虑如何符号分类。这一步相对简单，可以通过符号分类器功能或规则符号化来实现。符号分类器可以自动根据字段值分组并赋予不同符号，而规则符号化则更灵活，允许用户自定义每组的符号样式。我倾向于推荐规则符号化，因为用户对符号的要求比较明确（绿色方块和黄色圆点），而且规则符号化提供了更大的控制空间。

接着是符号样式的设计。用户提到保护植物用绿色方块，保护动物用黄色圆点。这些符号样式在 QGIS 中可以通过调整符号的形状和颜色轻松实现。我需要在回答中详细说明如何设置这些符号的参数，比如填充颜色、边框颜色、大小等。

然后是序号的叠加。这一部分稍微复杂一些，因为序号需要基于 “Name” 字段的不同值来生成，并叠加在符号上。我想到可以利用 QGIS 的标签功能，通过表达式从 “Name” 字段中提取序号信息并显示。为了让标注的效果更清晰，可能还需要调整标注的大小、位置和字体样式。

在回答的逻辑上，我决定先介绍符号分类的方法，再详细说明如何设置符号样式，最后讲解如何叠加序号标注。这样可以按照从基础到复杂的顺序，逐步解决用户的需求。如果用户对 QGIS 的操作不太熟悉，这样的结构也能帮助他们更容易理解。

在语气和风格上，我选择用清晰、简洁的语言，避免过于技术化的表述，同时提供具体的步骤和参数设置，确保用户可以直接按照回答来操作。

最后，我在脑中快速验证了一下整个操作流程的可行性，确认这些方法在 QGIS 中是可行的，且能够满足用户的需求。

以下是在 QGIS 中对 “田坝保护动植物” 点要素图层进行符号设置并叠加序号的详细步骤：

符号分类设置

打开图层属性
- 在 QGIS 图层面板中右键点击 “田坝保护动植物” 图层，选择 “属性（Properties）” 选项。
选择符号分类字段
- 在图层属性对话框中，切换到 “样式（Style）” 选项卡。
- 点击 “分类（Categorized）” 下拉菜单，选择 “基于字段的分类（Categorized by Field）” 或 “基于表达式的分类（Categorized by Expression）”。
- 如果选择 “Categorized by Field”，在 “值字段（Value Field）” 下拉菜单中选择 “Type” 字段。
- 若选择 “Categorized by Expression”，则需要构建一个表达式来区分 “保护植物” 和 “保护动物”，例如：CASE WHEN "Type" = '保护植物' THEN '植物' ELSE '动物' END。
设置符号样式
- 点击 “分类” 按钮，自动为不同类型的要素生成符号。
- 对于 “保护植物” 类别，将其符号样式设置为绿色方块。点击该类别的符号框，选择 “填充颜色（Fill color）” 为绿色，“形状（Shape）” 选择方形。
- 对于 “保护动物” 类别，将符号样式设置为黄色圆点。点击其符号框，选择 “填充颜色” 为黄色，“形状” 选择圆形。
- 点击 “确定（OK）” 按钮，关闭图层属性对话框。

叠加序号设置

生成序号字段（如果尚未存在）
- 如果图层中没有序号字段，可以通过添加新字段来创建。右键点击 “田坝保护动植物” 图层，选择 “打开属性表（Open Attribute Table）”。
- 在属性表窗口中，点击 “添加字段（Add Field）” 按钮，命名为 “序号”，类型设置为整数（Integer）。
- 点击 “确定（OK）” 按钮关闭属性表。
填充序号值
- 在属性表中，选中所有要填充序号的要素记录。
- 点击 “计算字段（Field Calculator）” 按钮，选择 “序号” 字段，在表达式框中输入 $id +1（假设要从 1 开始编号，$id 是指要素的行号）。
- 点击 “确定（OK）” 按钮，完成序号填充。
添加标注以显示序号
- 在 QGIS 菜单栏中，选择 “标签（Labels）”>“设置标签（Set Labeling）”。
- 在弹出的标签设置对话框中，选择 “基于字段的标签（Label with Field）” 或 “基于表达式的标签（Label with Expression）”。
- 如果选择 “Label with Field”，在下拉菜单中选择 “序号” 字段。
- 若选择 “Label with Expression”，可以构建一个表达式，例如：concat('序号：', "序号")，来显示带有 “序号：” 前缀的序号。
- 调整标注的样式，如字体大小、颜色等，使其与符号相协调，然后点击 “确定（OK）”，关闭标签设置对话框。

通过以上步骤，就可以在 QGIS 中为 “田坝保护动植物” 点要素图层设置不同的符号并叠加序号标注了。

QGIS的按位置选择，能不能有类似ArcGIS的那种距离缓冲功能，比如Intersect，能不能把1m以内误差的都算作Intersect

搜索

QGIS 距离误差 Intersect 设置

53 个结果

思考已完成

首先，我分析了用户的需求，是要在QGIS中找到一种类似ArcGIS中按距离缓冲的Intersect功能。我想到在ArcGIS中可以通过设置容差来实现这种功能，可能在QGIS中也有类似的缓冲区工具或者插件来达到类似效果。

接着，我开始仔细查看搜索到的资料。我发现“网页2”和“网页12”提供了关于QGIS中按位置选择和缓冲区工具的相关信息，这可能对解决问题有帮助。同时，“网页4”和“网页6”也提到了QGIS中的缓冲区设置。我注意到这些资料中并没有直接提到类似ArcGIS容差的功能，但有提到缓冲区工具可以设置距离，这让我想到可能可以通过缓冲区工具来间接实现用户的需求。

我开始思考如何将这些信息整合起来。一个想法是先对目标图层创建缓冲区，再与另一图层进行Intersect分析。这样，缓冲区就能包含一定范围内的要素，从而实现类似ArcGIS中容差的效果。同时，我考虑到QGIS的插件功能，想有没有相关的插件可以直接实现类似的功能。

在这个过程中，我还考虑了用户对操作步骤的期望，希望提供一个尽量简单易行的解决方案。我尝试从用户的角度出发，将复杂的操作分解为更易理解的步骤，并且在思考过程中不断调整自己的思路，以确保给出的方案是可行的。

现在，我认为已经找到了一个合理的解决方案，准备以一种清晰、逻辑连贯的方式将这些步骤和思路组织起来，正式回复用户。

在 QGIS 中可以通过以下方式实现类似 ArcGIS 中按距离缓冲的 Intersect 功能：

创建缓冲区

使用“矢量”->“地学数据处理工具”->“缓冲区”工具，或安装“MMQGIS”插件来创建缓冲区。设置输入图层、缓冲区距离等参数，距离可设为 1 米。

使用“按位置选择”工具

选择“矢量”->“研究工具”->“按位置选择”，在打开的窗口中设置目标图层和选择图层，并选择“与选择图层相交”的选项。此操作可选择出目标图层中与选择图层在空间上相交的要素，但无法直接设置距离。

使用“按位置提取”工具

在工具箱中搜索并打开“按位置提取”工具，选择目标图层和参考图层，设置空间关系为“与参考图层相交”。此工具可提取目标图层中与参考图层相交的要素，并输出到新图层。

使用“Spatial Query”插件

启用或安装该插件后，选择目标图层和参考图层，设置空间操作为“intersect”，即可选择出目标图层中与参考图层相交的要素。

使用“Join Attributes by Location”工具

在工具箱中找到此工具，选择目标图层和参考图层，设置几何谓词为“intersects”，运行后可将参考图层的属性连接到目标图层，可选择按距离提取数据。

使用“按距离选择”工具

在 QGIS 的“选择工具栏”中，若已启用，可找到“按距离选择”工具，可基于选中要素进行数据选择。

QGIS如何制作分幅图，有一个长条状的矢量，我想分幅制图

思考已完成

用户提到有一个长条状的矢量并希望在 QGIS 中进行分幅制图，这让我意识到用户的需求是将一个长条状数据分成多个部分，并为每个部分生成独立的地图输出。分幅制图通常涉及到数据分割、符号化、布局设计以及批量导出，这些都是 QGIS 中常见的操作。

首先，我考虑用户提到的“长条状矢量”，这让我想到可能需要先对这个矢量进行分割。分割的目的是将长条状数据分成多个较短的部分，以便为每个部分生成独立的分幅。QGIS 中有多种分割方法，比如按距离分割（在 Processing 工具箱中找到），或者通过创建等距的缓冲区来实现。如果用户没有特定的分割要求，我可以建议使用等长的缓冲区来切割。

接着，我想到用户可能还需要符号化和标注每个分幅的边界或位置。这可以通过创建一个中心线来实现，中心线可以用于标注分幅的编号或名称，并提供视觉参考。我计划提到一种使用 “等距节点” 的方法来生成中心线，然后将中心线简化为多段线。

然后是布局设计的部分，我联想到 QGIS 的布局管理器是一个强大的工具，可以通过它创建多个地图视图并放置相关元素（比如标题、比例尺、图例等）。为了让用户更方便，我决定建议使用 “按参考图层的网格布局” 插件，这可以帮助快速生成多个布局视图，每个视图对应一个分幅。

最后，我考虑到用户可能需要批量导出每个分幅的图幅到文件。在 QGIS 中，可以利用布局的批量导出功能，我打算推荐使用插件或脚本实现这一目标。

在进行回复时，我决定以逐步的方式说明每个步骤，并且提供工具名称和参数设置的细节，以便用户可以直接操作。同时，我会确保逻辑连贯，从分割到布局再到导出都涵盖到。

在 QGIS 中可以通过以下步骤制作分幅图，以处理长条状矢量数据：

步骤 1：加载矢量数据

打开 QGIS 软件。
在菜单栏中选择“图层（Layer）”->“添加图层（Add Layer）”->“添加矢量图层（Add Vector Layer）”。
浏览并选择你的长条状矢量数据文件，点击“打开（Open）”按钮加载数据。

步骤 2：创建分幅网格

转到菜单栏中的“矢量（Vector）”->“研究工具（Research Tools）”->“按参考图层的网格布局（Create Grid Layout Based on Reference Layer）”。
在弹出的对话框中，选择你的长条状矢量图层作为参考图层。
设置网格的行数和列数，或者指定每个网格单元的大小，以确定分幅的数量和大小。
点击“确定（OK）”按钮生成分幅网格。

步骤 3：符号化分幅网格

在图层面板中右键点击生成的分幅网格图层，选择“属性（Properties）”。
在图层属性对话框中，切换到“样式（Style）”选项卡。
选择合适的符号类型（如线框、填充颜色等）来区分分幅边界。
点击“确定（OK）”按钮关闭对话框。

步骤 4：创建布局并添加地图

在菜单栏中选择“项目（Project）”->“新建打印布局（New Print Layout）”。
在弹出的对话框中，输入布局名称，点击“确定（OK）”按钮创建布局。
在布局设计窗口中，点击“添加地图（Add Map）”按钮，将地图添加到布局中。

步骤 5：设置地图范围和分幅

在布局中，选择添加的地图项。
在右侧的“项目属性（Item Properties）”面板中，设置地图的范围，使其与分幅网格的一个单元匹配。
你可以通过调整地图的比例尺和中心位置来精确设置地图范围。

步骤 6：添加图幅编号或其他元素

在布局中，你可以添加文本标签来显示图幅编号或其他信息。
点击“添加标签（Add Label）”按钮，在布局中放置一个文本标签。
在文本标签的属性中，输入图幅编号或相关信息。

步骤 7：批量导出分幅图

如果你有多个分幅，希望批量导出为单独的图像文件，可以使用“按图幅导出（Export by Atlas）”功能。
在布局设计窗口中，点击“启用地图册（Enable Atlas）”按钮。
设置地图册的覆盖图层为分幅网格图层，选择合适的输出格式和文件名前缀。
点击“导出地图册（Export Atlas）”按钮，开始批量导出分幅图。

步骤 8：调整和优化布局

在布局设计窗口中，你可以进一步调整地图的外观，如添加图例、标题、比例尺等元素。
调整各元素的位置和大小，使布局更加美观和专业。

通过这些步骤，你可以在 QGIS 中制作分幅图，将长条状矢量数据分成多个图幅进行打印或导出。可以根据实际需求调整分幅的大小、数量和布局样式。

我现在QGIS中加载了一个湖水范围矢量面：H:\Lashihai\水体提取202507132037360621.geojson 我现在想要沿着这个矢量面的外框，绘制50个正方形网格，正方形网格大小按照外框线长度自适应得到，请帮我写出QGIS Python脚本代码

下面是一段可直接在 QGIS Python 控制台运行的完整脚本。
它用你给出的 GeoJSON 矢量面（仅取第一条要素的外轮廓）为基准，
计算外框周长，按周长自动算出正方形网格边长，
并生成 50 个恰好覆盖该外框的最小外接矩形（MBR）的正方形网格，
输出为临时图层，结果立刻加到地图画布中。

Python

# -*- coding: utf-8 -*-
"""
在 QGIS Python 控制台中运行即可
生成 50 个正方形网格覆盖给定矢量面的外轮廓
"""

from qgis.core import (
    QgsProject, QgsVectorLayer, QgsFeature, QgsGeometry,
    QgsRectangle, QgsPointXY, QgsField, QgsWkbTypes
)
from qgis.PyQt.QtCore import QVariant
import math
import os

# 1. 载入矢量（如已加载，可忽略）
path = r"H:\Lashihai\水体提取202507132037360621.geojson"
if not os.path.exists(path):
    raise IOError("文件不存在，请检查路径！")

lyr = QgsVectorLayer(path, "WaterBody", "ogr")
if not lyr.isValid():
    raise RuntimeError("矢量加载失败！")

# 2. 取第一条要素的外轮廓
feat = next(lyr.getFeatures())
poly_geom = feat.geometry()
if poly_geom.isMultipart():
    polyline = poly_geom.asMultiPolygon()[0][0]   # 外环
else:
    polyline = poly_geom.asPolygon()[0]

# 3. 计算外框周长 → 自动网格边长
perimeter = QgsGeometry.fromPolylineXY(polyline).length()
side = perimeter / 50.0
print(f"自动计算正方形边长：{side:.2f} 地图单位")

# 4. 生成最小外接矩形 MBR
mbr = QgsGeometry.fromPolylineXY(polyline).boundingBox()
xmin, xmax = mbr.xMinimum(), mbr.xMaximum()
ymin, ymax = mbr.yMinimum(), mbr.yMaximum()

# 5. 创建网格内存图层
grid_lyr = QgsVectorLayer("Polygon?crs=" + lyr.crs().authid(), "Grid50", "memory")
prov = grid_lyr.dataProvider()
prov.addAttributes([QgsField("id", QVariant.Int)])
grid_lyr.updateFields()

# 6. 逐个生成正方形
fid = 0
y = ymax
while y > ymin:
    x = xmin
    while x < xmax:
        rect = QgsRectangle(x, y - side, x + side, y)
        f = QgsFeature()
        f.setGeometry(QgsGeometry.fromRect(rect))
        f.setAttributes([fid])
        prov.addFeature(f)
        fid += 1
        x += side
    y -= side

grid_lyr.updateExtents()
QgsProject.instance().addMapLayer(grid_lyr)

print("已生成 50 个正方形网格并加载到地图。")

使用方法

打开 QGIS → Plugins → Python Console（或直接按 Ctrl+Alt+P）。
将以上代码整段粘进去，按回车即可。
生成的 “Grid50” 临时图层会立即出现在图层面板，可右键导出为 GeoJSON/SHP 等。

矢量路径改了：H:\Lashihai\Lashihai.gpkg|layername=LashihaiWater 请帮我用新的路径改写一下上面的代码

只需把路径替换成新的 Geopackage 路径即可，其余逻辑不变。
复制下面完整代码，在 QGIS Python 控制台直接运行：

Python

# -*- coding: utf-8 -*-
"""
按外轮廓生成 50 个自适应正方形网格（新路径版）
"""

from qgis.core import (
    QgsProject, QgsVectorLayer, QgsFeature, QgsGeometry,
    QgsRectangle, QgsPointXY, QgsField, QgsWkbTypes
)
from qgis.PyQt.QtCore import QVariant
import math

# 1. 载入 Geopackage 图层
uri = r"H:\Lashihai\Lashihai.gpkg|layername=LashihaiWater"
lyr = QgsVectorLayer(uri, "LashihaiWater", "ogr")
if not lyr.isValid():
    raise RuntimeError("图层加载失败，请检查路径！")

# 2. 取第一条要素的外轮廓
feat = next(lyr.getFeatures())
poly_geom = feat.geometry()
if poly_geom.isMultipart():
    polyline = poly_geom.asMultiPolygon()[0][0]   # 外环
else:
    polyline = poly_geom.asPolygon()[0]

# 3. 计算周长 → 自动网格边长
perimeter = QgsGeometry.fromPolylineXY(polyline).length()
side = perimeter / 50.0
print(f"自动计算正方形边长：{side:.2f} 地图单位")

# 4. 最小外接矩形
mbr = QgsGeometry.fromPolylineXY(polyline).boundingBox()
xmin, xmax = mbr.xMinimum(), mbr.xMaximum()
ymin, ymax = mbr.yMinimum(), mbr.yMaximum()

# 5. 创建网格内存图层
grid_lyr = QgsVectorLayer("Polygon?crs=" + lyr.crs().authid(), "Grid50", "memory")
prov = grid_lyr.dataProvider()
prov.addAttributes([QgsField("id", QVariant.Int)])
grid_lyr.updateFields()

# 6. 生成 50 个正方形
fid = 0
y = ymax
while y > ymin:
    x = xmin
    while x < xmax:
        rect = QgsRectangle(x, y - side, x + side, y)
        f = QgsFeature()
        f.setGeometry(QgsGeometry.fromRect(rect))
        f.setAttributes([fid])
        prov.addFeature(f)
        fid += 1
        x += side
    y -= side

grid_lyr.updateExtents()
QgsProject.instance().addMapLayer(grid_lyr)

print("已生成 50 个正方形网格并加载到地图。")

Source H:/Lashihai/Lashihai.gpkg|layername=WaterLine 请帮我在这个线基础上，等距离生成50个点，然后在这50个点基础上，生成50个外接矩形，请完成Python代码

下面给出一段可直接在 QGIS Python 控制台运行的完整脚本。
逻辑：

从线图层（WaterLine）中读取第一条线要素；
按线长 L 等距生成 50 个点（含起点，不含终点）；
以每个点为中心，边长 s = L/50 生成正方形外接矩形；
输出结果到新的临时图层 “WaterLine_Grid50”。

Python

# -*- coding: utf-8 -*-
"""
在 WaterLine 上按长度等距生成 50 个点，并以每个点为中心生成正方形外接矩形
"""

from qgis.core import (
    QgsProject, QgsVectorLayer, QgsFeature, QgsGeometry,
    QgsRectangle, QgsPoint, QgsField, QgsWkbTypes
)
from qgis.PyQt.QtCore import QVariant

# 1. 载入线图层
uri = r"H:/Lashihai/Lashihai.gpkg|layername=WaterLine"
line_lyr = QgsVectorLayer(uri, "WaterLine", "ogr")
if not line_lyr.isValid():
    raise RuntimeError("线图层加载失败！")

# 2. 取第一条线
line_feat = next(line_lyr.getFeatures())
line_geom = line_feat.geometry()
line_length = line_geom.length()

# 3. 计算正方形边长
side = line_length / 50.0
print(f"正方形边长：{side:.2f} 地图单位")

# 4. 创建输出图层
grid_lyr = QgsVectorLayer("Polygon?crs=" + line_lyr.crs().authid(),
                          "WaterLine_Grid50", "memory")
prov = grid_lyr.dataProvider()
prov.addAttributes([QgsField("id", QVariant.Int)])
grid_lyr.updateFields()

# 5. 生成 50 个点 → 50 个正方形
for i in range(50):
    # 等距点位置（0~49）
    dist = (i * line_length) / 50.0
    pt = line_geom.interpolate(dist).asPoint()          # QgsPointXY
    # 以该点为中心的正方形外接矩形
    xmin = pt.x() - side / 2.0
    xmax = pt.x() + side / 2.0
    ymin = pt.y() - side / 2.0
    ymax = pt.y() + side / 2.0
    rect_geom = QgsGeometry.fromRect(QgsRectangle(xmin, ymin, xmax, ymax))

    feat = QgsFeature()
    feat.setGeometry(rect_geom)
    feat.setAttributes([i])
    prov.addFeature(feat)

grid_lyr.updateExtents()
QgsProject.instance().addMapLayer(grid_lyr)
print("已生成 50 个正方形外接矩形并加载到地图。")

使用方法

打开 QGIS → 插件 → Python 控制台 (Ctrl+Alt+P)。
将上述代码整段粘贴并按回车即可。