I'm still a little confused about your end goal (polygonize all raster cells within a certain distance of the purple polygon?), but I think the Extraction toolbox in Spatial Analyst might be what you are looking for. You can either use the purple polygon as is, or buffer it (ex. "purple polygon + 100m around it"), and then use it as the input to Extract by ...
You can download the .style file at this Git repo: https://github.com/ISUEOGTP/GISTaskSheet mentioned in these PDF instructions: https://indicator.extension.iastate.edu/projects/tasksheets/reports/PM2082-14e.pdf
These instructions also explain how to import it once you have it.
You shouldn't need to use "Add Geometry Attributes" but just need to add a numeric field and calculate the value with "Calculate Geometry" directly in the table. For the 000229 error you might want to try to disable background geoprocessing as suggested in this ESRI technical article 000012177.
It sounds like (from your description) that you have some ...
In this case, in the Calculate Value tool...
The 'Code Block' field should be empty
The 'Expression' field should be: "%value%"[0:6]
(The expression should include the quotation marks.)
The 'Expression' is a simple one-liner Python expression that returns the result.
The 'Code Block' is an optional more complex Python script that will be ...
The default location for package install/download is:
If you browse to that directory in your file explorer you should find the attachments inside of the layer package sub-directory.
Example: If your layer package name was TestPackage the additional files would be located here: "C:\Users\xxxx\Documents\ArcGIS\...
I tried the below script with help of @Hornbydd script and it works in my case:
ws = arcpy.env.workspace = r'D:\BRB Snow Cover\Sample\Terra'
Terra_Rasters = arcpy.ListRasters("MOD*")
matches = 1
for item1 in Terra_Rasters:
for item2 in Terra_Rasters:
item1_num = int(item1[31:38])
item2_num = int(item2[31:38])
That is expected behavior for annotation. Since you have zoomed your data frame out from the annotation's reference scale, the fonts should be smaller. If you zoom the data frame into 1:15,000, they will be huge. You can try overriding the symbology on the annotation layer's properties using the "Substitute individual symbols in the symbol collection" ...
This is not a pure GIS question it's more about dealing with python in a GIS context. The following code should help you develop your code further:
r1 = r"c:\temp\2010001.tif"
r2 = r"c:\temp\2010002.tif"
r1_num = int(r1[12:15]) # returns 1
r2_num = int(r2[12:15]) # returns 2
if (r2_num - 1) == r1_num:
print "r2 is 1 day ahead of r1 do something"
If you want to create a new layer based on a subset of an existing layer you need to export after doing your selection (right click on your selected layer then Data>Export Data... and choose export selected data) this will create a new layer with a copy of the selected data. in this case you end up with two copy of your data and any edit made in one copy ...
The tools aren't trying to talk to a URL on the internet.
When you launch the tools, ArcMap is actually constructing tool pages in HTML (this is an incredibly simplified explanation). Anyways, the point is it makes use of ActiveX and HTML to show these tool dialogues, as such your Internet Settings on the machine come into play.
There are a lot of Esri ...
Give this a try (from a GeoNet post)
Open internet explorer and got to the tools, click on internet options, then Connections then LAN Settings. Uncheck the option that says Automatically detect settings. Basically the problem was that the basemap was searching for the host site through all my local networks and unchecking this option bypasses that.
Had exactly the same problem as you are having. It was as easy as clearing the image cache.
Customise > ArcMap Options > Display Cache > Clear Cache
In the last Windows upgrade the month it renamed your user name in AppData to XXXXX.000 which hides the added layer files from ESRI.
Well I had to go back and keep searching...I had the old external hard drive I have saved over the years...
I looked at the chrome bookmard I save and found this ...
This is what I have been searching since I move here from New Mexico.
The link to that one is : https://www.esri.com/arcgis-blog/products/product/mapping/can-you-read-me-now/
You've pretty much answered your own question. There are a couple different types of centroids: the geometric centroid (which is not necessarily within the polygon) and the center of mass. The spatial join and select by location tool are only capable of using the former. The only way to get the center of mass centroid, using the default tools, is the ...
It is critical that your map/datalayers are in a projected coordinate system when calculating area or distances. If the same projection is used consistently, you should get consistent answers when calculating geometry.
If you are going to discard the area that is flooded within the pond, you can use the Erase tool instead of the Union. Erase needs an ...
If you have follow the steps mention in Agent Analyst tool interface will not open on ArcGIS for Desktop? and then also you are getting same problem. This is due to Refresh.exe.config file. Your system will be taking as a text file. Download this file and paste it in C:\Repast 3\Agent Analyst\Refresh location. Open this file with notepad and change the ...
I found the same thing! The equation in the publication simply expresses the standard cost surface function, so I couldn't see how they converted it to depth. My maths isn't the best, but I remembered that the height of a right angled triangle could be calculated from its slope and distance (SOHCAHTOA and all that).
Height = cell resolution * Tangent of ...
The equation to linearly scale 100-0 to 1-10 is:
-0.09x + 10
The equation to linearly scale 200-500 to 1-10 is:
0.03x - 5
You could use nested Con functions in the raster calculator to apply these equations:
Con("distance" <= 100,
-0.09 * "distance" + 10,
Con("distance" < 200,
Con("distance" <= 500,
Multiple solutions possible to get your a very close result. If you can create 0.7 contour of raster, points on it are your best candidates.
Alternatively express your area in raster pixels and run focal statistics using relevant size window, e.g. 50*50 (not necessarily a square or rectangle) if your area is equivalent of 2500 cells. Use raster calculator:
Assuming I understand correctly, I would use a buffer instead using half the distance and dissolving adjacent buffers (if you buffer each polygon by half the distance when two or more polygons are near another the sum will equal the full distance). Then, use a spatial join on your buffer using the input polygons to get the count of the number of input ...
Some term this is as girth and the University of Connecticut, Center for Land Use Education & Research, has developed a free geoprocessing tool (Shape Metrics) to calculate this among many other shape metrics. This is mentioned in the link as an index but it is quite easy to find the radius from the formulae given. You can download it from the link, ...
As explained in comments, and demonstrated in the other answer (but not fully explained there), the solution was to make sure that each line is created with the same spatial reference as the polygons.
(I had actually tried this earlier, and got an error stating that spatial_reference is a read-only attribute. Turns out I had specified it incorrectly, and ...
Generate Points Along Line geoprocessing tool:
Creates point features along lines or polygons at fixed intervals or by percentage
Construct Points tool on the editing toolbar:
Creates new point features at intervals along a selected line.
Another option (thanks @mkennedy) is Create Fishnet geoprocessing tool with the label points option set to true:
Method "within" works fine, because this script
import itertools as itt
with arcpy.da.SearchCursor(pgons,"Shape@") as cursor:
for row in cursor:
I've developed a native iOS app using GDAL to manage the vector data, and OpenLayers (in a UIWebView) to display the layers (both raster and vector).
For vector layers, the app uses GeoJSON to communicate between OpenLayers and GDAL. It can use WMTS/TMS/WMS for raster layers.
There is a bit of a learning curve for both GDAL and for OpenLayers. But the ...
QGIS is out, it doesn't support raster attribute tables at all.
I don't know about GRASS.
GDAL would require you to (manually or in python code) handcraft VRTs with lookup tables from the attribute table fields for each "band", then composite the VRTs.
In ArcGIS, you can only use Lookup and Composite Bands. However, look at the versions of those functions ...
You can use your field calculator for this. Add a new text field to your table that will be your key field. Right-click on its field name in the attribute table and choose Field Calculator...
Use [ColA] & [ColB] & [ColC] for your field calculation.
Map packages are just zip files with a different extension. Use 7zip or whatever to unpack the package somewhere, delete the internal .gdb that hold the in-common data across all MXD, then repack the package.
When the MPK is delivered and opened the MXD will have broken layers; use the standard "Repair data sources" tools/methods to point to the shared in-...
More information about the comment made by @Tim Sexton can be found here:
This is an Esri technical support article regarding:
How To: Clear or delete trusted storage for concurrent-use and single-use licenses
As commented by @mkennedy a guess might be:
north pole-centered stereographic with a central meridian/longitude of origin around -20
Otherwise I think you will be best to try and contact the author(s) of the maps.
Sounds like the output of the Table To Excel tool is set to Intermediate.
Edit your model and right click on the green output bubble and verify it is not set to Intermediate. If it is, then un-tick it as ModelBuilder will delete it at the end of the model run, which would be the expected behaviour as it is set as Intermediate.
Unfortunately, I do not have an Advanced ArcMap license, so I could not use the Angle option in the Generate Near Table tool. However, I did find a work around, which (sorry for the major delay) I'm posting here in case someone else finds themselves stuck in the same situation. Here are the steps:
Use the Bearing Distance to Line tool to create a polyline ...
You can accomplish this using a combination of a Grid Index and Data Driven Pages to divide your country into a grid, and a Multiple Ring Buffer to provide a white background outside the country's border.
The following assumes that you have a polygon feature describing the country's border. The one I've used here is not a 100% perfect fit, but it will give ...
Try raster calculator, for example:
("RGBraster@1" > "RGBraster@2")*("RGBraster@1" > "RGBraster@3")*1
("RGBraster@2" > "RGBraster@1")*("RGBraster@2" > "RGBraster@3")*2
("RGBraster@3" > "RGBraster@1")*("RGBraster@3" > "RGBraster@2")*3
Should evaluate to 1 for red, 2 for blue and 3 for green. (When a condition is true it will return 1, ...
I would also like to see an answer to this questions using either QGIS or ArcGIS GUI's. Meanwhile I am posting a Python approach to this problem because both QGIS and ArcGIS provide a way to read a raster into a NumPy array.
If you are using QGIS, you already have GDAL installed. With the following code you can read the raster as an array:
Use the 'Dissolve' tool.
The Merge tool is for merging multiple feature classes into one new feature class, but still keeping the existing features separate.
The Dissolve tool is for merging features within a feature class into either a single feature, or multiple features based on same attribute values. If the features to be merged are not touching or ...
In ArcGIS Desktop 10.x you use arcpy.mapping.Layer
lyr = arcpy.mapping.Layer(r"path\to\file.lyr")
In ArcGIS Pro you use arcpy.mp.LayerFile:
lyrFile = arcpy.mp.LayerFile(r"path\to\file.lyrx")
for lyr in lyrFile.listLayers():
It seems to me that you are trying to load an AddIn compiled for ArcGIS 10.5 in ArcGIS 10.4. This is not supported. You'll need to update your ArcGIS version to ArcGIS 10.5 or higher. Higher ArcGIS versions normaly can handle older AddIn versions.
Heavy handed solution, no scripting:
Spatial join of counties one to many to your culvert lines,
i.e.pages will result in 1800 polygons, each having page name.
Add polygons to reference dataframe and apply DDP definition query
(match) on this layer. I called it MANY_POLYGONS.
Use properties of reference dataframe as follows:
Works like a charm.
Intersect the points and fishnet polygons. Then one way is to use python and the pandas module to count unique point ids per fishnet id.
import pandas as pd
#Change these four lines to match your data
inter_output = r'C:\data.gdb\intersect123'
fishnet_id_field = 'FID_fishnet'
vessel_id_field = 'VesselID'
adjp_field = 'AdjPounds'
#Create a ...
Never tried it before, but I guess it should be feasible this way:
Get the index layer's row object for the active or current page with pageRow property of DataDrivenPages object
Get the county name in the row (see here to understand what is a row and how to use it)
Select the county in the couty layer using the extracted couty name from the index layer's ...
After converting the LAS Dataset to a point feature class, then add x, y and z coordinates to the TOC with Add Geometry Attributes tool:
Adds new attribute fields to the input features representing the spatial or geometric characteristics and location of each feature, such as length or area and x-, y-, z-, and m-coordinates.
Another option is to use the ...
Other software creating polygons (so outside the ArcGIS system) cannot guarantee to conform to Shapefile specifications. What you are observing is almost certainly the result of malformed geometry, for example the vertices order is in the wrong direction. As @user2856 suggests run a repair geometry tool on the data first, this should fix it.
You can use the Spatial Join tool with the points as the target features and the polygons as the join features. You should leave the keep all target features option checked if you want to include all points in the output, even if some fall outside of any of the polygons. If you want the output to include the polygon object ID in a JoinFID field you should ...
You should be able to use the Conditional tool, either standalone or in Raster Calculator. Set your conditional statement so that wherever Mask=Water will be 1, and everywhere else will be 0. Something like this:
Con(Water=Mask, 1, 0)
This will produce a binary raster marking the waterline. If you need to know the height of the waterline you could then use ...