def main():
images = iglob('*.IMG')
isis.maptemplate(map='D_mdis_wac_eqr.map', projection='equirectangular',
clon=0.0, clat=0.0, resopt='mpp', resolution=1000, rngopt='user',
minlat=-40.0, maxlat=40.0, minlon=0.0, maxlon=360.0)
for img_name in images:
(cub_name, cal_name, pho_name, trim_name, proj_name) = file_variations(img_name,
['.cub', '.cal.cub', '.pho.cub', '.trim.cub', '.proj.cub'])
#basic_img_proc(img_name)
isis.photomet(from_=cal_name, to=pho_name, PHTNAME='HAPKEHEN',
Wh=0.215984749, Hh=0.075, B0=2.3, Theta=15.78892162,
HG1=0.206649235, HG2=0.811417942, zerob0standard='false',
NORMNAME='Albedo', Incref=30.0, Thresh=10E30, Incmat=0.0, Albedo=1.0)
isis.trim(from_=pho_name, to=trim_name, left=10)
isis.cam2map(from_=trim_name, to=proj_name, pixres='map',
map='D_mdis_wac_eqr.map', defaultrange='map')
# sort by resolution: output is .proj.cub
images = iglob('*.trim.cub')
sorted_images = sort_images(images, get_pixel_scale, reverse=True)
ext_len = -len('.trim.cub')
sorted_images = [img_name[:ext_len] + '.proj.cub'
for img_name in sorted_images]
write_file_list('mosaic.list', sorted_images)
isis.automos(fromlist='mosaic.list', mosaic=mosaic_name)
def create_mapfile(img_name):
""" Create a mapfile in isis from a LROC NAC
Args: img_name
"""
(cub_name, trim_name) = file_variations(img_name,
['.cub', '.trim.cub'])
isis.lronac2isis(from_=img_name, to=cub_name)
isis.spiceinit(from_=cub_name)
isis.trim(from_=cal_name, to=trim_name, left=45, right=45)
write_file_list('map.lis', trim_name)
isis.mosrange(fromlist='map.lis', to=img_name+'.map', precision=2,
projection='equirectangular')
