Latest Contributions by XanderBakker

[ATTACH=CONFIG]27251[/ATTACH][ATTACH=CONFIG]27252[/ATTACH][ATTACH=CONFIG]27253[/ATTACH] Hi, i got the footprint,lidar data and aerial photo which i need to use for constructing urban 3d model. i saw some tutorials online, but there is one problem for me. i can open lidar dataset in arcmap (as showd in picture 1 - ground lidar dataset), there are many points in 2d. but when i open it in arcscene 10.1, first several seconds there is a just a box (picture 2), after that nothing shows in the arcscene. but i saw online tutorial, they all can open las dataset in acrscene like cloudy points. Hi Yang, Maybe the amount of data causes this effect. You could try to reduce the amount of data to be displayed by using the "Area of Interest tool" in ArcScene. More info here: http://resources.arcgis.com/en/help/main/10.1/index.html#//00q800000144000000 Kind regards, Xander ... View more

Hello there, I'm trying to convert from raster to vector using the raster to polygon tool but I keep getting an error (000864). It says that the input is not in the defined domain and I'm not sure what that means. Do I have to do another conversion before this one? thank you, im extremely new to arcmap. Hi Carla, As Milan indicates you should have a look at the error description, which you can access through the Menu "Geoprocessing", option "Results". This will open the Results Window which lists the Geoprocessing tools and any errors that may have occurred. Expand the tool that generated the error and below expand the "Messages". See an example below: [ATTACH=CONFIG]27293[/ATTACH] This will reveal the error and provides a link to the Help on that error. In your case (in 10.2) 000864 will point to: http://resources.arcgis.com/en/help/main/10.2/index.html#//00vp00000012000864 000864 : <value>: The input is not within the defined domain. <value>. Description The specified input is not valid for the parameter. You might have tried to use a value for the parameter that is not one of the defined options or is not a valid type or perhaps made a typing error. The command syntax will identify the valid options for each function. For example, the options for the Filter type parameter are only LOW or HIGH. No other text string will be accepted. Solution Check the syntax and documentation to make sure you are specifying the command correctly and that the input types are valid. In order to help you, we will need to know what data you used for this tool. Kind regards, Xander ... View more

Hola a todos, no se si alguien me podria dar una ayuda, tengo un archivo en formato .mxd, y la verdad quiero graficarl en 3D Polyline, y me duele la cabeza buscando de que manera, es una area de un parque nacional, que me gustaria proyectalo en Espacio 3D para asi oder apreciar mucho mas todas sus caracteristicas que posee dicha area, tanto altura, profundidades, caminos, rios, pendientes y demas detalles que me son utiles pa asi monitorear mejor. Por favor denme una ayuda pa lograrlo. Muchas gracias, se lo agradecere de corazon. FYI; included link to Spanish Forum and some introductory topics in Spanish in answer below. Hola DanielAlexandar, Te recomendaría que colocas tu pregunta en el foro en Español. Lo encuentras aquí: http://forums.arcgis.com/forums/50-Espa%C3%B1ol Para contestar tu pregunta un poco; el archivo.mxd es un documento de ArcMap. Normalmente contiene referencias a capas de datos (layers) y su forma de visualización. Lo que entiendo de tu pregunta, es que quieres visualizarla en 3D. Para este fin ArcGIS tiene las aplicaciones ArcScene y ArcGlobe (que requieren la extensión 3D Analyst). Para conseguir más información te recomiendo comenzar por un recorrido rápido por ArcMap: http://help.arcgis.com/es/arcgisdesktop/10.0/help/index.html#/na/006600000332000000/ .. y luego leer el tutorial de ArcGIS 3D Analyst: http://help.arcgis.com/es/arcgisdesktop/10.0/help/index.html#//00q800000093000000 Cordialmente, Xander ... View more

