Latest Contributions by XanderBakker

‎10-08-2013

Hi Xander Thanks for your feed back. Yes I realised the necessity for a uniqueID but to perform the hotspot (getis-ord) it requires a count field. According to Lauren first the incident points have to be integrated and then have to be collected. But once these transformations have been done the time field of the original incidents are by default removed. when I create a spatial weights matrix using the original incidence which have the time field and later add a uniqueID to it I still find hard how to create a hotspot for temporal data without a count field I have attached an Image to show the problem I am facing hope you can please help me in this matter[ATTACH=CONFIG]28054[/ATTACH] Many Thanks Thivanka Hi Thivanka, If I look at image with the dialog for the 'Hot Spot Analysis (Getis-Ord Gi*)' calculation and the listbox for the Input Field that remains empty, I think this is caused due to not having a numeric field in your featureclass 'conftimetestsl_CopyFeatures_1'. In the Help the description is: The numeric field (number of victims, crime rate, test scores, and so on) to be evaluated. If you don't have this field (since each record is a single incident) you could create this field by adding a numeric field and fill it with the value 1. Not sure if that's the answer you're looking for... Kind regards, Xander Bakker

‎10-03-2013

In the below video What's New in Spatial Statistics: Space-Time Cluster Analysis at minutes 6.33 where it is shown how the spatial weights matrix is generated for time window, when the UniqueID is selected an ICOUNT field is also visible. My problem is once I integrate my data and collect the events the only selectable field is only ICOUNT. Please help me to create an associating UniqueID while having the relationships to the time fields in my original point incidents. Many Thanks Thivanka Hi Thivanka, In the video Lauren is using a field called "UniqID". This is a custom field she created. The point is that the field needs to be integer and normally an ID field like ObjectID is of type Long. You can create the integer field yourself and set it equal to the ObjectID, if the ID's are in the integer range... Kind regards, Xander

‎10-03-2013

Hi Xander, along these same lines, do you know how to run a least cost path between multiple points that results in a pairwise least cost path for all points? So, if you have 5 points, you would have 10 paths? Thanks, Stephanie Garvis Hi Stephanie, If you would like to calculate from a set of 5 points the least cost path between each pair of points, you shoud create cost distance rasters for each point. So we're talking about 1 calculation for each point (5 in total). Then for point 1 you take the cost distance raster and the 4 other points to calculate the 4 paths. For point number 2 you take points 3, 4 and 5 and calculate the paths. For point number 3 take points 4 and 5 en for point 4 take point 5. I do not see a faster way since to determine the path to a source you will need the cost distance raster specifically for that source. Kind regards, Xander

‎10-02-2013

Hi, 1. What is the difference between the layer file with �??end to end blue arrows�?? icon and the the layer file with the yellow diamond icon? 2. What is the format of the layer files with �??end to end blue arrows�?? icon? 3. Is the format of the layer files with a yellow diamond icon .lyr? 4. Should I select the layer file (e.g. elevation) with �??end to end blue arrows�?? icon or layer file with yellow diamond icon? Data variables are shown with the Model Variable icon (what you refer to as 'end to end blue arrows'), and layers from the ArcMap table of contents are shown with the Layer icon ('yellow diamond icon'). In your case they both refer to the same data source, but in order to avoid adding the same layer twice as a different variable, you should use the Model Variable icon ('end to end blue arrows'). You can read more about this in the topic Adding and connecting data and tools (lower part) 5. I would like to calculate the least cost distances (LCD) between more than two points, however the �??Spatial Analyst Tutorial�?? provides the workflow for only two points (the final site and destination in Stowe.gdb). Can you tell me the complete workflow for calculating the LCDs among more than two points (see, among the five schools in the tutorial)? If there are multiple cells or zones as input destinations, the least-cost path can be calculated from each cell (resulting in multiple paths, one path for each cell), each zone (one path from each zone), or by layer (only one path, the cheapest from any zone). A parameter in the tool specifies which process should be utilized. Source: http://resources.arcgis.com/en/help/main/10.2/index.html#//009z00000021000000 Your Cost Path is colored white since it cannot be executed. It requires at least three parameters and only two are specified. In your case the "in_destination_data" is missing. See Help topic: http://resources.arcgis.com/en/help/main/10.2/index.html#//009z00000019000000 So in short: The Cost Distance will create a accessibility from 1 or more sources The Cost Path will calculate the path from 1 or more destinations back to the best accessible source Kind regards, Xander

