|
34
|
1 Interfacing with Matlab and Octave
|
|
|
2 ==================================
|
|
|
3
|
|
|
4 Thanks to the REST API of Orthanc, it is easy to access DICOM images
|
|
|
5 from Matlab or Octave, as depicted in the following sample image:
|
|
|
6
|
|
|
7 .. image:: ../images/Matlab.png
|
|
|
8 :align: center
|
|
|
9 :width: 470px
|
|
|
10
|
|
|
11 Both Matlab and Octave have access to HTTP servers thanks to their
|
|
|
12 built-in `urlread() function
|
|
|
13 <http://nl.mathworks.com/help/matlab/ref/urlread.html>`__. Once must
|
|
|
14 simply install a Matlab/Octave library to decode JSON files. The
|
|
66
|
15 `JSONLab toolkit <https://github.com/fangq/jsonlab>`__ works perfectly
|
|
|
16 to this end.
|
|
34
|
17
|
|
|
18 .. highlight:: matlab
|
|
|
19
|
|
|
20 Using JSONlab, the following code will download and display a DICOM image::
|
|
|
21
|
|
|
22 SERIES = 'ae164c84-e5bd0366-ba937a6d-65414092-f294d6b6';
|
|
66
|
23 URL = 'http://demo.orthanc-server.com/';
|
|
34
|
24
|
|
|
25 # Get information about the instances in this DICOM series
|
|
|
26 instances = loadjson(urlread([ URL '/series/' SERIES '/instances' ]));
|
|
|
27
|
|
|
28 # Select one slice from the series
|
|
66
|
29 instance = instances(1,1).ID
|
|
34
|
30
|
|
|
31 # Decode the slice with Orthanc thanks to the "/matlab" URI
|
|
|
32 slice = eval(urlread([ URL '/instances/' instance '/matlab' ]));
|
|
|
33
|
|
|
34 # Compute the maximum value in this slice
|
|
|
35 max(max(slice))
|
|
|
36
|
|
|
37 # Display the slice
|
|
|
38 imagesc(slice)
|
|
|
39
|
|
|
40 # Annotate the graph with the patient name and ID
|
|
|
41 tags = loadjson(urlread([ URL '/instances/' instance '/tags?simplify' ]));
|
|
|
42 title([ 'This is a slice from patient ' tags.PatientID ' (' tags.PatientName ')' ])
|
