array(T).astype('float')
T = array(T).astype('float')
T.shape
figure()
plot(T[:,0],T[:,1])
clf()
plot(T[:,0],T[:,1],'s')
subplot(221)
plot(T[:,0],T[:,1],'s')
subplot(222)
plot(T[:,2],T[:,1],'s')
subplot(223)
plot(T[:,2]*T[:,0],T[:,1],'s')
644*14*5.3**2
import pyregion
r = pyregion.read_region('../crop_regions.reg')
r = pyregion.open('../crop_regions.reg')
r
r[0]
R = r[0]
R.coord_list
get_ipython().magic(u"pwd ")
3e5/50
0.02/5.
1/(0.02/5.)
get_ipython().system(u"ls -F /var")
np.fft.fftfreq([3,3])
help(np.fft.fft2)
help(np.fft.fftfreq)
np.fft.fftfreq(5)
np.fft.fftfreq(5,7.2)
return Header.fromtextfile("test.header")
def print_region(r):
for i, l in enumerate(r):
print("[region %d]" % (i + 1))
print()
print("%s; %s(%s)" % (l.coord_format, l.name, ", ".join([str(s) for s in l.coord_list])))
print(l.attr[0])
print(", ".join(["%s=%s" % (k, v.strip()) for k, v in list(l.attr[1].items())]))
print()
if __name__ == "__main__":
print("** coordinate in FK5 **")
print()
region_name = "test01_print.reg"
# region_name = "test_text.reg"
# region_name = "test01.reg"
r = read_region(open(region_name).read())
print_region(r)
print()
print()
print("** coordinate in image **")
print()
header = test_header()
r = read_region_as_imagecoord(open(region_name).read(), header)
print_region(r)
for i, l in enumerate(r):
print("[region %d]" % (i + 1))
print()
print(
"%s; %s(%s)" %
(l.coord_format, l.name, ", ".join([str(s)
for s in l.coord_list])))
print(l.attr[0])
print(", ".join(
["%s=%s" % (k, v.strip()) for k, v in list(l.attr[1].items())]))
print()
if __name__ == "__main__":
print("** coordinate in FK5 **")
print()
region_name = "test01_print.reg"
#region_name = "test_text.reg"
#region_name = "test01.reg"
r = read_region(open(region_name).read())
print_region(r)
print()
print()
print("** coordinate in image **")
print()
header = test_header()
r = read_region_as_imagecoord(open(region_name).read(), header)
print_region(r)