‎10-01-2013

Hi Yoshi, I'm glad I could be of any assistance. If you think it was helpful you can use the "arrow" button in order to help other members find useful information: [ATTACH=CONFIG]27933[/ATTACH] More info here: http://resources.arcgis.com/en/help/forums-mvp/ Kind regards, Xander

‎10-01-2013

Hi Josseline, I'm glad I could be of any assistance. If you think it was helpful you can use the "check mark" and "arrow" buttons in order to help other members find useful information: [ATTACH=CONFIG]27931[/ATTACH] More info here: http://resources.arcgis.com/en/help/forums-mvp/ Kind regards, Xander

‎10-01-2013

Hello Xander Using Phyton may be a good idea to go around this problem, but the problem is I am a newone for Phyton. What I want to do is as the followings: 1) use raster calcultor to evaluate "raster01 > 5000", 2) in the resulting calculation layer (calc01), resympolize the "0" values to be "No Color" or transparent. I checked ArcGIS Desktop Help page and tried to use Phyton but might have some mistakes. Any help or suggestion will be welcomed. >>> import arcpy >>> from arcpy.sa import * >>> myRas = arcpy.Raster('raster01') >>> from arcpy import env >>> env.workspace="C:\GISdata\workspace" >>> inRaster="raster01" >>> "raster01>5000" 'raster01>5000' >>> arcpy.CheckOutExtension("Spatial") u'CheckedOut' >>> Calcuation.save("C:\GISdata\workspace") Runtime error Traceback (most recent call last): File "<string>", line 1, in <module> NameError: name 'Calcuation' is not defined Hi Yoshi, You're almost there. A few things should be changed. import arcpy from arcpy import env env.workspace = r'C:\GISdata\workspace' # arcpy.CheckOutExtension("Spatial") myRas = arcpy.Raster('raster01') myResult = myRas > 5000 myResult.save('myOutputName') A few comments on the differences: When setting your workspace, use forward slashes, double backslashes or just put an 'r' before the path to avoid errors: env.workspace = r'C:\GISdata\workspace' Switching on the Spatial Analyst extension and setting the workspace (Geoprocessing/Environments) can also be done through the menu. The Python window will inherit these settings. I assume your TOC holds your input raster "raster01" When performing a calculation you need to specify an output (like "myResult"): myResult = myRas > 5000. At this point the SA license should be available. When saving the raster you specify the name you want to assign to the raster (example "myOutputName") The raster will be stored in the current workspace (in this case "C:\GISdata\workspace") The result is a raster with 0 (raster01 <= 5000) and 1 (raster01 > 5000) and NoData cells where the input has Nodata cells. Kind regards, Xander

‎09-30-2013

Hi Yoshi, I should have read your question more carefully. In this case I think it is best to contact Esri support and provide all the information about the system you're installing on, the version you're using. Just checked it on my systems and in 10.2 and 10.1 (SP1) I have a working expression window. The alternative at this moment is to use the Python window in ArcGIS: What is the Python window? Example: import arcpy from arcpy.sa import * # myRas = arcpy.Raster('the name of the raster in the TOC or path to raster') myRas = arcpy.Raster('J000_CumCO2_kgha') # now you are able to perform any type of raster calculator analysis; like for instance: myResult = Con(IsNull(myRas),FocalStatistics(myRas,NbrRectangle(5,5,"CELL"),"MEAN","DATA"),myRas) Some more reading: Executing tools in the Python window Kind regards, Xander

‎09-29-2013

