AO开发：将TIn 转换为栅格Grid需要的接口是那些？

楼主^#

更多发布于：2006-05-19 16:04

现在我需要将Tin生成Slope，在ArcMap里面当然可以直接生成了，但是用AO进行二次开发就是找不到需要的接口，希望知道的能指点一下，如果能列出AO里的samoles那跟更好了，呵呵。还有就是对生成的SLope栅格，表数据我查询了一下，有一个字段是“strech value”，另外一个是“slope value”，怎么获取这些点的值呢？希望大家能指点迷津，谢谢了。

喜欢0 评分0

举报回复

gis

管理员

注册日期2003-07-16
发帖数15951
QQ
铜币25345枚
威望15368点
贡献值0点
银元0个

加关注写私信

1楼^#

发布于：2006-05-19 17:34

这个函数可以转,

 ' Supported pixel types limited to float and long because output currently limited to native ESRI Grid ' This routine handles cancel tracking so passed TIN should not have its CancelTracker set. Public Function TinToRaster(ByRef pTin As ESRI.ArcGIS.Geodatabase.ITinAdvanced, ByRef eRastConvType As ESRI.ArcGIS.Geodatabase.esriRasterizationType, ByRef sDir As String, ByRef sName As String, ByRef ePixelType As ESRI.ArcGIS.Geodatabase.rstPixelType, ByRef cellsize As Double, ByRef pExtent As ESRI.ArcGIS.Geometry.IEnvelope, ByRef bPerm As Boolean) As ESRI.ArcGIS.Geodatabase.IRasterDataset
 ' The origin used by CreateRasterDataset is the lower left cell corner. ' The extent passed is that of the TIN's. ' Define the raster origin and number of rows and columns so that the raster ' is of sufficient extent to capture area defined by passed envelope. The cell ' center is located at the origin. Dim pOrigin As ESRI.ArcGIS.Geometry.IPoint pOrigin = pExtent.LowerLeft pOrigin.X = pOrigin.X - (cellsize * 0.5) pOrigin.Y = pOrigin.Y - (cellsize * 0.5)
 Dim nCol, nRow As Integer nCol = System.Math.Round(pExtent.Width / cellsize) + 1 nRow = System.Math.Round(pExtent.Height / cellsize) + 1
 Dim pGDS As ESRI.ArcGIS.Geodatabase.IGeoDataset pGDS = pTin Dim pSR As ESRI.ArcGIS.Geometry.ISpatialReference2 pSR = pGDS.SpatialReference
 Dim pRDS As ESRI.ArcGIS.Geodatabase.IRasterDataset pRDS = dbUtil.CreateRasterSurf(sDir, sName, "GRID", pOrigin, nCol, nRow, cellsize, cellsize, ePixelType, pSR, bPerm)
 System.Windows.Forms.Application.DoEvents()
 Dim pRawPixels As ESRI.ArcGIS.DataSourcesRaster.IRawPixels pRawPixels = dbUtil.GetRawPixels(pRDS, 0)
 Dim pCache As stdole.IUnknown pCache = pRawPixels.AcquireCache
 Dim pTinSurf As ESRI.ArcGIS.Geodatabase.ITinSurface pTinSurf = pTin
 Dim pRasterProps As ESRI.ArcGIS.DataSourcesRaster.IRasterProps pRasterProps = pRawPixels
 Dim nodataFloat As Single Dim nodataInt As Integer
 Dim dZMin As Double
 dZMin = pTin.Extent.ZMin
 Dim vNoData As Object If (ePixelType = ESRI.ArcGIS.Geodatabase.rstPixelType.PT_FLOAT) Then
 vNoData = CSng(dZMin - 1) Else
 vNoData = CInt(dZMin - 1) End If
 ' pRasterProps.NoDataValue = -9999 'vNoData
 Dim pOffset As ESRI.ArcGIS.Geodatabase.IPnt pOffset = New ESRI.ArcGIS.DataSourcesRaster.DblPnt
 ' Set blocksize. Restrict how large it is as not to consume too much memory for ' big output datasets. Dim lMaxBlockX As Integer lMaxBlockX = 2048 If (nCol < lMaxBlockX) Then lMaxBlockX = nCol End If
 Dim lMaxBlockY As Integer lMaxBlockY = 2048 If (nRow < lMaxBlockY) Then lMaxBlockY = nRow End If
 Dim pBlockSize As ESRI.ArcGIS.Geodatabase.IPnt pBlockSize = New ESRI.ArcGIS.DataSourcesRaster.DblPnt pBlockSize.X = lMaxBlockX pBlockSize.Y = lMaxBlockY
 Dim pPixelBlock As ESRI.ArcGIS.DataSourcesRaster.IPixelBlock3 pPixelBlock = pRawPixels.CreatePixelBlock(pBlockSize)
 Dim blockArray As Object blockArray = pPixelBlock.PixelDataByRef(0)
 Dim lBlockCount As Integer lBlockCount = System.Math.Round((nCol / lMaxBlockX) + 0.49) * System.Math.Round((nRow / lMaxBlockY) + 0.49)
 Dim pBlockOrigin As ESRI.ArcGIS.Geometry.IPoint pBlockOrigin = New ESRI.ArcGIS.Geometry.Point
 Dim lColOffset As Integer Dim lRowOffset As Integer
 ' Left to right, top to bottom, iteration of pixel blocks. Dim bReset As Boolean For lRowOffset = 0 To (nRow - 1) Step lMaxBlockY
 For lColOffset = 0 To (nCol - 1) Step lMaxBlockX
 ' See if pixelblock needs to be resized in X for last column chunk. ' RawPixel.Write will clip the pixelblock if it's too big, so the resize ' isn't absolutely necessary, but resizing will eliminate unecessary ' effort for TIN's QueryPixelBlock. If ((nCol - lColOffset) < lMaxBlockX) Then pBlockSize.X = (nCol - lColOffset) pPixelBlock = pRawPixels.CreatePixelBlock(pBlockSize) blockArray = pPixelBlock.PixelDataByRef(0) End If
 ' QueryPixelBlock takes an origin representing the upper left cell center. ' Calculate that cell center's position here. Calculate it based on the ' raster's origin (lower left) and current row/col offset. pBlockOrigin.X = pOrigin.X + (lColOffset * cellsize) + (cellsize * 0.5) pBlockOrigin.Y = pOrigin.Y + ((nRow - lRowOffset) * cellsize) - (cellsize * 0.5)
 pTinSurf.QueryPixelBlock(pBlockOrigin.X, pBlockOrigin.Y, cellsize, cellsize, eRastConvType, vNoData, blockArray)
 pOffset.X = lColOffset pOffset.Y = lRowOffset
 ' The offset for 'write' is the upper left of the pixel block by col/row number. ' Base is 0.
 pRawPixels.Write(pOffset, pPixelBlock)

 Next lColOffset
 ' See if pixelblock size needs to be reset for columns bReset = False If (pBlockSize.X <> lMaxBlockX) Then pBlockSize.X = lMaxBlockX bReset = True End If
 ' See if pixelblock size needs to be reset for rows If ((nRow - lRowOffset) < lMaxBlockY) Then pBlockSize.Y = (nRow - lRowOffset) bReset = True End If
 If (bReset) Then pPixelBlock = pRawPixels.CreatePixelBlock(pBlockSize) blockArray = pPixelBlock.PixelDataByRef(0) End If
 Next lRowOffset
 pRawPixels.ReturnCache(pCache)
 pCache = Nothing

 pRDS = Nothing
 pRawPixels = Nothing
 pPixelBlock = Nothing
 pRasterProps = Nothing
 blockArray = 0 pRDS = dbUtil.OpenRasterDataset(sDir, sName)
 If (lBlockCount = 1) Then
 pTin.TrackCancel = Nothing End If
 TinToRaster = pRDS Exit Function '
 End Function

GIS麦田守望者，期待与您交流。

举报回复喜欢评分

ybtang 路人甲注册日期2005-05-04 发帖数15 QQ 铜币175枚威望0点贡献值0点银元0个加关注写私信	2楼^# 发布于：2006-05-20 00:10 居然有这个函数，我找了Tin开头和raster开头得借口没有发现，真是太感谢了，明天去试试。
	举报回复(0) 喜欢(0) 评分

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

AO开发：将TIn 转换为栅格Grid需要的接口是那些？

最新喜欢：