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')