yuduyue

路人甲

注册日期2004-07-08
发帖数11
QQ101
铜币222枚
威望0点
贡献值0点
银元0个

加关注写私信

阅读：2790回复：5

AE开发TinToRaster的问题

楼主^#

更多发布于：2007-03-26 20:19

<DIV >在帮助系统里面有这样一段vb谢的代码,据说在vb环境下运行的结果是perfect但是我把它改成.net平台下用c@写就没有结果显示,在toc控件中显示value的值巨大,应该是double float的范围值,不知道有没有人做过类似的东东,何故没有结果显示呢 example中的示例代码(VB) ' 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(pTin As ITinAdvanced, eRastConvType As esriRasterizationType, _ sDir As String, sName As String, ePixelType As rstPixelType, cellsize As Double, pExtent As IEnvelope, _ bPerm As Boolean) As 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 IPoint Set pOrigin = pExtent.LowerLeft pOrigin.x = pOrigin.x - (cellsize * 0.5) pOrigin.y = pOrigin.y - (cellsize * 0.5) Dim nCol As Long, nRow As Long nCol = Round(pExtent.Width / cellsize) + 1 nRow = Round(pExtent.Height / cellsize) + 1 Dim pGDS As IGeoDataset Set pGDS = pTin Dim pSR As ISpatialReference2 Set pSR = pGDS.SpatialReference Dim pRDS As IRasterDataset Set pRDS = dbUtil.CreateRasterSurf(sDir, sName, "GRID", pOrigin, nCol, nRow, cellsize, cellsize, ePixelType, pSR, bPerm) DoEvents Dim pRawPixels As IRawPixels Set pRawPixels = dbUtil.GetRawPixels(pRDS, 0) Dim pCache As stdole.IUnknown Set pCache = pRawPixels.AcquireCache Dim pTinSurf As ITinSurface Set pTinSurf = pTin Dim pRasterProps As IRasterProps Set pRasterProps = pRawPixels Dim nodataFloat As Single Dim nodataInt As Long Dim dZMin As Double dZMin = pTin.Extent.zmin Dim vNoData As Variant If (ePixelType = PT_FLOAT) Then vNoData = CSng(dZMin - 1) Else vNoData = CLng(dZMin - 1) End If pRasterProps.NoDataValue = vNoData Dim pOffset As IPnt Set pOffset = New DblPnt ' Set blocksize. Restrict how large it is as not to consume too much memory for ' big output datasets. Dim lMaxBlockX As Long lMaxBlockX = 2048 If (nCol < lMaxBlockX) Then lMaxBlockX = nCol End If Dim lMaxBlockY As Long lMaxBlockY = 2048 If (nRow < lMaxBlockY) Then lMaxBlockY = nRow End If Dim pBlockSize As IPnt Set pBlockSize = New DblPnt pBlockSize.x = lMaxBlockX pBlockSize.y = lMaxBlockY Dim pPixelBlock As IPixelBlock3 Set pPixelBlock = pRawPixels.CreatePixelBlock(pBlockSize) Dim blockArray As Variant blockArray = pPixelBlock.PixelDataByRef(0) ' Set up cancel tracking and progress bar Dim pCancel As ITrackCancel Set pCancel = New CancelTracker pCancel.CancelOnClick = False pCancel.CancelOnKeyPress = True Dim pApp As IApplication Set pApp = New AppRef Dim pProg As IStepProgressor Set pProg = pApp.StatusBar.ProgressBar pCancel.Progressor = pProg Dim lBlockCount As Long lBlockCount = Round((nCol / lMaxBlockX) + 0.49) * Round((nRow / lMaxBlockY) + 0.49) pProg.Message = "Rasterizing. Press ESC to cancel..." pProg.Position = 0 If (lBlockCount = 1) Then ' tin querypixelblock can do the tracking/progressing with 1 block pProg.Show Set pTin.TrackCancel = pCancel Else ' more than 1 block requires this routine, rather than tin function, to track/progress pProg.MinRange = 0 pProg.MaxRange = lBlockCount pProg.StepValue = 1 pProg.Show End If DoEvents ' make sure the bar and the text get updated on screen Dim pBlockOrigin As IPoint Set pBlockOrigin = New Point Dim lColOffset As Long Dim lRowOffset As Long ' Left to right, top to bottom, iteration of pixel blocks. 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) Set 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 If (lBlockCount > 1) Then If (Not pCancel.Continue) Then GoTo Cancel Else If (pTin.ProcessCancelled) Then GoTo Cancel End If Next lColOffset ' See if pixelblock size needs to be reset for columns Dim bReset As Boolean 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 Set pPixelBlock = pRawPixels.CreatePixelBlock(pBlockSize) blockArray = pPixelBlock.PixelDataByRef(0) End If Next lRowOffset 'pProg.Message = "Returning cache..." pRawPixels.ReturnCache pCache Set pCache = Nothing ' need this for some reason with temporary integer grids 'If (Not bPerm) And (ePixelType = PT_LONG) Then ' pProg.Message = "Stats..." ' Dim pBand As iRasterBand ' Set pBand = pRawPixels ' Dim pStats As IRasterStatistics ' Set pStats = pBand.Statistics ' pStats.Recalculate 'End If 'If (bPerm) Then ' flush edits to disk by freeing all pointers 'pProg.Message = "Freeing and opening..." Set pRDS = Nothing Set pRawPixels = Nothing Set pPixelBlock = Nothing Set pRasterProps = Nothing blockArray = 0 Set pRDS = dbUtil.OpenRasterDataset(sDir, sName) 'End If pApp.StatusBar.HideProgressBar If (lBlockCount = 1) Then Set pTin.TrackCancel = Nothing End If Set TinToRaster = pRDS Exit Function Cancel: pApp.StatusBar.HideProgressBar Set TinToRaster = Nothing End Function Public Function CreateRasterSurf(sDir As String, sName As String, sFormat As String, _ pOrigin As IPoint, nCol As Long, nRow As Long, cellsizeX As Double, cellsizeY As Double, _ ePixelType As rstPixelType, pSR As ISpatialReference2, bPerm As Boolean) As IRasterDataset Dim rWksFac As IWorkspaceFactory Set rWksFac = New RasterWorkspaceFactory Dim wks As IWorkspace Set wks = rWksFac.OpenFromFile(sDir, 0) Dim rWks As IRasterWorkspace2 Set rWks = wks Dim numbands As Long numbands = 1 Dim pRDS As IRasterDataset Set pRDS = rWks.CreateRasterDataset(sName, sFormat, pOrigin, nCol, nRow, cellsizeX, cellsizeY, numbands, ePixelType, pSR, bPerm) Set CreateRasterSurf = pRDS End Function Public Function GetRawPixels(pRDS As IRasterDataset, band As Long) As IRawPixels Dim pBandCollection As IRasterBandCollection Set pBandCollection = pRDS Dim pRasterBand As iRasterBand Set pRasterBand = pBandCollection.Item(band) Set GetRawPixels = pRasterBand End Function 我改写的C# 的代码 有没有大哥知道究竟错在哪里, 我看esri online 的discussion上似乎这个QueryPixelBlock用在.net上有问题 而在vb里面似乎可以正常运行出结果</DIV>

