ruby-gr 0.0.23 → 0.58.1.0

data/lib/gr/ffi.rb

@@ -15,7 +15,7 @@ module GR
       raise LoadError, 'Could not find GR Framework'
     end
 
-    extend GRCommons::Extern
+    extend GRCommons::TryExtern
 
     # https://github.com/sciapp/gr/blob/master/lib/gr/gr.h
     # keep same order
@@ -41,6 +41,8 @@ module GR
                'int, int, int, int, int, int, int *)'
     try_extern 'void gr_polarcellarray(double, double, double, double, double, double, ' \
                'int, int, int, int, int, int, int *)'
+    try_extern 'void gr_nonuniformpolarcellarray(double, double, double *, double *, ' \
+               'int, int, int, int, int, int, int *);'
     try_extern 'void gr_gdp(int, double *, double *, int, int, int *)'
     try_extern 'void gr_spline(int, double *, double *, int, int)'
     try_extern 'void gr_gridit(int, double *, double *, double *, int, int, ' \
@@ -115,6 +117,7 @@ module GR
     try_extern 'void gr_setcolormap(int)'
     try_extern 'void gr_inqcolormap(int *)'
     try_extern 'void gr_setcolormapfromrgb(int n, double *r, double *g, double *b, double *x)'
+    try_extern 'void gr_inqcolormapinds(int *, int *)'
     try_extern 'void gr_colorbar(void)'
     try_extern 'void gr_inqcolor(int, int *)'
     try_extern 'int gr_inqcolorfromrgb(double, double, double)'
@@ -174,7 +177,7 @@ module GR
     try_extern 'void gr_shadepoints(int, double *, double *, int, int, int)'
     try_extern 'void gr_shadelines(int, double *, double *, int, int, int)'
     try_extern 'void gr_panzoom(double, double, double, double, double *, double *, double *, double *)'
-    # try_extern 'int gr_findboundary(int, double *, double *, double, double (*)(double, double), int, int *)'
+    try_extern 'int gr_findboundary(int, double *, double *, double, double (*)(double, double), int, int *)'
     try_extern 'void gr_setresamplemethod(unsigned int flag)'
     try_extern 'void gr_inqresamplemethod(unsigned int *flag)'
     try_extern 'void gr_path(int, double *, double *, const char *)'
@@ -182,6 +185,8 @@ module GR
     try_extern 'void gr_inqborderwidth(double *)'
     try_extern 'void gr_setbordercolorind(int)'
     try_extern 'void gr_inqbordercolorind(int *)'
+    try_extern 'void gr_selectclipxform(int)'
+    try_extern 'void gr_inqclipxform(int *);'
     try_extern 'void gr_setprojectiontype(int)'
     try_extern 'void gr_inqprojectiontype(int *)'
     try_extern 'void gr_setperspectiveprojection(double, double, double)'
@@ -202,5 +207,7 @@ module GR
     try_extern 'void gr_inqtext3d(double, double, double, char *, int axis, double *, double *)'
     try_extern 'void gr_settextencoding(int)'
     try_extern 'void gr_inqtextencoding(int *)'
+    try_extern 'void gr_loadfont(char *, int *)'
+    # gr_setcallback(char *(*)(const char *));
   end
 end

data/lib/gr/plot.rb

