2013-11-01 4 views
0

У меня есть данные по ряду экологических переменных, связанных с пространственными точками. Каждая точка имеет x & y координаты относительно ограничительной рамки, однако точки представляют собой круговые области различного диаметра. То, что я пытаюсь достичь, - проецировать площадь, занимаемую каждой точкой, на окно наблюдения, чтобы мы могли впоследствии пикселизовать область и получить степень перекрытия области каждой точки с каждым пикселем (ячейкой сетки). Раньше я мог достичь этого с помощью данных транца, преобразовывая в объект линии psp, а затем используя функцию pixellate в пакете spartstat, но я не уверен, как продолжить работу с этими круговыми областями. Мне кажется, что я должен использовать классы полигонов, но снова я не уверен, как их определять. Любое предложение было бы весьма признательным.Использование координат точки и диаметра для расчета площади проекции точек

+0

Это похоже на более общую проблему, не относящуюся к R. Вы можете построить сетку и для каждой точки сетки рассчитать, находится ли точка в пределах любого круга (меньше радиального расстояния от точек). В качестве альтернативы вы также можете попробовать моделирование по сетке. –

+1

Я бы предложил использовать пространственные типы данных, реализованные в пакетах ** sp ** и ** растр **. Затем вы можете использовать (что-то типа) 'rgeos :: gBuffer()' для создания круговых 'SpatialPolygons' нужных радиусов, окружающих каждую точку (идеально сохраненных как' SpatialPolygonsDataFrame' с прикрепленными данными атрибутов), 'raster :: rasterize () 'для преобразования SPDF в растровый файл, зарегистрированный на вашем существующем растере, и' raster :: crosstab() 'или' raster :: zonal() 'для извлечения информации о областях перекрытия. (Если бы вы включили минимальный воспроизводимый пример, я мог бы показать вам, как это сделать;) –

+0

Спасибо, R - это то, что я обычно использую, поэтому я идеально ищут решение r. Я подумал о нескольких способах сделать это «вручную», но я надеялся избежать этого, используя пространственную функциональность многих пакетов r. Просто не могу понять первый шаг. –

ответ

0

В spatstat пакете, функция discs будет принимать местоположения (х, у) и радиусы R (или диаметры, области и т.д.) и генерировать либо многоугольные или пиксел маски представления окружностей, и возвращать их либо как отдельные объекты или как единый комбинированный объект.

Смежные вопросы