I am trying to automate a process to populate a field in a point feature class with the name of district (a polygon layer) that it falls inside using python. Ideally this script would only populate the points that have NULL district names. I'm thinking this may be done programmatically using the select by attribute, select by location, and calculate field functions? If someone could point me in the right direction, that would be much appreciated!