cosmo align hdr

function [map_x2y, map_y2x] = cosmo_align(x, y)