User:Klaus/KingOfTheGraticules

From Geohashing
< User:Klaus
Revision as of 22:19, 28 August 2023 by Klaus (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

This is just a stupid idea which I had. I'll dump it here for now, until I improve it.

What is the idea?

The idea is to create a (daily) world map where the whole world cube is claimed by geohashers. Each point of the world is claimed by one geohasher at a time and he/she is the "king"/"queen" of that area. If you want to improve your kingdom, you can either visit more hashes in your existing kingdom which makes it harder for others to claim your area or you can visit more hashes in foreign weak areas and expand your kingdom.

Someone in Siberia or on the Atlantic ocean might claim a huge amount of land by just doing one geohash in the middle of it. People living in high density geohashing areas might not claim any land at all if they live close to someone with 100+ hashes.

I might only count "recent" hashes, too.

Things I need to decide

  • Should geohashes fade out linearly over time, i.e. after 10 years, a geohash does not count anymore?
  • What kind of weight function to use? Something exponential like "exp(-dist/C)" or something quadratic like "1/(1+dist)^2"? I'll probably decide as soon as I have some maps.
  • Count only reached hashes? I think so...
  • How fine? 100 pixels per Graticule or more?

Script for now

# !/usr/bin/python

import json
from math import radians, cos, sin, asin, sqrt, exp
from copy import deepcopy

# ugly parsing of https://fippe.de/alldata.js
data=open("alldata.js","r").read()
HEADER = "var expeditions = "
FOOTER = ",\n];\n"
assert(data.startswith(HEADER))
assert(data.endswith(FOOTER))
data=data[len(HEADER):]
data=data[:-len(FOOTER)]
data=data+"]"
data = json.loads(data)

# now we have a list of geohashes like
# ["2023-08-16_52_7",52.05395,7.73412,["Fippe"],true,1,["Bicycle_geohash_achievement", "Land_geohash_achievement"]],
# ["2010-01-12_global",33.67653,-81.27146,["NWoodruff","LuxMundi"],true,12,["Globalhash_achievement"]],

# filter only reached hashes
data = [i for i in data if i[4] == True]

# get set of users
users = set()
for i in data:
    users |= set(i[3])
print(users)

# https://stackoverflow.com/questions/4913349/haversine-formula-in-python-bearing-and-distance-between-two-gps-points/4913653#4913653
def haversine(lon1, lat1, lon2, lat2):
    """
    Calculate the great circle distance between two points 
    on the earth (specified in decimal degrees)
    """
    # convert decimal degrees to radians 
    lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])

    # haversine formula 
    dlon = lon2 - lon1 
    dlat = lat2 - lat1 
    a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
    c = 2 * asin(sqrt(a)) 
    r = 6371 # Radius of earth in kilometers. Use 3956 for miles
    return c * r


LONRANGE = [i for i in range(-90,91)]
LATRANGE = [i for i in range(-180,181)]
LONRANGE = [i for i in range(47,55+1)]
LATRANGE = [i for i in range(5,15+1)]


king = [ [ "" for lat in LATRANGE ] for lon in LONRANGE] 
print("init density done")


def weight(dist):
    assert(dist >= 0)
    return 1. / (dist+1)**2

for lonindex, lon in enumerate(LONRANGE):
    for latindex, lat in enumerate(LATRANGE):
        tmp = deepcopy(dict.fromkeys(users, 0.))
        for geohash in data:
            geolon, geolat = geohash[1], geohash[2]
            for geouser in geohash[3]:
                tmp[geouser] += weight(haversine(lon,lat,geolon,geolat))
        max_user  = max(tmp, key=tmp.get)
        #print(tmp)
        print(f"maximum for {lon},{lat} is {max_user} with weight {tmp[max_user]}")
        king[lonindex][latindex] = (max_user, tmp[max_user])

example output for Germany