@@ -33,17 +33,53 @@ module GR
     # Plot kinds conform to GR.jl
     PLOT_KIND = %i[line step scatter stem hist contour contourf hexbin heatmap
                    nonuniformheatmap wireframe surface plot3 scatter3 imshow
-                   isosurface polar polarhist polarheatmap trisurf tricont shade
-                   volume].freeze # the name might be changed in the future.
+                   isosurface polar polarhist polarheatmap nonuniformpolarheatmap
+                   trisurf tricont shade volume].freeze
 
     # Keyword options conform to GR.jl.
     KW_ARGS = %i[accelerate algorithm alpha ax backgroundcolor barwidth baseline
-                 clabels clear clim color colormap crange figsize grid horizontal
-                 isovalue kind label labels levels location nbins ratio rotation
-                 scale size spec subplot tilt title update xaxis xflip xform
-                 xlabel xlim xlog xrange xticks yaxis yflip ylabel ylim ylog
-                 zflip yrange yticks viewport vp where window zaxis zlabel zlim
-                 zlog zrange zticks].freeze
+                 clabels clear clim color colormap crange figsize font grid
+                 horizontal isovalue kind label labels levels linewidth location
+                 nbins ratio rotation scale size spec subplot tilt title update
+                 xaxis xflip xform xlabel xlim xlog xrange xticks yaxis yflip
+                 ylabel ylim ylog zflip yrange yticks viewport vp where window
+                 zaxis zlabel zlim zlog zrange zticks].freeze
+
+    FONTS = {
+      times_roman: 101,
+      times_italic: 102,
+      times_bold: 103,
+      times_bolditalic: 104,
+      helvetica_regular: 105,
+      helvetica_oblique: 106,
+      helvetica_bold: 107,
+      helvetica_boldoblique: 108,
+      courier_regular: 109,
+      courier_oblique: 110,
+      courier_bold: 111,
+      courier_boldoblique: 112,
+      symbol: 113,
+      bookman_light: 114,
+      bookman_lightitalic: 115,
+      bookman_demi: 116,
+      bookman_demiitalic: 117,
+      newcenturyschlbk_roman: 118,
+      newcenturyschlbk_italic: 119,
+      newcenturyschlbk_bold: 120,
+      newcenturyschlbk_bolditalic: 121,
+      avantgarde_book: 122,
+      avantgarde_bookoblique: 123,
+      avantgarde_demi: 124,
+      avantgarde_demioblique: 125,
+      palatino_roman: 126,
+      palatino_italic: 127,
+      palatino_bold: 128,
+      palatino_bolditalic: 129,
+      zapfchancery_mediumitalic: 130,
+      zapfdingbats: 131,
+      cmuserif_math: 232, # original: cmuserif-math
+      dejavusans: 233
+    }.freeze
 
     @last_plot = nil
     class << self
@@ -116,15 +152,20 @@ module GR
         vp1, vp2, vp3, vp4 = vp
       end
 
-      viewport = [vp1 + 0.125 * (vp2 - vp1),
-                  vp1 + 0.925 * (vp2 - vp1),
-                  vp3 + 0.125 * (vp4 - vp3),
-                  vp3 + 0.925 * (vp4 - vp3)]
-
-      if %i[contour contourf hexbin heatmap nonuniformheatmap polarheatmap
-            surface trisurf volume].include?(kind)
-        viewport[1] -= 0.1
-      end
+      left_margin = kvs.has_key?(:ylabel) ? 0.05 : 0
+      right_margin = if %i[contour contourf hexbin heatmap nonuniformheatmap polarheatmap
+                           nonuniformpolarheatmap surface trisurf volume].include?(kind)
+                       (vp2 - vp1) * 0.1
+                     else
+                       0
+                     end
+      bottom_margin = kvs.has_key?(:xlabel) ? 0.05 : 0
+      top_margin = kvs.has_key?(:title) ? 0.075 : 0
+
+      viewport = [vp1 + (0.075 + left_margin) * (vp2 - vp1),
+                  vp1 + (0.95 - right_margin) * (vp2 - vp1),
+                  vp3 + (0.075 + bottom_margin) * (vp4 - vp3),
+                  vp3 + (0.975 - top_margin) * (vp4 - vp3)]
 
       if %i[line step scatter stem].include?(kind) && kvs[:labels]
         location = kvs[:location] || 1
@@ -156,7 +197,7 @@ module GR
         GR.restorestate
       end
 
-      if %i[polar polarhist polarheatmap].include? kind
+      if %i[polar polarhist polarheatmap nonuniformpolarheatmap].include? kind
         xmin, xmax, ymin, ymax = viewport
         xcenter = 0.5 * (xmin + xmax)
         ycenter = 0.5 * (ymin + ymax)
@@ -167,7 +208,7 @@ module GR
 
     def set_window(kind)
       scale = 0
-      unless %i[polar polarhist polarheatmap].include?(kind)
+      unless %i[polar polarhist polarheatmap nonuniformpolarheatmap].include?(kind)
         scale |= GR::OPTION_X_LOG  if kvs[:xlog]
         scale |= GR::OPTION_Y_LOG  if kvs[:ylog]
         scale |= GR::OPTION_Z_LOG  if kvs[:zlog]
