Tag Info

New answers tagged

0

JTS is the gold-standard topology library for Java. The Geometry class has a full set of methods for relating geometries, including intersection.


0

ESRI has a GIS SDK on Android ..It has samples on spatial relationships.


1

If you don't want to zoom before around Heidelberg and if you are pretty sure that only one object is called Umweltzone Heidelberg, you can set a quick query like this : key name value Umweltzone Heidelberg in : empty value. This will create only an attribute query, not a geographic query. This kind of query is very quick. <osm-script output="xml" ...


3

I guess your query area Germany does not exist with that name. Using Deutschland and disabeling the timeout returns the data. I zoomed around Heidelberg and ran the query successfully and much quicker, querying the extent of the map canvas.


2

Use a colour scale, getting darker for each additional person buried in a plot. ColorBrewer is a great place to look for good colour schemes. I imagine most would just be one person, but ones that are darker indicate the presence of multiple people. If your problem is how you associate cemetary plots and a table of people buried therein, then you need to ...


2

As you have specified that the line from the point can pass in/out of the polygon then I think all you need to do is extract the Envelope of the polygon. You also state in the comment above you are look for an arcpy solution and you have a Basic license level. You can get most of the way with model builder and a sneaky use of the fishnet tool to extract the ...


0

Convert the given polygon into a set of line segments.Find the nearest segment from the given point feature. You can also use ICurve interface to query the farthest edge.


0

After seeing the raster image of the field, I though of this process - I used Arcgis 10.2.1. (Note: Just realised that I have got the length and widths around the wrong way when compared to the image, but If I was to recreate the answer I would use the same process below but amend the lengths and widths) Using the dimensions provided I worked out what the ...


3

You can try the Sort Tool. It can be used to sort by attribute, or more importantly for you, spatially. It allows you to pick which corner (i.e. upper right) to give priority when sorting.


2

I have never worked in QGIS. The R equivalent would be to use gUnion in the rgeos package: http://www.inside-r.org/packages/cran/rgeos/docs/gUnion Here is a reproducible example to demonstrate: library(sp) library(rgeos) #create polygon r1 <- cbind(c(641777, 642290, 642276, 641794), c(7036885, 7036743, 7036154, 7036146)) r2 <- cbind(c(642320, ...


2

