У меня есть ковариативные данные о количестве осадков. У меня также есть точечный образец, который представляет населенные пункты на карте в определенной области. У меня проблема в том, что в части моей области моего региона есть море, однако на море нет населенных пунктов, но данные о дождях действительно относятся к морской области. Я хочу установить значения дождя для ближайшей точки на суше.
есть идеи как это сделать? Возможно, чтобы найти, какие координаты осадков x и y представляют море, я хочу создать фрейм данных, содержащий центральные точки всех квадратных километров, а затем импортировать значения дождя для местоположений, которые совпадают, и посмотреть, какие местоположения не имеют значений. Но я застрял. Вот мой код:
> window<- data.frame(Lon=c(-1.560367, -1.078330 ), Lat=c( 50.576342, 51.243823))
> coordinates(window) <- ~Lon + Lat
>
>
> proj4string(window) = CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")
>
> proj4string(winch2) <- latlong
>
> window <- spTransform(window, bng)
>
> (floor (coordinates (window) / 1000) + 0.5) * 1000
>
> W2<- owin(c(431500,464500), c( 75500, 149500))
>
> Region<-Settlements[W2] ###Settlements is my data
>
> rain_im[W2] ###rain is my covariate as a pixel image
>
> as.data.frame(rain_im[W2]) ###Converted this into a pixel image
>
> `summary(Region)`
Marked planar point pattern: 308 points
Average intensity 1.261261e-07 points per square unit
Coordinates are given to 2 decimal places
i.e. rounded to the nearest multiple of 0.01 units
marks are numeric, of type ‘double’
Summary:
Min. 1st Qu. Median Mean 3rd Qu. Max.
2.375 19.000 47.500 103.029 88.364 5500.000
Window: rectangle = [431500, 464500] x [75500, 149500] units
Window area = 2.442e+09 square units
Пожалуйста, опишите ваши данные более четко. Вы говорите, что у вас есть «ковариативные данные о количестве осадков»: в каком формате эти данные? Является ли это пиксельным изображением (значения пикселей которого представляют собой количество дождя) или это фрейм данных с местоположениями датчиков дождя и записанными количествами дождя? У вас есть некоторые значения ковариаты дождя в пикселях, лежащих над морем: что вы хотите делать с такими значениями? Какие у вас есть данные, которые говорят вам, какие места находятся в море, а какие на суше?
Если бы у вас были полигональные данные, определяющие береговую линию, вы могли бы преобразовать их owin
, скажем, в оконный объект Land
. Затем, если ваши данные о дождевых осадках представляют собой пиксельное изображение, Rain
вы можете просто ClipRain <- Rain[Land, drop=FALSE]
ограничить пиксельное изображение областью суши (присвоить NA
любому пикселю над морем).
Если у вас есть данные об осадках в некоторых местах, скажем, точечный узор RainRecords
с метками, которые являются значениями осадков, то Smooth(RainRecords, sigma)
это пиксельное изображение, содержащее сглаженную ядром интерполяцию осадков, и nnmark(RainRecords)
пиксельное изображение, в котором каждое значение пикселя представляет собой количество осадков в ближайшем записанное местоположение. Обрезав или отобрав эти изображения, вы, вероятно, сможете получить то, что вам нужно.
Эта статья взята из Интернета, укажите источник при перепечатке.
Если есть какие-либо нарушения, пожалуйста, свяжитесь с[email protected] Удалить.
я говорю два предложения