Hello, I'm tryng to populate a field with this VBA code but it doesn't work. I really appreciate any help you can provide. Dim codigo as string If [CAMPO]="Captura" OR [CAMPO]="Replanteo" OR [CAMPO]="GPS" Then        codigo = "1"     ElseIf [CAMPO]="Plano" OR [CAMPO]="Ortofoto" Then         codigo = "2"     ElseIf [CAMPO]="Catastro" OR [CAMPO]="Ajuste parcelario" Then         codigo = "3"     Else codigo = [CAMPO] EndIf __esri_field_calculator_splitter__ codigo Hi alaska@pozo, The only thing I see in your code it the "Else codigo = [CAMPO]" being on 1 line, but changing it doesn't solve the problem. What you can do is change your code and use a Select Case. This simplifies the code and seems to work (at least for me in 10.2): Select Case [CAMPO] Case "Captura", "Replanteo", "GPS"   codigo = "1" Case "Plano", "Ortofoto"   codigo = "2" Case "Catastro","Ajuste parcelario"   codigo = "3" Case Else   codigo = [CAMPO] End Select __esri_field_calculator_splitter__ codigo Hope this works for you too. Kind regards, Xander BTW: since this is the Python forum, you could consider using the Python syntax: def ClfyCampo(campo):     if campo == "Captura" or campo == "Replanteo" or campo == "GPS":         return "1"     elif campo == "Plano" or campo == "Ortofoto":         return "2"     elif campo == "Catastro" or campo == "Ajuste parcelario":         return "3"     else:         return campo __esri_field_calculator_splitter__ ClfyCampo( !CAMPO! ) ... View more

Good morning, I am not sure if this is possible, and I am still fairly new at 3D Analyst. I am wondering if there is a way to convert a 3D shapefile into a TIN so that an image can be draped over it? Or if there is a simpler way to do this? Thanks for your time! Hi Maddy, In case your data consists of points and/or lines, it is better to create a TIN from it. The following links provide some help on that topic: http://resources.arcgis.com/en/help/main/10.2/index.html#//006000000010000000 http://resources.arcgis.com/en/help/main/10.2/index.html#/Create_TIN/00q90000001v000000/ If you have 3D polygons (Polygon Z) you could use it as the floating layer and drape the other layer on top of it. Some explanation of that can be found in the following help topics: http://resources.arcgis.com/en/help/main/10.2/#/Defining_data_as_a_draped_layer_in_ArcScene/00q8000000t5000000/ http://resources.arcgis.com/en/help/main/10.2/index.html#/Understanding_feature_based_heights_in_3D/00q80000013v000000/ http://resources.arcgis.com/en/help/main/10.2/index.html#//00q8000000zp000000 Since you say you're fairly new to 3D Analyst, it may be a good idea to follow an introductory course on the topic. Here you find a list of training provided by Esri: http://training.esri.com/gateway/index.cfm?fa=search.results&searchterm=&search=Search&orderby=&softwareid=20000 Kind regards, Xander ... View more

Hello, I am trying create a raster that shows road density (m/ha) within 3km radius of the cell.  The line density tool in Spatial Analyst appears to be the tool to do this as it allows for input of a line file (roads), a search radius (3000m), output raster cell size (10m), and area unit (ha). However, upon running the tool I getting results that I don't believe to be correct. For example the maximum value I am getting is 0.27 m/ha, which is way to low of a value considering to the high road density in the area. To double check I converted a portion of the raster to points and buffered those points by 3 km and found the area of road within this buffer. Using this information the road density values come out to around 9 or 10 m/ha which is much more reasonable for the density of roads in the area. The dataframe and road layer is in UTM Zone 12N so the units should be in meters. Any ideas why these results are not correct? I wonder if it has to do with the population field being set to none as in the tool help it indicates the population field is a: Numeric field denoting population values (the number of times the line should be counted) for each polyline. Values in the population field may be integer or floating point. The options and default behaviours for the field are listed below. Use None if no item or special value will be used and each feature will be counted once. Therefore I'm wondering if it is only counting each road line feature the first time it intersects a cell search radius instead of counting it within the search radius of each cell. Any advice or suggestions would be greatly appreciated. Below is my inputs if that helps. Thanks in advance. Hi Karl, If you look at the explanation in the Help you will see that the density is calculated as follows: Density = ((L1 * V1) + (L2 * V2)) / (area_of_circle) Lines L1 and L2 represent the length of the portion of each line that falls within the circle. The corresponding population field values are V1 and V2. Source: http://resources.arcgis.com/en/help/main/10.2/index.html#//009z00000012000000 In your case, you are not applying a population field, so it is reduced to: Density = (L1 + L2) / (area_of_circle) Since you are using a search radius of 3000m, the area of the circle is ± 28274333m² (or 2827ha or 28.27km²). This causes the density to result in very low values. In Help you will also find this: A default area unit is  selected based on the linear unit of the projection definition of the  input polyline feature data or as otherwise specified in the  output coordinate system environment setting. For line density, when an area unit factor is specified, it converts the units of both length and area. For  example, if the linear unit is meters, the output area units will  default to SQUARE_KILOMETERS and the resulting line density units will  convert to kilometers per square kilometer. The end result, comparing an  area scale factor of meters to kilometers, will be the density values  being different by a multiplier of 1,000. (... so not 1,000,000) Hope this helps you to interpret the results, although it may not be what you are looking for... Kind regards, Xander ... View more