Thanks for your feedback. I solved this by creating a second polygon around the entire raster extent, to which I assigned a value of zero. Previously, the shapefile only contained a single polygon, hence the cells outside of it had no value. Thus, the new shapefile and raster contain a value of 7.5 at the polygon and 0 elsewhere. (PS I am relatively new to ...


0

^Hi `Andrea, The Union will create duplicate polygons on top of each other (older one + newer one) at areas where they overlap (the selected polygons below have 2 overlapping). You could deduplicate these by running Find Identical (use the [SHAPE] field and check "Output only duplicate records". Then Join your union polygons based on [OBJECTID] and ...


0

I would try something like this assuming we're interested only in those geometries that lie inside of some other geometry - these are considered as higher valued. UPDATE polygons -- adjust geometry SET wkb_geometry = new_polygon FROM ( SELECT p1.id, -- take the difference of the geometries ST_Difference(p1.wkb_geometry, ...


2

As requested, comment is appended below as an answer: "Also, is Feature To Polygon designed to work on a layer selection? Perhaps export the layer as a standalone feature class with Copy Features (in memory perhaps) and try running the tool that way. If that doesn't solve it, there are some threads here and here that use arcpy.Geometry.cut() for those who ...


3

According to this QGIS Change Log entry for the latest 2.10 "Pisa" release of QGIS, support for true curve features, including circles, was only just added in this latest release. See the "New geometry engine" topic. It requires a datasource like PostGIS that supports true curved geometry types. However, the documentation there is slightly ambiguous. It ...


2

Take a look at the plug in CADDigitize. It may be helpful. edit : This tool allows the user to draw circles, squares, resctangels, ellipses from many ways.


1

To be honest, a python script (if you are familiar with Python) would be the most efficient method of creating your polygons in your use-case. However, if you are not familiar with scripting, then there are a variety of tools that you should investigate in QGIS (e.g. Vector->Research Tools->Vector Grid). You will have to use a crazy combination of ...


1

The following python script (which requires ArcGIS 10.1 or later) utilizes arcpy.da to take a shapefile as an input and create a spreadsheet with an entry for each vertex in each polygon present in the .shp (and I believe it works with lower level arcgis licenses). Object and sequence id's tie the points back to a specific position in a specific polygon. ...


2

Feature To Polygon will slice polygons by line feature classes when you use the two feature classes as its input. Inputs: Results: If you only want certain polygons sliced by certain lines, use an SQL expression and feature layers to limit your inputs.


1

You could use Zonal Statistics as Table to generate a table with a Unique ID field that can be used to join to the polygon buffer attribute table. The tool takes 1) a dataset that defines the zone (your buffer polygons), 2) a zone field (some unique ID field in the buffer polygons), 3) a raster that contains the values on which to calculate a statistic ...


2

A point inside a polygon canvas = qgis.utils.iface.mapCanvas() aLayer = canvas.layer(0) # or code to select the layer of interest for f in aLayer.getFeatures(): geom = f.geometry() p = geom.pointOnSurface() print p.asPoint() Get Inner Rings canvas = qgis.utils.iface.mapCanvas() aLayer = canvas.layer(0) # or code to select the layer of ...


0

It is possible with the LAStools Production toolbox, which can also be used from the ArcGIS toolboxes. According to rapidlasso GmbH the new toolbox has the advantage of batch automating LiDAR tasks across folders of las/laz files. A summary of Martin Isenburg's answer in the LAStools GoogleGroups is: a) Index the las files: lasindex -i folder/*.las ...


0

Others users might be interested to know that I managed to solve this. In order to extend the lines outside of the polygon (problem a): # Save centre point co-ordinates separately A.x <- Centre@coords[1] A.y <- Centre@coords[2] Using the great answer posted here: http://stackoverflow.com/a/7741655/4312169 I applied the code: # Generate points ...


0

relation //Name of country with native language // Use ["name:en"="France"] for english ["name"="France"] [boundary=administrative] [admin_level=2]; (._;>;); out body qt; It will return you Polygons, Lines and Points, so filter features by type, when import to GIS.


0

select desired polygons you are working with and from selection make layer. export layer to CAD & explode polygons and just keep lines that make up the edge of the road, add any lines needed to close lines into polygon and join all lines in cad to make it a polygon. add CAD drawing to map. select CAD polygon copy & paste into layer as polygon to ...


0

There is a tool in the linear referencing tools called Locate Features along routes. It sounds like you might have already tried using it, but if you increase the search radius to a distance greater than the furthest distance from your route you should be able to get not only measures, but also offset distance. I think the distance option only works with ...


0

Ok, here is the method I came up with and works for me. Union feature class by itself Use Find Identical on the unioned file Join the unioned file to the table produced in step 2. by the ObjectID and IN_FID fields. Dissolve by "FEAT_SEQ" field, putting the sum of the score in the statistics field.


2

Intersect will give you polygons representing the overlap. Then you can use Dissolve with a statistic of sum on your value field. Finally, you can Merge/Append back to the original dataset:


3

First you have to convert your polygon to lines. Then you can use Generate Near table. With this tool you can create table with few parameters (coordinates, angles) of nearest objects to your input objects. You can get this coordinates and create point layer if you need to visualise that.


1

Shapefiles should be converted to file geodatabases for performance and precision and for many options not available with shapefiles. One of the benefits of converting a shapefile to a feature class in a file geodatabase is that you get an automatically calculated area and perimeter length field for your polygons (length for polylines). If you keep it as a ...


1

Sounds like cleaning the geometry up will get you on the road to your non-raster solution ... herewith something of a kludge which does help with fixing bad geometry: # Load the library and problematic data library(rgeos) load("oneImage2_spdf.Rdata") >gIsValid(polysData, reason = T) Error in RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td, ...


1

Since you put GRASS into the tags, here's a solution based on GRASS: First, you need to know exactly what coordinate system the original data is in (as always with GRASS). I see that the *.prj file contains "TRANSVERSE MERCATOR" but it's not one of the standard UTM zones. Since you have not mentioned the EPSG code of this data, here's the proj4 string for ...


4

Without your original data, I can't be sure this will work, but I thought it might help you out. I didn't bring it all the way there, this solution still likely needs some level of automation, but might give you a general way forward First, I create some spatial polygons polypoints1 <- matrix(c(1,2,2,1,1,2,2,1,1,2),ncol=2) polypoints2 <- ...


0

easily done. the return geoJson will be multiple polygons in a GeoJson Response..this can be easily rendered in GoogleMap,etc. Very simple solution to this. sign up here: https://www.mashape.com/vanitysoft/boundaries-io then in your code: https://vanitysoft-boundaries-io-v1.p.mashape.com/reaperfire/rest/v1/public/boundary?zipcode=20037,20002 The results ...


1

I had the same problem and i tried to solve it with the answers above but it did not work for me. What actually worked is this var polygonLayer = new OpenLayers.Layer.Vector("Polygon"); var polygon = new OpenLayers.Control.DrawFeature(polygonLayer, OpenLayers.Handler.Polygon); map.addControl(polygon); polygon.handler.callbacks.point = function(point){ ...


2

I can think of two things off the top of my head. If you want to stick with an Add-In, you can use the onLinemethod and simply iterate through the line geometry the user creates and pass those points into an array and convert that to a polygon. The can add segments and finishes by double clicking. Or, another option is to make a script tool. With a ...


2

Yes, it's possible to fill a polygon with symbols on an angle. It cannot be done with Marker symbol type, only with Picture symbol type. Even with Picture it applies to the symbol as a whole, meaning you'd need to symbolize on your rotation field and set each one to the angle you want - the angle cannot be read automatically from a field (so far as I ...


0

You might want to try using the SpatialRelater and the 'intersects' test. This will find (and associate) all the features in one layer that touch another layer.


3

It seems that you're misunderstanding the way the AreaOnAreaOverlayer behaves. It actually intersects the polygons and creates new ones if they intersect. The documentation states that the tool: Performs an area-on-area overlay so that all input areas are intersected against each other and resultant area features are created and output. The ...


0

Without more information about your data sets and workspace its difficult to know why AreaOnAreaOverlayer isn't working as expected. Perhaps try side-stepping the need for the AreaOnAreaOverlayer in your workspace by putting one of your polygons through the CenterOfGravityReplacer and then using the resulting points along with your other polygon as the two ...


0

If you want to locate Shapes from A within B, I think you could use a "tester" transformer to filter according "_overlaps" attribute that is created on each of the output elements. Something like this: I hope it serves this idea.


0

Assuming you want to use it in Model Builder, I suggest you a "Calculate Value" script as follows: def LargestEcoregion(tabulate_intersection_result): rows=arcpy.SearchCursor(tabulate_intersection_result,sort_fields="PERCENTAGE D") row=rows.next() arcpy.AddMessage("%s--%s" ...



Top 50 recent answers are included