全国高分辨率土地利用数据服务 土地利用数据服务 土地覆盖数据服务 坡度数据服务 土壤侵蚀数据服务 全国各省市DEM数据服务 耕地资源空间分布数据服务 草地资源空间分布数据服务 林地资源空间分布数据服务 水域资源空间分布数据服务 建设用地空间分布数据服务 地形、地貌、土壤数据服务 分坡度耕地数据服务 全国大宗农作物种植范围空间分布数据服务
多种卫星遥感数据反演植被覆盖度数据服务 地表反照率数据服务 比辐射率数据服务 地表温度数据服务 地表蒸腾与蒸散数据服务 归一化植被指数数据服务 叶面积指数数据服务 净初级生产力数据服务 净生态系统生产力数据服务 生态系统总初级生产力数据服务 生态系统类型分布数据服务 土壤类型质地养分数据服务 生态系统空间分布数据服务 增强型植被指数数据服务
多年平均气温空间分布数据服务 多年平均降水量空间分布数据服务 湿润指数数据服务 大于0℃积温空间分布数据服务 光合有效辐射分量数据服务 显热/潜热信息数据服务 波文比信息数据服务 地表净辐射通量数据服务 光合有效辐射数据服务 温度带分区数据服务 山区小气候因子精细数据服务
全国夜间灯光指数数据服务 全国GDP公里格网数据服务 全国建筑物总面积公里格网数据服务 全国人口密度数据服务 全国县级医院分布数据服务 人口调查空间分布数据服务 收入统计空间分布数据服务 矿山面积统计及分布数据服务 载畜量及空间分布数据服务 农作物种植面积统计数据服务 农田分类面积统计数据服务 农作物长势遥感监测数据服务 医疗资源统计数据服务 教育资源统计数据服务 行政辖区信息数据服务
Landsat 8 高分二号 高分一号 SPOT-6卫星影像 法国Pleiades高分卫星 资源三号卫星 风云3号 中巴资源卫星 NOAA/AVHRR MODIS Landsat TM 环境小卫星 Landsat MSS 天绘一号卫星影像
先来一个效果图:
创建符合要求的表
//输出结果为一个张表,这张表有3个字段,其中面ID为面要素数据的FID
public ITable CreateTable (string _TablePath, string _TableName)
{
IWorkspaceFactory pWks = new ShapefileWorkspaceFactoryClass ();
IFeatureWorkspace pFwk = pWks.OpenFromFile (_TablePath, 0) as IFeatureWorkspace;
//用于记录面中的ID;
IField pFieldID = new FieldClass ();
IFieldEdit pFieldIID = pFieldID as IFieldEdit;
pFieldIID.Type_2 = esriFieldType.esriFieldTypeInteger;
pFieldIID.Name_2 = "面ID";
//用于记录个数的;
IField pFieldCount = new FieldClass ();
IFieldEdit pFieldICount = pFieldCount as IFieldEdit;
pFieldICount.Type_2 = esriFieldType.esriFieldTypeInteger;
pFieldICount.Name_2 = "个数";
//用于添加表中的必要字段
ESRI.ArcGIS.Geodatabase.IObjectClassDescription objectClassDescription = new ESRI.ArcGIS.Geodatabase.ObjectClassDescriptionClass ();
IFields pTableFields = objectClassDescription.RequiredFields;
IFieldsEdit pTableFieldsEdit = pTableFields as IFieldsEdit;
pTableFieldsEdit.AddField (pFieldID);
pTableFieldsEdit.AddField (pFieldCount);
ITable pTable = pFwk.CreateTable (_TableName, pTableFields, null, null,"");
return pTable;
}
统计需要的数据
// 第一个参数为面数据,第二个参数为点数据,第三个为输出的表
public void StatisticPointCount (IFeatureClass _pPolygonFClass, IFeatureClass _pPointFClass, ITable _pTable)
{
IFeatureCursor pPolyCursor = _pPolygonFClass.Search (null, false);
IFeature pPolyFeature = pPolyCursor.NextFeature ();
while (pPolyFeature != null)
{
IGeometry pPolGeo = pPolyFeature.Shape;
int Count = 0;
ISpatialFilter spatialFilter = new SpatialFilterClass ();
spatialFilter.Geometry = pPolGeo;
spatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelContains;
IFeatureCursor pPointCur = _pPointFClass.Search (spatialFilter,false);
if (pPointCur != null)
{
IFeature pPointFeature = pPointCur.NextFeature ();
while (pPointFeature != null)
{
pPointFeature = pPointCur.NextFeature ();
Count++;
}
}
if (Count != 0)
{
IRow pRow = _pTable.CreateRow ();
pRow.set_Value (1, pPolyFeature.get_Value (0));
pRow.set_Value (2, Count);
pRow.Store ();
}
pPolyFeature = pPolyCursor.NextFeature ();
}
}
效果如下:
上面这个例子只是用了空间过滤,没有用到属性过滤,我们将上面的代码稍微改动下,加上一句代码即可, 如下:
结果对照: