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')
def project(img_from, img_to, **options): with img_from.make_child('map') as mapfile: with img_from.make_child('lis') as fromlist: write_file_list(unicode(fromlist), [unicode(img_from)]) isis.mosrange( fromlist = fromlist, to = mapfile, precision = 2, projection = 'equirectangular' ) isis.cam2map( from_ = img_from, to = img_to, pixres = 'map', map = mapfile )
def process_image(image, title, caption): image = ImageName(splitext(image)[0]) isis.lronac2isis.print_cmd(from_=image.IMG, to=image.cub) isis.lronac2isis(from_=image.IMG, to=image.cub) isis.spiceinit(from_=image.cub) #lronaccal default is IOF calibration isis.lronaccal(from_=image.cub, to=image.cal.cub) isis.trim(from_=image.cal.cub, to=image.trim.cub, left=45, right=45) write_file_list(image.map.lis, [image.trim.cub]) isis.mosrange(fromlist=image.map.lis, to=image.nac_eqr.map, precision=2, projection='equirectangular') isis.cam2map(from_=image.trim.cub, to=image.proj.cub, pixres='map', map=image.nac_eqr.map) create_yml(image, title) create_caption(image, caption)
def process_images(images, title, caption): """ Gets passed "images", with multiple "image" the extensions are called "ext_set" the img_set has the imagename + ext_set look at named tuples for future reference """ images = [ImageName(splitext(image)[0]) for image in images] with IsisPool() as isis_pool: for image in images: isis_pool.lrocnac2isis(from_=image.IMG, to=image.cub) with IsisPool() as isis_pool: for image in images: isis_pool.spiceinit(from_=image.cub) with IsisPool() as isis_pool: for image in images: #lronaccal default is IOF calibration isis_pool.lronaccal(from_=image.cub, to=image.cal.cub) with IsisPool() as isis_pool: for image in images: isis_pool.trim(from_=image.cal.cub, to=image.trim.cub, left=45, right=45) with IsisPool() as isis_pool: for image in images: write_file_list(image.map.lis, [image.trim.cub]) isis_pool.mosrange(fromlist=image.map.lis, to=image.nac_eqr.map, precision=2, projection='equirectangular') with IsisPool() as isis_pool: for image in images: isis_pool.cam2map(from_=image.trim.cub, to=image.proj.cub, pixres='map', map=image.nac_eqr.map) create_yml(image, title) create_caption(image, caption)