Tag Info

Hot answers tagged

17

Try these steps with ArcMap 10: Buffer your point feature (ArcToolbox > Analysis Tools > Proximity > Buffer). Make sure to select the correct distance in the Linear unit box. Input your newly created buffers into the Feature Envelope to Polygon tool (Data Management Tools > Features > Feature Envelope to Polygon). Make sure to select the "Create multpart ...


17

The Python GDAL/OGR Cookbook has some sample code to Buffer a Geometry. from osgeo import ogr wkt = "POINT (1198054.34 648493.09)" pt = ogr.CreateGeometryFromWkt(wkt) bufferDistance = 500 poly = pt.Buffer(bufferDistance) print "%s buffered by %d is %s" % (pt.ExportToWkt(), bufferDistance, poly.ExportToWkt()) and to Calculate intersection between two ...


15

The buffer size is always applied in the layer CRS units. Therefore, the layer CRS has to use meters if you want to buffer in meters. You don't need ftools to change the CRS. Open the original layer in WGS84 CRS. Right-click in layer list and select "Save as ...". (DON'T change the CRS in layer options!) Set the target CRS to NAD83/Maryland and save. ...


15

Use the Erase (Analysis) Tool:


10

Proper one-sided buffers were supposed to have landed in 1.5, but it looks to me that while the styles did land, sidedness didn't make it in. There is however a current patchset which exposes GEOSSingleSidedBuffer and performs the one-sided buffer as expected, under the name ST_OffsetCurve; see further background in ticket #413. In use: select ...


10

ArcGIS is utterly woeful for dissolving/merging. We had to do a buffer/merge for 3 million points recently and soon gave up on using ArcGIS -- their help desk didn't have much clue either. Postgres did it in less than an hour using the st_union function. see http://blog.cleverelephant.ca/2009/01/must-faster-unions-in-postgis-14.html


10

You can use the "Intersect" tool on the buffer and land use layers. This should create a third layer where buffers are "cut out of" the land use polygons. Then you can calculate the area of the resulting polygons using Field Calculator. Village buffer before and after the Intersect operation and calculation of area using Field Calculator: Divide by ...


10

Summary This answer places the question into a larger context, describes an efficient algorithm applicable to the shapefile representation of features (as "vectors" or "linestrings" of points), shows some examples of its application, and gives working code for using or porting into a GIS environment. Background This is an example of a morphological ...


10

Let's break this down into simple pieces. By doing so, all the work is accomplished in just a half dozen lines of easily tested code. First, you will need to compute distances. Because the data are in geographic coordinates, here is a function to compute distances on a spherical datum (using the Haversine formula): # # Spherical distance. # `x` and `y` ...


9

The steps to do this are: Select your layer by clicking on it From the Editor toolbar, select Start Editing On the Editor Menu, select Buffer Write a Negative Distance Amount to create an inside buffer... Without Inner Buffer: With Inner Buffer:


9

Here is my list of Python geoprocessing software. Shapely, python OGR, python QGIS, pyqgis, python SagaGIS, python Grass, python spatialite, pyspatialite, python PostreSQL/PostGIS, Psycopg, python R Project, rpy2, python Whitebox GAT, python GeoScript, jython


9

I made a custom Create Buffer Interval toolbox for you: http://ianbroad.com/downloads/CreateBufferIntervals.tbx In ArcMap open ArcToolbox, right click in the whitespace and click Add Toolbox. Browse to the one I made and run the Create Buffer at Interval tool. Here's a screenshot of the parameters: It should be pretty straightforward, but let me know ...


8

As Underdark comments ST_Dwithin is the recommended way of finding geometries at a distance. In many other gis-systems the buffer method is the only way but building buffers is quite costly. But there is other use cases when you need to buffer. One reason can be to visualize a buffer. If you for instance wants to show the area closer than 100 meters from a ...


8

