Latest Contributions by RichardFairhurst

Is is possible to create and apply a definition query in ModelBuilder? The model would be the same as in the attached screenshot except that the tool would be "Definition Query" not "Select Layer By Attribute". [ATTACH=CONFIG]26204[/ATTACH] I feel like this should be easy to do. Any help is appreciated. Thanks. pbh Use the Make Feature Layer tool with the query option specified.  That is the definition query of the layer as far as I know.  A layer without a query can feed a layer with one, but not the other way around as far as I know. ... View more

Does the shape file have a spatial reference projection?  If that is corrupt or missing you won't be able to calculate geometry outputs.   Does it have a Geographic (angular units like long/lat degrees) or Projected (linear coordinate units like meters or kilometers) spatial reference? Let us know whether or not the normal field calculator is disabled as well.  If it is disabled then there is a lock of some kind on the file.  Could anyone else have the file open or could you have two copies of it open in separate maps?  Also have you started an edit session or not?  If you started an edit session on a different workspace that might affect the geometry calcuator (although I don't recall that situation preventing it from working).  You are using a shapefile and not an SDE database I am assuming. ... View more

Dear all, In the attribute table of a dbf file I added a field "Area" and wanted to calculate the area value of the features with function "Calculate Geometry". But this function is grey and does not work. What could be the reason for that and what do I have to change? Thanks in advance and kind regards, Matthias Area is a reserved word and field name for geodatabases.  It is therefore rendered as a non-editable field.  You may not even be able to removed that field by deleting it now either.  So add a field like AREA_SQFT and that will work. See this post and specifically this help. ... View more

If I understood correctly, how would including a text file with a path be different than simply hard coding paths in the logic script? The issue at hand is to read files outside the folder where addin is placed. So it would be something like this: say addin is in following folder: 'C:/proj/addin/addin.esriaddin' and I want addin to recognise files in: 'C:/proj/data/raster1' or 'raster2.....' In above instance, data and addin are placed in C drive but I want addin to be configured in a way so that it can read data even if it is placed in D or E or F drive without the need to hard code the root directory path in logic script. This way, user can simply place updated addin in the designated folder irrespective of what drive it is residing in. Please let me know if my answer isn't clear. What you are describing makes no sense.  You have to program a path somewhere that is hard coded for at least the root if you are not allowed to deduce the path from the add-in location. Think about it.  Lets say I want to make phone call to someone in my own area code.  I don't have to tell anyone what the area code is because we can all deduce it from my location. Now lets say I want to dial the same number and not give the area code, but I want the phone to connect with someone in a different area code altogether.  How do you expect the phone company to program their phone system to do that?  Randomly try every area code for me until the call goes through where I wanted?  They could try that approach, but I won't be providing any confirmation that they have connected me to the correct party.  If more than one connection is possible, how will they know they got the right one if I won't give them any feedback?  Do I have a right to expect my phone call to go through under these conditions?  What could do this other than a magical mind reading app?  If I could successfully develop that kind of app I would already have a Nobel prize. This is how you have described your program requirements so far.  Describe any information I actually was expected to provide to the phone company or the kinds of feedback loops and confirmations I was suppose to give them in the scenario above that would make my phone call possible. Fianlly, you do realize that if you are wanting the user to access data on a network drive that you should not use a drive letter at all for the path ever.  Drive letters are dependent on the local computer.  What you must specify is the domain URL to the network server (which is what the drive letter actually points to).  The path:  r'//myserver/files/data/...' will go through on any computer in my network with access permissions to myserver in my domain.  A domain name server can be used to abstract the path further, making the path independent of an actual specific server, so that the domain administrator can transfer data and repoint all users to a new server keeping the same path for his end users without touching any local machines. ... View more

