Subscribe RSS

Posts Tagged ‘python’

ImageJ in a Sage notebook, an example of Python calling Java

November 10th, 2009 by fmn | 11 Comments | Filed in Enseignement, Research

Sage is a wonderful system: a big set of mathematical languages and libraries glued with Python. This is quite remarquable, but i use frequently Java for my image processing and analysis applications. In particular, the ImageJ framework allows to write plugins, making easy the application distribution. The aim is thus to easily access a Java library from a Python code. The complete method is described in this post, allowing to use ImageJ from a Sage notebook.

As Sage rely on Python, the key is to find a way to call Java code from Python. I found several solutions (excluding do-it-yourself stuff with RMI):

  • Javaclass : can’t make it work. However, seems to need work in order to pass data between Java and Python worlds.
  • JCC : a little bit clumsy, it particular it needs to add very large options when launching Python.
  • JPype : the simplest. A Python module with an easy usage.

JPype installation

It is quite regular:

unzip JPype-0.5.4.1.zip
cd JPype-0.5.4.1
sage -python setup.py build

Here a bug usually appears. JPype rely on the sets module, which is neither included in recent Python. There is not a lot of modification. You can download a modified version here. Beware, you must still follow the README to provide the path to your JDK installation. If the compilation is OK (you can send me a mail if not), the module is installed with:

sage -python setup.py install

Now we can work inside a Sage notebook (downloadable here).

JPype use

Quite simple:

  1. import the module,
  2. launch a java virtual machine,
  3. import one or more Java packages.

from jpype import *
startJVM(’/usr/lib/jvm/java-6-openjdk/jre/lib/i386/client/libjvm.so’, ‘-Djava.class.path=/home/fredn/.libjar/ij.jar’)
ij = JPackage(’ij’)

The first parameter given to startJVM() is the path to the dynamical library of the JVM. To find it, you can do:

locate libjvm.so

The second parameter allows to extend the classpath in order to add libraries, here Image.

At this point, all the methods of Java-package ij are available, directly. We can by example load an image and display it in the notebook:

im = ij.IJ.openImage(”http://rsb.info.nih.gov/ij/images/lena.jpg”) ij.IJ.save(im, os.getcwd()+’/tmp.png’)
JPype allows to convert Java objects in Python objects (and vis versa), very easily. Thus a Python object is constructed from a Java one, all the methods available:
lena = ij.IJ.openImage(”http://rsb.info.nih.gov/ij/images/lena.jpg”).getProcessor().convertToByte(True) lena
<jpype._jclass.ij.process.ByteProcessor object at 0xb7ad52c>
<jpype._jclass.ij.process.ByteProcessor object at 0xb7ad52c>

We can now find Lena edges by calling a ByteProcessor method:

lenaedges = lena.duplicate()
lenaedges.findEdges()
lenaedges
<jpype._jclass.ij.process.ByteProcessor object at 0xb313652c>
<jpype._jclass.ij.process.ByteProcessor object at 0xb313652c>

To display side by side Lena and its edges, i need to get the pixel array of each ImageProcessor, by calling getIntArrray():

lenaedges.getIntArray()
<jpype._jarray.int[][] object at 0xb313a20c>
<jpype._jarray.int[][] object at 0xb313a20c>

This object is Python indexable, as a regular array. I can give it directly to the Python functions that display images. Here, the implicit conversion between primitive types is realized by JPype:

import pylab as pl
pl.gray()
pl.figure()
pl.subplot(1,2,1)
pl.imshow(lena.getIntArray())
pl.subplot(1,2,2)
pl.imshow(lenaedges.getIntArray())
pl.savefig(’tmp.png’)

(more…)

Tags: , , , , , , ,

Détection de motifs : exploitation de la phase (suite de l’inter-corrélation)

December 18th, 2008 by fmn | 5 Comments | Filed in Popularization

Translate original post with Google Translate

Dans un billet précédent, nous avons vu comment (en théorie) une inter-corrélation pouvait permettre de localiser un objet dans une image. La conclusion était que l’inter-corrélation n’est pas très efficace car sa valeur dépend énormément du niveau de gris des images et assez peu de leur information spatiale. Nous allons voir comment corriger cela.

Le notebook sage de ce billet est : detecteur_de_motifs_phase.sws

(more…)

Tags: , , , , ,

Détection de motifs par intercorrélation

November 28th, 2008 by fmn | 7 Comments | Filed in Popularization

Translate original post with Google Translate

Avec cet article, je commence une petite série destinée à expliquer quelques méthodes permettant de trouver des objets dans une image. Toute les méthodes seront accompagnées d’illustrations reproductible sous sage. Les codes sources seront également téléchargeables sous la forme d’un notebook sage.

Notebook

Vous pouvez télécharger le notebook sage contenant le code complet présenté ici accompagné des images de test : detecteur_de_motifs_base__sur_une_intercorrelation.sws

Objectif

Pour ce premier article, imaginons que j’ai l’image suivante (que j’appelle image reférence)  :
Image de référence

Je pense que vous avez tous remarqués le mignon petit ourson au centre de l’image. Tentons de le retrouver. Il faut d’abord posséder une image contenant un exemplaire de l’objet à chercher, en voici une que j’appelle image motif :
ourson
(more…)

Tags: , , , , , ,