Latest Contributions by NeilSang

Hi, I keep getting an error message with the following code. The error message is just generic but I suspect the problem is with the output feature class but I can't figure out exacltly what it is wanting for input (blank FClass, FClass with all the fields prepared?) Any ideas / example code much appreciated :                        Dim TargetLayer As IFeatureLayer = New FeatureLayer             TargetLayer.FeatureClass = pFClass             Dim JoinLayer As IFeatureLayer = New FeatureLayer             JoinLayer.FeatureClass = hchange_dFClass             'SpatialJoin             Dim JoinFeat As IFeatureClass = ws.CreateFeatureClass("JoinFeat", _fields, _uid, Nothing,  esriFeatureType.esriFTSimple, "SHAPE", "")             'Dim JoinFeatObj As Object = New Object             Dim ISJ As New SpatialJoin             ISJ.join_features = JoinLayer             ISJ.target_features = TargetLayer             ISJ.out_feature_class = JoinFeat             gptool.Execute(ISJ, Nothing) ... View more

Hi, I am trying to simply add a layer in programatically. The layer appears in ArcMap, but there is no data there in. If I try to look at layer properties, ArcMap crashes. I've tried both shape and Access workspacefactory , and various shape files in various paths (no spaces etc..) but the same problem. Refreshing the data doesn't make any difference. I can see from the featurelayer properties that the correct data has been found, with the correct number of fields. What am I doing wrong? Thanks Neil     Private Sub AddShapefile(ByRef m_Application)         Dim mapDoc As IMapDocument = m_Application.Document         ' Create a ShapefileWorkspaceFactory object and         ' open a shapefile folder. The path works with a standard 9.3 installation.         Dim AccworkspaceFactory As IWorkspaceFactory = New AccessWorkspaceFactoryClass()         Dim featureWorkspace As IFeatureWorkspace = TryCast(AccworkspaceFactory.OpenFromFile("C:\WindaWorkSpace\Data\WindaPGdb.mdb", 0), IFeatureWorkspace)         ' Create a FeatureLayer and assign a shapefile to it.         Dim featureLayer As IFeatureLayer = New FeatureLayerClass         featureLayer.FeatureClass = featureWorkspace.OpenFeatureClass("winda")         Dim layer As ILayer = featureLayer         layer.Name = featureLayer.FeatureClass.AliasName         ' Add the layer to the focus map.         Dim pMap As ESRI.ArcGIS.Carto.IMap = mapDoc.ActiveView.FocusMap         pMap.AddLayer(featureLayer)         mapDoc.ActiveView.Refresh()         Dim t As Integer         t = 1     End Sub EDIT : Once again Neil Clemmens saves the day.. When automating from an external application you need to use the IObjectsFactory : http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#//0001000001nn000000 Given that most people must now be using VS 2010 or similar to develop, this should be written in big red letters the first page of the esri user help pages! ... View more

Hi, I have a GWR model with good r2 and AICC. I've tried as far as possible to remove co-linearity from the model, the explanatories for which were selected first by longitudinal MLR of the aggregated data (to avoid spatial auto-correllation), and afew more weeded out through OLS. I know there is still some spatial co-linearity there but as the explanatories vary in significance over space none of the remaining parameters may be considered redundant to the global model. I would like now to predict the dependent variable's map for new scenarious by plugging in different parameter values to one or more of the co-efficients, just as one might in MLR (I appreciate that for GWR this relies on the substantial assumption that the co-efficients will be stable over time accross space). But I'm concerned that local co-linearity may be inflating the effect in some locations for some parameters. Any suggests as to how to model this? Is it legitmate in GWR to attempt to look at the effect of idividual variables in this way? Thanks, Neil p.s. The data is already normalised to a percentage change from each variables mean (i.e. the regression is built on the variance not the absolute values). ... View more

Hello, I have a script which identifies the locations of intersects between n-polygon boundaries and a line (a vector) using shape files without any spatial index. At present it generates a polyline and uses a spatial filter to identify the relevant polygons, then loops through the vertices of each to find the intersection points. It is rather slow as it has to be repeated thousands of times. Implementing a spatial index will be lengthy as it will mean recoding all the references to shp files to access a gdb throughout the software and as I didn't write it I'd rather restrict my tinkering to the one, easily restored, class. My question is what alterative functions users might suggest that could be quicker. I am considering afew : iRay to identify the interesect points Or LocateFeaturesAlongRoutes_lr? Or convert everything to lines and clip to a buffer of the vector first before doing intersect testing.. I am sure there are otherways to skin the cat, but which is likely to see a time saving? Any suggestions or references to a good source of explainations about the algorithms behind the classes would be welcome also, for example I assume spatial filter deploys some sort of prior minimum bounding geometry test but I'd like to know for sure. Thanks for your time, Neil ... View more