Tag Info

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 ...

13

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 ...

12

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. ...

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

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 ...

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

8

i dont know is this work for you, but you can do: . 1.Select your layer with clicking on it 2.From Editor toolbar, select Start Editing 3.in Editor Menu, select Buffer 4.Write Negative Distance Amount for inner buffer... . Without Inner Buffer: With Inner Buffer: manual inner buffer solution which i use. i hope it helps you...

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).

8

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 ...

7

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 ...

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

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 ...

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

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

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

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

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

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 ...

6

To simplify, Shapely: manual allows all geometry processing of PostGIS in Python. The first premise of Shapely is that Python programmers should be able to perform PostGIS type geometry operations outside of an RDBMS... The first example of PolyGeo from shapely.geometry import Point, LineString, Polygon, mapping from shapely.wkt import loads pt = ...

5

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 ...

5

You can perform operations with WKT you return from consultations in postgis SELECT ST_ASTEXT(ST_Buffer(ST_GeomFromText('POINT(100 100)'), 500)); you recieve WKT and use Openlayer to convert the wkt to geometry and feature. for example : var elemento = New OpenLayers.Format.WKT().read('textWktFromRecordColumn'); I hope it helps. greetings

5

ArcGIS is the software suite name, whereas ArcView, ArcEditor, & ArcInfo are licensing levels for ArcGIS. ArcView has limited buffer options, and ArcInfo has the full, advanced buffer options (see bullet point #4 under "Usage Tips"). http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=buffer_%28analysis%29

5

I am assuming you are using ArcGIS 10. Let me know if you are using a different version and I will change this for that version. First, make sure that your result shapefile from the intersection is in a projected coordinate system. It should be with the actions you have done already, but it is a good idea to check. Run Dissolve (Data Management > ...

5

OK, how I did this, is loop through creating ring buffers on objects, and saving the output as a new featureclass in a standlone datasets collection. I then merged the featureclasses, like so: tmp_merge = "path to area"\+ "tmp_merge" final_out_path_agg = final_staging + "_Agg" final_out_path_simp = final_staging + "_Simp" #check if files exist and ...

5

Select by Location with negative 1000m buffer, switch selection, select by location "Within" again from reversed selection set Step 1: Selection method: Select Features From (default) Target Layer: Houses Source Layer: Big Polygon Spatial Selection Method: Target layers(s) features intersect the Source layer feature (default) Apply a search distance = ...

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