parser.add_argument('--helper','-H',nargs=1,default=[0],type=int,help='Number of helper contours to add between each contour') parser.add_argument('--datum','-d',nargs=1,default=[0.0],type=float,help='Height datum offset for contours') parser.add_argument('--idatum','-D',nargs=1,default=[0],type=int,help='Height datum offset for index contours as an integer number of standard contours') parser.add_argument('--index','-x',nargs=1,default=[0],type=int,help='Number of contours between index contours') parser.add_argument('--interval','-c',nargs=1,default=[5.0],type=float,help='Contour interval') parser.add_argument('--add_unfiltered','-u',action='store_true',default=False,help='Add unfiltered contours') parser.add_argument('--map_bounds','-b',nargs=4,type=float,help='Map bounds in UTM coordinates (W S E N)') parser.add_argument('--UTMZone','-z',nargs=1,default=[30],type=int,help='UTM Zone of map (e.g 30)') cargs = parser.parse_args() files = os.listdir(cargs.input[0]) maptiles = [] bounds = [+np.inf,+np.inf,-np.inf,-np.inf] #Initialise output cOMAP = omap.omap(map_type='contour') cOMAP.UTMZone = cargs.UTMZone[0] cOMAP.setGMLCRS(init=DTMCRScode) cOMAP.setOMAPCRS(proj='utm',zone=cargs.UTMZone[0],ellps='WGS84') for file in files: m = re.match('([a-z])([a-z])(\d{2})(\d{2})_DTM_2M\.asc',file) if m: with open(os.path.join(cargs.input[0],file)) as mapfile: #Find metadata for each data file metadata = {'filename':file} while True: line = mapfile.readline().strip().split() if len(line) > 2: #Break when we finish the header and reach the data break metadata[line[0]] = int(line[1])
rc = [] for node in nodelist: if node.nodeType == node.TEXT_NODE: rc.append(node.data) return ''.join(rc) if len(sys.argv) > 1: gmlname = sys.argv[1] else: gmlname = 'test.gml' print('Reading features from ' + gmlname) unrecognisedFeatures = {} #misfitOMAPs = {}; mainOMAP = omap.omap() mainOMAP.setGMLCRS(init=GMLCRScode) #mainOMAP.setOMAPCRS(OMAPCRScode) objectCount = 0 #Setup streaming XML reader srcstrm = pulldom.parse(gmlname) #Find a coordinate (any coordinate) to choose the UTM Zone for event, node in srcstrm: if event == pulldom.START_ELEMENT and node.tagName == 'gml:coordinates': srcstrm.expandNode(node) mainOMAP.setUTMZoneFromCoord(node) break #Setup streaming XML reader again
type=float, help='Map bounds in UTM coordinates (W S E N)') parser.add_argument('--UTMZone', '-z', nargs=1, default=[30], type=int, help='UTM Zone of map (e.g 30)') cargs = parser.parse_args() files = os.listdir(cargs.input[0]) maptiles = [] bounds = [+np.inf, +np.inf, -np.inf, -np.inf] #Initialise output cOMAP = omap.omap(map_type='contour') cOMAP.UTMZone = cargs.UTMZone[0] cOMAP.setGMLCRS(init=DTMCRScode) cOMAP.setOMAPCRS(proj='utm', zone=cargs.UTMZone[0], ellps='WGS84') for file in files: m = re.match('([a-z])([a-z])(\d{2})(\d{2})_DTM_2M\.asc', file) if m: with open(os.path.join( cargs.input[0], file)) as mapfile: #Find metadata for each data file metadata = {'filename': file} while True: line = mapfile.readline().strip().split() if len( line
def getText(nodelist): rc = [] for node in nodelist: if node.nodeType == node.TEXT_NODE: rc.append(node.data) return ''.join(rc) if len(sys.argv) > 1: gmlname = sys.argv[1] else: gmlname = 'test.gml' print('Reading features from '+gmlname) unrecognisedFeatures = {}; #misfitOMAPs = {}; mainOMAP = omap.omap(); mainOMAP.setGMLCRS(init=GMLCRScode) #mainOMAP.setOMAPCRS(OMAPCRScode) objectCount = 0 #Setup streaming XML reader srcstrm = pulldom.parse(gmlname) #Find a coordinate (any coordinate) to choose the UTM Zone for event,node in srcstrm: if event == pulldom.START_ELEMENT and node.tagName == 'gml:coordinates': srcstrm.expandNode(node) mainOMAP.setUTMZoneFromCoord(node) break