geofeatureclassm1.cpp

arcgis geodatabase Access 数据库创建

// GeoFeatureClassM.cpp: implementation of the CGeoFeatureClassM class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "TestGeoDatabase.h"
#include "GeoFeatureClassM.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

CGeoFeatureClassM::CGeoFeatureClassM()
: INT_MAX_VALUE(2147473648.0)
{
	//pFeatureWorkspace As IFeatureWorkspace, className As String, GeomType As esriGeometryType, pFeatureClass As IFeatureClass, progress As Long
	
	m_bSetZDomain = FALSE;
	m_bSetMDomain = FALSE;
	m_dZMin   = 0.0;
	m_dZMax   = 0.0;
	m_dMMin   = 0.0;
	m_dMMax   = 1000.0;
	m_dMinX   = 0.0;
	m_dMaxX   = 1000.0;
	m_dMinY   = 0.0;
	m_dMaxY   = 1000.0;
	
}

CGeoFeatureClassM::~CGeoFeatureClassM()
{
	
}

BOOL CGeoFeatureClassM::CreateGeoFeatureClass( CString strLayerName,esriFeatureType  FeatureType,IFeatureClassPtr &ipFClass)
{
    BSTR bstrName=_com_util::ConvertStringToBSTR(strLayerName);
	CComBSTR bstrLocalName(bstrName);
    if (bstrLocalName == L"")
    return E_POINTER;

	IUIDPtr pCLSID;

	if (0==pCLSID) {
		HRESULT result =pCLSID.CreateInstance(CLSID_UID);
		if(FALSE(result))
		{
			AfxMessageBox("无法创建IUID接口，新建图层失败。");
			return FALSE;
		}
		this->SetGeometryType(FeatureType,pCLSID);
	}
	
	////////////////////////////////////////////////
	 // establish a fields collection
  IFieldPtr pField;
  IFieldsPtr pFields;
  if (0 == pFields)
  {
    pFields.CreateInstance(CLSID_Fields);
    IFieldsEditPtr ipFieldsEdit(pFields);

    // create the geometry field
    IGeometryDefPtr ipGeomDef(CLSID_GeometryDef);
    IGeometryDefEditPtr ipGeomDefEdit(ipGeomDef);

    // assign the geometry defintion properties.
    ipGeomDefEdit->put_GeometryType(m_GeometryType);
    ipGeomDefEdit->put_GridCount(1);
    ipGeomDefEdit->put_GridSize(0, 10);
    ipGeomDefEdit->put_AvgNumPoints(2);
    ipGeomDefEdit->put_HasM(VARIANT_FALSE);
    ipGeomDefEdit->put_HasM(VARIANT_FALSE);
	// assign the spatial reference
    pField.CreateInstance(CLSID_Field);
    IFieldEditPtr ipFieldEdit(pField);
	////////////////////////////////////////////////////
    ipFieldEdit->put_Name(CComBSTR(L"shape"));
    ipFieldEdit->put_AliasName(CComBSTR(L"geometry"));
    ipFieldEdit->put_Type(esriFieldTypeGeometry);

//	ISpatialReferencePtr ipUnkSpatial(CLSID_UnknownCoordinateSystem);
//     ipUnkSpatial->SetFalseOriginAndUnits(0,0,100000);
//     ipGeomDefEdit->putref_SpatialReference(ipUnkSpatial);
     ipFieldEdit->putref_GeometryDef(ipGeomDef);
	 ipFieldsEdit->AddField(pField);

    //create the object id field
    pField.CreateInstance(CLSID_Field);
    ipFieldEdit = pField;
    ipFieldEdit->put_Name(CComBSTR(L"OBJECTID"));
    ipFieldEdit->put_AliasName(CComBSTR(L"object identifier"));
    ipFieldEdit->put_Type(esriFieldTypeOID);
	
    ipFieldsEdit->AddField(pField);
	///////////////////////////////////////////////

  }
  
  this->CheckFields(pFields,pField);
  
  IUIDPtr pCLSEXT;
  HRESULT hr;
  hr = m_ipFDataset->CreateFeatureClass(bstrName, pFields, m_pCLSID, 0, FeatureType,m_bstrShapeFld, CComBSTR(L""), &ipFClass);
 /////////////////////////////////////////
  return hr;

}

