def main(): author = 'Milena Veneziani' # make a geometric features object that knows about the geometric data # cache up a couple of directories gf = GeometricFeatures(cacheLocation='../../geometric_data') fc = FeatureCollection() # ********* OSNAP array West ********* # Read in OSNAP West lon,lat OSNAPwestLonLat = csv.reader(open('./OSNAParrayWest20210322.csv', 'r')) # Skip 4 header lines next(OSNAPwestLonLat, None) next(OSNAPwestLonLat, None) next(OSNAPwestLonLat, None) next(OSNAPwestLonLat, None) coords = [] for line in OSNAPwestLonLat: coords.append([float(line[0]), float(line[1])]) feature = {} feature['type'] = 'Feature' feature['properties'] = {} feature['properties']['name'] = 'OSNAP section West' feature['properties']['tags'] = 'arctic_sections' feature['properties']['object'] = 'transect' feature['properties']['component'] = 'ocean' feature['properties']['author'] = author feature['geometry'] = {} feature['geometry']['type'] = 'LineString' feature['geometry']['coordinates'] = coords fcOW = FeatureCollection([feature]) fcOW.plot(projection='northpole') fc.merge(fcOW) # ********* OSNAP array East ********* # Read in OSNAP East lon,lat OSNAPeastLonLat = csv.reader(open('./OSNAParrayEast20210322.csv', 'r')) # Skip 4 header lines next(OSNAPeastLonLat, None) next(OSNAPeastLonLat, None) next(OSNAPeastLonLat, None) next(OSNAPeastLonLat, None) coords = [] for line in OSNAPeastLonLat: coords.append([float(line[0]), float(line[1])]) feature = {} feature['type'] = 'Feature' feature['properties'] = {} feature['properties']['name'] = 'OSNAP section East' feature['properties']['tags'] = 'arctic_sections' feature['properties']['object'] = 'transect' feature['properties']['component'] = 'ocean' feature['properties']['author'] = author feature['geometry'] = {} feature['geometry']['type'] = 'LineString' feature['geometry']['coordinates'] = coords fcOE = FeatureCollection([feature]) fcOE.plot(projection='northpole') fc.merge(fcOE) # "split" these features into individual files in the geometric data cache gf.split(fc) # update the database of feature names and tags write_feature_names_and_tags(gf.cacheLocation) # move the resulting file into place shutil.copyfile('features_and_tags.json', '../../geometric_features/features_and_tags.json') plt.show()
def main(): author = 'Xylar Asay-Davis, Alice Barthel, Nicolas Jourdain' tags = 'Antarctic;ISMIP6' # make a geometric features object that knows about the geometric data # cache up a couple of directories gf = GeometricFeatures('../../geometric_data') bedmap2_bin_to_netcdf('bedmap2.nc') fcContour1500 = get_longest_contour(contourValue=-1500., author=author) fc = FeatureCollection() lons = [-65., -25., -25., -65.] lats = [-80., -80., -77., -71.] fc.merge( shelf_polygon(lons, lats, name='ISMIP6 Weddell Sea', author=author, tags=tags, fcContour=fcContour1500)) lons = [-128., -128., -90., -90.] lats = [-76., -69., -69., -76.] fc.merge( shelf_polygon(lons, lats, name='ISMIP6 Amundsen Sea', author=author, tags=tags, fcContour=fcContour1500)) lons = [45., 45., 90., 90.] lats = [-70., -60., -60., -70.] fc.merge( shelf_polygon(lons, lats, name='ISMIP6 Amery Sector', author=author, tags=tags, fcContour=fcContour1500)) lons = [-22.5, -22.5, 22.5, 22.5] lats = [-75., -65., -65., -75.] fc.merge( shelf_polygon(lons, lats, name='ISMIP6 Dronning Maud Land', author=author, tags=tags, fcContour=fcContour1500)) lons = [110., 110., 130., 130.] lats = [-70., -60., -60., -70.] fc.merge( shelf_polygon(lons, lats, name='ISMIP6 Totten Region', author=author, tags=tags, fcContour=fcContour1500)) lons = [165., 165., 180., 180.] lats = [-80., -71., -73., -80.] fc_ross = shelf_polygon(lons, lats, name='ISMIP6 Western Ross Sea', author=author, tags=tags, fcContour=fcContour1500) lons = [-180., -180., -150., -150.] lats = [-80., -73., -77., -80.] fc_ross.merge( shelf_polygon(lons, lats, name='ISMIP6 Eastern Ross Sea', author=author, tags=tags, fcContour=fcContour1500)) old_props = fc_ross.features[0]['properties'] fc_ross = fc_ross.combine('ISMIP6 Ross Sea') props = fc_ross.features[0]['properties'] for prop in ['tags', 'zmin', 'zmax']: props[prop] = old_props[prop] fc.merge(fc_ross) fc.plot(projection='southpole') fc.to_geojson('ismip6_antarctic_ocean_regions.geojson') # "split" these features into individual files in the geometric data cache gf.split(fc) # update the database of feature names and tags write_feature_names_and_tags(gf.cacheLocation) # move the resulting file into place shutil.copyfile('features_and_tags.json', '../../geometric_features/features_and_tags.json') plt.show()
# create a GeometricFeatures object that points to a local cache of geometric # data and knows which branch of geometric_feature to use to download # missing data gf = GeometricFeatures('./geometric_data') # create a FeatureCollection containing all iceshelf regions wtih one of the # 27 IMBIE basin tags tags fc = FeatureCollection() for basin in range(1, 28): print('Adding feature from IMBIE basin {:d}'.format(basin)) basinName = 'Antarctica_IMBIE{:d}'.format(basin) tags = [basinName] # load the iceshelf regions for one IMBIE basin fcBasin = gf.read(componentName='iceshelves', objectType='region', tags=tags) # combine all regions in the basin into a single feature fcBasin = fcBasin.combine(featureName=basinName) # merge the feature for the basin into the collection of all basins fc.merge(fcBasin) # save the feature collection to a geojson file fc.to_geojson('Extended_Antarctic_Basins.geojson') if plot: fc.plot(projection='southpole') plt.show()
def main(): author = 'Xylar Asay-Davis' timTags = 'Antarctic;Timmermann' orsiTags = 'Antarctic;Orsi' # make a geometric fieatures object that knows about the geometric data # cache up a couple of directories gf = GeometricFeatures('../../geometric_data') bedmap2_bin_to_netcdf('bedmap2.nc') fcContour700 = get_longest_contour(contourValue=-700., author=author) fcContour800 = get_longest_contour(contourValue=-800., author=author) fc = FeatureCollection() fcWeddell = split_rectangle(lon0=-63., lon1=0., lat0=-80., lat1=-65., name='Weddell Sea', author=author, tags=timTags, fcContour=fcContour800) # get rid of the Weddell Sea because we're not that happy with this # definition, but keep the deep/shelf ones fcWeddell.features = fcWeddell.features[1:] fc.merge(fcWeddell) fcEW = split_rectangle(lon0=-20., lon1=25., lat0=-80., lat1=-55., name='Eastern Weddell Sea', author=author, tags=orsiTags, fcContour=fcContour800) fc.merge(fcEW) fcWW = split_rectangle(lon0=-63., lon1=-20., lat0=-80., lat1=-60., name='Western Weddell Sea', author=author, tags=orsiTags, fcContour=fcContour800) fc.merge(fcWW) # The Weddell feature is the sum of the Eastern and Western features before # splitting into shelf/deep fcWeddell = FeatureCollection() fcWeddell.features.append(fcEW.features[0]) fcWeddell.features.append(fcWW.features[0]) # now combine these into a single feature fcWeddell = fcWeddell.combine('Weddell Sea') props = fcWeddell.features[0]['properties'] props['tags'] = orsiTags props['zmin'] = -1000. props['zmax'] = -400. fc.merge(fcWeddell) # add the Weddell Sea back as the sum of Eastern and Western fc.merge( make_rectangle(lon0=-63., lon1=45., lat0=-80., lat1=-58., name='Weddell Sea', author=author, tags=orsiTags)) fc.merge( split_rectangle(lon0=-100., lon1=-63., lat0=-80., lat1=-67., name='Bellingshausen Sea', author=author, tags=timTags, fcContour=fcContour700)) fc.merge( split_rectangle(lon0=-140., lon1=-100., lat0=-80., lat1=-67., name='Amundsen Sea', author=author, tags=timTags, fcContour=fcContour800)) fc.merge( split_rectangle(lon0=-180., lon1=-140., lat0=-80., lat1=-67., name='Eastern Ross Sea', author=author, tags=timTags, fcContour=fcContour700)) fc.merge( split_rectangle(lon0=160., lon1=180., lat0=-80., lat1=-67., name='Western Ross Sea', author=author, tags=timTags, fcContour=fcContour700)) fc.merge( split_rectangle(lon0=25., lon1=160., lat0=-80., lat1=-62., name='East Antarctic Seas', author=author, tags=orsiTags, fcContour=fcContour800)) fc.merge( make_rectangle(lon0=-180., lon1=180., lat0=-80., lat1=-60., name='Southern Ocean 60S', author=author, tags=timTags)) fcSO = gf.read('ocean', 'region', ['Southern Ocean']) props = fcSO.features[0]['properties'] props['zmin'] = -1000. props['zmax'] = -400. fc.merge(fcSO) fc.plot(projection='southpole') fc.to_geojson('antarctic_ocean_regions.geojson') # "split" these features into individual files in the geometric data cache gf.split(fc) # update the database of feature names and tags write_feature_names_and_tags(gf.cacheLocation) # move the resulting file into place shutil.copyfile('features_and_tags.json', '../../geometric_features/features_and_tags.json') plt.show()