@@ -182,11 +223,11 @@ module GR
         kvs[:xrange] = [xmin, xmax]
         kvs[:yrange] = [ymin, ymax]
       else
-        minmax
+        minmax(kind)
       end
 
       major_count = if %i[wireframe surface plot3 scatter3 polar polarhist
-                          polarheatmap trisurf volume].include?(kind)
+                          polarheatmap nonuniformpolarheatmap trisurf volume].include?(kind)
                       2
                     else
                       5
@@ -202,15 +243,15 @@ module GR
         xmax += 0.5
       end
       xtick, majorx = if (scale & GR::OPTION_X_LOG) == 0
-                        unless %i[heatmap polarheatmap].include?(kind)
-                          unless kvs.has_key?(:xlim)
-                            xmin, xmax = GR.adjustlimits(xmin, xmax) unless kvs[:panzoom]
-                          end
+                        if !%i[heatmap polarheatmap].include?(kind) &&
+                           !kvs.has_key?(:xlim) &&
+                           !kvs[:panzoom]
+                          xmin, xmax = GR.adjustlimits(xmin, xmax)
                         end
                         if kvs.has_key?(:xticks)
                           kvs[:xticks]
                         else
-                          [GR.tick(xmin, xmax) / major_count, major_count]
+                          [auto_tick(xmin, xmax) / major_count, major_count]
                         end
                       else
                         [1, 1]
@@ -231,15 +272,15 @@ module GR
         end
       end
       ytick, majory = if (scale & GR::OPTION_Y_LOG) == 0
-                        unless %i[heatmap polarheatmap].include?(kind)
-                          unless kvs.has_key?(:ylim)
-                            ymin, ymax = GR.adjustlimits(ymin, ymax) unless kvs[:panzoom]
-                          end
+                        if !%i[heatmap polarheatmap].include?(kind) &&
+                           !kvs.has_key?(:ylim) &&
+                           !kvs[:panzoom]
+                          ymin, ymax = GR.adjustlimits(ymin, ymax)
                         end
                         if kvs.has_key?(:yticks)
                           kvs[:yticks]
                         else
-                          [GR.tick(ymin, ymax) / major_count, major_count]
+                          [auto_tick(ymin, ymax) / major_count, major_count]
                         end
                       else
                         [1, 1]
@@ -254,7 +295,7 @@ module GR
                           if kvs.has_key?(:zticks)
                             kvs[:zticks]
                           else
-                            [GR.tick(zmin, zmax) / major_count, major_count]
+                            [auto_tick(zmin, zmax) / major_count, major_count]
                           end
                         else
                           [1, 1]
@@ -264,15 +305,16 @@ module GR
       end
 
       kvs[:window] = xmin, xmax, ymin, ymax
-      if %i[polar polarhist polarheatmap].include?(kind)
+      if %i[polar polarhist polarheatmap nonuniformpolarheatmap trisurf].include?(kind)
         GR.setwindow(-1, 1, -1, 1)
       else
         GR.setwindow(xmin, xmax, ymin, ymax)
       end
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
         rotation = kvs[:rotation] || 40
-        tilt     = kvs[:tilt]     || 70
-        GR.setspace(zmin, zmax, rotation, tilt)
+        tilt     = kvs[:tilt]     || 60
+        GR.setwindow3d(xmin, xmax, ymin, ymax, zmin, zmax)
+        GR.setspace3d(-rotation, tilt, 30, 0)
       end
 
       kvs[:scale] = scale
@@ -291,10 +333,10 @@ module GR
       GR.setlinecolorind(1)
       diag = Math.sqrt((viewport[1] - viewport[0])**2 + (viewport[3] - viewport[2])**2)
       GR.setlinewidth(1)
-      charheight = [0.018 * diag, 0.012].max
-      GR.setcharheight(charheight)
       ticksize = 0.0075 * diag
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
+        charheight = [0.024 * diag, 0.012].max
+        GR.setcharheight(charheight)
         ztick, zorg, majorz = kvs[:zaxis]
         if pass == 1 && drawgrid
           GR.grid3d(xtick, 0, ztick, xorg[0], yorg[1], zorg[0], 2, 0, 2)
@@ -304,6 +346,8 @@ module GR
           GR.axes3d(0, ytick, 0, xorg[1], yorg[0], zorg[0], 0, majory, 0, ticksize)
         end
       else