void CGeoFeatureClassM::SetGeometryType(esriFeatureType FeatureType,IUIDPtr pCLSID)
{
	   m_pCLSID=pCLSID;
	   switch (FeatureType)
	   {
       case esriFTSimple:
		   m_GeometryType=esriGeometryPolygon;
		   m_pCLSID->put_Value(CComVariant(L"{52353152-891a-11d0-bec6-00805f7c4268}")); // IFeature's GUID
		   break;
       case esriFTSimpleJunction:
		   m_GeometryType = esriGeometryPolygon;
		   m_pCLSID->put_Value(CComVariant(L"{cee8d6b8-55fe-11d1-ae55-0000f80372b4}")); // ISimpleJunctionFeature's GUID
		   break;
       case esriFTComplexJunction:
		   m_pCLSID->put_Value(CComVariant(L"{df9d71f4-da32-11d1-aeba-0000f80372b4}")); // IComplexJunctionFeature's GUID
		   break;
       case esriFTSimpleEdge:
		   m_GeometryType = esriGeometryPolyline;
		   m_pCLSID->put_Value(CComVariant(L"{e7031c90-55fe-11d1-ae55-0000f80372b4}")); // ISimpleEdgeFeature's GUID
		   break;
       case esriFTComplexEdge:
		   m_GeometryType = esriGeometryPolygon;
		   m_pCLSID->put_Value(CComVariant(L"{A30E8A2A-C50B-11D1-AEA9-0000F80372B4}"));//a30e8a2a-c50b-11d1-aea9-0000f80372b4}")); // IComplexEdgeFeature's GUID
		   break;
       case esriFTAnnotation://注释注解		   
		   m_pCLSID->put_Value(CComVariant(L"{e3676993-c682-11d2-8a2a-006097aff44e}"));//a30e8a2a-c50b-11d1-aea9-0000f80372b4}")); // IComplexEdgeFeature's GUID
		   break;
		   return ;
       default:
		   return ;
	   }
}

void CGeoFeatureClassM::CheckFields(IFieldsEditPtr ipFields,IFieldPtr pField)
{
	esriFieldType fieldType;
	long lNumFields;
	ipFields->get_FieldCount(&lNumFields);
	for (int i = 0; i<lNumFields; ++i)
	{
		ipFields->get_Field(i,&pField);
		pField->get_Type(&fieldType);
		if (esriFieldTypeGeometry == fieldType)
		{
			pField->get_Name(&m_bstrShapeFld);
		}
	}
	
}




void CGeoFeatureClassM::AddExtendField()
{
	
}

BOOL CGeoFeatureClassM::CreateSpatialReference(ISpatialReferencePtr &ipSReference)
{
	// 范围与精度的关系：（最大坐标 - 最小坐标）* 精度 = 2147473648
	// 其中，2147473648 是整数的最大值
	HRESULT result = ipSReference.CreateInstance(CLSID_UnknownCoordinateSystem);
	
	if (FAILED(result)) 
	{
		AfxMessageBox("无法创建空间参考接口。");
		return FALSE;
	}
	double dMaxX = m_dMaxX - m_dMinX;
	double dMaxY = m_dMaxY - m_dMinY;
	double dMax = dMaxX > dMaxY ? dMaxX : dMaxY;
	ipSReference->SetDomain(m_dMinX, m_dMaxX, m_dMinY, m_dMaxY);
	ipSReference->SetFalseOriginAndUnits(m_dMinX, m_dMinY, INT_MAX_VALUE / dMax);
		
	if (m_bSetZDomain) 
	{
		ipSReference->SetZDomain(m_dZMin, m_dZMax);
		ipSReference->SetZFalseOriginAndUnits(m_dZMin, INT_MAX_VALUE / (m_dZMax - m_dZMin));
	}
	
	if (m_bSetMDomain)
	{
		ipSReference->SetMDomain(m_dMMin, m_dMMax);
		ipSReference->SetMFalseOriginAndUnits(m_dMMin, INT_MAX_VALUE / (m_dMMax - m_dMMin));
	}
	return TRUE;
	//IGeometryCollection
	
}
BOOL CGeoFeatureClassM::CreateDataset(CGeoDatabaseM *geo,CString strDatasetName, IFeatureDatasetPtr &ipFDataset)
{

  HRESULT hr;
  CComPtr<ISpatialReference> ipSRef;
  BSTR name=_com_util::ConvertStringToBSTR(strDatasetName);
  if (!ipSRef)
  {
    ipSRef.CoCreateInstance(CLSID_UnknownCoordinateSystem);
    ipSRef->SetDomain(0, 21474.83645, 0, 21474.83645);
    ipSRef->SetFalseOriginAndUnits(0,0,100000);
    ipSRef->SetMDomain(0, 21474.83645);
    ipSRef->SetZDomain(0, 21474.83645);
    ipSRef->SetZFalseOriginAndUnits(0, 100000);
  }
	
  CComPtr<IFeatureWorkspace> ipFeatWksp;
  if (FAILED(hr = geo->m_ipWorkspace->QueryInterface(&ipFeatWksp)))
    return hr;

   hr=ipFeatWksp->CreateFeatureDataset(name, ipSRef, &ipFDataset);
   m_ipFDataset=ipFDataset;

  return hr;
	/*
	ISpatialReferencePtr ipSRef;
	ipSRef.CreateInstance(CLSID_UnknownCoordinateSystem);
	 ipSRef->SetDomain( 0, 21474.83645, 0, 21474.83645);
     ipSRef->SetFalseOriginAndUnits(0,0,100000);
	if (!CreateSpatialReference(ipSRef)) 
		return FALSE;
	BSTR name=_com_util::ConvertStringToBSTR(strDatasetName);
	IFeatureWorkspacePtr ipFeatWksp =geo->m_ipWorkspace; 
	HRESULT hr = ipFeatWksp->CreateFeatureDataset(name, ipSRef, &ipFDataset);
	if(SUCCEEDED(hr)) 
		m_ipFDataset=ipFDataset;
	return (SUCCEEDED(hr) && NULL != ipFDataset); */
}