A possible solution would be to create your "normal" round buffers using the standard ESRI buffer tool with whatever radius you would like and then performing a Feature Envelope To Polygon on that resulting feature class of buffers. This creates a square envelope feature around the extent of each feature. Feature Envelope to Polygon is located within Data ...


8

Since Erase (as @Jens linked) only is available with an Advanced license, you can download ET Geowizards. It can be installed as an Arcmap toolbox. Although you have to pay for the full suite, there's a free part of the program and the Erase function is included there (Overlay group).


7

You don't say which software you're using, but the thing you're looking for is Voronoi polygons (AKA Theissen polygons). This is the set of polygons such that any point within a polygon is nearest to its seed point. You will find that the polygons tessellate, which might be a problem if your offices have a maximum distance of responsibility. If that is the ...


7

If you are referring to the tools in ArcGIS, then one difference is that using the tool provides an output dataset, whereas Select by Location does not.


7

Defining a line's side is rather straight forward on an orientable surface, which a 2D plane in a GIS is. If you define a start point and an end point for a line, you can unambiguously define "left" and "right" sides. This is the case even if the line crosses itself. From a more practical standpoint, simple workflow for creating a one-sided buffer includes ...


7

From the menu choose: Vector --> Geoprocessing tools --> Buffer(s) To aggregate, check box 'Dissolve buffer results'.


7

For this application, I would use an Azimuthal Equidistant projection centered in the middle of your source points. This projection has the nice feature of all radial distances around the center of the projection being accurate. That particular projection is not part of QGIS standard projections. You can define your own using Settings/Custom CRS with the ...


7

It sounds like you are entering a value in miles or km (ie 10) but the projection is in geographic. The buffer tool will interpolate that as a 10 degree buffer. You will need to do one of two things. 1) translate you value to equal the buffer size in degrees, or, 2) reproject your data into a projection that is not geographic


7

Shapely gives python access to GEOS which can do buffers/intersects/etc. GEOS is the library most OSGeo programs use to perform those operations.


6

Use the Buffer Tool from the Analysis toolbox and enter a negative value for the buffer distance. Pick OUTSIDE_ONLY for the side type and then dissolve on the relevant attributes. Requires Arcinfo.


6

Hallo Here I think is a little fun way of doing it in PostGIS. This I think could be extended so the expansion follows some linestring representing the current. But now it just expands in one direction. It iterates 50 times and for each iteration takes the polygon from last iteration, moves it, bufferes it (simplifies it to make things run smother) and ...


6

If you buffer a point, the result will be a circle only as long as you don't re-project the layer. I guess you have on-thy-fly re-projection enabled. Try turning it off.


6

Forgive me if I'm missing something, but it seems that all you need here is your point shapefile and the buffer tool which is located at: Analysis Tools -> Proximity -> Buffer Specify your input points and a buffer distance. Use the "Show Help" pane in the tool if you need guidance on some of the tool options.


6

Because this question is buried in comments, let's assemble them in one place to see what's the issue: (1) I have created shapefiles that identify the location of operational incinerators and areas that have refused incinerators planning permission and wish to compare the buffers of each. (2) I also wish to ... see if there is any correlation ...


6

There are 2 ways you could approach this; server side or client side. The approach you take will depend on whether you want the processing overhead on the client or on the server (or even if you have access to a backend server like GeoServer). Method 1: Client Side There is a pure JavaScript port of the Java Topology Suite called JSTS which contains ...


6

As MappaGnosis indicated, you could write a little script for this. Here's one called polygonbuffer which takes three arguments: The output file name, the radius of your buffer, and the number of corners of the polygons. Open a Python console in QGIS, paste the script and press enter to define the function, then call using something like ...


6

Background map and your data fit together because you have "on-the-fly"-reprojection enabled. This makes sense if you want to see WGS-degree-data on projected map backgrounds. But If you want to create buffers, both layers must have the same CRS. For buffers in metres, you have to change the degree data to a projected CRS (like RGF93 Lanmbert). So you have ...



Only top voted, non community-wiki answers of a minimum length are eligible