+        charheight = [0.018 * diag, 0.012].max
+        GR.setcharheight(charheight)
         if %i[heatmap nonuniformheatmap shade].include?(kind)
           ticksize = -ticksize
         elsif drawgrid
@@ -345,29 +389,29 @@ module GR
         GR.axes(xtick, ytick, xorg[1], yorg[1], -majorx, -majory, -ticksize)
       end
 
-      if kvs[:title]
+      if kvs.has_key?(:title)
         GR.savestate
         GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_TOP)
-        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title])
+        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title].to_s)
         GR.restorestate
       end
       if %i[wireframe surface plot3 scatter3 trisurf volume].include?(kind)
-        xlabel = kvs[:xlabel] || ''
-        ylabel = kvs[:ylabel] || ''
-        zlabel = kvs[:zlabel] || ''
+        xlabel = (kvs[:xlabel] || '').to_s
+        ylabel = (kvs[:ylabel] || '').to_s
+        zlabel = (kvs[:zlabel] || '').to_s
         GR.titles3d(xlabel, ylabel, zlabel)
       else
         if kvs.has_key?(:xlabel)
           GR.savestate
           GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_BOTTOM)
-          text(0.5 * (viewport[0] + viewport[1]), vp[2] + 0.5 * charheight, kvs[:xlabel])
+          text(0.5 * (viewport[0] + viewport[1]), vp[2] + 0.5 * charheight, kvs[:xlabel].to_s)
           GR.restorestate
         end
         if kvs.has_key?(:ylabel)
           GR.savestate
           GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_TOP)
           GR.setcharup(-1, 0)
-          text(vp[0] + 0.5 * charheight, 0.5 * (viewport[2] + viewport[3]), kvs[:ylabel])
+          text(vp[0] + 0.5 * charheight, 0.5 * (viewport[2] + viewport[3]), kvs[:ylabel].to_s)
           GR.restorestate
         end
       end
@@ -386,7 +430,7 @@ module GR
       GR.setcharheight(charheight)
       GR.setlinetype(GR::LINETYPE_SOLID)
 
-      tick = 0.5 * GR.tick(rmin, rmax)
+      tick = auto_tick(rmin, rmax)
       n = ((rmax - rmin) / tick + 0.5).round
       (n + 1).times do |i|
         r = i.to_f / n
@@ -475,7 +519,7 @@ module GR
       if kvs.has_key?(:title)
         GR.savestate
         GR.settextalign(GR::TEXT_HALIGN_CENTER, GR::TEXT_VALIGN_TOP)
-        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title])
+        text(0.5 * (viewport[0] + viewport[1]), vp[3], kvs[:title].to_s)
         GR.restorestate
       end
       GR.selntran(1)
@@ -548,7 +592,7 @@ module GR
       charheight = [0.016 * diag, 0.012].max
       GR.setcharheight(charheight)
       if kvs[:scale] & GR::OPTION_Z_LOG == 0
-        ztick = 0.5 * GR.tick(zmin, zmax)
+        ztick = auto_tick(zmin, zmax)
         GR.axes(0, ztick, 1, zmin, 0, 1, 0.005)
       else
         GR.setscale(GR::OPTION_Y_LOG)
@@ -577,16 +621,34 @@ module GR
         # Not yet.
       end
 
-      # The following fonts are the default in GR.jl
-      # Japanese, Chinese, Korean, etc. cannot be displayed.
-      # GR.settextfontprec(232, 3) # CM Serif Roman
+      if kvs.has_key?(:font)
+        name = kvs[:font]
+        # 'Cmuserif-Math' => :cmuserif_math
+        sym_name = name.to_s.gsub('-', '_').downcase.to_sym
+        if FONTS.include?(sym_name)
+          font = FONTS[sym_name]
+          GR.settextfontprec(font, font > 200 ? 3 : 0)
+        else
+          font = GR.loadfont(name)
+          if font >= 0
+            GR.settextfontprec(font, 3)
+          else
+            warn "Unknown font name: #{name}"
+          end
+        end
+      else
+        # The following fonts are the default in GR.jl
+        # Japanese, Chinese, Korean, etc. cannot be displayed.
+
+        # GR.settextfontprec(232, 3) # CM Serif Roman
+      end
 
       set_viewport(kind, kvs[:subplot])
       unless kvs[:ax]
         set_window(kind)
         if %i[polar polarhist].include?(kind)
           draw_polar_axes
