各类图形库创建图层的方式示例

不同的图形库和地理信息系统（GIS）软件提供了不同的方式来创建图层。以下是一些常见的图形库和GIS软件的创建图层的方式和示例代码：

1. **使用GDAL/OGR库**：
   GDAL（Geospatial Data Abstraction Library）和OGR（Simple Features Library）是用于读写地理空间数据的库。你可以使用它们来创建图层并写入地理空间数据。

   示例代码（C++）：
   ```cpp

  #include "gdal/ogr_api.h"

   int main() {
       // 创建一个源
       OGRRegisterAll();
       OGRDataSource *poDS = OGRSFDriverRegistrar::Open("data_source_name", TRUE);

       // 创建一个图层
       OGRLayer *poLayer = poDS->CreateLayer("layer_name", NULL, wkbPoint, NULL);

       // 添加
       OGRFieldDefn oField("field_name", OFTString);
       oField.SetWidth(32);
       poLayer->CreateField(&oField);

       // 创建要素并写入图
       OGRFeature *poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());
       poFeature->SetField("field_name", "field_value");
       OGRPoint pt;
       pt.setX(1.0);
       pt.setY(1.0);
       poFeature->SetGeometry(&pt);
       poLayer->CreateFeature(poFeature);

       // 释放资源       OGRFeature::DestroyFeature(poFeature);
       OGRDataSource::DestroyDataSource(poDS);
       
       return 0;
   }
   ```

2. **使用Python的Fiona库**：
   Fiona是Python中用于读写地理空间数据的库，它可以与GDAL/OGR一起使用来创建图层并写入地理空间数据。

   示例代码（Python）：
   ```python  

import fiona
   from fiona.crs import from_epsg

   # 创建图层
   schema = {
       'geometry': 'Point',
       'properties': {'name': 'str'}
   }
   with fiona.open('output.shp', 'w', 'ESRI Shapefile', schema, crs=from_epsg(4326)) as layer:
       # 添加要素
       point = {
           'type': 'Point',
           'coordinates': (1.0, 1.0)
       }
       layer.write({
           'geometry': point,
           'properties': {'name': 'example'}
       })
   ```

3. **使用QGIS软件**：
   在QGIS软件中，你可以通过界面操作或者使用Python脚本来创建图层。

   示例代码（Python）：
   ```python

  from qgis.core import QgsVectorLayer, QgsField, QgsPoint, QgsFeature, QgsGeometry

   # 创建图层
   layer = QgsVectorLayer("Point", "layer_name", "memory")
   pr = layer.dataProvider()
   pr.addAttributes([QgsField("name",  QVariant.String)])
   layer.updateFields()

   # 添加素
   feature = QgsFeature()
   feature.setGeometry(QgsGeometry.fromPointXY(QgsPoint(1.0, 1.0)))
   feature.setAttributes(["example"])
   pr.addFeature(feature)

   # 刷新地图
   layer.updateExtents()
   QgsProject.instance().addMapLayer(layer)
   ```

未完待续……