AnsweredAssumed Answered

Buffer geocoded Address

Question asked by ahmadabusaleh on Oct 10, 2018
Latest reply on Oct 18, 2018 by ahmadabusaleh
  • Comment • 5

Hello,

I am using Python API to geocode a list of address, then buffer the geocoded address by 1 mile, what I have so far:

 

addlist=['adress1',
         'address2',
         'address3']

bgeocode=batch_geocode(addlist)


for address in bgeocode:
   

 poly=buffer(geometries=[Geometry(address['location'])],in_sr=3857,
            distances=[1],unit='Miles')[0]
    
 m1.draw(poly)



The result:

 

 

Now my question is why the result buffer looks like an ellipse not a circle? also, the circles are more than 1 mile ? is there another method to buffer the address and show them on the map, because here I am using the buffer for each address?

 

Thank you,

Ahmad 

simo
simo Oct 15, 2018 9:55 PM
Correct Answer

The buffer function is a bit confusing and clunky, here is what I mean:

 

When you specify the in_sr value, the unit parameter is not relevant any more. for example, if you specify in_sr=4326, the unit will be "degree" regardless whatever value you specify for the unit parameter. if you set in_sr=3857, the input values will be "Meters".

 

In regard to which projection to choose, for a 1 Mile radius circle, you should not see big difference between SRID 3857 and a relative accurate local projection --- the distortion caused by projection become more visible over long distances or large areas.

 

I would recommend to use use_proximity.create_buffers() for your task.

 

 

addresses = ['Adelaide, South Australia',
            'Mount Barker, South Australia',
            'Brookside St, Oakden,South Australia']
results = batch_geocode(addresses)

fs=[Feature(geometry=Geometry(res['location'])) for res in results]
fset = FeatureSet(fs)

sms = {
    "color": [255,0,0,150],
    "size": 10,
    "angle": 0,
    "xoffset": 0,
    "yoffset": 0,
    "type": "simple-marker",
    "outline": {
        "color": [
            255,
            255,
            255,
            255
        ],
        "width": 1,
        "type": "esriSLS",
        "style": "esriSLSSolid"
    }
}

for res in results:
    map1.draw(res['location'],symbol=sms)

prox=use_proximity.create_buffers(fset.to_dict(),distances=[1],units='Miles')
map1.add_layer(prox)

See the reply in context
No one else had this question

Outcomes