-        elsif !%i[imshow isosurface polarheatmap].include?(kind)
+        elsif !%i[imshow isosurface polarheatmap nonuniformpolarheatmap].include?(kind)
           draw_axes(kind)
         end
       end
@@ -607,8 +669,12 @@ module GR
 
         when :line
           mask = GR.uselinespec(spec)
-          GR.polyline(x, y) if hasline(mask)
-          GR.polymarker(x, y) if hasmarker(mask)
+          linewidth = kvs[:linewidth]
+          # Slightly different from Julia,
+          # Because the implementation of polyline and polymarker is different.
+          z ||= linewidth # FIXME
+          GR.polyline(x, y, z, c) if hasline(mask)
+          GR.polymarker(x, y, z, c) if hasmarker(mask)
 
         when :step
           mask = GR.uselinespec(spec)
@@ -644,20 +710,9 @@ module GR
 
         when :scatter
           GR.setmarkertype(GR::MARKERTYPE_SOLID_CIRCLE)
-          if z || c
-            if c
-              cmin, cmax = kvs[:crange]
-              c = c.map { |i| normalize_color(i, cmin, cmax) }
-              cind = c.map { |i| (1000 + i * 255).round }
-            end
-            x.length.times do |i|
-              GR.setmarkersize(z[i] / 100.0) if z
-              GR.setmarkercolorind(cind[i]) if c
-              GR.polymarker([x[i]], [y[i]])
-            end
-          else
-            GR.polymarker(x, y)
-          end
+          z = z&.map { |i| i * 0.01 }
+          c = c&.map { |i| normalize_color(i, *kvs[:crange]) }
+          GR.polymarker(x, y, z, c)
 
         when :stem
           GR.setlinecolorind(1)
@@ -707,7 +762,7 @@ module GR
             GR.fillarc(-ρ[i], ρ[i], -ρ[i], ρ[i], θ[i - 1], θ[i])
           end
 
-        when :polarheatmap
+        when :polarheatmap, :nonuniformpolarheatmap
           w, h = z.shape
           cmap = colormap
           cmin, cmax = kvs[:zrange]
@@ -716,7 +771,15 @@ module GR
           data.reverse(axis: 1) if kvs[:yflip]
           colors = data * 255 + 1000
           colors = colors.transpose # Julia is column major
-          GR.polarcellarray(0, 0, 0, 360, 0, 1, w, h, colors)
+          case kind
+          when :polarheatmap
+            GR.polarcellarray(0, 0, 0, 360, 0, 1, w, h, colors)
+          when :nonuniformpolarheatmap
+            ymax = y.max.to_f
+            ρ = y.map { |i| i / ymax }
+            θ = x.map { |i| i * 180 / Math::PI }
+            GR.nonuniformpolarcellarray(θ, ρ, w, h, colors)
+          end
           draw_polar_axes
           kvs[:zrange] = [cmin, cmax]
           colorbar
@@ -738,7 +801,7 @@ module GR
             zmin = kvs[:zlim].first if kvs[:zlim].first
             zmax = kvs[:zlim].last if kvs[:zlim].last
           end
-
+          GR.setprojectiontype(0)
           GR.setspace(zmin, zmax, 0, 90)
           levels = kvs[:levels] || 0
           clabels = kvs[:clabels] || false
@@ -748,9 +811,10 @@ module GR
           else
             h = levels
           end
-          if kind == :contour
+          case kind
+          when :contour
             GR._contour_(x, y, h, z, clabels ? 1 : 1000)
-          elsif kind == :contourf
+          when :contourf
             GR._contourf_(x, y, h, z, clabels ? 1 : 0)
           end
           colorbar(0, h.length)
@@ -972,7 +1036,7 @@ module GR
 
     def to_svg
       ## Need IRuby improvemend.
-      GR.show(false) if ENV['GKSwstype'] == 'svg'
+      GR.show(false) if ENV['GKS_WSTYPE'] == 'svg'
     end
 
     private
@@ -1045,8 +1109,8 @@ module GR
 
     # Normalize a color c with the range [cmin, cmax]
     #   0 <= normalize_color(c, cmin, cmax) <= 1
