Changeset 22291

Timestamp:

12/05/17 16:18:31 (7 years ago)

Author:

jbondzio

Message:

ADD:adding shaded option for use with googlemaps plotting

File:

• issm/trunk-jpl/src/m/plot/plot_googlemaps.m (modified) (1 diff)

issm/trunk-jpl/src/m/plot/plot_googlemaps.m

@@ -85 +85 @@
 image_rgb = ind2rgb(uint16((data_grid - data_min)*(length(colorm)/(data_max-data_min))),colorm);
 
+if exist(options,'shaded'),
+        a    = -45;
+        scut = 0.2;
+        c    = 1;
+        % computes lighting from elevation gradient
+        [fx,fy] = gradient(data_grid,x_m,y_m);
+        fxy = -fx*sind(a) - fy*cosd(a);
+        clear fx fy % free some memory...
+        fxy(isnan(fxy)) = 0;
+
+        % computes maximum absolute gradient (median-style), normalizes, saturates and duplicates in 3-D matrix
+        r = repmat(max(min(fxy/nmedian(abs(fxy),1 - scut/100),1),-1),[1,1,3]);
+
+        % applies contrast using exponent
+        rp = (1 - abs(r)).^c;
+        image_rgb = image_rgb.*rp;
+
+        % lighter for positive gradient
+        k = find(r > 0);
+        image_rgb(k) = image_rgb(k) + (1 - rp(k));
+end
+
 alpha=ones(size(data_grid));
 alpha(find(~data_nan))=transparency;