Hello, I have a table with Expiration Dates.  I want to use the select by attributes to tell me what dates are between now and exactly one year.  Thanks in advance for any help. Mark I am going to assume you are using a File geodatabase, since I am most familiar with that database and you did not say what database you actually are using.  All SQL help is dependent on the database you are working with, so until you provide that information you really can't get help that deals with your specific needs. This expression will select all records based on the current date and within the last year when time of day is not a factor (accounting for leap year): (EXTRACT(MONTH FROM CURRENT_DATE) = EXTRACT(MONTH FROM CURRENT_DATE - 365) AND EXTRACT(DAY FROM CURRENT_DATE) = EXTRACT(DAY FROM CURRENT_DATE - 365) AND "CREATED" <= CURRENT_DATE AND "CREATED" > CURRENT_DATE - 365) OR (EXTRACT(MONTH FROM CURRENT_DATE) = EXTRACT(MONTH FROM CURRENT_DATE - 366) AND EXTRACT(DAY FROM CURRENT_DATE) = EXTRACT(DAY FROM CURRENT_DATE - 366) AND "CREATED" <= CURRENT_DATE AND "CREATED" > CURRENT_DATE - 366) OR (EXTRACT(MONTH FROM CURRENT_DATE) = 2 AND EXTRACT(DAY FROM CURRENT_DATE) = 29 AND "CREATED" <= CURRENT_DATE AND "CREATED" > CURRENT_DATE - 366) This expression will select all records based on the current date and time stamp to immediately after the exact same time on the same date 1 year ago (accounting for leap year): (EXTRACT(MONTH FROM CURRENT_TIMESTAMP) = EXTRACT(MONTH FROM CURRENT_TIMESTAMP - 365) AND EXTRACT(DAY FROM CURRENT_TIMESTAMP) = EXTRACT(DAY FROM CURRENT_TIMESTAMP - 365) AND "CREATED" <= CURRENT_TIMESTAMP AND "CREATED" > CURRENT_TIMESTAMP - 365) OR (EXTRACT(MONTH FROM CURRENT_TIMESTAMP) = EXTRACT(MONTH FROM CURRENT_TIMESTAMP - 366) AND EXTRACT(DAY FROM CURRENT_TIMESTAMP) = EXTRACT(DAY FROM CURRENT_TIMESTAMP - 366) AND "CREATED" <= CURRENT_TIMESTAMP AND "CREATED" > CURRENT_TIMESTAMP - 366) OR (EXTRACT(MONTH FROM CURRENT_TIMESTAMP) = 2 AND EXTRACT(DAY FROM CURRENT_TIMESTAMP) = 29 AND "CREATED" <= CURRENT_TIMESTAMP AND "CREATED" > CURRENT_TIMESTAMP - 366) I am using matching date logic for leap years rather than end of the month logic so that 2/28/2013 looks back to 2/28/2012.  I have not worked out how to get end of the month logic to work so that you go from 2/28/2013 back to 2/29/2012. Speaking to the programmer for the file geodatabase at the ESRI UC this year I found out that file geodatabase SQL complies with all by a few SQL 92 specifications, so you can look up help on the internet for the SQL 92 syntax and normally apply it to the file geodatabase.  That is how I found out that CURRENT_TIMESTAMP works with a file geodatabase. ... View more

Is there a way to have a feature's symbol change based on the scale of the map? For example, at 1:250,000 I would prefer for my sewer main layer to be symbolized on main-type only, but when the user adjusts the scale to say 1:10,000, it would be great if the direction of flow would be added to the symbol (but not before). The only work around I have found to adjusting symbol detail is to add the layer in more than one time and adjust the symbol with scale ranges. This gets the job done, but is a bit messy in the TOC. Can Symbol Level drawing be used to control this? Unlike Label classes, the symbology dialog does not provide any way to define several standard layer symbology classes within a single layer that are scale dependent.  I believe that Representations can have rules set up to do behaviors that may resemble the behavior you want, but I have not used representations enough to be absolutely certain.  Representations have to be built with a particular target scale in mind, so they do not have the same level of flexibility that standard symbols have.  I have always created multiple copies of the same layer to make different scale dependent standard symbol sets and typically included them all in the same group layer to simplify their TOC presentation. ... View more

I just started working with models so I don't have much experience. I'm trying to take a database with multiple line feature classes and merge them into one final feature class.  I made this model but it isn't working.  I'm ending up with a new feature class with only 1 line in it instead of about 6 lines. [ATTACH=CONFIG]25978[/ATTACH] With the merge tool you don't use an iterator.  You just include all 5 feature classes as inputs to the tool assuming they are all the same geometry type and run it.  Every feature from those 5 inputs will be combined into a single output feature classes. ... View more

I need to end up with a series of line segments that are equal to the original line feature. I was looking to split the line at the midpoint between every two vertices(vertices every 0.1 feet) If that isn't possible I was wondering if I could split the line every 0.05 feet. I would then merge every 2 segments into 1 segment to get the result I need. Split the line at every vertex then apply the method I specified to get the midpoints of each of these split segments exactly midway between the original vertices. Extract those points as I described previously. Then you can split the original line on those points and except for the end vertices of the line (or any 2 point line) a single vertex of the original line will now be somewhere in the center of these newly split lines. That seems to be what you are describing as the result of the splitting and merging that you want. The ends of the line will be connected to split line segments that are 1/2 the distance to the original vertex that was closest to them. So if the original line looked like this (where red periods are line ends and black periods are pseudo-node vertices). .______.____.______.____. The end result of the process would look like this (assume all periods have no actual width so the set of 5 lines below are in fact the exact same total length as the single line above): .___.___.__.__.___.___.__.__. If this is not the result you want, use the above method of displaying lines and vertices to depict how the first line shown should look at the end of the process you really want to accomplish on that line. ... View more

I am wanting to be able to, in an automated fashion through a Model or other means, rename a field in an attribute table once the dataset is created.  I have read in earlier versions of ArcMap, the only way to "rename" a field would be to add the new field, copy over the data, and delete the old field, and this functionality was added to 10.x whereby you can manually rename a field using ArcCatalog.  However, I'm wanting to rename a field within a Model. There's an older article at http://ideas.arcgis.com/apex/ideaView?id=0873000000087wgAAA which indicates this functionality at one time did not exist, but now that you can do it manually, did Esri implement the ability to do it using some tool or other means? This is the only help I see on the subject and it only applies to manually renaming a field through ArcCatalog properties. ... View more

