Latest Contributions by PriyankaMehta

http://forums.esri.com/Thread.asp?c=93&f=993&t=176065&mc=2#msgid518770 Regards, Priyanka ... View more

Hi, Here is a VB code for it. You will have to write it on load event. Private Sub Tool_Load(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles MyBase.Load         ''''' Populates layer combo box. Iterates through TOC and populates the combo box with available layers         On Error GoTo EH      Dim pMxDoc As IMxDocument     Dim pMap As IMap         pMxDoc = My.ArcMap.Application.Document         pMap = pMxDoc.FocusMap         Dim pLayer As ILayer         Dim i As Integer         If pMap.LayerCount = 0 Then             Exit Sub         End If         For i = 0 To pMap.LayerCount - 1             pLayer = pMap.Layer(i)             cmbLayer.Items.Add(pMap.Layer(i).Name)         Next i         Exit Sub EH:         Debug.Print(Err.Number & " " & Err.Description)     End Sub Please try searching for solution on forums and ESRI help documents before posting a thread. You will find most of your problems are already asked and resolved. Below is the link to archived forum http://forums.esri.com/Search.asp?c=93 Its very useful especially in the beginning stage. Hope it helps ! Regards, Priyanka ... View more

Hi,    I am not too sure of what problem you are exactly facing, but I guess this might help. This code joins layer to layer         Dim pFeatLayer As IFeatureLayer         Dim pFeatLayerjn As IFeatureLayer         Dim pDispTable As IDisplayTable         Dim pFCLayer As IFeatureClass         Dim pTLayer As ITable         Dim pActiveView As IActiveView         pMxDoc = My.ArcMap.Document         pMap = pMxDoc.FocusMap         pActiveView = pMap         pMap.ClearSelection()         pActiveView.Refresh()         Dim pTLayer2 As ITable         If pMap.LayerCount < 2 Then             MsgBox("Must have at least two layer", MsgBoxStyle.Information)             Exit Function         End If         ''''''''''''''''*****************************''''''''''''''''''''''''''         pFeatLayerjn = FindLayerByName(pMap, LayerName)          If pFeatLayerjn Is Nothing Then             MsgBox("Could not find Layer " & LayerName & "for joining", MsgBoxStyle.Exclamation)             Exit Function         End If         pDispTable = pFeatLayerjn         pFCLayer = pDispTable.DisplayTable         pTLayer = pFCLayer         pFeatLayer = pMap.Layer(0)           pDispTable = pFeatLayer         pFCLayer = pDispTable.DisplayTable         pTLayer2 = pFCLayer                Dim strJnField As String         Dim JnField2 As String         strJnField = joinfield1           ''''''' FID         JnField2 = joinfield2             ''''''' PntName         ' Create virtual relate         Dim pMemRelFact As IMemoryRelationshipClassFactory         Dim pRelClass As IRelationshipClass         pMemRelFact = New MemoryRelationshipClassFactory         pRelClass = pMemRelFact.Open("LayertoLayer", pTLayer2, strJnField, pTLayer, JnField2, "forward", "backward", esriRelCardinality.esriRelCardinalityOneToOne)         ' use Relate to perform a join         Dim pDispRC As IDisplayRelationshipClass pDispRC = pFeatLayer pDispRC.DisplayRelationshipClass(pRelClass, esriJoinType.esriLeftOuterJoin)           MsgBox("Fields Joined", MsgBoxStyle.Information) Referring this one might also help http://forums.esri.com/Thread.asp?c=93&f=992&t=97996&mc=3#msgid274992 Regards, Priyanka ... View more

Hi,    Are you by any chance decreasing the resolution of your raster i.e. increasing the pixel size from original. I tired the code on a raster with a pixel size of 5 originally. If I resample it with cell size 7 I will definitely loose information. Even the size of the raster on hard disk will decrease. On the other hand, if I change the size to 2, then I obviously wont loose any information and there are no round offs. Moreover, in your case, since your pixels are not square originally, that might also be a reason for lost information somewhere. Try the same code on a standard raster and see if there is any difference. Does it work fine when you do it manually by using 'Data' --> 'Export Data'. ... View more

Hi, I guess this should work. You can exclude killifexistsraster function if you want, because it deletes the output raster without warning when you try it again. Dim pMxDoc As IMxDocument         Dim pMap As IMap         Dim pRasterLayer As IRasterLayer         Dim pRaster As IRaster         Dim pGeometryProc As IRasterGeometryProc         pMxDoc = My.ArcMap.Application.Document         pMap = pMxDoc.FocusMap         pRasterLayer = pMap.Layer(0)         pRaster = pRasterLayer.Raster                pGeometryProc = New RasterGeometryProc         pGeometryProc.Resample(rstResamplingTypes.RSP_NearestNeighbor, 7, pRaster)              Dim pDataset As IDataset         pDataset = pRasterLayer         Dim pWS As IWorkspace         pWS = pDataset.Workspace         Dim str As String         str = pWS.PathName         killifexistsRaster("try", str)         Dim RasterDataset As IRasterDataset         Dim pRasBandCol As IRasterBandCollection         pRasBandCol = pRaster         'pDS = pRasBandCol.SaveAs("try3", pWS, "GRID")         RasterDataset = pRasBandCol.SaveAs("try", pWS, "GRID")         '***********************************         Dim pRasterLy As IRasterLayer         pRasterLy = New RasterLayer         ' pRasterLayer.CreateFromRaster(pRaster)         pRasterLy.CreateFromDataset(RasterDataset)         pMap.AddLayer(pRasterLy)         pMxDoc.ActiveView.Refresh()         '***********************************         SaveAsWithoutBuildingPyramids(pRaster)     End Sub     Sub SaveAsWithoutBuildingPyramids(ByVal pRaster As IRaster) 'Or as IRasterDataset)         killifexistsRaster("tryexp", "C:\temp")         Dim pSaveAs2 As ISaveAs2         pSaveAs2 = pRaster         Dim pStorageDef As IRasterStorageDef2         pStorageDef = New RasterStorageDef         pStorageDef.PyramidLevel = 0         Dim pWs As IWorkspace         pWs = SetRasterWorkspace("c:\temp")         pSaveAs2.SaveAs("try_exp", pWs, "GRID")     End Sub     Public Function SetRasterWorkspace(ByVal sPath As String) As IRasterWorkspace         Dim pWsFact As IWorkspaceFactory         pWsFact = New RasterWorkspaceFactory         SetRasterWorkspace = pWsFact.OpenFromFile(sPath, 0)     End Function     Public Function killifexistsRaster(ByVal RasterName As String, ByVal strFolder As String)         '''''' Delete raster from disk         Dim pRWSF As IWorkspaceFactory         pRWSF = New RasterWorkspaceFactory         Dim pRWS As IRasterWorkspace         pRWS = pRWSF.OpenFromFile(strFolder, 0)         Dim pRasterDS As IRasterDataset         On Error Resume Next         pRasterDS = pRWS.OpenRasterDataset(RasterName)         If Not pRasterDS Is Nothing Then             Dim pDS As IDataset             pDS = pRasterDS             If pDS.CanDelete Then pDS.Delete()         End If         System.IO.Directory.Delete(strFolder & "\" & RasterName, True)         DeleteFolder(strFolder & "\" & RasterName)     End Function Public Function DeleteFolder(ByVal strPath As String) As Boolean         Dim objDir As New IO.DirectoryInfo(strPath)         Try             objDir.Delete(True)             Return True         Catch             Return False         End Try     End Function ... View more

So now what I understand is, that you are being able to standardize it in terms of pixel size but it becomes an 8 bit raster because it keeps forcing the renderer. Do you want to share your code developed till now .. may b we can catch something ... View more

What exactly are you attempting to explain when you say 'non standard raster' ? If you want to change the pixel width, have you tried resampling(because your final raster has x=y cellsize) before exporting. Regards, Priyanka ... View more

Hi Neil,   Thanks for the reply. I still have some doubts. Create a tool that allows the user to click the map. Earlier I was using 9.3 there I used to create a UI Tool control for such purposes. Now that I am using version 10, I am not really being able to find / create a tool that would allow the user to click the map !! So basically, the main problem I am facing is "clicking on screen/map" since I cant find UI Tool control or something similar to that in version 10. Regards, Priyanka ... View more

Thanks a lot Sandhya! The codes I was trying available from the old forums were running but were too slow. I found this tool http://code.google.com/p/typeconvert/ which also has 'remove duplicates' and runs pretty fast and also creates 2 shapefiles : (1) With no duplicate polylines (2) the duplicate polylines that were deleted So now I use (2) to spatially locate and select polylines exactly overlapping (1) .. export these selected polylines . Join fields and do the calculation.. seems to be a bit of round about but is working for me for now.. Any suggestions to make it a simpler process are most welcome thanks and regards, Priyanka ... View more

Thanks a lot again Brian, The code looks neat to me. I am yet to test it though.. The reason I want to do it to determine the sky view factor. It is a value that determines how much sky is obstructed. The raster pixel values are nothing but the height of the buildings. If i am standing in the midst of tall buildings the amount of sky i would see is the sky view factor. The first tallest building that comes in my sight is the one that would be taken into consideration (Reason for the first maximum pixel) If there are consecutive pixels of the same elevation, then they might just be the breadth of that building It is calculated by a formula and the output from above code is one of the input of the formula. It is an interesting thing to study. Usually done for temperature modeling Thanks again, Priyanka ... View more

Found one nice VBA code. http://forums.esri.com/Thread.asp?c=93&f=987&t=172882&mc=6#msgid508922 Thanks, Priyanka ... View more