|
POST
|
Just to add a little comment to the solution Dan Patterson provided. There is a difference between integer and floating rasters. The integer rasters have an attribute table (so it has attribute fields and allows you to construct a condition) while the floating (decimal point) rasters do not. Hence the empty query builder screen. Way back in the older versions of ArcGIS, one needed to add the spaces around floating rasters and the operators when composing map algebra. I think this is an example of what has not been corrected (yet).
... View more
11-29-2014
09:42 AM
|
1
|
1
|
2495
|
|
POST
|
Vince Angelo is absolutely right! (+1 for that) It is better to get things right, create a structure that is durable and enables analysis and statistics rather than start with all types of customization to get results you could obtain without customization if your schema was set up correctly. Invest now in setting up your data correctly, the ROI will be high.
... View more
11-29-2014
08:10 AM
|
0
|
0
|
2828
|
|
POST
|
The video that I included in this thread (question): Use Briefing Book offline on an iPad includes a map with the time slider. I wonder if this remains a prototype or is coming to us soon. Kind regards, Xander
... View more
11-28-2014
05:05 PM
|
0
|
2
|
1316
|
|
POST
|
About a year ago Esri APL (Esri Applications Prototype Lab) published this video: ... showing a demo at 8:28 of a native app consuming a briefing book "iPad Native Application consuming the newly Authored book". Is this possible with the Briefing Book application that is available now or is this still a prototype? Thanx in advance, Xander
... View more
11-28-2014
05:02 PM
|
0
|
0
|
3157
|
|
POST
|
Hi Gary Crowley, Thanx for the virtual beer! For now I think it's best to mark the answer as correct so other people may find the solution more easy.
... View more
11-28-2014
04:21 PM
|
1
|
0
|
3229
|
|
POST
|
Thanks Dan Patterson for explaining the reason for checking the pnt object.
... View more
11-28-2014
04:19 PM
|
0
|
0
|
3229
|
|
POST
|
If you have a look at: Spatial Reference List -- Spatial Reference you will see that the Well Know ID does not exist. I think the list is not correct. When I search for the description in ArcGIS Desktop I get WKID = 103129. Give it a try.
... View more
11-28-2014
12:28 PM
|
0
|
0
|
1101
|
|
POST
|
You could try something like this: import arcpy
import os
# specify your folder with the .sde connection files
sde_folder = r"C:\Users\xbakker\AppData\Roaming\ESRI\Desktop10.2\ArcCatalog"
arcpy.env.workspace = sde_folder
# list the sde files
sdes = arcpy.ListWorkspaces("*", "SDE")
# define featureclass name and field name
fcname = "state_forest_area"
fldname = "total_acres"
flds = (fldname)
total_acres = 0
# loop through the sde files
for sde in sdes:
try:
sde_path, sde_name = os.path.split(sde)
fc = os.path.join(sde, fcname) # assuming data is in sde root and not in a feature dataset
lst_acres = [r[0] for r in arcpy.da.SearchCursor(fc, flds)]
county_acres = sum(lst_acres)
total_acres += county_acres
print "County '{0}' has {1} acres '{2}'".format(sde_name[:-4], county_acres, fcname)
except Exception as e:
print "Error: {0}".format(e)
print "check existence of the featureclass and field in the current sde connection"
print "Total acres: {0}".format(total_acres) Kind regards, Xander
... View more
11-28-2014
10:17 AM
|
0
|
0
|
2828
|
|
POST
|
I'm afraid there is not much support for listing domains without the data access (da) module. You may have to use one of the other suggestions (like the one by Jake) directly on the database.
... View more
11-28-2014
07:18 AM
|
1
|
0
|
1436
|
|
POST
|
Statistically the amount of a class changed can be derived from the count fields in both attributes tables shown in your screen shots. If you want to see which pixels have changed you can use the "Not Equal" tool (in Spatial Analyst Tools\Math\Logical). This will generate a result where all pixels that have changed will have value 1 and the ones that remained the same will have value 0. Maybe you are interested to see the change in each class (which amount of non-forest changed to upland-forest or mangrove). In that case you can use the "Combine" tool (in Spatial Analyst Tools\Local) which will produce a raster that contains the a unique value for each combination of input and output values. Using the attribute table of the result will provide you the information on switched between classes and the raster itself allows you to identify each combination spatially. Kind regards, Xander
... View more
11-28-2014
06:51 AM
|
0
|
0
|
815
|
|
POST
|
Just to emphasize... it is technically possible, but from a database design point of few, you shouldn't use this.. unless the application you use it in (AGOL and its related products) requires such a strategy since you're trying to do something that is not supported by the product, but functional from a users perspective.
... View more
11-27-2014
08:29 PM
|
0
|
0
|
1700
|
|
POST
|
Since the line is long and the width limited here is what I was referring to (check the how the statement ends): get_dtype="[('YEAR', '<f8'), ... , ('RANKPRODSUBC', '<f8')]"
... View more
11-27-2014
08:24 PM
|
0
|
0
|
2053
|
|
POST
|
With standard tools I agree with Dan Patterson. You could achieve this with some Python or ArcObjects programming. It would for instance involve a Near function to determine for each point which point is nearest. If you would connect each point with its nearest neighbor, you could intersect that line with the boundary of the polygon. The intersection would be the point you are looking for.
... View more
11-27-2014
08:16 PM
|
2
|
0
|
1517
|
|
POST
|
OK, to throw in an example using 10.0 syntax (although in run it on 10.2.2) would be: import arcpy
fc = r"C:\Forum\EditVertices\test.gdb\polygons"
oid = 2 # objectid to change geometry
# Identify the required fields
desc = arcpy.Describe(fc)
shapefieldname = desc.ShapeFieldName
oidfieldname= desc.OIDFieldName
# create a where clause
where = "{0} = {1}".format(arcpy.AddFieldDelimiters(fc, oidfieldname), oid)
# Create update cursor
rows = arcpy.UpdateCursor(fc, where_clause=where)
# Enter for loop for each feature/row
for row in rows:
# Create the geometry object
feat = row.getValue(shapefieldname)
# create an array for the updated feature
arr_feat = arcpy.Array()
# Step through each part of the feature
partnum = 0
for part in feat:
# create an array for the part
arr_part = arcpy.Array()
# Step through each vertex in the feature
for pnt in feat.getPart(partnum):
if pnt:
x = pnt.X
y = pnt.Y
# do something with the point coordinates
x += 5
y += 5
# add the point to the part array
arr_part.add(arcpy.Point(x, y))
else:
# If pnt is None, this represents an interior ring
pass
# add part to feature array
arr_feat.add(arr_part)
partnum += 1
# create polygon
polygon = arcpy.Polygon(arr_feat)
# set polygon to shapefield and update row
row.setValue(shapefieldname, polygon)
rows.updateRow(row) What happens in the code is that the polygon with OBJECTID is moved (x+5 and y+5). See image below: If you want to edit vertices of a polygon individually, make sure that the first and last point coincide to force closure of the polygon. Kind regards, Xander
... View more
11-27-2014
08:04 PM
|
2
|
0
|
3229
|
|
POST
|
LiDAR data is a pointcloud, from which you can derive raster data. LiDAR data is normally provided as LAS files. You do not want to store the point clouds in your SDE database. In ArcGIS you work with LAS data using LAS datasets. Some links: What is lidar: ArcGIS Help (10.2, 10.2.1, and 10.2.2) A quick tour of lidar in ArcGIS: ArcGIS Help (10.2, 10.2.1, and 10.2.2) Creating a LAS Dataset: ArcGIS Help (10.2, 10.2.1, and 10.2.2)
... View more
11-27-2014
09:41 AM
|
0
|
1
|
2624
|
| Title | Kudos | Posted |
|---|---|---|
| 1 | 01-09-2020 09:26 AM | |
| 6 | 12-20-2019 08:41 AM | |
| 1 | 01-21-2020 07:21 AM | |
| 2 | 01-30-2020 12:46 PM | |
| 1 | 05-30-2019 08:24 AM |
| Online Status |
Offline
|
| Date Last Visited |
yesterday
|