I�??m trying to make a DEM of an old mining site. I have point data in a 10*10 meter grid, some lines representing cliffs, and I also have a lake boundary polygon, representing the flooded area of the cavity. When I try to create the DEM with the create TIN or Topo to raster functions, it appears with peaks in the measured height points place. The only way I managed to produce a smooth surface was with the spline tool, but it does not allow the use of barriers. Is there another way to create a smooth surface with those constraints? Here�??s a printscreen of the two surfaces (converted to TIN) in the arcscene (left: spline, right: topo to raster with constrains): [ATTACH=CONFIG]27850[/ATTACH] Thank you and apologies for my bad English. Hi Anton, What you could do is create a raster that contains you constraints and has NoData where there are none. Let's call this raster "Constraints". In the example below DEM is the DEM created with constraints, in your image on the right. You could perform the calculation below: Con ( IsNull ( Constraints ) , FocalStatistics ( DEM , NbrRectangle(3,3,"CELL") , "Mean" , "DATA" ) , DEM ) This will detect those areas with out constraint and calculate the mean elevation in a 3x3 cell neighborhood. The areas with constraints will maintain their value and will not be affected. Some more reading: Focal Statistics (Spatial Analyst) How Focal Statistics works Conditional evaluation with Con Is Null (Spatial Analyst) Kind regards, Xander

‎09-29-2013

Hi Yoshi, The Raster Calculator is still there, but has moved from the Spatial Analyst toolbar to the Toolbox: Spatial Analyst Tools\Map Algebra\Raster Calculator [ATTACH=CONFIG]27864[/ATTACH] Kind regards, Xander

‎09-29-2013

Hi all, I have a question about using the DEM created by Topo to Raster in ArcScene (10.2). When I try to set the base heights so that the DEM displays as a 3D surface, I don't have the option to "Use elevation values in the layer's feature." Any suggestions on how to fix this, or what I may have done wrong? I generated the DEM from a georeferenced contour line file. I should mention that I have only recently started getting back into GIS after several years away, and it's entirely possible that I have forgotten a crucial step. Thanks, Meg [ATTACH=CONFIG]27817[/ATTACH] Hi Meg, In the top of the dialog you'll see "Elevation from surface" you should change the option to "Floating on a custom surface". Then your raster demtest3 will be used for elevation. The option "Use elevation values in the layer's feature" is only for vector data and your DEM is raster. Kind regards, Xander

‎09-29-2013

Hello, I recently acquired LIDAR files which are in .asc format. My objective is to create an automatic model using modelbuilder so any LIDAR files of a watershed area can be processed in order to extract paramters such as geographical extend... But I'm stuck at the first step, I have the raw LIDAR data in .asc format and when I open it in notepad++, it shows only 3 columns and about 3 millions rows of datas, the matter is that I have no spatial reference, I don't have any information such as : NCOLS xxx NROWS xxx XLLCORNER xxx YLLCORNER xxx CELLSIZE xxx NODATA_VALUE xxx I can specify myself the number of rows, columns and the nodata value but I can't specify the cellsize or the geographical extend of the dataset. I can't create a Raster using the ASCII to Raster Tool (which is strange by the way considering a raster doesn't need spatial references to exist). I'd like to know how I could convert such files into Raster data (preferably TIFF) which I can then specify spatial references. Thanks in advance, Josseline Hi Josseline, The data you received is an XYZ ASCII point cloud (vector) which is quite different from the raster format you're trying to obtain. Please read these help topics for more information on both formats. Esri ASCII raster format and What is lidar data? LiDAR data is normally a lot richer than XYZ only. In your case it has been stripped down to XYZ. ArcGIS has several ways of working with LiDAR data, but those are normally based on LAS files which are imported to a LAS dataset: What is a LAS dataset? At this moment Esri still does not provide a simple straightforward way of converting XYZ (ASCII) files to LAS. If you also feel that this is a necessary tool you can support the idea here (it's under consideration): Ability to convert xyz point data to LAS format In the mean time there are several other ways of doing this. I found a tool on the ArcScripts site (although I have no idea what version of ArcGIS it supports): ArcScripts: Convert XYZ,ASCII Raster files to ESRI GRID-Batch Another possibility is using the ASCII 3D To Feature Class (requires 3D Analyst) which results in a PointZ featureclass followed by a Point to raster. You can read more about this in this Esri Blog: Lidar Solutions in ArcGIS_part2: Creating raster DEMs and DSMs from large lidar point collections. A very powerful toolset is LAStools: http://www.cs.unc.edu/~isenburg/lastools/ and http://rapidlasso.com/lastools/ Please bear in mind that LAStools is not free for commercial applications: http://rapidlasso.com/pricing/ Some more reading on Creating raster DEMs and DSMs from large lidar point collections Kind regards, Xander