Awesome idea Xander.  Thank you.  Although I only have 8 polygons to extrude over the whole USA.  So mine won't look as cool as yours.  Plus some of my polygons with small values are behind some with much larger values.  It would cover the stack. I could create a new polygon shapefile with just squares at the center of each large region, then just join the region $ data to each of these new polygons. Then extrude just this layer so as not to have anything overlapping.  And maybe extrude the regional polygons negatively so that it doesn't cover up the bar.  Or I could export the tilted map with the regions shaded to a 2D png/jpg and insert it into a PowerPoint slide and insert each 3D bar graph over the regions with the original excel file. Its easy for my purposes and looks pretty cool for not having an expensive graphics program.  I'll try both and see what goes over better with the boss.  Thanks again. Hi Jessica, You have enough creative ideas to make this work. Maybe you could share the result on this forum to inspire others ... Kind regards, Xander ... View more

Is this in the process of being fixed?  It would be nice to make maps like this http://www.researchgate.net/post/How_do_you_create_this_visual_pseudo-3d_map_representation without having to buy and learn a graphics program. Hi Jessica (and Jonathan), Not sure if this issue was fixed, but with the current functionality you can do some interesting visualizations too. For a project I wanted to show the spatial distribution of ammoniac emission by farm animals in the Netherlands and I choose to do this in ArcScene (10.1 at the time). This was the result (extracted from a presentation): [ATTACH=CONFIG]27190[/ATTACH] What I had was 4 attribute columns with the emission values per square (2.5x2.5 km²). I added this layer 4 times to ArcScene (or copy and paste the layer in the TOC). For the first emission (species) set the Base Height to 0 and the Extrusion to the attribute value itself. For the second emission (species), set the Base Height to the first emission value and the Extrusion to the attribute value itself and so on. The last one starts (Base Height) at the sum of the other emissions... BTW; I just tried to use the stack graph in 10.2 and the data still seems to disappear. Kind regards, Xander ... View more

Hi, I�??d like to normalize map data in raster format to the data range using Raster Calculator. I'm not so much concerned with the accuracy of the method. e.g. v�??(i) = (v(i) - minA) / (maxA - minA) What is the syntax I need to do this? I am getting problems using the following script: ("raster" - min("raster")) / (max("raster") - min("raster")) Thanks! Hi, If I could convince you to do this in the Python Window, it would be a lot easier. In that case you can make use of the Raster object in arcpy. Just have a look at the code below: import arcpy my_raster = arcpy.Raster(r'obj\objecten') result = (my_raster - my_raster.minimum) / (my_raster.maximum - my_raster.minimum) Basically you load the arcpy module, you define a variable "my_raster" to which you will assign the raster you want to process. In this case I used a raster called "objecten" which resides in a grouplayer called "obj", but you can just specify the name of the raster as it is called in your table of content. The last line performs the calculation (make sure you have the Spatial Analyst extension switched on). As you can see you can access the minimum and maximum values of a raster, since they are properties of the raster object. See also: http://resources.arcgis.com/en/help/main/10.1/index.html#//002100000017000000 http://resources.arcgis.com/en/help/main/10.1/index.html#//018z00000051000000 To save the result, you can add another line: result.save('yourNameGoesHere') It will save it in your current workspace. Hope this helps you. Kind regards, Xander ... View more

