views:

155

answers:

5

Does Matlab have any support for hashtables?


Some background

I am working on a problem in Matlab that requires a scale-space representation of an image. To do this I create a 2-D Gaussian filter with variance sigma*s^k for k in some range., and then I use each one in turn to filter the image. Now, I want some sort of mapping from k to the filtered image.

If k were always an integer, I'd simply create a 3D array such that:

arr[k] = <image filtered with k-th guassian>

However, k is not necessarily an integer, so I can't do this. What I thought of doing was keeping an array of ks such that:

arr[find(array_of_ks_ = k)] = <image filtered with k-th guassian>

Which seems pretty good at first thought, except I will be doing this lookup potentially a few thousand times with about 20 or 30 values of k, and I fear that this will hurt performance.

I wonder if I wouldn't be better served doing this with a hash table of some sort so that I would have a lookup time that is O(1) instead of O(n).


Now, I know that I shouldn't optimize prematurely, and I may not have this problem at all, but remember, this is just the background, and there may be cases where this is really the best solution, regardless of whether it is the best solution for my problem.

+3  A: 

Matlab does not have support for hashtables. EDIT Until r2010a, that is; see @Amro's answer.

To speed up your look-ups, you can drop the find, and use LOGICAL INDEXING.

arr{array_of_ks==k} = <image filtered with k-th Gaussian>

or

arr(:,:,array_of_ks==k) = <image filtered with k-th Gaussian>

However, in all my experience with Matlab, I've never had a lookup be a bottleneck.


To speed up your specific problem, I suggest to either use incremental filtering

arr{i} = GaussFilter(arr{i-1},sigma*s^(array_of_ks(i)) - sigma*s^(array_of_ks(i-1)))

assuming array_of_ks is sorted in ascending order, and GaussFilter calculates the filter mask size based on the variance (and uses, 2 1D filters, of course), or you can filter in Fourier Space, which is especially useful for large images and if the variances are spaced evenly (which they most likely aren't unfortunately).

Jonas
+2  A: 

You could use java for it.

In matlab:

dict = java.util.Hashtable;
dict.put('a', 1);
dict.put('b', 2);
dict.put('c', 3);
dict.get('b')

But you would have to do some profiling to see if it gives you a speed gain I guess...

tauran
+8  A: 

Containers.Map:

k = {'Jan', 'Feb', 'Mar', 'Apr', 'May', 'Jun', ...
  'Jul', 'Aug', 'Sep', 'Oct', 'Nov', 'Dec', 'Annual'};

v = {327.2, 368.2, 197.6, 178.4, 100.0,  69.9, ...
  32.3,  37.3,  19.0,  37.0,  73.2, 110.9, 1551.0};

rainfallMap = containers.Map(k, v)
rainfallMap = 
  containers.Map handle
  Package: containers

  Properties:
        Count: 13
      KeyType: 'char'
    ValueType: 'double'
  Methods, Events, Superclasses
Amro
Wow, I didn't know that! +1. Do you know whether they're much faster than logical indexing?
Jonas
Containers.Map was added in MATLAB 7.7 (R2008b) see http://www.mathworks.com/access/helpdesk/help/techdoc/rn/brqyzax-1.html . New in R2010a is a constructor to specify key type as well as value type. M = containers.Map('KeyType', kType, 'ValueType', vType)
zellus
@Jonas: I haven't used them extensively, it would be interesting to see how they compare against using logical indexing for lookup..
Amro
@zellus: my bad, I guess I only noticed them recently :)
Amro
@zellus,@amro: Isn't it annoying how there is no history of the commands in Matlab?
Jonas
@Jonas: Good point, a command history might be interesting. Overall the documentation is outstanding, compared with other products. I was waiting for the new constructor, therefore I knew that the maps where available before R2010a.
zellus
+1  A: 

It's a little clugey, but I'm surprised nobody has suggested using structs. You can access any struct field by variable name as struct.(var) where var can be any variable and will resolve appropriately.

dict.a = 1;
dict.b = 2;

var = 'a';

display( dict.(var) ); % prints 1
Mark E
It would break if you use a number as a fieldname: `dict.('2')`: http://www.mathworks.com/access/helpdesk/help/techdoc/matlab_prog/br04bw6-38.html#br1v5bu-1
Amro
@Mark: Also, the variables have to be integers: `dict.(['k',num2str(1)])` works, but `dict.(['k',num2str(1.1)])` fails, and if the values are integers, you can use them to index directly. It's a nice idea otherwise.
Jonas
@Amro, @Jonas, fair points, if the the keys were integers you wouldn't need to use this trick (an array would make more sense)...if the keys are arbitrary floats this is a little more challenging, but I'd prefix with a letter and replace the `.` with an `_`.
Mark E
+7  A: 

Matlab R2008b (7.7)’s new containers.Map class is a scaled-down Matlab version of the java.util.Map interface. It has the added benefit of seamless integration with all Matlab types (Java Maps cannot handle Matlab structs for example) as well as the ability since Matlab 7.10 (R2010a) to specify data types.

Serious Matlab implementations requiring key-value maps/dictionaries should still use Java’s Map classes (java.util.EnumMap, HashMap, TreeMap, LinkedHashMap or Hashtable) to gain access to their larger functionality if not performance. Matlab versions earlier than R2008b have no real alternative in any case and must use the Java classes.

A potential limitation of using Java Collections is their inability to contain non-primitive Matlab types such as structs. To overcome this, either down-convert the types (e.g., using struct2cell or programmatically), or create a separate Java object that will hold your information and store this object in the Java Collection.

You may also be interested to examine a pure-Matlab object-oriented (class-based) Hashtable implementation, which is available on the File Exchange.

Yair Altman
Another Matlab class-based implementation posted today: http://www.mathworks.com/matlabcentral/fileexchange/28586
Yair Altman

related questions