projekt-fredrika/islands.py at main · Cqsi/projekt-fredrika · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
import pandas as pd
import requests
import re

from math import cos, asin, sqrt, pi

def distance(lat1, lon1, lat2, lon2):
    p = pi/180
    a = 0.5 - cos((lat2-lat1)*p)/2 + cos(lat1*p) * cos(lat2*p) * (1-cos((lon2-lon1)*p))/2
    return 12742 * asin(sqrt(a)) #2*R*asin...

excel_file = "excel\\nagu_testing.xlsx"

api_key = "e99a7c7b-40c7-43e6-a062-b25cdbe7cd52"

# stop = 20
# counter = 0

dist = []
hits = []

try:
    df = pd.read_excel(excel_file, engine='openpyxl')
    print("Excel read.")
except Exception as e:
    print("Something went wrong. Check the error:")
    print(e)

def filter_islands(island):
    if (island["properties"]["label:placeType"] == "Saari tai luoto" or island["properties"]["label:placeType"] == "Saari- tai luotoryhmä") and island["properties"]["label:municipality"] == "Parainen":
        return True
    else:
        return False


for index in df.index:

    # counter+=1

    # if counter == stop:
    #     break

    cur = df.loc[index, "itemLabel"]

    response = requests.get("https://avoin-paikkatieto.maanmittauslaitos.fi/geocoding/v1/pelias/search?sources=geographic-names&text=" + cur + "&api-key=" + api_key)

    if response.status_code == 200:

        islands = response.json()["features"]

        islands_parainen = list(filter(filter_islands, islands))
        length = len(islands_parainen)

        coords = ""

        if length==0:
            print("Island not found.")
            print()
        elif length>1:
            print("More than one result.")
            print()
        else:
            prop = islands_parainen[0]["properties"]

            label = prop["label"]
            placeId = prop["placeId"]
            placeElevation = prop["placeElevation"]
            municipality = prop["label:municipality"]
            coordinates = islands_parainen[0]["geometry"]["coordinates"]

            # Printing out values in terminal
            print(label)
            print(placeId)
            print(placeElevation)
            print(municipality)
            print(coordinates)

            wd_coords = re.findall(r'\d+\.\d+', df.loc[index, "coords"])

            coords = round(distance(float(wd_coords[0]), float(wd_coords[1]), coordinates[0], coordinates[1])*1000, 2)

            print(str(coords) + " meters")
            print()


            df.loc[index, "MML_label"] = label
            df.loc[index, "MML_placeId"] = placeId
            df.loc[index, "MML_coordinates"] = str(coordinates[0]) + ", " + str(coordinates[1])
            df.loc[index, "placeElevation"] = placeElevation
            df.loc[index, "MML_label:municipality"] = municipality

        hits.append(length)
        dist.append(coords)


    else:
        print("Something went wrong. Status code " + response.status_code)


df["distance"] = dist
df["träffar"] = hits

# Export dataframe as excel
df.to_excel("excel\\edited_nagu.xlsx", index = False)