水文分析模块在VC中的调用问题

<P>我在VC中调用IHydrologyOp接口,实现小流域的水文分析,进行河网提取,编译没什么问题,就是得不到结果,代码如下,请高手指教,欢迎大家讨论:</P>
<P>HRESULT hr;<BR> IRasterLayerPtr pRasterLy(CLSID_RasterLayer);<BR> hr=pRasterLy->CreateFromFilePath(_bstr_t(m_inputpath));<BR> IRasterPtr ipRaster;<BR>    pRasterLy->get_Raster(;ipRaster);<BR>    IWorkspaceFactoryPtr wf(CLSID_ShapefileWorkspaceFactory);<BR> IWorkspacePtr ipWorkspace;<BR>    hr = wf->OpenFromFile((CComBSTR)m_OutputPath, 0, ;ipWorkspace);<BR> if (FAILED(hr) || ipWorkspace == 0)<BR> {<BR>  return ;<BR> }<BR>    IRasterAnalysisEnvironmentPtr ipRastAnalEnv(CLSID_RasterHydrologyOp);<BR>    ipRastAnalEnv->putref_OutWorkspace(ipWorkspace);<BR>    IHydrologyOpPtr ipHydroOp(ipRastAnalEnv);<BR>    IGeoDatasetPtr surfaceRasterDs(ipRaster);<BR>    IRasterDatasetPtr surfaceRaster(surfaceRasterDs);<BR>    ISpatialReferencePtr ipSpaRef;<BR>    surfaceRasterDs->get_SpatialReference(;ipSpaRef);<BR>    ipRastAnalEnv->putref_OutSpatialReference(ipSpaRef);<BR>    IGeoDatasetPtr directionResult;<BR>    hr = ipHydroOp->FlowDirection(surfaceRasterDs, FALSE, FALSE, ;directionResult);</P>

<P>难道没人搞过这方面的开发吗？我怀疑我组件是不是有问题，但在ArcGis8.3该功能可以用</P>