喜欢0 评分0

asdfdf

举报回复

yuduyue

路人甲

注册日期2004-07-08
发帖数11
QQ101
铜币222枚
威望0点
贡献值0点
银元0个

加关注写私信

1楼^#

发布于：2007-03-26 20:28

<DIV 12px">我写的C#的代码 public IRasterDataset TinToRaster(ITinAdvanced pTin,esriRasterizationType eRastConvType,string sDir,string sName,rstPixelType ePixelType,double cellsize, IEnvelope pExtent,bool bPerm) { try { IPoint pOrigin; pOrigin = pExtent.LowerLeft; pOrigin.X = pOrigin.X - (cellsize * 0.5); pOrigin.Y = pOrigin.Y - (cellsize * 0.5); int nCol; int nRow; nCol = (int)Math.Round(pExtent.Width/cellsize) + 1; nRow = (int)Math.Round(pExtent.Height/cellsize) + 1; IGeoDataset pGDS; pGDS = pTin as IGeoDataset; ISpatialReference2 pSR; pSR = pGDS.SpatialReference as ISpatialReference2; IRasterDataset pRDS; pRDS = CreateRasterSurf(sDir, sName, "GRID", pOrigin, nCol, nRow, cellsize,cellsize, ePixelType, pSR, bPerm); IRawPixels pRawPixels; pRawPixels = GetRawPixels(pRDS, 0); IRasterProps pRasterProps; pRasterProps = pRawPixels as IRasterProps; stdole.IUnknown pCache; pCache = pRawPixels.AcquireCache() as stdole.IUnknown; ITinSurface pTinSurf; pTinSurf = pTin as ITinSurface; float nodataFloat; long nodataInt; double dZMin; dZMin = pTin.Extent.ZMin; object vNoData; if(ePixelType == rstPixelType.PT_FLOAT) { vNoData = Convert.ToSingle(dZMin - 1); } else { vNoData = Convert.ToInt64(dZMin - 1); } pRasterProps.NoDataValue = vNoData; IPnt pOffset=null; pOffset = new DblPnt(); long lMaxBlockX; lMaxBlockX = 2048; if (nCol < lMaxBlockX) { lMaxBlockX = nCol; } long lMaxBlockY; lMaxBlockY = 2048; if (nRow < lMaxBlockY) { lMaxBlockY = nRow; } IPnt pBlockSize=null; pBlockSize = new DblPnt(); pBlockSize.X = lMaxBlockX; pBlockSize.Y = lMaxBlockY; IPixelBlock3 pPixelBlock3; pPixelBlock3 = pRawPixels.CreatePixelBlock(pBlockSize) as IPixelBlock3; // Get pixeldata array System.Array blockArray; blockArray=(System.Array)pPixelBlock3.get_PixelDataByRef(0); / // Create pixelblock IPnt pBlockOrigin; pBlockOrigin = new DblPntClass(); //定义栅格化时高程的内插方法--自然邻居内插方法 pTinSurf.RasterInterpolationMethod=esriSurfaceInterpolationType.esriNaturalNeighborInterpolation; // Write the pixeldata back System.Object cachePointer; cachePointer = pRawPixels.AcquireCache(); long lColOffset; long lRowOffset; // 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. for(lRowOffset = 0; lRowOffset<=(nRow - 1); lRowOffset=lRowOffset+lMaxBlockY) { for(lColOffset = 0; lColOffset<=(nCol - 1); lColOffset=lColOffset+lMaxBlockX) { if ((nCol - lColOffset) < lMaxBlockX) { pBlockSize.X = (nCol - lColOffset); //See if pixelblock size needs to be reset for columns bool bReset; bReset = false; if (pBlockSize.X != lMaxBlockX) { pBlockSize.X = lMaxBlockX; bReset = true; } //See if pixelblock size needs to be reset for rows if ((nRow - lRowOffset) < lMaxBlockY) { pBlockSize.Y = (nRow - lRowOffset); bReset = true; } if (bReset) { pPixelBlock3 = pRawPixels.CreatePixelBlock(pBlockSize) as IPixelBlock3; blockArray = (System.Array)pPixelBlock3.get_PixelDataByRef(0); } } pPixelBlock3 = pRawPixels.CreatePixelBlock(pBlockSize) as IPixelBlock3; blockArray = (System.Array)pPixelBlock3.get_PixelDataByRef(0); //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; pRawPixels.Write(pOffset, (IPixelBlock)pPixelBlock3); } pRawPixels.ReturnCache(cachePointer); return pRDS; } catch(Exception ex) { System.Diagnostics.Debug.WriteLine(ex.Message); return null; } } //*************************************** //*目的：返回IRawPixels //*输入：Raster数据集 //*************************************** public IRawPixels GetRawPixels(IRasterDataset pRDS, int band) { IRawPixels pRP; IRasterBandCollection pBandCollection; pBandCollection = pRDS as IRasterBandCollection; IRasterBand pRasterBand; pRasterBand = pBandCollection.Item(band); pRP = pRasterBand as IRawPixels; return pRP; } //*************************************** //*目的：创建Raster数据集 //*输入：Raster数据 //*************************************** public IRasterDataset CreateRasterSurf(string sDir, string sName, string sFormat, IPoint pOrigin,int nCol, int nRow, double cellsizeX,double cellsizeY, rstPixelType ePixelType, ISpatialReference2 pSR, bool bPerm) { IWorkspaceFactory rWksFac; rWksFac = new RasterWorkspaceFactory(); IWorkspace wks; wks = rWksFac.OpenFromFile(sDir, 0); IRasterWorkspace2 rWks; rWks = wks as IRasterWorkspace2; int numbands; numbands = 1; IRasterDataset pRDS; pRDS = rWks.CreateRasterDataset(sName, sFormat, pOrigin, nCol, nRow, cellsizeX, cellsizeY, numbands, ePixelType, pSR, bPerm); return pRDS; } 有没有大哥知道究竟错在哪里, 急啊 我看esri online 的discussion上似乎这个QueryPixelBlock用在.net上有问题 而在vb里面似乎可以正常运行出结果</DIV>