maximum for 47,5 is Baarde with weight 0.0023012356385766286
maximum for 47,6 is Crox with weight 0.007699984689657429
maximum for 47,7 is Crox with weight 0.038897544181984746
maximum for 47,8 is Tungmar with weight 0.018528739383917605
maximum for 47,9 is MykaDragonBlue with weight 0.040467624501727084
maximum for 47,10 is Ekorren with weight 0.00479617522599326
maximum for 47,11 is LeslieHapablap with weight 0.010826449652301917
maximum for 47,12 is LeslieHapablap with weight 0.004651845155012959
maximum for 47,13 is Reinhard with weight 0.00181486196596779
maximum for 47,14 is B2c with weight 0.0013414195123689567
maximum for 47,15 is B2c with weight 0.0017957337500081071
maximum for 48,5 is GeorgDerReisende with weight 0.0018293951759897626
maximum for 48,6 is Baarde with weight 0.0036729989744967197
maximum for 48,7 is Baarde with weight 0.03905682877442473
maximum for 48,8 is RecentlyChanged with weight 0.0634133129194444
maximum for 48,9 is Ekorren with weight 0.0573629377999951
maximum for 48,10 is Ekorren with weight 0.013606021617061224
maximum for 48,11 is TheOneRing with weight 0.10490586164127187
maximum for 48,12 is Harleydavidsonandy with weight 0.049155557595685906
maximum for 48,13 is Robert with weight 0.014065845915421692
maximum for 48,14 is B2c with weight 0.008400864577314968
maximum for 48,15 is B2c with weight 0.0048942445828206275
maximum for 49,5 is GeorgDerReisende with weight 0.002298096872817964
maximum for 49,6 is Baarde with weight 0.004815191763979033
maximum for 49,7 is Baarde with weight 0.010452285862118691
maximum for 49,8 is Q-Owl with weight 0.1308363622024999
maximum for 49,9 is DODO with weight 0.11608691966236
maximum for 49,10 is DODO with weight 0.16226044004570395
maximum for 49,11 is Danatar with weight 0.015067194959247333
maximum for 49,12 is dawidi with weight 0.9275373146653234
maximum for 49,13 is dawidi with weight 0.054531989434608054
maximum for 49,14 is Danatar with weight 0.0026655669753878207
maximum for 49,15 is B2c with weight 0.003035046716508224
maximum for 50,5 is GeorgDerReisende with weight 0.0028013838715444327
maximum for 50,6 is GeorgDerReisende with weight 0.004819725875603006
maximum for 50,7 is Gefrierbrand with weight 0.011305720533205395
maximum for 50,8 is Mampfred with weight 0.2493283500677317
maximum for 50,9 is GeorgDerReisende with weight 0.10955262573524727
maximum for 50,10 is Danatar with weight 0.5574659857259744
maximum for 50,11 is Klaus with weight 0.03246377476871615
maximum for 50,12 is Danatar with weight 0.04510605156118596
maximum for 50,13 is Danatar with weight 0.006932465340946476
maximum for 50,14 is Thepiguy with weight 0.016236373964000697
maximum for 50,15 is redguy with weight 0.0036676901158184134
maximum for 51,5 is NeThuS with weight 0.02278942467026077
maximum for 51,6 is GeorgDerReisende with weight 0.007758063094770804
maximum for 51,7 is GeorgDerReisende with weight 0.046127187525615246
maximum for 51,8 is GeorgDerReisende with weight 0.031176024530840462
maximum for 51,9 is GeorgDerReisende with weight 0.11182640443754852
maximum for 51,10 is GeorgDerReisende with weight 0.20197117796069514
maximum for 51,11 is Rincewind with weight 0.3068679952657682
maximum for 51,12 is Paintedhell with weight 0.18043674734153983
maximum for 51,13 is Danatar with weight 0.042979847300641985
maximum for 51,14 is Danatar with weight 0.10304354931677923
maximum for 51,15 is Danatar with weight 0.05158105699417946
maximum for 52,5 is JwB with weight 0.08384123652599894
maximum for 52,6 is FelixTheCat with weight 0.03751767255243819
maximum for 52,7 is Arvid with weight 0.10008613866271686
maximum for 52,8 is Fippe with weight 0.04333370440046326
maximum for 52,9 is Fippe with weight 0.11099314064914968
maximum for 52,10 is Fippe with weight 0.23921887962904775
maximum for 52,11 is Fippe with weight 0.028824683331931167
maximum for 52,12 is Reinhard with weight 0.012544702986250447
maximum for 52,13 is relet with weight 0.18398722671708406
maximum for 52,14 is relet with weight 0.012282745304953786
maximum for 52,15 is Danatar with weight 0.003753119935988901
maximum for 53,5 is Eupeodes with weight 0.014214693402230935
maximum for 53,6 is Shevek with weight 0.006978796178056222
maximum for 53,7 is Blauwe_BIC with weight 0.15384901950770558
maximum for 53,8 is Fippe with weight 0.05434795167890141
maximum for 53,9 is Fippe with weight 0.09637262165602019
maximum for 53,10 is Fippe with weight 0.06578415936139503
maximum for 53,11 is Fippe with weight 0.029673455220319597
maximum for 53,12 is relet with weight 0.02162159157388358
maximum for 53,13 is relet with weight 0.0200411966585085
maximum for 53,14 is relet with weight 0.01004751712891905
maximum for 53,15 is relet with weight 0.002251959454330422
maximum for 54,5 is Sourcerer with weight 0.002593321050567368
maximum for 54,6 is GeorgDerReisende with weight 0.0030226125726690372
maximum for 54,7 is Fippe with weight 0.005984561424176756
maximum for 54,8 is Fippe with weight 0.01124357033856347
maximum for 54,9 is π_π_π with weight 0.02349812462248916
maximum for 54,10 is π_π_π with weight 0.407150149215557
maximum for 54,11 is π_π_π with weight 0.012193396496538186
maximum for 54,12 is Trurl with weight 0.010269151954365657
maximum for 54,13 is SastRe.O with weight 0.021989678217166495
maximum for 54,14 is π_π_π with weight 0.07916964180279705
maximum for 54,15 is JRK with weight 0.0026736658022132624
maximum for 55,5 is Sourcerer with weight 0.0020892046174319605
maximum for 55,6 is GeorgDerReisende with weight 0.0022102755980369965
maximum for 55,7 is Fippe with weight 0.002978486381969255
maximum for 55,8 is Fippe with weight 0.004453126739758758
maximum for 55,9 is π_π_π with weight 0.005982335660953738
maximum for 55,10 is π_π_π with weight 0.009470347531698415
maximum for 55,11 is π_π_π with weight 0.005449596705163217
maximum for 55,12 is Pastori with weight 0.004622361573023937
maximum for 55,13 is Llavids with weight 0.00900023902201388
maximum for 55,14 is Llavids with weight 0.008172263901904994
maximum for 55,15 is Llavids with weight 0.0029021091702394843