资源环境信息系统培训资料.ppt

基于VisualC 和ArcEngine的二次开发 系统功能需求 简单图层控制功能的实现数据的转换和可视化表达 一 简单图层控制功能的实现 图层属性表的添加和显示图层的添加和移除图层间关系的调整 1图层属性表的添加和显示 步骤1 设计图层控制的菜单 响应ContextMenuStrip菜单控件 为图层控件添加菜单项 显示属性表菜单和删除图层菜单 该菜单是在鼠标右键点击图层时弹出的快捷菜单 添加的菜单项 修改菜单name属性 步骤2 实现图层与图层控制菜单间的关联 响应TOCControl控件的OnMouseDown事件 添加相关代码 usingESRI ArcGIS Carto 因为程序中用到ILayer和 IBasicMap接口 所以需要添加命名空间ILayerlayer newFeatureLayerClass privatevoidaxTOCControl1 OnMouseDown objectsender ESRI ArcGIS Controls ITOCControlEvents OnMouseDownEvente this axTOCControl1 ContextMenuStrip null IBasicMapmap newMapClass System Objectother null System Objectindex null ESRI ArcGIS Controls esriTOCControlItemitem ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemNone this axTOCControl1 HitTest e x e y refitem refmap reflayer refother refindex if item ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemLayer 运行程序显示结果如下 步骤3 设计一个属性表窗体 1 添加一个属性表窗体 2 为属性表窗体添加DataGridView控件 并修改控件Anchor属性为 Top Bottom Left Right 3 为属性表窗体编写相关代码 实现图层属性信息往属性表中的添加 publicpartialclass属性表 Form privateILayerlayer public属性表 ILayerlyr 对属性表窗口进行初始化 InitializeComponent layer lyr this Text layer Name 属性表 privatevoid属性表 Load objectsender EventArgse try ITablelyrtable ITable layer DataTabletable newDataTable IFieldfield for inti 0 i lyrtable Fields FieldCount i field lyrtable Fields get Field i table Columns Add field Name object values newobject lyrtable Fields FieldCount IQueryFilterqueryFilter newQueryFilterClass ICursorcursor lyrtable Search queryFilter true IRowrow while row cursor NextRow null for intj 0 j lyrtable Fields FieldCount j objectob row get Value j values j ob table Rows Add values this dataGridView1 DataSource table catch Exceptione1 MessageBox Show 无法显示属性表 this Close 由于该段代码中用到了ILayer接口 所以要在属性表 cs开头添加命名空间ESRI ArcGIS Carto 同时也用到ITable IQueryFilter ICursor IRow接口 所以要添加命名空间ESRI ArcGIS Geodatabase 即添加如下代码 usingESRI ArcGIS Carto usingESRI ArcGIS Geodatabase 但此时运行程序会报错 命名空间ESRI ArcGIS中不存在类型名或命名空间名称Geodatabase 是缺少程序集引用吗 这时需要在资源管理器中将该程序集添加进来 右键点击引用 从弹出的菜单中点击 添加引用 菜单项 再在弹出的对话框中将ESRI ArcGIS Geodatabase命名空间添加进来 步骤4 建立图层控制菜单与属性表的连接 实现属性表界面的显示 双击图层控制菜单上的 显示属性表 菜单项 在响应的Click事件中添加如下代码 privatevoid显示属性表ToolStripMenuItem Click objectsender EventArgse 属性表form new属性表 layer form Show 步骤5 运行程序 查看图层属性 右键点击该图层 从弹出快捷菜单中选择 显示属性表 菜单项 弹出该图层的属性表界面 步骤6 这里实现一个通过属性表中属性查找图形要素的小功能 响应datagridView1控件的SelectionChanged事件 先在主窗口中编写一个公有函数 以获取主窗口的MapControl控件 因为该控件是Private类型 在其它窗体中不能直接调用 publicESRI ArcGIS Controls AxMapControlgetMapControl returnaxMapControl1 再在SelectionChanged事件中添加如下代码 privatevoiddataGridView1 SelectionChanged objectsender EventArgse DataGridViewSelectedRowCollectionSelRows this dataGridView1 SelectedRows DataGridViewRowrow 主窗口form 主窗口 Application OpenForms 0 IMapm form getMapControl Map m ClearSelection for inti 0 i SelRows Count i row SelRows i intID Convert ToInt32 row Cells FID Value ToString IFeatureLayerflyr IFeatureLayer layer IFeatureClassfeaturecls flyr FeatureClass IFeaturefeature featurecls GetFeature ID m SelectFeature layer feature 获取属性表中选中行对应的图形要素 form getMapControl Refresh 运行程序 显示结果如下 属性表中选择一个点要素时 将其高亮显示 属性表中拖动选择多个点要素时 将其高亮显示 2图层的添加和移除 步骤1 在ToolBarControl控件中直接添加 图层添加 按钮 即可以实现图层添加功能 搭建的框架中已有该按钮 步骤2 前面建立的图层控制菜单中已经包含了删除图层的菜单 此时只需要响应该菜单 编辑相关程序代码 即可以实现图层删除功能 双击该菜单项 响应Click事件 在响应的Click事件中添加相关代码 privatevoid删除此图层ToolStripMenuItem Click objectsender EventArgse for inti 0 i this axMapControl1 Map LayerCount i if this axMapControl1 Map get Layer i layer this axMapControl1 DeleteLayer i 运行程序 显示如下 删除biaozhidian图层前 删除biaozhidian图层后 3图层间关系的调整 在ArcMap中 要调整图层的现实顺序 只需要简单的鼠标拖拽就可以完成 既简单又直观 但是TOCControl控件并未实现此功能 本实例即实现此功能 在本实例中主要利用TOCControl控件的OnMouseDown OnMouseUp两个时间和HitTest Update 方法 OnMouseDown事件方法的主要功能是获取源图层pMovelayer OnMouseUp事件方法的主要功能是 1 获取目标图层及其索引号 2 利用IMap提供的Movelayer pMoveLayer Toindex 将源图层pMoveLayer移到目标图层 Toindex 的下方 3 使用ITOCControl提供的Update 方法更新TOCControl控件中显示的内容 步骤1 在主窗口 cs中申明两个全局变量 并响应TOCControl控件的OnMouseDown事件 前面已响应此事件 添加相关代码 申明全局变量ILayerpMoveLayer intToindex privatevoidaxTOCControl1 OnMouseDown objectsender ESRI ArcGIS Controls ITOCControlEvents OnMouseDownEvente 前面已添加代码 这里不需要再添加this axTOCControl1 ContextMenuStrip null IBasicMapmap newMapClass System Objectother null System Objectindex null ESRI ArcGIS Controls esriTOCControlItemitem ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemNone this axTOCControl1 HitTest e x e y refitem refmap reflayer refother refindex if item ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemLayer layer是前面申明的一个全局变量 步骤2 响应TOCControl控件的OnMouseUp事件 并添加相关代码 双击OnMouseUp事件 privatevoidaxTOCControl1 OnMouseUp objectsender ESRI ArcGIS Controls ITOCControlEvents OnMouseUpEvente if e button 1 ESRI ArcGIS Controls esriTOCControlItemitem ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemNone IBasicMapmap null layer null objectother null objectindex null this axTOCControl1 HitTest e x e y refitem refmap reflayer refother refindex HitTest是返回特定坐标下 的TOCControl控件的itemIMappMap axMapControl1 ActiveView FocusMap 获取当前窗口中地图 if item ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemLayer layer null if pMoveLayer layer ILayerpTempLayer for inti 0 i pMap LayerCount i pTempLayer pMap get Layer i if pTempLayer layer Toindex i pMap MoveLayer pMoveLayer Toindex axMapControl1 ActiveView Refresh this axTOCControl1 Update 这里可以添加命名空间usingESRI ArcGIS Controls 这样就可以在编写一些代码时简化程序 如 ESRI ArcGIS Controls esriTOCControlItemitem ESRI ArcGIS Controls esriTOCControlItem esriTOCControlItemNone 中 ESRI ArcGIS Controls 就可以省略 步骤3 运行程序 显示结果 图层位置调整前 图层位置调整后 二 数据的转换和可视化表达 根据环保局 规划局等部门的需求分析 我们往往需要能在地形图上直观查看一些污染源信息 即要实现重点污染源 一般污染源数信息 数值数据 的可视化表达 图形数据 本实例中即需要实现重点污染源信息的可视化表达 这是一个典型的将数值型数据转换为图形数据的过程 功能演示 点击 查看重点污染源 菜单后 G101表中的重点污染源信息 主要是由经纬度标志的位置信息 生成点图层显示在图形显示窗口中 新生成的点图层 新生成点图层在地图上的显示 功能实现流程 1添加一个查看重点污染源的菜单 2点击该菜单时 生成一个 重点污染源 点图层 并在地图显示区域显示出来 1 从数据源 G101表 中获取重点污染源的代码 名称 中心经度 度 中心经度 分 中心经度 秒 中心纬度 度 中心纬度 分 中心纬度 秒 字段信息 2 利用坐标转换公式将重点污染源的地理坐标 经纬度坐标 转换为高斯坐标 这里是北京54投影坐标系 3 将转换坐标后的数据生成点图层 显示在图层显示区 并为该图层添加属性表信息 步骤1 添加数据显示菜单以及查看重点污染源菜单项 步骤2 双击 查看重点污染源 菜单项 响应该菜单命令 并编辑相关代码 获取G101表中重点污染源的经纬度信息 通过相关坐标转换公式将经纬度转换为高斯坐标 并生成点图层在地图显示区显示出来 在主窗口 cs开头添加如下命名空间 其中 ESRI ArcGIS Geodatabase ESRI ArcGIS DataSourcesGDB等需要在引用中 添加进来usingESRI ArcGIS Geodatabase usingSystem Data OleDb usingESRI ArcGIS Geometry usingESRI ArcGIS esriSystem usingESRI ArcGIS DataSourcesGDB 一个创建图层的函数publicIFeatureLayerCreateFeatureLayerInmemeory stringDataSetName stringAliaseName ISpatialReferenceSpatialRef esriGeometryTypeGeometryType IFieldsPropertyFields IWorkspaceFactoryworkspaceFactory newInMemoryWorkspaceFactoryClass ESRI ArcGIS Geodatabase IWorkspaceNameworkspaceName workspaceFactory Create MyWorkspace null 0 ESRI ArcGIS esriSystem INamename IName workspaceName ESRI ArcGIS Geodatabase IWorkspaceinmemWor IWorkspace name Open IFieldoField newFieldClass IFieldsoFields newFieldsClass IFieldsEditoFieldsEdit null IFieldEditoFieldEdit null IFeatureClassoFeatureClass null IFeatureLayeroFeatureLayer null try oFieldsEdit oFieldsasIFieldsEdit oFieldEdit oFieldasIFieldEdit for inti 0 i PropertyFields FieldCount i oFieldsEdit AddField PropertyFields get Field i IGeometryDefgeometryDef newGeometryDefClass IGeometryDefEditgeometryDefEdit IGeometryDefEdit geometryDef geometryDefEdit AvgNumPoints 2 5 geometryDefEdit GeometryType 2 GeometryType geometryDefEdit GridCount 2 1 geometryDefEdit HasM 2 false geometryDefEdit HasZ 2 false geometryDefEdit SpatialReference 2 SpatialRef 创建Shape字段oFieldEdit Name 2 SHAPE oFieldEdit Type 2 esriFieldType esriFieldTypeGeometry oFieldEdit GeometryDef 2 geometryDef oFieldEdit IsNullable 2 true oFieldEdit Required 2 true oFieldsEdit AddField oField oFeatureClass inmemWorasIFeatureWorkspace CreateFeatureClass DataSetName oFields null null esriFeatureType esriFTSimple SHAPE oFeatureClassasIDataset BrowseName DataSetName oFeatureLayer newFeatureLayerClass oFeatureLayer Name AliaseName oFeatureLayer FeatureClass oFeatureClass catch finally try System Runtime InteropServices Marshal ReleaseComObject oField System Runtime InteropServices Marshal ReleaseComObject oFields System Runtime InteropServices Marshal ReleaseComObject oFieldsEdit System Runtime InteropServices Marshal ReleaseComObject oFieldEdit System Runtime InteropServices Marshal ReleaseComObject name System Runtime InteropServices Marshal ReleaseComObject workspaceFactory System Runtime InteropServices Marshal ReleaseComObject workspaceName System Runtime InteropServices Marshal ReleaseComObject inmemWor System Runtime InteropServices Marshal ReleaseComObject oFeatureClass catch GC Collect returnoFeatureLayer 响应查看重点污染源菜单项的Click事件privatevoid查看重点污染源ToolStripMenuItem Click objectsender EventArgse 创建图层IFieldsPropertyFields newFieldsClass IFieldsEditFieldsEdit IFieldsEdit PropertyFields 创建属性表 添加属性字段并编辑列名IFieldField newFieldClass IFieldEditFieldEdit IFieldEdit Field FieldEdit Name 2 FID FieldEdit AliasName 2 FID FieldEdit Type 2 esriFieldType esriFieldTypeOID FieldsEdit AddField FieldEdit Field newFieldClass FieldEdit IFieldEdit Field FieldEdit Name 2 代码 FieldEdit AliasName 2 代码 FieldEdit Type 2 esriFieldType esriFieldTypeString FieldsEdit AddField FieldEdit Field newFieldClass FieldEdit IFieldEdit Field FieldEdit Name 2 x FieldEdit AliasName 2 x FieldEdit Type 2 esriFieldType esriFieldTypeDouble FieldsEdit AddField FieldEdit Field newFieldClass FieldEdit IFieldEdit Field FieldEdit Name 2 y FieldEdit AliasName 2 y FieldEdit Type 2 esriFieldType esriFieldTypeDouble FieldsEdit AddField FieldEdit 添加点特征图层IFeatureLayerPointFeatureLyr CreateFeatureLayerInmemeory 重点污染源 重点污染源 this axMapControl1 Map SpatialReference esriGeometryType esriGeometryPoint PropertyFields ITabletable ITable PointFeatureLyr 图层属性表中各个属性的索引intindex 名称 table FindField 名称 intindex 代码 table FindField 代码 intindex x table FindField x intindex y table FindField y intindex shape table FindField shape 建立与数据库的连接 通过SQL获取需要属性字段值stringsql SELECTG101 代码 G101 名称 G101 中心经度 度 G101 中心经度 分 G101 中心经度 秒 G101 中心纬度 度 G101 中心纬度 分 G101 中心纬度 秒 FROMG101 OleDbConnectionconn newOleDbConnection Provider Microsoft Jet OLEDB 4 0 DataSource Application StartupPath Database mdb OleDbCommandcmd newOleDbCommand sql conn intindex code 0 index name 0 index jd D 0 index jd F 0 index jd M 0 index wd D 0 index wd F 0 index wd M 0 intjd D 0 jd F 0 jd M 0 wd D 0 wd F 0 wd M 0 doublejd exp M 经度用秒表示doublewd exp M 纬度用秒表示doubleresult x 0 result y 0 转换之后的北京54投影坐标stringfield Name try conn Open OleDbDataReaderreader cmd ExecuteReader intcols reader FieldCount object values newobject cols for intindex Field 0 index Field cols index Field field Name reader GetName index Field switch field Name case 代码 index code index Field break case 名称 index name index Field break case 中心经度 度 index jd D index Field break case 中心经度 分 index jd F index Field break case 中心经度 秒 index jd M index Field break case 中心纬度 度 index wd D index Field break case 中心纬度 分 index wd F index Field break case 中心纬度 秒 index wd M index Field break IPointpt while reader Read reader GetValues values IRowrow table CreateRow row set Value index 名称 values index name row set Value index 代码 values index code try jd D Convert ToInt32 values index jd D ToString catch Exceptione1 jd D 0 try jd F Convert ToInt32 values index jd F ToString catch Exceptione1 jd F 0 try jd M Convert ToInt32 values index jd M ToString catch Exceptione1 jd M 0 try wd D Convert ToInt32 values index wd D ToString catch Exceptione1 wd D 0 try wd F Convert ToInt32 values index wd F ToString catch Exceptione1 wd F 0 try wd M Convert ToInt32 values index wd M ToString catch Exceptione1 wd M 0 jd exp M jd D 3600 jd F 60 jd M wd exp M wd D 3600 wd F 60 wd M GeoToGauss jd exp M wd exp M 20 6 refresult y refresult x row set Value index x result x row set Value index y result y pt newPointClass pt X result y pt Y result x row set Value index shape pt row Store reader Close catch Exceptionex MessageBox Show ex Message MessageBoxButtons OK MessageBoxIcon Error finally conn Close this axMapControl1 AddLayer ILayer PointFeatureLyr GeoToGauss函数功能说明 1 将地理坐标 wd jd 转换成绝对的高斯坐标 y x 2 本函数支持基于六度带 或三度带 克拉索夫斯基椭球进行转换 适用范围 本函数适用于将地球东半球中北半球 即东经0度到东经180度 北纬0度至90度 范围内所有地理坐标到高斯坐标的转换 使用说明 调用本函数后返回的结果应在满足精度的条件下进行四舍五入 doublejd 输入参数 地理坐标的经度 以秒为单位 doublewd 输入参数 地理坐标的纬度 以秒为单位 shortDH 输入参数 三度带或六度带的带号 六度带 三度带 的带号是这样得到的 从东经0度到东经180度自西向东按每6度 3度 顺序编号 编号从1开始 这个顺序编号就称为六度带 三度带 的带号 因此 六度带的带号的范围是1 30 三度带的带号的范围是1 60 如果一个点在图号为TH的图幅中 那麽该点所处的六度带的带号就可以这样得到 将该图号的第3 4位组成的字符串先转换成数字 再减去30 例如某点在图幅06490701中 该点所在的带号就是49 30 即19 如果调用本函数去进行一般的从地理坐标到基于六度带高斯坐标的变换 非邻