Как разделить один график smoothScatter на другой?

Question

Я хочу знать, насколько часто встречаются определенные географические названия.Из национального справочника я сделал два smoothScatter() в R, один со всеми местами, другой с местами, названия которых мне интересны.

Все места:

Места с определенными названиями:

Теперь, как я могу разделить второе на первое,получить плотность имен, представляющих интерес по всем именам?Это может быть решение R или ImageMagick, GIMP ...

Answer 1

В Imagemagick 6 вы можете разделить два изображения следующим образом:

convert second.png first.png +swap -compose divide -composite division.png

Если вы используете Imagemagick 7, используйтемагия, а не конвертация.

Imagemagick 7 по умолчанию компилируется в режиме HDRI;тогда как Imagemagick 6 не компилируется в режиме HDRI по умолчанию.

Это означает, что результаты могут быть вне нормального 16-битного диапазона.Таким образом, деление может привести к значениям больше 1. (См. Комментарий ниже о нормализованных значениях, используемых в Imagemagick).Однако, как правило, они фиксируются на черно-белых крайностях, чтобы иметь возможность сохранять изображения с обычными типами изображений, такими как PNG.

Но можно отключить ограничение и получить результат в таком формате, как PFM, который может сохранять значения вне нормального диапазона.Или мы можем разделить на достаточно, чтобы диапазон все еще находился в пределах диапазона, соответствующего PNG, или можно растянуть результат до полного диапазона черного и белого.Так, например:

magick second.png first.png +swap -define compose:clamp=off -compose divide -composite -evaluate divide 10 division2.png

magick second.png first.png +swap -define compose:clamp=off -compose divide -composite -auto-level division3.png

Это дает больше деталей, чем первый подход выше.