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
