arcengine栅格数据使用总结.doc

arcengine栅格数据使用总结 两个星期以来一直与栅格数据打交道，对ao的栅格部分应该有了一定的理解，下面是自己的一点体会，希望高手指教：）1、栅格数据的存储类型栅格数据一般可以存储为esri grid（由一系列文件组成），tiff格式（包括一个tif文件和一个aux文件），imagine image格式 在ae中一般调用isaveas接口来保存栅格数据2、栅格数据集和栅格编目的区别一个栅格数据集由一个或者多个波段（rasterband）的数据组成，一个波段就是一个数据矩阵。对于格网数据（dem数据）和单波段的影像数据，表现为仅仅只有一个波段数据的栅格数据集，而对于多光谱影像数据则表现为具有多个波段的栅格数据集栅格编目（rastercatalog）用于显示某个研究区域内各种相邻的栅格数据，这些相邻的栅格数据没有经过拼接处理合成一副大的影像图3、irasterworkspaceex与irasterworkspace ,irsterworkspace2的区别1).irasteworkspaceex接口主要是用来读取geodatabase中的栅格数据集和栅格编目2) . irasterworkspace ,irsterworkspace2主要是用来读取以文件格式存储在本地的栅格数据4、加载栅格数据（以存储在本地的栅格数据文件为例）1.直接用irasterlayer接口打开一个栅格文件并加载到地图控件irasterlayer rasterlayer = new rasterlayerclass();rasterlayer.createfromfilepath(filename); / filename指存本地的栅格文件路径axmapcontrol1.addlayer(rasterlayer, 0);2. 用irasterdataset接口打开一个栅格数据集iworkspacefactory workspacefactory = new rasterworkspacefactory();iworkspace workspace;workspace = workspacefactory.openfromfile(inpath, 0); /inpath栅格数据存储路径if (workspace = null)console.writeline(could not open the workspace.);return;irasterworkspace rastwork = (irasterworkspace)workspace;irasterdataset rastdataset;rastdataset= rastwork.openrasterdataset(inname);/inname栅格文件名if (rastdataset = null)console.writeline(could not open the raster dataset.);return;5、如何读取栅格数据的属性和遍历栅格数据栅格数据的属性包括栅格大小，行数，列数，投影信息，栅格范围等等，见下面代码（假设当前加载的栅格文件栅格值存储方式为：ushort类型）irasterprops rasterprops = (irasterprops)clipraster;int dheight = rasterprops.height;/当前栅格数据集的行数int dwidth = rasterprops.width; /当前栅格数据集的列数double dx = rasterprops.meancellsize().x; /栅格的宽度double dy = rasterprops.meancellsize().y; /栅格的高度ienvelope extent=rasterprops.extent; /当前栅格数据集的范围rstpixeltype pixeltype=rasterprops.pixeltype; /当前栅格像素类型ipnt pntsize = new pntclass();pntsize.setcoords(dx, dy);ipixelblock pixelblock = clipraster.createpixelblock(pntsize);ipnt pnt = new pntclass();for (int i = 0; i dheight; i+)for (int j = 0; j dwidth; j+)pnt.setcoords(i, j);clipraster.read(pnt, pixelblock);if (pixelblock != null)object obj = pixelblock.getval(0, 0, 0);messagebox.show( convert.touint32(obj).tostring();6、如何提取指定的范围的栅格数据提取指定范围内的栅格数据通常用两种方法irasterlayerexport(esricarto), iextractionop, iextractionop2 (url=esrispatialanalyst/url)，irasterlayerexport接口提供的栅格数据提取功能有限，只能以矩形范围作为提取范围，而iextractionop接口提供了多边形，圆，属性，矩形等几种形式作为提取栅格数据.1).irasterlayerexport接口irasterlayerexport rlayerexport = new rasterlayerexportclass();rlayerexport.rasterlayer = rasterlayer;/ rasterlayer指当前加载的栅格图层rlayerexport.extent = clipextent;/clipextent指提取栅格数据的范围if (prospatialref != null)rlayerexport.spatialreference = prospatialref;/ prospatialref当前栅格数据的投影信息iworkspacefactory pwf = new rasterworkspacefactoryclass();tryiworkspace prasterworkspace = pwf.openfromfile(_folder, 0);/ _folder指栅格文件保存路径irasterdataset outgeodataset = rlayerexport.export(prasterworkspace, code, strrastertype);/调用isaveas接口将导出的数据集保存.catch(exception ex)throw new argumention(ex.message);2iextractionop接口(调用此接口前，应该先检查空间许可)iextractionop extraction = new rasterextractionopclass();tryigeodataset geodataset = extraction.rectangle(igeodataset)clipraster, clipextent, true);iraster raster = geodataset as iraster;if (raster != null)iworkspacefactory wf = new rasterworkspacefactoryclass();iworkspace rasterworkspace = wf.openfromfile(_folder, 0);isaveas saveas = (isaveas)raster;saveas.saveas(“result.tif”, rasterworkspace, tiff);catch (exception ex)messagebox.show(ex.message);7栅格数据重采样栅格数据的重采样主要基于三种方法：最邻近采样（nearest），双线性ilinear）和三次卷积采样（cubic）。（1）.最邻近采样：它用输入栅格数据中最临近栅格值作为输出值。因此，在重采样后的输出栅格中的每个栅格值, 都是输入栅格数据中真实存在而未加任何改变的值。这种方法简单易用，计算量小，重采样的速度最快。（2）.双线性采样：此重采样法取待采样点（x，y）点周围四个邻点，在y方向（或x方向）内插两次，再在x方向（或y方向）内插一次，得到（x，y）点的栅格值。（3）.三次卷积采样：这是进一步提高内插精度的一种方法。它的基本思想是增加邻点来获得最佳插值函数。取待计算点周围相邻的16个点，与双线性采样类似，可先在某一方向上内插，如先在x方向上，每四个值依次内插四次，再根据四次的计算结果在y方上内插，最终得到内插结果代码示例：采用双线性采样irastergeometryproc rastergeometryproc = new rastergeometryprocclass();rastergeometryproc.resample(rstresamplingtypes.rsp_cubicconvolution, newcellsize, clipraster);tag标签: 栅格数据 arcengine标签：栅格 重分类 分类：ae二次开发 public static irasterlayer setviewshedrenderer(iraster pinraster,string sfield,string spath) irasterdescriptor prd = new rasterdescriptorclass(); prd.create(pinraster, new queryfilterclass(), sfield); ireclassop preclassop = new rasterreclassopclass(); igeodataset pgeodataset=pinraster as igeodataset; irasteranalysisenvironment penv = preclassop as irasteranalysisenvironment; iworkspacefactory pwsf=new rasterworkspacefactoryclass(); iworkspace pws = pwsf.openfromfile(spath, 0); penv.outworkspace = pws; object objsnap = null; object objextent = pgeodataset.extent; penv.setextent(esrirasterenvsettingenum.esrirasterenvvalue, ref objextent, ref objsnap); penv.outspatialreference = pgeodataset.spatialreference; irasterlayer prlayer = new rasterlayerclass(); irasterbandcollection prsbandcol = pgeodataset as irasterbandcollection; /定义波段集irasterband prasterband = prsbandcol.item(0); prasterband.computestatsandhist(); irasterstatistics prasterstatistic = prasterband.statistics; /获取像元统计信息double dmaxvalue = prasterstatistic.maximum ; double dminvalue = prasterstatistic.minimum ; inumberremap pnumremap = new numberremapclass(); /定义inumberremap设定阈值pnumremap.maprange(dminvalue, 0, 0); /设置值区间，输出值pnumremap.maprange(0, dmaxvalue, 1); iremap premap = pnumremap as iremap; /转换成iremapiraster poutraster = preclassop.reclassbyremap(pgeodataset, premap, false) as iraster ; /调用reclassbyremap方法prlayer.createfromraster(poutraster); return prlayer; 栅格图层和矢量图层的属性表浏览 if (plyr is ifeaturelayer) datatable ptable = new datatable(); ifeaturelayer pfealyr = plyr as ifeaturelayer; ifeatureclass pfcls = pfealyr.featureclass; string shape = ; if (pfcls.shapetype = esrigeometrytype.esrigeometrypoint) shape = point; else if (pfcls.shapetype = esrigeometrytype.esrigeometrypolyline) shape = polyline; else if (pfcls.shapetype = esrigeometrytype.esrigeometrypolygon) shape = polygon; for (int i = 0; i pfcls.fields.fieldcount; i+) ptable.columns.add(pfcls.fields.get_field(i).name); ifeaturecursor pcursor = pfcls.search(null, false); int ishape = pfcls.fields.findfield(shape); ifeature pfea = pcursor.nextfeature(); while (pfea != null) datarow prow = ptable.newrow(); for (int i = 0; i pfcls.fields.fieldcount; i+) if (i = ishape) prow = shape; continue; prow = pfea.get_value(i).tostring(); ptable.rows.add(prow); pfea = pcursor.nextfeature(); datagridview1.datasource = ptable; else if (plyr is irasterlayer) irasterlayer prlyr = plyr as irasterlayer; iraster praster = prlyr.raster; irasterprops pprop = praster as irasterprops; pprop.pixeltype = rstpixeltype.pt_long; if (pprop.pixeltype = rstpixeltype.pt_long) irasterbandcollection pbcol = praster as irasterbandcollection; irasterband pband = pbcol.item(0); itable prtable = pband.attributetable; datatable ptable = new datatable(); for (int i = 0; i prtable.fields.fieldcount; i+) ptable.columns.add(prtable.fields.get_field(i).name); icursor pcursor= prtable.search(null, false); irow prrow= pcursor.nextrow(); while (prrow != null) datarow prow = ptable.newrow(); for (int i =0 ;iprrow .fields .fieldcount ;i+) prow = prrow.get_value(i).tostring () ; ptable.rows.add(prow); prrow = pcursor.nextrow(); datagridview1.datasource = ptable; irasterworkspace2 irasterdataset createrasterdataset c# public irasterdataset createfilerasterdataset(string directoryname, string filename) / this function creates a new img file in the given workspace / and then assigns pixel values try irasterdataset rasterdataset = null; ipoint originpoint = new pointclass(); originpoint.putcoords(0, 0); / create the dataset irasterworkspace2 rasterworkspace2 = null; rasterworkspace2 = createrasterworkspace(directoryname); rasterdataset = rasterworkspace2.createrasterdataset(filename, imagine image, originpoint, 200, 100, 1, 1, 1, rstpixeltype.pt_uchar, new unknowncoordinatesystemclass(), true); irawpixels rawpixels = null; ipixelblock3 pixelblock3 = null; ipnt pixelblockorigin = null; ipnt pixelblocksize = null; irasterbandcollection rasterbandcollection; irasterprops rasterprops; / qi for irawpixels and irasterprops rasterbandcollection = (irasterbandcollection)rasterdataset; rawpixels = (irawpixels)rasterbandcollection.item(0); rasterprops = (irasterprops)rawpixels; / create pixelblock pixelblockorigin = new dblpntclass(); pixelblockorigin.setcoords(0, 0); pixelblocksize = new dblpntclass(); pixelblocksize.setcoords(rasterprops.width, rasterprops.height); pixelblock3 = (ipixelblock3)rawpixels.createpixelblock(pixelblocksize); / read pixelblock rawpixels.read(pixelblockorigin, (ipixelblock)pixelblock3); / get pixeldata array system.object, pixeldata; pixeldata = (system.object,)pixelblock3.get_pixeldatabyref(0); / loop through all the pixels and assign value for (int i = 0; i rasterprops.width; i+) for (int j = 0; j rasterprops.height; j+) pixeldatai, j = (i * j) % 255; / write the pixeldata back system.object cachepointer; cachepointer = rawpixels.acquirecache(); rawpixels.write(pixelblockorigin, (ipixelblock)pixelblock3); rawpixels.returncache(cachepointer); / return raster dataset return rasterdataset; catch (exception ex) system.diagnostics.debug.writeline(ex.message); return null; public irasterworkspace2 createrasterworkspace(string pathname) / create rasterworkspace iworkspacefactory workspacefactory = new rasterworkspacefactoryclass(); return workspacefactory.openfromfile(pathname, 0) as irasterworkspace2; #4楼楼主 2009-02-09 15:38 | 尤文之鹤 public irasterdataset tin2raster(string tempbathytin,string geopath, string gridname) string tinfolder = system.io.path.getdirectoryname(tempbathytin); string tinname = system.io.path.getfilename(tempbathytin); irasterdataset rasterdataset = new rasterdatasetclass(); try string rasterpath = system.io.path.getdirectoryname(geopath); iworkspacefactory tinwf = new tinworkspacefactory(); itinworkspace tinwk = tinwf.openfromfile(tinfolder,0)as itinworkspace; itinadvanced2 tinad = tinwk.opentin(tinname) as itinadvanced2; ienvelope extent = tinad.extent; ipoint origin = extent.lowerleft; origin.x = origin.x - (5 * 0.5); origin.y = origin.y - (5 * 0.5); int ncol = (int)math.round(extent.width / 5) + 1; int nrow = (int)math.round(extent.height / 5) +1; ispatialreference2 spatialref = (ispatialreference2)extent.spatialreference; iworkspacefactory rasterwf = new rasterworkspacefactoryclass(); irasterworkspace2 workspace = (irasterworkspace2)rasterwf.openfromfile(rasterpath,0); rasterdataset = workspace.createrasterdataset(gridname, grid, origin,ncol,nrow,5,5,1,esri.arcgis.geodatabase.rstpixeltype.pt_float, spatialref,true); irasterbandcollection bandcoll = (irasterbandcollection) rasterdataset; irasterband rasterband = bandcoll.item(0); irawpixels rawpixels = (irawpixels)rasterband; ipnt blocksize = new dblpntclass(); blocksize.x = ncol; blocksize.y = nrow; ipixelblock3 pixelblock = (ipixelblock3)rawpixels.createpixelblock(blocksize); itinsurface tinsurface = (itinsurface)tinad; irasterprops rasterprops = (irasterprops)rawpixels; object nodatafloat; /long nodataint; object val = pixelblock.get_pixeldatabyref(0); messagebox.show(val.tostring(); double cellsize = 5; origin.x = origin.x + (5 * 0.5); origin.y = origin.y + (5 * nrow) - (5 * 0.5); nodatafloat = convert.todouble(rasterprops.nodatavalue.tostring(); tinsurface.querypixelblock(origin.x,origin.y,cellsize,cellsize,esrirasterizationtype.esrielevationasraster,nodatafloat,val); ipnt offset = new dblpntclass(); offset.x = 0; offset.y = 0; rawpixels.write(offset,pixelblock as ipixelblock); catch(exception ex) messagebox.show(ex.tostring(); return rasterdataset; ilayer player = paxmapcontrol.get_layer(this.cmblayer.selectedindex); irasterlayer pr