Hi Yang, I suppose you have experience with ArcGIS 9.x and are used to find the functionality of 3D analyst in the 3D Analyst toolbar. In 9.x you would see something like this: [ATTACH=CONFIG]27167[/ATTACH] All of these tools (and more) can be accessed now through the ArcToolbox: [ATTACH=CONFIG]27168[/ATTACH] This screenshot is from 10.2, but you will see something similar at 10.1. Hope this is what you are looking for... Kind regards, Xander ... View more

Hi Zuoqi Chen, Having a dataset locked can be a difficult thing to solve. IMHO I thing the locking mechanism in ArcGIS has been implemented to drastic, causing more problems than it avoids... I assume you are writing your raster datasets to a file geodatabase. It is possible to test the scema locking before deleting and place it in a "try except" (as shown below), but this doesn't solve anything if you have a lock on your dataset. Always be sure to delete any references to your dataset and delete cursors, rows, etc... if arcpy.TestSchemaLock(data):     try:         arcpy.Delete(data)     except:         # catch error here else:      print "Unable to acquire schema lock" source: http://resources.arcgis.com/en/help/main/10.2/index.html#/TestSchemaLock/018v0000002m000000/ What I have also come across is using the garbage collector. This may sometimes solve the problem. See snippet below:     # at beginning of script     import gc # Garbage Collector     gc.collect()     # at end of script     collected = gc.collect()     print "Garbage collector: collected %d objects." % (collected)     # try to delete now... sources: http://gis.stackexchange.com/questions/19408/arcgis-10-0-python-searchcursor-file-locking http://docs.python.org/2/library/gc.html http://pymotw.com/2/gc/ If this doesn't work, you could go back to the gold old Esri grid data format (write raster to folder), which is probably one of the last formats for which the locking mechanism hasn't been implemented. 😄 This did solve some nasty locking problems I had... ... View more

Hi yang li, Judging from your first screenshot, the 3D extension is already switched on (otherwise the ArcScene and ArcGlobe icons would have been grayed out). To enable the icons (functionality) on the 3D Analyst toolbar, you need to add data to your table of content. When you add a valid raster you will see this: [ATTACH=CONFIG]27159[/ATTACH] Most options are enabled, only the Profile Graph option remains disabled, since it requires a selection of 3D lines. Kind regards, Xander ... View more

Hi Peter, I'm glad that worked out for you. Could you click on the green check mark? See also the note below: Did you post and someone replied with a good answer? Help the community  to find these great posts with answers by clicking the green check mark  to indicate that it has been answered for you. Kind regards, Xander ... View more

Hi Stuart, Below you will find a snippet of Python code that can be used to calculate the maximum M for each feature in a featureclass. It is a standalone script that can be executed from your Python IDE of choice. Before you do, you should change the location/name of the featureclass used and the name of the (existing) field where the result will be written. In short it will: access your featureclass loop through the features calculate per feature the maximum M, using GetMmax() update the feature (row) def main():     import arcpy     # settings, replace with your own dataset and existing column name:     FClassOrTable = r'C:\Project\_tmp\test.gdb\Mawares'     outputField = 'MMax'     # .. or make variable     # FClassOrTable = arcpy.GetParameterAsText(0)     # outputField = arcpy.GetParameterAsText(1)     # update featureclass     fields = ['SHAPE@',outputField]     with arcpy.da.UpdateCursor(FClassOrTable, fields) as cursor:         for row in cursor:             row[1] = GetMmax(row[0])             cursor.updateRow(row)         del row         del cursor def GetMmax(geom):     # get the max(M) from a geometry     mmax = -9999     for part in geom.getPart():         for point in part:             m = point.M             if m > mmax:                 mmax = m     return mmax if __name__ == '__main__':     main() You can also add the script to a toolbox. The tool can eventually be used in modelbuilder in case you want to combine it with additional tools. In case you want to create a tool from it, you should uncomment the "arcpy.GetParameterAsText()" lines and define the parameters properly. first: the location/name of the featureclass (Type: Feature Layer) second: the field that will store the max M (Type: Field, obtained from the featureclass, and apply a fields filter to only show the double or floating fields). Links to additional information: http://resources.arcgis.com/en/help/main/10.1/index.html#/UpdateCursor/018w00000014000000/ Cheers, Xander ... View more