页面载入中....  
页面载入中....
时 间 记 忆
页面载入中....
最 新 评 论
页面载入中....
专 题 分 类
页面载入中....
最 新 日 志
页面载入中....
最 新 留 言
页面载入中....
搜 索
用 户 登 录
页面载入中....
友 情 连 接
博 客 信 息
页面载入中....


 
转载:AE创建泰森多边形
[ 2008-5-5 16:07:00 | By: luyuan-cyj ]
 

AE创建泰森多边形 
 

ArcEngin提供了TIN组件,可以利用等高线或离散的高程值生成TIN的方法,其实等高线和离散点都可以看作是要素类—feature class
AE的组件TIN提供了一个非常重要的接口ITinEdit利用这个接口,可以方便地生成三维的TIN,
思路:1找到一个数据集,作为生成TIN的数据源,我们用shp文件来做,它必须有一个字段里含用高程值,2.利用ITinEdit接口生成TIN

重要的是这两个方法:InitNew(int pExtent,,IEnvelope)它需要一个IEnvelope接口的对象

以及方法:SetSpatialReference(int pSpatialReference,ISPatialReference);因此我们首先创建一个ISpatialReference接口对象,利用get_SpatialReference属性得到要素类的空间引用并将空间引用赋给ISpatialReference接口对象,再创建一个IEnvelope接口对象,利用putref_ SpatialReference属性将前面创建的ISpatialReference接口对象的空间引用赋给IEnvelope接口对象,再创建一个TIN对象并得到其ITinEdit接口,利用此接口的InitNew方法创建TIN,而前面创建的IEnvelope接口对象作为该方法的输入对象,这样就可以生成TIN了.
3.上面已能生成TIN,但生成的TIN并没有高程信息,要获得高程字段,再生成三维TIN,这要用到IFields接口和IField的指针,用其get_Fields方法获取高程信息.

代码如下:
//创建TIN数据,同时根据以有的要素类来创建泰森多边形
 
Public Sub CreateTin()
  Dim TinSurface As ITin
  Set TinSurface = New Tin
  Dim FeatClass As IFeatureClass
  Dim Pdocment As IMxDocument
  Set Pdocment = ThisDocument
  Dim pFeatureLayer As IFeatureLayer
  Set pFeatureLayer = Pdocment.ActiveView.FocusMap.Layer(0)
 
  Dim pEnv As IEnvelope
  Set pEnv = pFeatureLayer.AreaOfInterest.Envelope
 
  Dim pTinEdit As ITinEdit
  Set pTinEdit = New Tin
  pTinEdit.InitNew pEnv
  pTinEdit.SaveAs "H:\New Folder\tin\a3", False '创建放的路径
  pTinEdit.StartEditing
   
 
  Set FeatClass = pFeatureLayer.FeatureClass
 
  Dim pTagFeild As IField
  Set pTagFeild = New Field
 
 
  Dim pHightFeild As IField
  Set pHightFeild = FeatClass.Fields.Field(0)
  'MsgBox pHightFeild.Name
 
'//其中FeatClass.Fields.Field(2), 是要素类中的某一字段  作为Tin的高程信息
  pTinEdit.AddFromFeatureClass FeatClass, Nothing, FeatClass.Fields.Field(2), FeatClass.Fields.Field(0), esriTinMassPoint
  pTinEdit.StopEditing (True)
  pTinEdit.Refresh
 
  Dim pTinNodeCollection As ITinNodeCollection
  Set pTinNodeCollection = pTinEdit
  MsgBox pTinNodeCollection.NodeCount
  Dim ptin As ITin
  Set ptin = pTinEdit
  MsgBox ptin.DataNodeCount
  Dim pNewfeatureclass As IFeatureClass
  Set pNewfeatureclass = OpenFeatureClass_Example
 
  MsgBox ptin.DataEdgeCount
     ///创建泰森多边形  OpenFeatureClass_Example是 找到的要素类图层
  pTinNodeCollection.ConvertToVoronoiRegions OpenFeatureClass_Example, Nothing, Nothing, "NodeIndex", "asdf"
 
End Sub
 
//添加要素类的图层函数
Public Function OpenFeatureClass_Example() As IFeatureClass
  Dim pWorkspaceFactory As IWorkspaceFactory
  Set pWorkspaceFactory = New ShapefileWorkspaceFactory
  Dim pFeatureWorkspace As IFeatureWorkspace
  Set pFeatureWorkspace = pWorkspaceFactory.OpenFromFile("H:\New Folder", 0)
 
'  Dim pWorkspaceEdit  As IWorkspaceEdit
'  Set pWorkspaceEdit = pFeatureWorkspace
'  pWorkspaceEdit.StartEditing True
 
  Dim pFeatureClass As IFeatureClass
  Set pFeatureClass = pFeatureWorkspace.OpenFeatureClass("Vor3")
  Set OpenFeatureClass_Example = pFeatureClass
 end function

 
发表评论:
页面载入中....
Powered by Oblog.