Vjylku

ESPP--any reason not to go all in?

What would be the most expensive material to an intergalactic society?

Why restrict private health insurance?

(Codewars) Linked Lists-Sorted Insert

How do spaceships determine each other's mass in space?

Converting from "matrix" data into "coordinate" data

Is it possible to clone a polymorphic object without manually adding overridden clone method into each derived class in C++?

How to install round brake pads

-1 to the power of a irrational number

Was it really inappropriate to write a pull request for the company I interviewed with?

What is better: yes / no radio, or simple checkbox?

Logistic regression BIC: what's the right N?

Locked Away- What am I?

How do you make a gun that shoots melee weapons and/or swords?

Movie: boy escapes the real world and goes to a fantasy world with big furry trolls

Boss Telling direct supervisor I snitched

Is it appropriate to ask a former professor to order a book for me through an inter-library loan?

If sound is a longitudinal wave, why can we hear it if our ears aren't aligned with the propagation direction?

Yet another question on sums of the reciprocals of the primes

How can I portion out frozen cookie dough?

Help! My Character is too much for her story!

What does *dead* mean in *What do you mean, dead?*?

Should we avoid writing fiction about historical events without extensive research?

How do we create new idioms and use them in a novel?

Creating square buffers around point and merging them if they intersect using shapely?

Creating square buffer around point feature using ArcGIS for Desktop?point in polygon starting with csv and shp in pythonHow to create square buffers around points in QGIS with Python?Creating polygons for points with different radius in QGISmedA easy way to create buffer in QGIS?Missing 1.6% of expected area after buffering SpatialPoints from vector of widthsMerge contiguous lines unless they intersect with a pointCreate a Square Buffer around a Raster ( .vrt ) file in QGIS Python ConsoleCreate buffers around points and merge with SpatialPolygonsDataFrame to create a list of matched polygons and the area of intersectionHow to create buffers around points in ArcMap when there is a lack of OID

0

I am trying to create square buffers around given points, I am able to create circular buffers but not square ones.

from shapely.ops import transform

from shapely.geometry import Point

    local_azimuthal_projection = "+proj=aeqd +R=6371000 +units=m +lat_0={} +lon_0={}".format(lat, lon)

                    wgs84_to_aeqd = partial(

                        pyproj.transform,

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                        pyproj.Proj(local_azimuthal_projection),

                    )

                    aeqd_to_wgs84 = partial(

                        pyproj.transform,

                        pyproj.Proj(local_azimuthal_projection),

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                    )



                    point_transformed = transform(wgs84_to_aeqd, Point(float(lon), float(lat)))

                    buffer = point_transformed.buffer(0.5*1000)

                    buffered_geom = transform(aeqd_to_wgs84, buffer).exterior.coords[:]

Also, if I have some buffers that overlap, how can I merge them?

python geometry coordinates buffer

share|improve this question

edited 1 min ago

PolyGeo♦

53.6k1781243

asked 1 hour ago

AtihskaAtihska

1146

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  
  Please Edit the question to provide a picture of a "square buffer"
  
  – Vince
  1 hour ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  As per the Tour there should be only one question asked per question.
  
  – PolyGeo♦
  33 secs ago
  

  
  
  
  

add a comment | 

0

I am trying to create square buffers around given points, I am able to create circular buffers but not square ones.

from shapely.ops import transform

from shapely.geometry import Point

    local_azimuthal_projection = "+proj=aeqd +R=6371000 +units=m +lat_0={} +lon_0={}".format(lat, lon)

                    wgs84_to_aeqd = partial(

                        pyproj.transform,

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                        pyproj.Proj(local_azimuthal_projection),

                    )

                    aeqd_to_wgs84 = partial(

                        pyproj.transform,

                        pyproj.Proj(local_azimuthal_projection),

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                    )



                    point_transformed = transform(wgs84_to_aeqd, Point(float(lon), float(lat)))

                    buffer = point_transformed.buffer(0.5*1000)

                    buffered_geom = transform(aeqd_to_wgs84, buffer).exterior.coords[:]

Also, if I have some buffers that overlap, how can I merge them?

python geometry coordinates buffer

share|improve this question

edited 1 min ago

PolyGeo♦

53.6k1781243

asked 1 hour ago

AtihskaAtihska

1146

• 
  
  
  

  
  1
  

  
  
  
  
  
  

  
  
  Please Edit the question to provide a picture of a "square buffer"
  
  – Vince
  1 hour ago
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  As per the Tour there should be only one question asked per question.
  
  – PolyGeo♦
  33 secs ago
  

  
  
  
  

add a comment | 

I am trying to create square buffers around given points, I am able to create circular buffers but not square ones.

from shapely.ops import transform

from shapely.geometry import Point

    local_azimuthal_projection = "+proj=aeqd +R=6371000 +units=m +lat_0={} +lon_0={}".format(lat, lon)

                    wgs84_to_aeqd = partial(

                        pyproj.transform,

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                        pyproj.Proj(local_azimuthal_projection),

                    )

                    aeqd_to_wgs84 = partial(

                        pyproj.transform,

                        pyproj.Proj(local_azimuthal_projection),

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                    )



                    point_transformed = transform(wgs84_to_aeqd, Point(float(lon), float(lat)))

                    buffer = point_transformed.buffer(0.5*1000)

                    buffered_geom = transform(aeqd_to_wgs84, buffer).exterior.coords[:]

Also, if I have some buffers that overlap, how can I merge them?

python geometry coordinates buffer

share|improve this question

edited 1 min ago

PolyGeo♦

53.6k1781243

asked 1 hour ago

AtihskaAtihska

1146

from shapely.ops import transform

from shapely.geometry import Point

    local_azimuthal_projection = "+proj=aeqd +R=6371000 +units=m +lat_0={} +lon_0={}".format(lat, lon)

                    wgs84_to_aeqd = partial(

                        pyproj.transform,

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                        pyproj.Proj(local_azimuthal_projection),

                    )

                    aeqd_to_wgs84 = partial(

                        pyproj.transform,

                        pyproj.Proj(local_azimuthal_projection),

                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),

                    )



                    point_transformed = transform(wgs84_to_aeqd, Point(float(lon), float(lat)))

                    buffer = point_transformed.buffer(0.5*1000)

                    buffered_geom = transform(aeqd_to_wgs84, buffer).exterior.coords[:]

share|improve this question

edited 1 min ago

PolyGeo♦

53.6k1781243

asked 1 hour ago

AtihskaAtihska

1146

share|improve this question

edited 1 min ago

PolyGeo♦

53.6k1781243

asked 1 hour ago

AtihskaAtihska

1146

edited 1 min ago

PolyGeo♦

53.6k1781243

edited 1 min ago

PolyGeo♦

53.6k1781243

asked 1 hour ago

AtihskaAtihska

1146

asked 1 hour ago

AtihskaAtihska

1146