将shp转化为geodatabase

       将shp转化为geodatabase

       如何将shape文件转化成GeoDataBase(各种文件格式的转换)

       本例演示的是如何将shape文件转化成personal GeoDatabase文件，其它格式间的与此转换类似。主要用到IFeatureDataConverter接口的ConvertFeatureClass方法。

       l   要点

       首先，创建新的GeoDataBase数据库，并创建IFeatureDatasetName对象。创建定义两个IFeatureClassName接口对象分别引用输入表（shape文件）和输出表。

然后设置输出表的Shape字段的GeormetryDef属性。这一步非常关键，因为其中包含了数据库和shape文件的空间参考信息。

       最后调用IFeatureDataConverter.ConvertFeatureClass方法完成功能。

       l   程序说明

       过程UIBConvert_Click是实现模块，调用过程ConvertShapeToGeodatabase实现功能。

       sDataPath定义了数据与工程文件的相对路径。SHAPE_NAME描述了要转化的shape文件的文件名。MDB_NAME和F_DS_NAME分别描述了Access数据库名和库的数据集的名称。

       Option Explicit

       Private Sub UIBConvert_Click()

       Call ConvertShapeToGeodatabase

       End Sub

Private Sub ConvertShapeToGeodatabase()

    Dim pOutWorkspaceFactory    As IWorkspaceFactory

    Dim pOutWorkspaceName       As IWorkspaceName

    Dim pInWorkspaceName        As IWorkspaceName

    Dim pOutFeatureDSName       As IFeatureDatasetName

    Dim pOutDSName              As IDatasetName

    Dim pInFeatureClassName     As IFeatureClassName

    Dim pInDatasetName          As IDatasetName

    Dim pOutFeatureClassName    As IFeatureClassName

    Dim pOutDatasetName         As IDatasetName

    Dim iCounter                As Long

    Dim pOutFields              As IFields

    Dim pInFields               As IFields

    Dim pFieldChecker           As IFieldChecker

    Dim pGeoField               As IField

    Dim pOutGeometryDef         As IGeometryDef

    Dim pOutGeometryDefEdit     As IGeometryDefEdit

    Dim pName                   As IName

    Dim pInFeatureClass         As IFeatureClass

    Dim pShpToFeatClsConverter  As IFeatureDataConverter

    Dim pVBProject              As VBProject

    Dim sDataPath               As String

    Const SHAPE_NAME As String = "country"

    Const MDB_NAME As String = "countryDB"

    Const F_DS_NAME As String = "World"

    On Error GoTo ErrorHandler

    Set pVBProject = ThisDocument.VBProject

    sDataPath = pVBProject.FileName & "\..\..\..\..\data\"

    If Not "" = Dir(sDataPath & MDB_NAME & ".mdb") Then

        MsgBox MDB_NAME & ".mdb already exist"

        Exit Sub

    Else

        ' Create a new Access database

        Set pOutWorkspaceFactory = New AccessWorkspaceFactory

        Set pOutWorkspaceName = pOutWorkspaceFactory.Create(sDataPath, MDB_NAME, Nothing, 0)

        ' create a new feature datset name object for the output Access feature dataset, call

        ' it "World"

        Set pOutFeatureDSName = New FeatureDatasetName

        Set pOutDSName = pOutFeatureDSName

        Set pOutDSName.WorkspaceName = pOutWorkspaceName

        pOutDSName.Name = F_DS_NAME

        ' Get the name object for the input shapefile workspace

        Set pInWorkspaceName = New WorkspaceName

        pInWorkspaceName.PathName = sDataPath

        pInWorkspaceName.WorkspaceFactoryProgID = _

"esriCore.ShapefileWorkspaceFactory.1"

        Set pInFeatureClassName = New FeatureClassName

        Set pInDatasetName = pInFeatureClassName

        pInDatasetName.Name = SHAPE_NAME

        Set pInDatasetName.WorkspaceName = pInWorkspaceName

        ' Create the new output FeatureClass name object that will be passed

        '   into the conversion function

        Set pOutFeatureClassName = New FeatureClassName

        Set pOutDatasetName = pOutFeatureClassName

        ' Set the new FeatureClass name to be the same as the input FeatureClass name

        pOutDatasetName.Name = pInDatasetName.Name

        ' Open the input Shapefile FeatureClass object, so that we can get its fields

        Set pName = pInFeatureClassName

        Set pInFeatureClass = pName.Open

        ' Get the fields for the input feature class and run them through

        '   field checker to make sure there are no illegal or duplicate field names

        Set pInFields = pInFeatureClass.Fields

        Set pFieldChecker = New FieldChecker

        pFieldChecker.Validate pInFields, Nothing, pOutFields

        ' Loop through the output fields to find the geometry field

        For iCounter = 0 To pOutFields.FieldCount

            If pOutFields.Field(iCounter).Type = esriFieldTypeGeometry Then

                Set pGeoField = pOutFields.Field(iCounter)

                Exit For

            End If

        Next iCounter

        ' Get the geometry field's geometry definition

        Set pOutGeometryDef = pGeoField.GeometryDef

        ' Give the geometry definition a spatial index grid count and grid size

        Set pOutGeometryDefEdit = pOutGeometryDef

        pOutGeometryDefEdit.GridCount = 1

        pOutGeometryDefEdit.GridSize(0) = 1500000

        ' Now use IFeatureDataConverter::Convert to create the output FeatureDataset and

        '   FeatureClass.

        Set pShpToFeatClsConverter = New FeatureDataConverter

        pShpToFeatClsConverter.ConvertFeatureClass pInFeatureClassName, Nothing, _

                pOutFeatureDSName, pOutFeatureClassName, Nothing, pOutFields, "", 1000, 0

        MsgBox "Convert operation complete!", vbInformation

    End If

    Exit Sub

ErrorHandler:

    MsgBox Err.Description

End Sub