Skip to content

Splitting image to radial and non-radial components ImageJ macro

License

Notifications You must be signed in to change notification settings

UU-cellbiology/radialitymap

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

13 Commits
 
 
 
 
 
 

Repository files navigation

Radiality map

DOI
ImageJ/FIJI macro that splits image to radial and non-radial components using output of OrientationJ plugin. radialitymap logo

How to install and run macro

  1. You need to download and install ImageJ or FIJI on your computer first.
  • you will need ImageJ version 1.48h or higher for macro to work.
  1. Download OrientationJ_.jar and place it in the "plugins" folder of ImageJ/FIJI.
  2. Download macro file "radiality_script...ijm" and store it somewhere on your hard disk.
  3. Load your image to ImageJ/FIJI. Currently macro works only with 8-/16-/32-bit images.
  • regular RGB images can be converted in ImageJ using Image->Type->...bit command
    example image
  1. Choose the origin of coordinates (center with respect to what image considered radial) using "Point" selection tool: point tool example image with dot
  2. Go to Plugins->Macros->Runs.. menu and choose previously downloaded macro file ("radiality_script...ijm" from step 3).
  • If you try Plugins->Macros->Install.. and choose macro file, it will appear as a separate row in Plugins->Macros menu
    parameters window
  1. Choose values of parameters. For "gaussian window" a good initial guess is an approximate thickness of lines on your image. Try values 2 times higher or less and see result. The difference in "Map calculation method" is not so critical, pick the best one judging from final result.
  2. Done!

Output images

Macro/script generates five extra images. Here is description:

  1. Local orientation map (32-bit) generated by OrientationJ with local directionality angle in each pixel in radians:
    orientation map
  2. Local radiality map (32-bit) ranging from 0 to 1. It is absolute value of cosine of difference between the local orientation angle and the angle of vector drawn from the new origin of coordinates (marked by Point tool) to the current pixel position
    radiality map
  3. Multiplication of original image and radiality map (32-bit), i.e. "radial" intensity component. By default it uses "Fire" color lookup table, but can be anything else.
    radiality mult map
  4. Local non-radiality map (32-bit) ranging from 0 to 1. It is just one minus radiality map
    nonradiality map
  5. Multiplication of original image and non-radiality map (32-bit), i.e. "non-radial" intensity component. By default it uses "Fire" color lookup table, but can be anything else.
    nonradiality mult map

In general, you can combine radial/non-radial images (3 and 5) together with different LUTs using Image->Color->Merge Channels.. command of ImageJ (cyan=radial, magenta=non-radial):
composite map

Example application to microtubules shapes

Here is image of a cytoplast at round micropattern kindly provided by Manuel Thery (gray) and separation result (cyan=radial, magenta=non-radial):
cytoplast

It is easy to see some artefact at the center. The density of microtubules is high there, they overlap over each other and so it is impossible to distinguish individual filaments shapes. But line features can be enchanced using Hessian FeatureJ function (smallest eighen value, no absolute comparison), so shapes can be distinguished a bit better:
cytoplast hessian

How to cite this macro?

You can use Zenodo DOI for the link to the specific release (see DOI badge on the top of the page):

Katrukha E. 2019, RadialityMap macro for ImageJ, v0.2, Zenodo, doi:10.5281/zenodo.6817259

Updates history

2019.10.17 Thanks to Junnan Fang feedback macro is updated to the current version of OrientationJ


Developed in Cell Biology group of Utrecht University.
Email katpyxa @ gmail.com for any questions/comments/suggestions.