C#+AE分区统计像元数量

- 日理万妓 2024-04-18 13:22 37阅读 0赞

因为近期需要进行对某一区域内进行像元数量的统计，而Arcgis内置的分区统计功能输出的又是一个栅格数据感觉不是我想要的，所以就自己编了这样一个小功能。话不多说，开始进入正题。

## 一、设计思路 ##

按照老习惯，每次要编代码之前先梳理一下思路：  
1、获取栅格数据和矢量数据（栅格数据是要统计的对象，矢量数据是统计的区域）  
2、获取矢量数据中要素的包络范围（IEnvelope）并以此创建像素块（IPixelBlock）用来存储包络范围内栅格数据的像元值  
3、遍历矢量数据中的要素根据判断条件来统计像元个数  
4、显示输出结果

## 二、功能的实现 ##

1、定义全局变量

long pixCount = 0;   //记录像元的个数
    double xmin, ymin, xmax, ymax;  //记录地理坐标范围
    int RowMin, ColumnMin, RowMax, ColumnMax;    //记录像元坐标范围
    double GeoX, GeoY;   //地理坐标
    int red,green,blue;

2、获取输入数据

//获取栅格图层
    IRasterLayer rasterLayer = GetLayer(comboBox1_StaticFile.Text) as IRasterLayer;
    IRaster raster = rasterLayer.Raster;
    IRaster2 raster2 = raster as IRaster2;//创建栅格数据

//获取矢量图层
    IFeatureLayer featureLayer = GetLayer(comboBox2_StaticZone.Text) as IFeatureLayer;
    IFeatureCursor cursor = featureLayer.Search(null, false);   //矢量游标
    IFeature feature = cursor.NextFeature();
    IPoint point = new PointClass();

3、遍历矢量数据中的要素，并统计各个要素内满足条件的像元个数

while (feature != null)
    {
         IPolygon polygon = feature.Shape as IPolygon;
    
         //获取当前矢量范围IEnvelope行列号范围，ColMin，ColMin，ColMax，RowMax
         IEnvelope envelope = feature.Extent;
         envelope.QueryCoords(out xmin, out ymin, out xmax, out ymax);
         raster2.MapToPixel(xmin, ymin, out ColumnMin,out RowMax);
         raster2.MapToPixel(xmax, ymax, out ColumnMax,out RowMin);
    
         //定义像素块的尺寸
         IPnt size = new PntClass();
         size.X = Math.Abs(ColumnMax - ColumnMin);       //像素块中的每一个像元都对应着栅格图中的一个点
         size.Y = Math.Abs(RowMax - RowMin);             //像素块只需要足够大就可以达到目的
    
         //创建像素块
         IPixelBlock pb = raster.CreatePixelBlock(size);
    
         //定义像素块起始位置坐标
         IPnt orginPoint = new PntClass();
         orginPoint.SetCoords(ColumnMin, RowMin);
    
         //读取数据到像素块中
         raster.Read(orginPoint, pb);
         IPixelBlock3 pb3 = pb as IPixelBlock3;
    
         //判断点与面之间的关系器
         IRelationalOperator pRelOpt = polygon as IRelationalOperator;
    
         //遍历每一个像元
         for (int i = ColumnMin; i < ColumnMax; i++)
         {
             for (int j = RowMin; j < RowMax; j++)
             {
                 raster2.PixelToMap(i, j, out GeoX, out GeoY);
                 point.PutCoords(GeoX, GeoY);
                 //获取像元值
                 red = Convert.ToInt32(pb3.GetVal(0,i-ColumnMin,j-RowMin));
                 green = Convert.ToInt32(pb3.GetVal(1, i - ColumnMin, j - RowMin));
                 blue = Convert.ToInt32(pb3.GetVal(2, i - ColumnMin, j - RowMin));
    
                 if (pRelOpt.Contains(point) && green > 50) //判断条件
                 {
                     pixCount++;
                 }
             }
         }
    
         MessageBox.Show("像元个数：" + pixCount.ToString());
         pixCount = 0;
         feature = cursor.NextFeature();
     }

## 三、实验的效果和小结 ##

![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70]  
小结：既然经历了这个过程，就要总结一下在这中间所遇到的问题，要不然就白痛苦了哈哈。  
1、IPixelBlock的创建是通过栅格数据的CreatePixelBlock函数来创建的，像素块的尺寸通过IPnt来确定  
2、在将数据读取到像素块中起始的坐标也是通过IPnt来确定  
3、判断点线面的关系可以使用IRelationalOperator接口

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ_size_16_color_FFFFFF_t_70]: https://img-blog.csdnimg.cn/20190907115439955.PNG?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L2RheXVoYWl0YW5nMQ==,size_16,color_FFFFFF,t_70