【求助】点图层生成表面，总不成功

楼主^#

更多发布于：2006-08-23 10:56

在用插值方法生成表面时总是得不到要输出的结果，代码如下（delphi），其中蓝色变量为输出变量，但一直为空，请高手帮忙？
procedure TfrmMain.Button3Click(Sender: TObject); var pWFact: IWorkspaceFactory; pWks: IWorkspace; pFWks: IFeatureWorkspace; pInFCls: IFeatureClass; pFClsDes: IFeatureClassDescriptor; pRaster: IRaster; pRLayer: IRasterLayer; pLayer: ILayer; pInGeoDataset, pOutGeoDataset: IGeoDataset; pRadius: IRasterRadius; pInterpolatOp: IInterpolationOp; pEnv: IRasterAnalysisEnvironment; pPointNum: integer; pMaxDistance, pBarrier, CellSize: OleVariant;
 pFieldName, pFieldName1: widestring; i: integer; begin pPointNum := 12; pMaxDistance := 0; CellSize := 3; pFieldName := 'GJZ';
 pWFact := CoShapeFileWorkspaceFactory.Create as IWorkspaceFactory; pWFact.OpenFromFile('D:\Projects\ZLS\temp', 0, pWks); pFWks := pWks as IFeatureWorkspace; pFWks.OpenFeatureClass('BasicPoint.shp', pInFCls);
 //set fieldname pFClsDes := CoFeatureClassDescriptor.create as IFeatureClassDescriptor; pFClsDes.Create(pInFCls, nil, pFieldName); pInGeoDataset := pFClsDes as IGeoDataset; pFClsDes.Get_FieldName(pFieldName1);
 //Set cell size for output raster. pInterpolatOp := CoRasterInterpolationOp.Create as IInterpolationOp; pEnv := pInterpolatOp as IRasterAnalysisEnvironment; pEnv.SetCellSize(esriRasterEnvValue, CellSize);
 //Set Radius pRadius := CoRasterRadius.Create as IRasterRadius; pRadius.SetVariable(pPointNum, pMaxDistance);
 //Create GeoDataset for Raster pInterpolatop.Krige(pInGeoDataset, esriGeoAnalysisCircularSemiVariogram, pRadius, true, pBarrier, pOutGeoDataset);
 //Add the Raster to map pRaster := pOutGeoDataset as IRaster; pRLayer := CoRasterLayer.Create as IRasterLayer; pRLayer.CreateFromRaster(pRaster); pLayer := pRLayer as ILayer; pMap.AddLayer(pLayer, 0);
end;

[此贴子已经被作者于2006-8-23 10:59:53编辑过]

喜欢0 评分0

举报回复

cftao2008

路人甲

注册日期2005-03-09
发帖数141
QQ
铜币568枚
威望0点
贡献值0点
银元0个

加关注写私信

1楼^#

发布于：2006-08-23 12:44

下面是转换后的VB代码，错误内容是“there are no spatial analyst or 3Danalyst licenses currently available” ，请问怎样获得这两个扩展模块儿的授权？急。。。 
Dim pWFact As IWorkspaceFactory Dim pWks As IWorkspace Dim pFWks As IFeatureWorkspace Dim pInFCls As IFeatureClass Dim pFClsDes As IFeatureClassDescriptor Dim pRaster As IRaster Dim pRLayer As IRasterLayer Dim pInGeoDataset, pOutGeoDataset As IGeoDataset Dim pRadius As IRasterRadius Dim pInterpolatOp As IInterpolationOp Dim pEnv As IRasterAnalysisEnvironment Dim pPointNum As Integer Dim pMaxDistance, pBarrier, CellSize As Variant Dim pFieldName As String pPointNum = 12 pMaxDistance = 0 CellSize = 3 pFieldName = "GJZ" 'Open shape file Set pWFact = New ShapefileWorkspaceFactory Set pWks = pWFact.OpenFromFile("D:\Projects\ZLS_VB\temp", 0) Set pFWks = pWks Set pInFCls = pFWks.OpenFeatureClass("BasicPoint") 'Set filed Name Set pFClsDes = New FeatureClassDescriptor pFClsDes.Create pInFCls, Nothing, pFieldName 'Create RasterInterpolationOp object Set pInterpolatOp = New RasterInterpolationOp 'Set cell size for output raster Set pEnv = pInterpolatOp pEnv.SetCellSize esriRasterEnvValue, CellSize 'Create raset radius using variable distance Set pRadius = New RasterRadius pRadius.SetVariable 12 Set pRaster = pInterpolatOp.Krige(pFClsDes, esriGeoAnalysisCircularSemiVariogram, pRadius, True) Set pRLayer = New RasterLayer pRLayer.CreateFromRaster pRaster m_Map.AddLayer prlaeyr, 0

[此贴子已经被作者于2006-8-23 12:45:13编辑过]

举报回复喜欢评分

hsghxm

路人甲

注册日期2004-10-27
发帖数517
QQ4604052
铜币1枚
威望0点
贡献值0点
银元0个

加关注写私信

2楼^#

发布于：2006-08-24 09:13

MSN:hsghxm@163.com QQ:4604052 (很少用,最好别加) Email:hsghxm@163.com 我的BOLG：http://www.gisempire.com/blog/user1/864/index.htm

举报回复喜欢评分

您需要登录后才可以回帖，登录或者注册

【求助】点图层生成表面，总不成功

最新喜欢：