Update app.py

app.py

@@ -6,6 +6,9 @@ import pandas as pd
 from PIL import Image
 import numpy as np
 from streamlit_folium import st_folium
+from typing import Dict, List
+import requests
+from time import sleep
 
 def process_image(image_path):
     # Read the image
@@ -59,17 +62,49 @@ def create_map(location_name, latitude, longitude):
     return map
 
 def analyze_location(location_name, latitude, longitude):
-    # Dummy function to gather information about the location
-    # In a real application, this would query a geospatial database or API
     location_info = {
         'Location Name': location_name,
         'Latitude': latitude,
         'Longitude': longitude,
-        'Environmental Impact': 'Moderate',
-        'Geographical Features': 'Mountains, Rivers',
-        'Agricultural Impact': 'High',
-        'Urbanisation Impact': 'Low'
     }
+
+    # Get address data from Nominatim API
+    nominatim_url = f"https://nominatim.openstreetmap.org/reverse?format=json&lat={latitude}&lon={longitude}&zoom=10"
+    headers = {'User-Agent': 'DroneImageProcessor/1.0'}
+    response = requests.get(nominatim_url, headers=headers)
+    sleep(1)  # Respect rate limiting
+
+    if response.status_code == 200:
+        osm_data = response.json()
+        address = osm_data.get('address', {})
+
+        # Extract relevant information
+        location_info.update({
+            'Country': address.get('country', 'Unknown'),
+            'State/Region': address.get('state', address.get('region', 'Unknown')),
+            'City': address.get('city', address.get('town', address.get('village', 'Unknown'))),
+            'Land Use': address.get('landuse', 'Unknown'),
+        })
+
+    # Get elevation data from OpenTopoData API
+    elevation_url = f"https://api.opentopodata.org/v1/aster30m?locations={latitude},{longitude}"
+    response = requests.get(elevation_url)
+    sleep(1)  # Respect rate limiting
+
+    if response.status_code == 200:
+        elevation_data = response.json()
+        elevation = elevation_data['results'][0]['elevation']
+        location_info['Elevation'] = f"{elevation:.1f} meters"
+
+        # Basic terrain classification based on elevation
+        if elevation < 100:
+            terrain = "Lowland"
+        elif elevation < 500:
+            terrain = "Hills"
+        else:
+            terrain = "Mountains"
+        location_info['Terrain Type'] = terrain
+
     return location_info
 
 # Streamlit app
@@ -89,11 +124,11 @@ longitude = longitude_deg if longitude_dir == "E" else -longitude_deg
 # Enter button to find and display the location on the map
 if st.button("Enter"):
     map = create_map(location_name, latitude, longitude)
-    st_folium(map, width=700, height=500)
     location_info = analyze_location(location_name, latitude, longitude)
     st.write("### Location Information")
     for key, value in location_info.items():
         st.write(f"{key}: {value}")
+    st_folium(map, width=700, height=500)
 
 # File upload
 uploaded_files = st.file_uploader("Upload Images", type=["jpg", "jpeg", "png"], accept_multiple_files=True)
@@ -117,9 +152,8 @@ if uploaded_files and location_name and latitude and longitude:
 
     # Display images
     display_images(images)
-    # Display results
-    from typing import Dict, List
 
+    # Display results
     combined_results: Dict[str, List[Dict[str, str]]] = {
         'Geological Analysis': [],
         'Agricultural and Environmental Analysis': [],