‎09-25-2013

I have a set of x,y coordinates spanning the entire northern hemisphere of Mars, and I want to see how many of them fall within 10° of the edges of a given shapefile (polygon). Is there a way I can quickly edit the shapefile to create a version with the boundaries extended by 10°, or is there a way to calculate the intersection of my existing files to account for this 10°? Since the shapes are irregular it would be time consuming to manually extend the boundaries by creating new shapes. Hi, Although I don't know what support ArcGIS has for data on Mars, the normal approach would be to use the Select By Location. Extended Help topic is found here: Using Select By Location (10.2) Use the Intersect spatial selection method. Switch on the option "Apply search distance", enter "10" and use decimal degrees as units. In case you only want to only obtain those feature that are outside of your polygon, you can perform a second select by location on the same layers, but now choosing the "remove from currently selected features" option for the selection method (top of dialog) and switch off the option apply search distance. Kind regards, Xander

‎09-25-2013

Hello! Is it possible to perform Cluster Analysis (Need Hot Spot Analysis) for multiple critearia? I have a .shp file with polygons, each polygon has several fields: A,B,C,D,E. Is it possible to find clusters based of ALL fields A-E? How can I do it? If its not what is the best way to carry out this cluster analysis? I can perform Hot spot analysis for each field - have several layers - combine them? Thank you very much! Hi Jaroslav, I think the new "Grouping Analysis" (10.2) could be just what you are looking for: Help on Grouping Analysis (Spatial Statistics) 10.2 and How Grouping Analysis works Kind regards, Xander

‎09-24-2013

Hi David, Not sure if this has anything to do with it, but I notice in your original post that you obtain your "SubmitDate" through: SubmitDate = arcpy.GetParameterAsText(1) Your "SubmitDate" is a string. I don't know how you have defined the type of your field "SUBMITTALDATE", but if it is a date field, it may be better to provide an actual date. To convert a string to a datetime object see snippet below (but, you need to know the format how it's specified): from datetime import datetime submitDate = '24-9-2013' dateObj = datetime.strptime(submitDate,'%d-%m-%Y') submitDate = '9/24/2013' dateObj = datetime.strptime(submitDate,'%m/%d/%Y') More info on syntax: http://docs.python.org/2/library/datetime.html#strftime-and-strptime-behavior The alternative could be to read it as an object if it's provided as a date: SubmitDate = arcpy.GetParameter(1) See also the Help: http://resources.arcgis.com/en/help/main/10.1/index.html#//018v00000038000000 Kind regards, Xander

Online Status	Offline
Date Last Visited	3 weeks ago

My Ideas

Latest Contributions by XanderBakker

Re: How to create a UniqueID field together with ICOUNT- Space Time Analysis

Re: How to create a UniqueID field together with ICOUNT- Space Time Analysis

Re: The least cost distances between more than two points

Re: The least cost distances between more than two points

Re: raster calculator expression box missing - ArcGIS 10

Re: LIDAR asc files without spatial references processing ArcGIS 10.0

Re: raster calculator expression box missing - ArcGIS 10

Re: raster calculator expression box missing - ArcGIS 10

Re: Creating smooth DEM with constraints

Re: raster calculator expression box missing - ArcGIS 10

Re: Problems with Topo to Raster DEM

Re: LIDAR asc files without spatial references processing ArcGIS 10.0

Re: Calculate Intersection Within Range of Polygon Edge?

Re: Cluster Analysis for multiple criteria

Re: Problem with where clause and Update Cursor

Re: Dissolve/merge features without creating new f...

Conditional Field display with Arcade in Pop Ups (...

Re: Pop-up Arcade Expression- Show List of Field N...

Re: Arcade Text Formatting - Bold Text

Re: Why does the HTML pop-up strip internal spaces...

Esri Forestry Group (EFG)

Human Brain WebGIS

Python Snippets

LIDAR Resources