asdfdf

举报回复喜欢评分

gis

管理员

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

加关注写私信

2楼^#

发布于：2007-03-26 21:47

我看esri online 的discussion上似乎这个QueryPixelBlock用在.net上有问题 
是有这样的问题,不知道你是不是使用的9.2版本?

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

举报回复喜欢评分

yuduyue

路人甲

注册日期2004-07-08
发帖数11
QQ101
铜币222枚
威望0点
贡献值0点
银元0个

加关注写私信

3楼^#

发布于：2007-03-27 08:17

我用的是AE9.1 SP3 
不知道有没有解决的办法啊

asdfdf

举报回复喜欢评分

cl991036

管理员

注册日期2003-07-25
发帖数5917
QQ14265545
铜币29669枚
威望217点
贡献值0点
银元0个

加关注写私信

4楼^#

发布于：2007-03-27 09:02

vb可以自己转换数值的精度
c#就不回自己转换,要在你转换的数值前加上Convert.ToDouble辅助转换

没钱又丑，农村户口。头可断，发型一定不能乱。邮箱:gisempire@qq.com

举报回复喜欢评分

yuduyue

路人甲

注册日期2004-07-08
发帖数11
QQ101
铜币222枚
威望0点
贡献值0点
银元0个

加关注写私信

5楼^#

发布于：2007-03-27 09:59

谢谢阿
但是我在 if(ePixelType == rstPixelType.PT_FLOAT) { vNoData = Convert.ToSingle(dZMin - 1); }
已经把它转成了单精度浮点型阿,并且rstPixelType.PT_FLOAT就是Pixel values are single precision floating point.应该不需要转成Convert.ToDouble双精度吧,我刚才试了Convert.ToDouble,运行会报错.

asdfdf

举报回复喜欢评分

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

AE开发TinToRaster的问题

最新喜欢：