I am studying air pollution exposure estimation. Does anyone know the calculation principle of GeoDensityPlot? I don't know if it can replace the calculation of complex Kriging or Gaussian process regression?

data = Table[
   city -> AirTemperatureData[city], {city, 
    GeoNearest["City", 
     Entity["City", {"NewYork", "NewYork", "UnitedStates"}], {All, 
      Quantity[10, "Kilometers"]}]}];
Show[GeoListPlot[data[[All, 1]], ImageSize -> Large, 
            GeoBackground -> "StreetMap", 
            GeoRangePadding -> Quantity[2, "Kilometers"]],
           GeoDensityPlot[data, ColorFunction -> "TemperatureMap", 
            MaxPlotPoints -> Infinity, GeoBackground -> None, 
            PlotLegends -> Automatic]]