def match_multi_test(): shandong = gio.read_shp('../data/shape/shandong.shp') shandong = shandong.to_crs(3857) info = gutil.shp2box(shandong, (2048, 1536), 0.05, 1) paper = gnp.frombox(*info, dtype=np.int16) fs = glob('../data/modis/*.hdf') rasters = [gio.read_hdf(i, 0) for i in fs] gmt.match_multi(rasters, paper, out='in') plt.imshow(paper) plt.show()
def run(self, ips, imgs, para=None): table = TableManager.get(para['temp']).get_subtab() box = gutil.shp2box(table, para['scale'] or (para['width'], para['height']), para['mar']) chans = ['Channel %s' % i for i in range(ips.get_nchannels())] chans = [chans.index(i) for i in para['chans']] order = {'nearest': 0, 'linear': 1}[para['order']] rst = gmt.match_multi(imgs, box, chans, step=para['step'], order=order) IPy.show_img([rst], ips.title + '-merge')
def match_multi_test(): shandong = gio.read_shp('../data/shape/shandong.shp') shandong = shandong.to_crs(3857) info = gutil.shp2box(shandong, (2048, 1536), 0.05) fs = glob('../data/modis/*.hdf') rasters = [gio.read_hdf(i, 0) for i in fs] rst = gmt.match_multi(rasters, info, 0) plt.imshow(rst) plt.show()
def run(self, ips, imgs, para=None): ipst = ImageManager.get(para['temp']) chans = ['Channel %s' % i for i in range(ips.get_nchannels())] chans = [chans.index(i) for i in para['chans']] order = {'nearest': 0, 'linear': 1}[para['order']] rst = gmt.match_multi(imgs, ipst.img, chans, step=para['step'], order=order) ipst.update()
def run(self, ips, imgs, para=None): prj = None if para['wkt'] != '': with open(para['wkt']) as f: prj = f.read() crs = gutil.makecrs(prj or para['crs']) table = gmt.build_index(imgs).to_crs(crs) box = gutil.shp2box(table, para['scale'] or (para['width'], para['height']), para['mar']) chans = ['Channel %s' % i for i in range(ips.get_nchannels())] chans = [chans.index(i) for i in para['chans']] order = {'nearest': 0, 'linear': 1}[para['order']] rst = gmt.match_multi(imgs, box, chans, step=para['step'], order=order) IPy.show_img([rst], ips.title + '-merge')
def run(self, ips, imgs, para = None): temp = WindowsManager.get(para['temp']).ips rst = gmt.match_multi(imgs, temp.img.getbox()) IPy.show_img([rst], ips.title+'-match')