Господа! Кто-нить работал с Spatial Analy

Question

как вааще мона программно интерполировать точечную тему?


esriCore::IQueryFilterPtr pFilterLine(esriCore::CLSID_QueryFilter);
     pFilterLine->put_SubFields(_bstr_t("PV_1"));
     esriCore::IFeatureClassPtr pFeatureClass;
     pFeatureLayer->get_FeatureClass(&pFeatureClass);
     

     esriCore::IFeatureClassDescriptorPtr pFeatureClassDescriptor(esriCore::CLSID_FeatureClassDescriptor);
     BSTR field1 = _bstr_t("PV_1");
     pFeatureClassDescriptor->Create(pFeatureClass,pFilterLine,field1);

     esriSpatialShared::IInterpolationOpPtr pInterpolationOp(esriSpatialShared::CLSID_RasterInterpolationOp);
     esriCore::IRasterAnalysisEnvironmentPtr pEnv;
     pEnv = pInterpolationOp;
     pEnv->SetCellSize(esriCore::esriRasterEnvValue,&variant_t(edSize));
     pEnv->SetExtent(esriCore::esriRasterEnvValue);



     esriSpatialShared::IRasterRadiusPtr pRadius(esriSpatialShared::CLSID_RasterRadius);
     //pRadius->SetVariable(atol(aPoints));
     pRadius->SetFixed (1000);


     esriCore::IGeoDatasetPtr pGeoDatasetI;
     pGeoDatasetI = pFeatureClass;
     

     esriCore::IRasterPtr pOutRaster(esriCore::CLSID_Raster);

     esriCore::IGeoDatasetPtr pGeoDatasetO;
     //pGeoDatasetO = pOutRaster;

     pInterpolationOp->Trend(pGeoDatasetI,esriSpatialShared::esriGeoAnalysisLinearTrend,10,&pGeoDatasetO);

     pOutRaster = pGeoDatasetO;
     

Всегда возврящяется пустой geodataset

Answer 1

Интересует вААбще рабочий код на сях по работе с spatial analyst