-    # Note: narray.map{|i| normalize_color(i)} There's room for speedup.
     def normalize_color(c, cmin, cmax)
+      # NOTE: narray.map{|i| normalize_color(i)} There's room for speedup.
       c = c.to_f # if c is Integer
       c = c.clamp(cmin, cmax) - cmin
       c /= (cmax - cmin) if cmin != cmax
@@ -1054,6 +1118,7 @@ module GR
     end
 
     def inqtext(x, y, s)
+      s = s.to_s
       if s.length >= 2 && s[0] == '$' && s[-1] == '$'
         GR.inqmathtex(x, y, s[1..-2])
       elsif s.include?('\\') || s.include?('_') || s.include?('^')
@@ -1064,6 +1129,7 @@ module GR
     end
 
     def text(x, y, s)
+      s = s.to_s
       if s.length >= 2 && s[0] == '$' && s[-1] == '$'
         GR.mathtex(x, y, s[1..-2])
       elsif s.include?('\\') || s.include?('_') || s.include?('^')
@@ -1081,7 +1147,7 @@ module GR
       [a, b]
     end
 
-    def minmax
+    def minmax(kind)
       xmin = ymin = zmin = cmin = Float::INFINITY
       xmax = ymax = zmax = cmax = -Float::INFINITY
       scale = kvs[:scale]
@@ -1094,6 +1160,8 @@ module GR
           x0, x1 = x.minmax
           xmin = [x0, xmin].min
           xmax = [x1, xmax].max
+        elsif kind == :volume
+          xmin, xmax = -1, 1
         else
           xmin = 0
           xmax = 1
@@ -1105,6 +1173,8 @@ module GR
           y0, y1 = y.minmax
           ymin = [y0, ymin].min
           ymax = [y1, ymax].max
+        elsif kind == :volume
+          ymin, ymax = -1, 1
         else
           ymin = 0
           ymax = 1
@@ -1146,6 +1216,21 @@ module GR
       end
     end
 
+    def to_wc(wn)
+      xmin, ymin = GR.ndctowc(wn[0], wn[2])
+      xmax, ymax = GR.ndctowc(wn[1], wn[3])
+      [xmin, xmax, ymin, ymax]
+    end
+
+    def auto_tick(amin, amax)
+      scale = 10.0**Math.log10(amax - amin).truncate
+      tick_size = [5.0, 2.0, 1.0, 0.5, 0.2, 0.1, 0.05, 0.02, 0.01]
+      i = tick_size.find_index do |tsize|
+        ((amax - amin) / scale / tsize) > 7 # maximum number of tick marks
+      end
+      tick = tick_size[i - 1] * scale
+    end
+
     def legend_size
       scale = GR.inqscale
       GR.selntran(0)
@@ -1163,20 +1248,12 @@ module GR
       [w, h]
     end
 
-    # NOTE: duplicated definition (GRCommonUtils)
     def equal_length(*args)
-      lengths = args.map(&:length)
-      unless lengths.all? { |l| l == lengths[0] }
-        raise ArgumentError,
-              'Sequences must have same length.'
-      end
-
-      lengths[0]
+      GRCommons::GRCommonUtils.equal_length(*args)
     end
 
-    # NOTE: duplicated definition (GRCommonUtils)
     def narray?(data)
-      defined?(Numo::NArray) && data.is_a?(Numo::NArray)
+      GRCommons::GRCommonUtils.narray?(data)
     end
   end
 
@@ -1240,6 +1317,20 @@ module GR
       end
     end
 
+    # (Plot) Draw a nonuniformpolarheatmap.
+    def nonuniformpolarheatmap(*args)
+      # FIXME
+      args, kv = format_xyzc(*args)
+      _x, _y, z = args
+      ysize, xsize = z.shape
+      z = z.reshape(xsize, ysize)
+      create_plot(:nonuniformpolarheatmap, kv) do |plt|
+        plt.kvs[:xlim] ||= [0.5, xsize + 0.5]
+        plt.kvs[:ylim] ||= [0.5, ysize + 0.5]
+        plt.args = [[(1..xsize).to_a, (1..ysize).to_a, z, nil, '']]
+      end
+    end
+
     alias _contour_ contour
     # (Plot) Draw a contour plot.
     def contour(*args)