1

Can Python be used to extract flood hazard classifications from National Flood Hazard Layer (NFHL) site using any given lat/long point?

Specifically the following seven classifications:

  • 1% Annual Chance Flood Hazard
  • Regulatory Floodway
  • Special Floodway
  • Area of Undetermined Flood Hazard
  • 0.2% Annual Chance Flood Hazard
  • Future Conditions 1% Annual Chance Flood Hazard
  • Area with Reduced Risk Due to Levee

Flood Hazard Zones (28) appears to contain these seven classifications.

Additional information:

NFHL GIS Services user guide

Definitions of FEMA Flood Zone Designations

some notes:

nfhl_layers = fiona.open(county_nfhl_file_path) nfhl_layer_1 = shape(nfhl_layers['geometry'] lat = 47.661694 long= -122.192574

-- Given the lat/long, calculate the seven flood hazards for that point

This function returns true if that point is in the layer and false if it is not.

def is_lat_long_in_nfhl_layer(lat, long, nfhl_layer): pass print is_lat_long_in_nfhl_layer(lat,long, nfhl_layer_1) def get_nfhl_predicates(lat,long): ...... return {"1%_annual_chance" : percent_chance_1, "0.2%_annual_chance": percent_chance_0_point_2,}

closed as off-topic by Vince, whyzar, Dan C, John Powell, mkennedy May 25 '17 at 19:46

This question appears to be off-topic. The users who voted to close gave this specific reason:

  • "Questions seeking help to debug/write/improve code must include the desired behavior, a specific problem or error and the shortest code necessary to reproduce it in the question itself. Providing a clear problem statement and evidence of a code attempt will help others to help you. See: How to create a Minimal, Complete, and Verifiable example." – Vince, whyzar, Dan C
If this question can be reworded to fit the rules in the help center, please edit the question.

  • What have you tried? When asking for Python coding help, it is expected that you provide your best effort to solve the problem. – Vince May 25 '17 at 18:09
  • I am a Python novice... some notes: nfhl_layers = fiona.open(county_nfhl_file_path) nfhl_layer_1 = shape(nfhl_layers['geometry'] lat = 47.661694 long= -122.192574 -- GIven the lat/long, calculate the seven flood hazards for that point This function returns true if that point is in the layer and false if it is not. def is_lat_long_in_nfhl_layer(lat, long, nfhl_layer): pass print is_lat_long_in_nfhl_layer(lat,long, nfhl_layer_1) def get_nfhl_predicates(lat,long): ...... return {"1%_annual_chance" : percent_chance_1, "0.2%_annual_chance": percent_chance_0_point_2,} – sirgeo May 25 '17 at 18:21
  • 1
    Python can be used to do basically anything, seeing as it is Turing complete. Time to start learning it, as it will vastly enhance what you can do in GIS :-) – John Powell May 25 '17 at 19:21

Browse other questions tagged or ask your own question.