In [1]:

    

import json
import pandas as pd


    

In [2]:

    

with open('geoJson/all_atlases.geojson') as fp:
    geojson = json.load(fp)


    

In [15]:

    

urns = {
    "blaeu":"urn-3:FHCL:5347094",
    "colom":"urn-3:FHCL:6750170",
    "dewit":"",
    "dudleyV1":"urn-3:FHCL:9695784",
    "dudleyV3":"urn-3:FHCL:9695784",
    "goos":"urn-3:FHCL:5362897",
    "keulenV1":"urn-3:FHCL:6712457",
    "keulenV2":"urn-3:FHCL:6712457",
    "renard":"urn-3:FHCL:865460",
    "waghenaer":"urn-3:FHCL:8645097"
}


    

In [16]:

    

for f in geojson['features']:
    f['properties']['atlasURN'] = urns[f['properties']['collection']]


    

In [17]:

    

geojson['features'][50]


    

    
Out[17]:





{'geometry': {'coordinates': [[[0.2125836346515833, 60.029404824274394],
    [0.2125836346515833, 64.40985482427439],
    [12.188133634651585, 64.40985482427439],
    [12.188133634651585, 60.029404824274394]]],
  'type': 'Polygon'},
 'id': None,
 'properties': {'DRS_ID': 36673588,
  'HOLLIS': '001644159',
  'MAP': None,
  'OwnerSuppliedName': '001644159_0287',
  'SEQUENCE': 287,
  'UNIQUE_ID': '36673588',
  'URN': 'urn-3:FHCL:5347381',
  'atlasURN': 'urn-3:FHCL:5347094',
  'collection': 'blaeu',
  'degree_area': 52.45829799749996,
  'fname': '36673588_geotiff.tif',
  'geographic_scope': 'Western Coast of Norway',
  'lat_extent': 4.380449999999996,
  'lng_extent': 11.975550000000002,
  'maxLat': 64.40985482427439,
  'maxLong': 12.188133634651585,
  'maxZoom': 9,
  'minLat': 60.029404824274394,
  'minLong': 0.2125836346515833,
  'minZoom': 4,
  'pds_link': '\n<p><a href="http://pds.lib.harvard.edu/pds/view/36673588?n=287">View map in Page Delivery Service</a></p>\n',
  'pretty_collection': 'Blaeu',
  'title': 'Blaeu | Western Coast of Norway'},
 'type': 'Feature'}

In [19]:

    

with open("geoJson/all_atlases.geojson","w") as fp:
    json.dump(geojson,fp,indent=2,sort_keys=True)