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')
Exemple #3
0
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
        )
Exemple #4
0
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)
Exemple #5
0
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)