The lines that built expression and expression2 were messed up (missing escape character to make single quote around your address and street literal) and expression3 badly written.  They are corrected below. import arcpy import os import re basedir = r"C:\Test" fc = r"G:\GIS\ParcelsPublic.shp" field1 = "ADD1" field2 = "STREETNAME" try:     for fn in os.listdir(basedir):         add = re.sub(r"\D"," ", fn)         add = str.lstrip(add)         add = str.upper(add)         street = re.sub(r"[^a-zA-Z]"," ",fn)         street = str.lstrip(street)         street = str.upper(street)         print (add + " " + street)         expression = '"' + field1 + '" = ' + add         expression2 = '"' + field2 + '" LIKE \'%' + street + '%\''         expression3 = expression + " AND " + expression2         rows = arcpy.SearchCursor(fc, fields="STREETNAME; ARN; ADD1", where_clause=expression3) I assumed ADD1 is a numeric field?  If it is a string field rewrite the line that builds expression as:         expression = '"' + field1 + '" = \'' + add + '\'' And if you want the LIKE syntax applied to the address string use:         expression = '"' + field1 + '" = \'' + add + '%\'' ... View more

I'm looking for a way to split a line at the midpoint between two vertices.  Or if I have a series of line segments, is there anyway to go down the line and merge every 2 segments together?  I have a series of slope measurements on a road.  Just splitting the segment at each vertex does not accurately reflect the data which was collected at each vertex.  The vertex needs to be the midpoint of each line segment.  It can be a script or a program, etc as long as it does the job.  Thanks. Create two double fields called X_LINE_MID and Y_LINE_MID on the lines.  Use the Geometry calculator to calculate the X and Y MIdpoint values of the lines into those fields.  Open the tableview of the lines and export that to a table (not a feature class).  Use the Make X/Y Event Layer tool to create points from the Midpoint X and Y coordinates of the exported table.  Export the layer its context menu in the table of contents to permanently covert the event table layer into an actual point Feature Class.  That should be what you want and each point will have all of the attributes of the lines they came from. ... View more

Hi, I was hoping that someone might be able to help me with a problem I'm having.  I'm using ArcGIS10.1 for Desktop with an ArcView license. I have two point shapefile datasets (I've called them 'long' and 'short') of ship positions for the month.  Each shapefile has the same fields but different amount of data.  The fields are more complete in the long shapefile but this only has positional data for every 15 minutes.  The short shapefile has positional information every 3 minutes but with fewer of the fields complete.  I want to use the short shapefile as there are more points over time but I have found out that there is some missing data points in the short shapefile (database issues which I have no control over!) which are in the long shapefile.  My first thought was merging the data but this also incorporated points occupying the same spatial and temporal area (identical ships).  I also tried selecting by location but this threw up problems with two points occupying the same spatial location on different days (I need to keep this information). So, my question is; is it possible to select by location and attribute at the same time and for two shapefiles? I was thinking something along the lines of select when two points are within the same source layer feature and the shipnames are identical and the timstamps are identical. Any help would be greatly appreciated. Kind regards, Richard Use Spatial Join with an appropriate tolerance and then correct the result by selecting for matching attributes that are joined together.  The more dense point set should be the target and the less dense points set should be the Join.  Use the One to Many option to create all match up points and then filter where the attribute matches make sense.   Now you can do an attribute selection that also takes account of the spatial relationship of the points. After selecting the matches on attributes, export that set so that you hopefully have a one-to-one match back to the dense point set.  Do a standard join of the original dense point set with this spatial join matched and exported dense point set to transfer attributes picked up from the low density points that have correct attributes using the field calculator and then break the join.  That should be all you need to do. ... View more

My formula used the outdated 9.2 format, which does not work at 9.3 or later in the regular field calculator.  To update the calculation for python script and the field calculator the XY concatenation calculation to a point feature class should be: # Process: Calculate X_Y_Link Field (2)... arcpy.CalculateField_management(Output_Feature_Class, "X_Y_LINK", "Output(!SHAPE.CENTROID!)", "PYTHON_9.3", "def Output(FirstPoint):\\n FPX = round(float(FirstPoint.X), 8)\\n FPY = round(float(FirstPoint.Y), 8)\\n return \"{%(FX)013.8f}{%(FY)013.8f}\" % {'FX': FPX, 'FY': FPY}") In the normal 9.3 and later field calculator you would output to a 30 characer field with the calculation entered as shown below: Parser: Python Show Codeblock: Checked Pre-Logic Codeblock: def Output(FirstPoint):   FPX = round(float(FirstPoint.X), 8)   FPY = round(float(FirstPoint.Y), 8)   return "{%(FX)013.8f}{%(FY)013.8f}" % {'FX': FPX, 'FY': FPY} Expression: Output(!SHAPE.CENTROID!) ... View more