BOOL CGeoFeatureClassM::CheckParameter(CGeoDatabaseM *geodatabase)
{
	
	ASSERT(NULL !=geodatabase->m_ipWorkspace);
	// ASSERT(NULL != m_ipFields);
	ASSERT(m_dZMin <= m_dZMax);
	ASSERT(m_dMMin <= m_dMMax);
	ASSERT(m_dMinX <= m_dMaxX && m_dMinY <= m_dMaxY);
	
	if (NULL == geodatabase->m_ipWorkspace) return FALSE;
	//if (NULL == m_ipFields) return FALSE;
	if (m_dZMin > m_dZMax) return FALSE;
	if (m_dMMin > m_dMMax) return FALSE;
	if (m_dMinX > m_dMaxX || m_dMinY > m_dMaxY) return FALSE;
	
	return TRUE;
}

HRESULT CGeoFeatureClassM::OpenFeatureClass(BSTR bstrWorkspace, BSTR bstrFeatureClass, IFeatureClass **ppFeatureClass)
{
	
  // Create the workspace name object
  IWorkspaceNamePtr ipShpWorkspaceName(CLSID_WorkspaceName);
  ipShpWorkspaceName->put_PathName(bstrWorkspace);
  ipShpWorkspaceName->put_WorkspaceFactoryProgID(L"esriDataSourcesFile.shapefileworkspacefactory.1");

  // Create the feature class name object
  IDatasetNamePtr ipDatasetName(CLSID_FeatureClassName);
  ipDatasetName->put_Name(bstrFeatureClass);
  ipDatasetName->putref_WorkspaceName(ipShpWorkspaceName);

  // Open the feature class
  INamePtr ipName(ipDatasetName);
  IUnknownPtr ipUnknown;
  ipName->Open(&ipUnknown);
  IFeatureClassPtr ipFeatureClass(ipUnknown);

  *ppFeatureClass = ipFeatureClass;
  if (*ppFeatureClass)
    (*ppFeatureClass)->AddRef();
  else
    return E_FAIL;

  return S_OK;


}

HRESULT CGeoFeatureClassM::LoadFeatures()
{
// Open shape_file where new features will be written to // For simplicity, sample does not contain code to create a new shapefile
  IFeatureClassPtr ipOutFeatureClass;
  OpenFeatureClass(CComBSTR(L"d:\\data\\usa"), CComBSTR(L"test"), &ipOutFeatureClass);
  if (0 == ipOutFeatureClass)
    return E_FAIL;

  IFeatureCursorPtr ipInsertFeatureCursor;
  ipOutFeatureClass->Insert(VARIANT_TRUE, &ipInsertFeatureCursor);

  IFeatureBufferPtr ipInsertFeatureBuffer;
  ipOutFeatureClass->CreateFeatureBuffer(&ipInsertFeatureBuffer);

  // Open shape_file containing the features that will be copied
  IFeatureClassPtr ipInFeatureClass;
  OpenFeatureClass(CComBSTR(L"d:\\data\\usa"), CComBSTR(L"counties"), &ipInFeatureClass);
  if (0 == ipInFeatureClass)
    return E_FAIL;

  CComVariant varVariant;
  IGeometryPtr ipGeom;
  int nNewFeatureCount(0);
  // Loop through all the features in InFeatureClass
  IFeatureCursorPtr ipSearchFeatureCursor;
  ipInFeatureClass->Search(0, VARIANT_TRUE, &ipSearchFeatureCursor);
  IFeaturePtr ipFeature;
  ipSearchFeatureCursor->NextFeature(&ipFeature);
  while (ipFeature)
  {
    // Add the original feature's geometry to the feature buffer
    ipFeature->get_Shape(&ipGeom);
    ipInsertFeatureBuffer->putref_Shape(ipGeom);
    // Add all the original feature's fields to the feature buffer
 //   AddFields(ipInsertFeatureBuffer, ipFeature);
    // Insert the feature into the cursor
    ipInsertFeatureCursor->InsertFeature(ipInsertFeatureBuffer, &varVariant);
    nNewFeatureCount += 1;
    // Flush the feature cursor every 100 features // This is safer because you can write code to handle a flush error // If you don't flush the feature cursor it will automatically flush but // after all of your code executes at which time you have no control
    if (100 == nNewFeatureCount)
    {
      ipInsertFeatureCursor->Flush();
      nNewFeatureCount = 0;
    }
    ipSearchFeatureCursor->NextFeature(&ipFeature);
  }
  ipInsertFeatureCursor->Flush();   // Flush the cursor one last time
  return S_OK;

}