notebook.community

Edit and run



In [ ]:

    
from shapely.geometry.polygon import Polygon
from shapely.geometry import MultiPolygon
import pyugrid
import json



In [ ]:

    
#url = 'http://www.smast.umassd.edu:8080/thredds/dodsC/FVCOM/NECOFS/Forecasts/NECOFS_GOM3_FORECAST.nc'
url = 'http://www.smast.umassd.edu:8080/thredds/dodsC/FVCOM/NECOFS/Forecasts/NECOFS_FVCOM_OCEAN_MASSBAY_FORECAST.nc'



In [40]:

    
ug = pyugrid.UGrid.from_ncfile(url)



In [16]:

    
lon = ug.nodes[:,0]
lat = ug.nodes[:,1]
nv = ug.faces[:]



In [69]:

    
mp = MultiPolygon([Polygon(zip(lon[element],lat[element])) for element in nv])



In [47]:

    
with open('ugrid.json','w') as f:
    f.write(json.dumps(mp))









    Out[47]:





u'/home/usgs'