Python C/C++ IPMU IRAF

Size: px

Start display at page:

Download "Python C/C++ IPMU IRAF"

ふみなすわ
6 years ago
Views:

1 Python C/C++ IPMU IRAF

2 Python Swig Numpy array Image

3 Python swig numpy pyfits 2.3 pyds9 1.1 svn co hjp://svn.scipy.org/svn/numpy/tags/1.5.0/doc/swig swig/numpy.i /usr/local/share/swig/1.3.40/python swig

4 Python

5 Hello Python Python python/hello.py % python Python (r266:84292, Nov , 16:27:48) [GCC (Red Hat )] on linux2 Type help, copyright, credits or license for more information. >>> print Hello Python! Hello Python! >>> print >>> ^D (Crtl+D) % % python hello.py Hello Python! 8 %

6 python/loop.py a = [ abc, def, ghi, jkl ] for b in a: print b for for i in range(5): print i range(5) [0,1,2,3,4]

7 I/O python/io.py file = open( test.txt ) file = open( out.txt, w ) line = file.readline() file.write(line)

8 python/func.py def times(x, y): return x * y import func from func import * x = func.times(3.14, 4)

9 pyfits FITS hdulist = pyfits.open( test.fits ) hdulist[0].header[ OBJECT ] hdulist[0].data[y, x] pyfits/sample.py numpy array 0- indexed ds9 [y, x]

10 pyds9 ds9 X Public Access (XPA) XPA Access Points d = ds9.ds9() ds9 d.set( file test.fits ) pyds9/sample.py array d.set( array [xdim=2080,ydim=4100,bitpix=-32], data)

11 Python Numpy hjp:// pyfits/stats.py Python Numpy

12 Python C/C++ Python C/C++ Python C/C++

13 Python C/C++ boost.python swig ctypes swig

14 Swig C swig/example.*

15 C /* File : example.c */ swig/example.c #include <time.h> double My_variable = 3.0; int fact(int n) { if (n <= 1) return 1; else return n*fact(n-1); } int my_mod(int x, int y) { return (x%y); } char *get_time() { time_t ltime; time(&ltime); return ctime(&ltime); }

16 Swig swig/example.* C Swig

17 Swig swig/example.i /* example.i */ %module example Python %{ /* Put header files here or function declarations like below */ extern double My_variable; extern int fact(int n); extern int my_mod(int x, int y); extern char *get_time(); %} extern double My_variable; extern int fact(int n); extern int my_mod(int x, int y); extern char *get_time();

18 Swig C Swig swig % swig -python example.i example.py, example_wrap.c % gcc fpic c example.c example_wrap.c I/adc/local/include/ python2.6/ example.o, example_wrap.o shared object % gcc -shared example.o example_wrap.o -o _example.so _example.so swig/example.* compile.sh

19 import example swig/sample.py print example.cvar.my_variable print example.fact(4) print example.my_mod(10,3) print example.get_time() <>.cvar.< > Python

20 Swig output/example.* void add(double a, double b, double *result) %include typemaps.i %apply double *OUTPUT { double *result }; typemaps.i void mod(int a, int b, int *quo, int *rem) %include typemaps.i # %apply int *OUTPUT { int*quo, int *rem };

21 Swig import example output/sample.py a = example.add_0(3,4) print a b = example.add(3,4) print b quo, rem = example.mod(7,3) print quo, rem

22 Swig C++ vector double average(std::vector<int> v) std::vector<double> half(const std::vector<double>& v) %include std_vector.i %template(intvector) std::vector<int>; %template(doublevector) std::vector<double>; vector/example.*

23 Swig from example import * iv = IntVector(4) for i in range(0,4): iv[i] = i print average(iv) vector<int> vector/sample.py print average([0,1,2,3]) print half([1,2,3]) Swig

24 Numpy array C Numpy array rms C double rms(double *seq, int n) Python v = rms(array) numpy/rms.* numpy.itypemaps

25 Numpy array C %module rms %{ #define SWIG_FILE_WITH_INIT #include rms.h %} %include numpy.i %init %{ import_array(); %} numpy/rms.i %apply (double *IN_ARRAY1, int DIM1) {(double *seq, int n)}; %include rms.h 1 array numpy_swig.pdf

26 CNumpy array ezrange/ezrange.i Cnumpy array Python %apply (double *IN_ARRAY1, int DIM1) {(double *seq, int n)}; %apply (int *ARGOUT_ARRAY1, int DIM1) {(int *rangevec, int n)}; array Python

27 array 2 array(=) %apply (double *IN_ARRAY1, int DIM1) %apply (double *IN_ARRAY2, int DIM1, DIM2) C array double mean(double *array, int n, int m) numpy2/stats.* Numpy

28 Python C++

29 Image image0/image.* Numpy array Image Numpy array vector class Image { int npx, npy, npix; double *data; public: Image(); Image(double *array, int n, int m); ~Image(); int getnpx() const; int getnpy() const; int getnpix() const; void getnaxis(int *npx, int *npy); double getvalue(int x, int y) const; };

30 Image image0/image.* Image::Image(double *array, int n, int m) { data = new double[n*m]; // for (int j = 0; j < n; j++) { for (int i = 0; i < m; i++) { data[i+j*m] = array[i+j*m]; // } } npx = m; npy = n; npix = n*m; }

31 Image image0/sample.py pyfits Image img = Image.Image(hdulist[0].data) ds9 image01/image.*

32 median median image1/image.* medianmedian

33 Image::Image(int npx, int npy) void Image::setValue(int x, int y double v) void Image::getData(float *oarray, int n) numpy array Image& Image::operator/=(double a) a double Image::median() median static Image Image::getMedian(std::vector<Image>& v) vmedian image1/image.*

34 ds9 image2/image.* n x n

35 Image& Image::operator/=(const Image& a) image2/image.* Image Image::boxavg(int width) width x width std::vector<point> findpeaks(double thres) thres vector Point class Point { int x, y; public: Point(int x, int y); int getx(); int gety(); }

36 Image median Gaussian

37 Swig C/C++ Python Swig Python C++

untitled

untitled 1 IPCC1988 WMOUNEP 2 3 4 CO2 CO2 2 2 2 2 2 2 CO2 CO2 CO2 CO2 5 6 CO2 CO2 7 5 A. A-1 A-2 B. 1 6 (a)(b)(c) (d)(e)(f) (g)(h)(i) (j)(k)(l) B-1 3 B-2 C. 10 3 2 0 9 9 1 C. 5 8,2009,,