Search This Blog

Friday, 17 March 2017

Python program to print countryName, countryCode, minLat, maxLat, minLong,maxLong

from osgeo import ogr

shapefile = ogr.Open('D:/vector/TM_WORLD_BORDERS-0.3/TM_WORLD_BORDERS-0.3.shp')

layer = shapefile.GetLayer(0)
countries =[]
for i in range(layer.GetFeatureCount()):
    feature = layer.GetFeature(i)
    countryCode = feature.GetField('ISO3')
    countryName = feature.GetField('NAME')
    print '{}:{}'.format(countryCode, countryName)
    geometry = feature.GetGeometryRef()
    minLong, maxLong, minLat, maxLat = geometry.GetEnvelope()
    countries.append((countryName, countryCode, minLat, maxLat, minLong,maxLong))


countries.sort()
for name,code, minLat, maxLat, minLong, maxLong in countries:
    print '{} ({}) lat = {:.4f}..{:.4f}..Long={:.4f}..{:.4f}'.format(name, code, minLat, maxLat, minLong, maxLong)

No comments:

Post a Comment