Создание спектральных тепловых карт или карт интенсивности из данных CDIP с использованием Ruby

Я очень спешу, поэтому не могу закончить прямо сейчас, но, поскольку никто еще не ответил, вот первый подход:

Код в системе Mathematica (извините, как я уже сказал, сейчас нет времени):

a = Import["http://cdip.ucsd.edu/data_access/MEM_2dspectra.cdip?100", 
   "Table"];

Graphics@Flatten[Table[

    (*colors, dont mind*)
    {ColorData["CMYKColors"][(a[[r, t]] - .000007)/(.0003 - 0.000007)], 

    (*point size, dont mind*)
    PointSize[1/Sqrt[r]/10], 

    (*Coordinates for your points "a" is your data matrix *)
       Point[
            {(rr =Log[.025 + (.58 - .25)/64 r]) Cos@(tt = t 5 Degree), 
              rr Sin@tt}]
            }

     (*values for the iteration*)
     , {r, 7, 64}, {t, 1, 72}], 1] 

     (*Rotation, dont mind*)
     /. gg : Graphics[___] :> Rotate[gg, Pi/2]  

Я до сих пор не могу подобрать нужную цветовую гамму:

Чип, я знаю, не использует Ruby, но если предположить, что расчет точек Велисария в порядке, я бы использовал _ 1_ вместо этого, так как это намного проще для понимания и дает более ясную картину.

(* Read in data, the web address can be specified *)
a = Import[<url of cpid data>, "Table"];

Import остается в теге pre в первом подсписке и как подсписок с одним элементом в конце, это удаляет его

dat = a[[ ;; -2]][[All, -72;; ]];

сначала взяв все, кроме последнего элемента, а затем взяв последние 72 элемента из всех оставшихся подсписок.

ListContourPlot ожидает список точек в форме {{x, y, f}, ...}, поэтому нам нужно преобразовать dat в эту форму, как описано в в другом месте:

pts =Flatten[
    Table[ {
          (rr =Log[.025 + (.58 - .25)/64 r]) Cos@(tt = t 5 Degree), (* x-coord *)
          rr Sin@tt, (* y-coord *)
          dat[[r,t]] (* function value *)
           },
           {r, 7, 64}, {t, 1, 72}
         ],
    1 ]

а затем построить их

ListContourPlot[pts,
  ColorFunctionScaling -> False,
  ColorFunction -> (ColorData["CMYKColors"][(# - .000007)/(.0003 - 0.000007)]&)
  ]

Это дает:

что можно исправить, исправив обрезку и ColorFunction масштабирование. Также, поработав, можно добавить радиальные контуры.