开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > GIS编程 > 查看源码
DlgCreateRasterIDW.cpp
资源名称：GISCode.rar　[点击查看]
上传用户：juying163
上传日期：2014-09-24
资源大小：5867k
文件大小：9k
源码类别：

GIS编程

开发平台：

Visual C++

DlgCreateRasterIDW.cpp：源码内容
  1. // DlgCreateRasterIDW.cpp : implementation file
  2. //
  3. #include "stdafx.h"
  4. #include "3D0214.h"
  5. #include "DlgCreateRasterIDW.h"
  7. #ifdef _DEBUG
  8. #define new DEBUG_NEW
  9. #undef THIS_FILE
  10. static char THIS_FILE[] = __FILE__;
  11. #endif
  13. /////////////////////////////////////////////////////////////////////////////
  14. // CDlgCreateRasterIDW dialog
  17. CDlgCreateRasterIDW::CDlgCreateRasterIDW(CWnd* pParent /*=NULL*/)
  18. : CDialog(CDlgCreateRasterIDW::IDD, pParent)
  19. {
  20. //{{AFX_DATA_INIT(CDlgCreateRasterIDW)
  21. m_edit_IDWInputCellsize = 0.0;
  22. m_edit_IDWOutputRaster = _T("");
  23. m_edit_IDWPointInput = _T("");
  24. m_edit_IDWPower = 0.0;
  25. //}}AFX_DATA_INIT
  26. }
  29. void CDlgCreateRasterIDW::DoDataExchange(CDataExchange* pDX)
  30. {
  31. CDialog::DoDataExchange(pDX);
  32. //{{AFX_DATA_MAP(CDlgCreateRasterIDW)
  33. DDX_Control(pDX, IDC_EDIT_IDW_OUTPUT_RASTER, m_Control_IDWOutputRaster);
  34. DDX_Control(pDX, IDC_STATIC_Set2, m_static_settings2);
  35. DDX_Control(pDX, IDC_STATIC_Set1, m_static_settings1);
  36. DDX_Control(pDX, IDC_EDIT_IDW_RADIUSSETTINGS_2, m_edit_settings2);
  37. DDX_Control(pDX, IDC_EDIT_IDW_RADIUSSETTINGS_1, m_edit_settings1);
  38. DDX_Control(pDX, IDC_COMBO_IDW_SEARCH_RADIUS_TYPE, m_combo_IDWsearchRadiusType);
  39. DDX_Control(pDX, IDC_COMBO_IDW_ZVALUE_FIELD, m_combo_IDWZvalueField);
  40. DDX_Text(pDX, IDC_EDIT_IDW_INPUT_CELL_SIZE, m_edit_IDWInputCellsize);
  41. DDX_Text(pDX, IDC_EDIT_IDW_OUTPUT_RASTER, m_edit_IDWOutputRaster);
  42. DDX_Text(pDX, IDC_EDIT_IDW_POINT_INPUT, m_edit_IDWPointInput);
  43. DDX_Text(pDX, IDC_EDIT_IDW_WIGHT, m_edit_IDWPower);
  44. //}}AFX_DATA_MAP
  45. }
  48. BEGIN_MESSAGE_MAP(CDlgCreateRasterIDW, CDialog)
  49. //{{AFX_MSG_MAP(CDlgCreateRasterIDW)
  50. ON_CBN_SELCHANGE(IDC_COMBO_IDW_SEARCH_RADIUS_TYPE, OnSelchangeComboIdwSearchRadiusType)
  51. ON_BN_CLICKED(IDC_BUTTON_IDW_ADDPOINTLAYER, OnButtonIdwAddpointlayer)
  52. ON_BN_CLICKED(IDC_BUTTON_SELSAVERASTERPATH, OnButtonSelsaverasterpath)
  53. //}}AFX_MSG_MAP
  54. END_MESSAGE_MAP()
  56. /////////////////////////////////////////////////////////////////////////////
  57. // CDlgCreateRasterIDW message handlers
  59. BOOL CDlgCreateRasterIDW::OnInitDialog() 
  60. {
  61. CDialog::OnInitDialog();
  63. // TODO: Add extra initialization here
  64. //将搜索范围默认选项设置为"可变的
  65. m_combo_IDWsearchRadiusType.SetCurSel(0);
  66. m_edit_settings1.SetWindowText("12");
  67. m_selIndex=0;
  68. // m_combo_IDWZvalueField. AddString("可变的");
  70. return TRUE;  // return TRUE unless you set the focus to a control
  71.               // EXCEPTION: OCX Property Pages should return FALSE
  72. }
  74. void CDlgCreateRasterIDW::OnSelchangeComboIdwSearchRadiusType() 
  75. {
  76. // TODO: Add your control notification handler code here
  77. if(m_combo_IDWsearchRadiusType.GetCurSel()!=CB_ERR)
  78. {
  79. m_selIndex=m_combo_IDWsearchRadiusType.GetCurSel();
  80. if ( m_selIndex == 0 ) //**当选中"可变的"
  81. {
  82. m_static_settings1.SetWindowText("点的数目:");
  83. m_static_settings2.SetWindowText("最大距离:");
  84. m_edit_settings1.SetWindowText("12");
  85. m_edit_settings2.SetWindowText("");
  86. }
  87. if( m_selIndex ==1) //*当选中"固定的"
  88. {
  89. m_static_settings1.SetWindowText("距离:");
  90. m_static_settings2.SetWindowText("最少点数:");
  91. m_edit_settings1.SetWindowText("");
  92. m_edit_settings2.SetWindowText("0");
  93. }
  94. }
  95. }
  97. void CDlgCreateRasterIDW::OnButtonIdwAddpointlayer() 
  98. {
  99. // TODO: Add your control notification handler code here
  100. ::OleInitialize(NULL);//组件操作初始化
  101. //打开对话框
  102. IGxDialogPtr ipGxDialog(CLSID_GxDialog);
  104.  IGxObjectFilterPtr ipGxObjectFilter(CLSID_GxFilterPointFeatureClasses);   
  105. //注意使用了GxFilterPointFeatureClasses
  106. //可以过滤点特征集，满足功能只用点内插的需要
  107.     HRESULT hr=ipGxDialog->putref_ObjectFilter(ipGxObjectFilter);
  109. //HRESULT hr;
  111. if(FAILED(hr)) return;
  113. hr=ipGxDialog->put_Title(::SysAllocString(L"选择数据"));
  115. if(FAILED(hr)) return;
  116. IEnumGxObjectPtr ipEnumGxObject;
  117. VARIANT_BOOL bResult;
  118. hr=ipGxDialog->DoModalOpen(0,&ipEnumGxObject,&bResult);
  119. if(FAILED(hr)) return;
  120. if(bResult==VARIANT_FALSE) return;
  121. hr=ipEnumGxObject->Reset();///
  122. if(FAILED(hr)) return;
  123. IGxObjectPtr ipGxObject;
  125. while(!ipEnumGxObject->Next( &ipGxObject))
  126. {
  127. IGxDatasetPtr ipGxDataset( ipGxObject);
  128. if(ipGxDataset!=0)
  129. {
  130. IDatasetPtr ipDataset;
  131. hr=ipGxDataset->get_Dataset( &ipDataset);
  132. if(FAILED(hr)) return;
  133. IFeatureClassPtr ipFeatureClass(ipDataset);
  134. m_ipfeat=ipFeatureClass;
  135. CComBSTR bstrName ;//获得输人点层的全名
  136. ipGxObject->get_FullName(&bstrName);
  137. CString str(bstrName);
  138. m_edit_IDWPointInput=str;
  139. UpdateData(FALSE);
  140. //将点层ipFeatureClass的所有字段添加到。ombo box控件列表中
  141. m_combo_IDWZvalueField.ResetContent();//先清空
  142. HRESULT hr;
  143. long lFieldCount;
  144. IFieldsPtr pFClassPointFields;
  145. m_ipfeat->get_Fields(&pFClassPointFields);
  146. pFClassPointFields->get_FieldCount(&lFieldCount);
  147. for(int i=0;i<lFieldCount;i++)
  148. {    
  149. CComBSTR pname;    
  150. IFieldPtr pField;   
  151. hr=pFClassPointFields->get_Field(i,&pField);    
  152. if(FAILED(hr )) return;   
  153. hr=pField->get_Name(&pname);    
  154. if( FAILED(hr)) return;   
  155. CString ppname(pname);    
  156. m_combo_IDWZvalueField. AddString(ppname);
  158. }
  160. m_combo_IDWZvalueField.SetCurSel(0);
  161. //初始化显示。omboBox(3)
  162. }
  163. }
  164. }
  166. void CDlgCreateRasterIDW::OnOK() 
  167. {
  168. // TODO: Add extra validation here
  169. UpdateData(TRUE);
  170. if(m_ipfeat!=NULL)
  171. {
  172. IDWOp(m_ipfeat);
  173. }
  174. else
  175. {
  176. AfxMessageBox("没有输人点层!"); return;
  177. }
  178. CDialog::OnOK();
  179. }
  180. void CDlgCreateRasterIDW::IDWOp(IFeatureClass * ipfeatureclass) 
  181. {
  182. IFeatureClassDescriptorPtr pFDescr(CLSID_FeatureClassDescriptor);
  183. //获得选中字段的字符串
  184. int index=m_combo_IDWZvalueField.GetCurSel();
  185. CString sFieldName;
  186. m_combo_IDWZvalueField.GetLBText(index,sFieldName);
  187. CComBSTR s=(CComBSTR)sFieldName;
  188. HRESULT hr;
  189. hr=pFDescr->Create(ipfeatureclass,NULL,s);
  190. //创建RasterlnterpolationOp对象
  192. IInterpolationOpPtr pIntOp(CLSID_RasterInterpolationOp);
  193. //设置输出RASTER网格大小，默认的RASTER输出空间范围和输人数据相同
  195. IRasterAnalysisEnvironmentPtr pEnv(pIntOp);
  196. VARIANT cs;
  197. cs.vt=VT_R8;
  198. cs .dblVal=m_edit_IDWInputCellsize; //对话框输入cell 尺寸
  199. hr = pEnv->SetCellSize(esriRasterEnvValue,&cs);
  201. //搜索范围设置
  202. IRasterRadiusPtr pRadius(CLSID_RasterRadius);
  203.   if(m_selIndex==0)
  204. {
  205. CString strl;
  206. m_edit_settings1.GetWindowText(strl);
  207. long count=atoi(strl); CString str2;
  208. m_edit_settings2.GetWindowText(str2);
  209. double distance=atof ( str2);
  211. VARIANT dis;
  212. dis.vt=VT_R8;
  213. dis. dblVal=distance;    
  214. hr=pRadius->SetVariable(count,&dis);
  215. }
  216.   if(m_selIndex==1)
  217. {     
  218. CString str1,str2;    
  219. m_edit_settings1.GetWindowText(str1);    
  220. double distance=atof( str1);    
  221. m_edit_settings2.GetWindowText(str2);
  222.      long count=atoi(str2);    
  223. VARIANT cou;
  224. cou. vt=VT_I4;
  225. cou.lVal=count;
  226. hr=pRadius->SetFixed(distance,&cou); 
  227. }
  230. //使用FeatureClassDescriptor作为IInterpolationOp的输人对象进行内插操作
  232. IGeoDatasetPtr pFDescr2Geo(pFDescr);
  233. IGeoDatasetPtr pOutRaster2Geo;
  234. hr=pIntOp->IDW(pFDescr2Geo,m_edit_IDWPower,pRadius,NULL,&pOutRaster2Geo);
  235. if(FAILED("hr"))
  236. {
  237. AfxMessageBox ("内插操作失败!");
  238. return;
  239. }    
  240. IRasterPtr pOutRaster(pOutRaster2Geo);    
  241. IRasterBandCollectionPtr pRasBandC(pOutRaster);    
  242. //输出RASTER路径、名称及其格式设置    
  243. IWorkspaceFactoryPtr pWSF(CLSID_RasterWorkspaceFactory);    
  244. IWorkspacePtr pWS; 
  245. hr= pWSF->OpenFromFile((CComBSTR)m_edit_IDWOutputRaster,NULL,&pWS); 
  246.      if(FAILED(hr)) AfxMessageBox("failed");
  247. //RASTER格式可以为"GRID" , " TIFF"或者"IMAGINE Image"    
  248. IDatasetPtr prasterdataset;    
  249. hr=pRasBandC->SaveAs(CComBSTR("Raster"),pWS,CComBSTR("GRID"),&prasterdataset);    
  250. if(FAILED( hr))
  251.      {        
  252. AfxMessageBox ( " DEM保存失败!");    
  253. }
  255. }
  256. void CDlgCreateRasterIDW::OnButtonSelsaverasterpath() 
  257. {
  258. // TODO: Add your control notification handler code here
  259.     CString strFilter="File Folder(*.File Folder)|*.File Folder||";//设置打开文件的过滤条件
  260. CString strTitle="选择保存路径";
  261.         DWORD MAXFILE = 2412; // allocate enough memory space
  262.      CString strPath,strInit,strTemp;
  263.         CFileDialog pDlg(FALSE,0,0,OFN_ALLOWMULTISELECT| OFN_EXPLORER | OFN_LONGNAMES|
  264.                         OFN_FILEMUSTEXIST |OFN_HIDEREADONLY,strFilter, NULL);
  266.         //设置对话框               
  267.         pDlg.m_ofn.nMaxFile = MAXFILE; // set the buffer size        
  268.         pDlg.m_ofn.lpstrTitle=strTitle;
  269.         pDlg.m_ofn.lpstrInitialDir=strInit;
  270. char* buf = new char[MAXFILE];
  271. pDlg.m_ofn.lpstrFile = buf;
  272. pDlg.m_ofn.lpstrFile[0] = NULL;
  273. if(pDlg.DoModal()== IDOK)
  274. {
  275. //POSITION pos = pDlg.GetStartPosition();
  276. strPath=pDlg.GetPathName();
  277. //strPath=pDlg.GetFolderPath();
  278. m_Control_IDWOutputRaster.SetWindowText(strPath);
  279. AfxMessageBox(strPath,MB_YESNO);
  280. //AfxMessageBox("Do you really want to exit?",MB_YESNO);
  281. }
  282. }