Можете ли вы объяснить следующий код R?

Question

Можете ли вы объяснить этот код?

Какая польза от х1 и х2

Я хочу создать контур из этого кода.

объясните, пожалуйста, использование interp ().

как этот код может быть полезен при построении контура?

у нас есть текстовый файл, в котором у нас есть три отдельных столбца с длиной long. как использовать interp ()?

f=read.table("data2_input.txt",sep="\t",header=T)
attach(f)
library(akima)
library(reshape2)
c1=length(f[1,])
r1=length(f[,1])
lat=0,long=0,conc=0,x1=0 y1=0 midx=0 x2=0 y2=0 midy=0 
conc=0 r=r1                            
c=c1
lat=f[,1]
long=f[,2]
conc=f[,3]
ak=interp(long,lat,conc, xo=seq(min(long),max(long),length=200),yo=seq(min(lat),max(lat),length=75))
mak=melt(ak$z)
names(mak)=c('x','y','value')
mak$lon=ak$x[mak$x]
mak$lat=ak$y[mak$y]
lon_diff=round(ak$x[2]-ak$x[1],5)
lat_diff=round(ak$y[2]-ak$y[1],5)
mak1=mak[,3:5]
mak1=subset(mak1,mak1$value!='NA')
names(mak1)=c('Concentration','Longitude','Latitude')

/* why use x1 and x2*/
x1=round(mak1$Latitude-lat_diff/2,4)
y1=round(mak1$Longitude-lon_diff/2,4)
x2=round(mak1$Latitude+lat_diff/2,4)
y2=round(mak1$Longitude+lon_diff/2,4)
/*why use round function*/          
z1=as.data.frame(cbind(x1,y1,x2,y2,round(mak1$Latitude,4),round(mak1$Longitude,4),round(mak1$Concentration,4)))
names(z1)=c('latitude_x1','longitude_y1','latitude_x2','longitude_y2','midx','midy','dvalue')
write.csv(z1,'data2_output.csv',row.names=F,quote=F)

Answer 1

Добро пожаловать в Stack Overflow и спасибо за ваш вопрос (однако опубликованные комментарии по-прежнему актуальны:)

• x1, y1 и x2, y2 - векторы координат нижних правых и верхних левых вершин прямоугольников, которые вы собираетесь создать;
• ваши пространственные данные о концентрации распределены неравномерно на поверхности, поэтому, если вы смогли нарисовать их, вы увидите много пустых мест, сбоев и неровностей на карте. Ваш контур и основное изображение будут нечитаемыми и искаженными.
• Чтобы избежать этого, вы должны интерполировать данные (или поместить данные в обычную сетку), поэтому вот почему вам нужно использовать функцию interp пакета akima.
• крайний пример (однако для ggplot2 пакета) отсутствия интерполяции / сглаживания, который вы можете увидеть на Контурная диаграмма не заполнена полностью
• base графическое построение невозможно вызвать само по себе для нерегулярных данных, поскольку в качестве аргумента должна быть предоставлена ​​матрица регулярно распределенных z-значений.