38
You can use rgdal to access feature classes in Esri file geodatabases.
require(rgdal)
# The input file geodatabase
fgdb <- "C:/path/to/your/filegeodatabase.gdb"
# List all feature classes in a file geodatabase
subset(ogrDrivers(), grepl("GDB", name))
fc_list <- ogrListLayers(fgdb)
print(fc_list)
# Read the feature class
fc <- readOGR(dsn=fgdb,...
32
Vector data can have properties stored at two levels. File level or Layer level.
For feature class (a type of vector data), file level information is stored in a geodatabase. This consists of properties like feature class geometry, attribute table, assigned subtypes and domains, a coordinate system, etc:. All these file level properties make up the actual ...
28
This works without needing to remove the layer and recreate a new one as described above:
geojson_layer.eachLayer(function (layer) {
if(layer.feature.properties.NAME == 'feature 1') {
layer.setStyle({fillColor :'blue'})
}
});
It seems to be quite a bit more efficient than removing and recreating the geoJson layer.
From the docs, a GeoJSON ...
24
You just gotta use the LIKE operator.
For example,
"TITLE" LIKE '%MINERAL CLAIM%'
The % symbol acts like a wildcard.
LIKE is case-sensitive, whereas ILIKE is not.
19
I have written down a small code to style specific geojson feature using leaflet. you can try it on JSFiddle (Original, non-functional), Functional JSFiddle 2018-02-17, or use the following code test locally.
For this example i am using us-states.json files but it can be used for any geojson file.
I hope it will help.
Here is the code:
<!DOCTYPE html&...
17
I know this is old but I was looking for the same as I don't have ArcInfo for the FeatureVerticesToPoints tools. After using Search cursor solution above I went forward to simplify the code and found that using NumPy Arrays in the Data Access Module a simple and very quick script could be produced. I'm using this as a script tool.
Note: The key is the ...
16
The Feature Class to Feature Class tool is a Conversion tool. The idea is to convert a feature class in one format (e.g. shapefile) to another (e.g. geodatabase). It also allows the user to control the Field Map - i.e. which fields will be copied over.
The Copy Features tool simply creates an exact copy of the input features. Note - this tool can also be ...
16
Answering your two questions:
You can change your feature values from the layer object, no need to access the dataProvider().
Yes, you can use the iterator in a for loop.
Check the code below:
layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line')
layer = layers[0]
it = layer.getFeatures()
layer.startEditing()
for feat in it:
layer....
16
Here you find the map including the legend. The legend is in German, but I am German, so I tried to translate it as precise as possible.
I think (and please @Giancarlo correct me, if you disagree):
1 is "Reichsstraße, ausgebaut". That means a well-developed street.
2 is propably a sort of drainage system or stream (thank you @AndreJ). It is not named in ...
15
you can use the export (multiple).
In ArcCatalog...
Just right click on the database and choose "export to shapefile (multiple)"
If you do it at the database level
you will see everything in the database in the tool listed (there is a remove button if there are a few you don't want).
If you do it at the feature dataset level
you see everything in ...
15
The line where you create and bind your popup should have been included in the onEachFeature option of your L.geoJSON factory.
var layerGroup = L.geoJSON(data, {
onEachFeature: function (feature, layer) {
layer.bindPopup('<h1>'+feature.properties.f1+'</h1><p>name: '+feature.properties.f2+'</p>');
}
}).addTo(map);
Demo: http:/...
14
Which part of the script is actually taking up most of the time? There are about 5 other steps going on before you actually start deleting stuff.
You might want to break your script down into bite-sized tests. For example, instead of creating a temporary connection file, listing a bunch of datasets, listing their contents, counting their records, and then ...
14
Sure it is. Depending on what you need:
zoom to whole layer
right click on layer in Layers list and first item is "Zoom to layer"
zoom to selected feature(s)
select feature in map and click on icon "Zoom to selection" in the main toolbar
or right click on layer in layers list and choose "Open Attribute Table" (there is also icon in the main toolbar/...
14
The drawn feature are added to the layer named vector. See draw interactions constructors.
// Get the array of features
var features = vector.getSource().getFeatures();
// Go through this array and get coordinates of their geometry.
features.forEach(function(feature) {
console.log(feature.getGeometry().getCoordinates());
});
13
These terms specifically relate to the ESRI software and data schemas.
A feature class is a set of common geographic features, for example roads, or parcels, or land use, as they are stored in a geodatabase or shapefile, or other storage format. Here is the help file on Feature Class Basics.
A feature layer is the representation of a feature class after ...
13
OpenLayers 2:
You can unselect features with unselect and unselectAll methods of SelectFeature control:
selectControl = new OpenLayers.Control.SelectFeature(vectorLayer);
...
map.addControls([selectControl]);
selectControl.activate();
// unselect any specific feature...
selectControl.unselect(vectorLayer.features[0]);
// ...or all features
selectControl....
12
Notwithstanding whether there are non-geographical features, I think of a feature as being something with geometry (x,y etc) and attributes i.e. a feature is more than geometry. Then when a feature (as part of a feature class) is displayed on a map it also has symbology (points drawn in red dots of 10pts, etc) and can be referred to as part of a layer.
In ...
11
Deleting all but one feature from a set of features with the same ID, keeping one with largest area?
Here's how I would do it:
Make sure your data are in an appropriate projected coordinate system.
Add a new Double field to the attribute table called Area. Run field calculator on this new field with the expression: !shape.area! and Python as the parser.
Run Summary Statistics, performing COUNT on the Name field, and MAX on the Area field we added in step 2....
11
If you want a solution that does not involve any extra extensions or "high-grade" licences you can try this:
Union A1 and A2 to make A3.
Then select by location where A3 does not have its center in A2 (you may need to select A3 where it DOES have its centre in A2 and then switch selection). OR select by attributes as Union will append attributes from A1 ...
10
You can use the Erase Tool in ArcToolBox. According to ESRI website :
"The Erase tool is used to perform overlay analysis on feature classes. This tool creates a feature class from those features or portions of features outside the erase feature class."
Have a look at this webpage, they have really nice explanations.
To do an "Erase" using a Basic level ...
10
Officially, .shp and .dbf files are limited to 2GB. Using GDAL Shapefile driver this can be exceeded, but as stated on this site, this is not recommended due to compatibility issues (See point #Size Issues).
As @Gabriel C. already mentioned, geopackage is a good alternative. You could also use a SpatialLite or SQLite database.
9
In simple terms ...
I try to think of a feature class being a spatial dataset, i.e. not much more than x,y coords and attributes on disk, while a layer applies symbology (and many other lesser known properties) to a spatial dataset.
A layer does not store data, just symbology and a link to where the data is located on a disk somewhere.
Feature classes ...
9
If upgrading to ArcGIS 10.1 (now released) is an option then I just found this in the What's New in ArcGIS 10.1 PDF:
New tool to delete all rows from a table
The TruncateTable geoprocessing tool in the Tables toolset can be used
to delete all rows from a table or feature class. You should use the
Truncate Table tool instead of the Delete Rows ...
9
The most basic difference is the level at which each of the functions are working.
The Feature Class to Feature Class tool is working at the Layer level, meaning that regardless of the features selected, it treats them as a whole. A new feature class is going to be created. This means that you have control over the structure of the new featureclass, ...
9
I recommend using the Select_analysis tool which will accept a where clause.
arcpy.Select_analysis("C:/CGDP.shp", "C:/New.shp", ' "GDP_TOP10" = 1 ')
9
Features in the most general sense are analagous to records in computer science:
In computer science, records (also called tuples, structs, or
compound data) are among the simplest data structures. A record is a value that contains other values, typically in fixed number and
sequence and typically indexed by names. The elements of records are
...
9
You want to use the Append_Management tool.
If you know the attribute tables will match up (including data type), use schema_type NO_TEST anyway (even though the documentation would suggest TEST).
If the attribute tables do not match up you will have to deal with field mappings, which can be a huge pain in arcpy.
(If you are using NO_TEST and a subtype, ...
9
You can add a geometry to an existing feature which has no geometry:
Put the layer in edit mode
Select the feature from the attribute table
Use that "add part" tool from the "Advanced Digitizing Toolbar":
Draw the geometry you want to attach to the feature
8
You can use ITopologicalOperator.ConstructUnion.
The steps would be as follows:
Create an instance of the Polygon class and cast that to ITopologicalOperator. This polygon instance will contain the result.
Since ITopologicalOperator.ConstructUnion takes an IEnumGeometry to specify the geometries to union (or in your words, concatenate), create an instance ...
8
Yes, there is. Go to Options in the Settings menu and open the tab Map tools.
Then you can set the Identify Mode. You probably want the Top down, stop at first.
Only top voted, non community-wiki answers of a minimum length are eligible
