Sample Data: France Bastides

Source Notebook

Locations of bastides (walled cities) in France annotated with name and foundation year

Details

Locations of medieval bastides (walled cities) in the observation region in the South-West France, annotated with name and foundation year.

Examples

Basic Examples (1) 

In[1]:=
ResourceData[\!\(\*
TagBox["\"\<Sample Data: France Bastides\>\"",
#& ,
BoxID -> "ResourceTag-Sample Data: France Bastides-Input",
AutoDelete->True]\), "Data"]
Out[1]=

Summary of the spatial point data:

In[2]:=
ResourceData[\!\(\*
TagBox["\"\<Sample Data: France Bastides\>\"",
#& ,
BoxID -> "ResourceTag-Sample Data: France Bastides-Input",
AutoDelete->True]\), "Data"]["Summary"]
Out[2]=

Visualizations (2) 

Plot the spatial point data:

In[3]:=
GeoListPlot[ResourceData[\!\(\*
TagBox["\"\<Sample Data: France Bastides\>\"",
#& ,
BoxID -> "ResourceTag-Sample Data: France Bastides-Input",
AutoDelete->True]\), "Data"]]
Out[3]=

Visualize the varying intensity:

In[4]:=
GeoSmoothHistogram[ResourceData[\!\(\*
TagBox["\"\<Sample Data: France Bastides\>\"",
#& ,
BoxID -> "ResourceTag-Sample Data: France Bastides-Input",
AutoDelete->True]\), "Data"]["Points"]]
Out[4]=

Analysis (4) 

Compute probability of finding a point within given radius of an existing point - NearestNeighborG is the CDF of the nearest neighbor distribution:

In[5]:=
nnG = NearestNeighborG[ResourceData[\!\(\*
TagBox["\"\<Sample Data: France Bastides\>\"",
#& ,
BoxID -> "ResourceTag-Sample Data: France Bastides-Input",
AutoDelete->True]\), "Data"]]
Out[5]=
In[6]:=
maxR = nnG["MaxRadius"]
Out[6]=
In[7]:=
res = Table[nnG[r], {r, range = Range[maxR/100, maxR, maxR/100]}];
ListPlot[Transpose[{range, res}], AxesLabel -> {"radius", "probability"}]
Out[8]=

NearestNeighborG as the CDF of nearest neighbor distribution can be used to compute the mean distance between a typical point and its nearest neighbor - the mean of a positive support distribution can be approximated via a Riemann sum of 1- CDF. To use Riemann approximation create the partition of the support interval from 0 to maxR into 100 parts and compute the value of the NearestNeighborG at the middle of each subinterval:

In[9]:=
step = maxR/100;
middles = Subdivide[step/2, maxR - step/2, 99];
values = nnG[middles];

Now compute the Riemann sum to find the mean distance between a typical point and its nearest neighbor:

In[10]:=
Total[(1 - values)*step]
Out[10]=

Test for complete spacial randomness:

In[11]:=
SpatialRandomnessTest[ResourceData[\!\(\*
TagBox["\"\<Sample Data: France Bastides\>\"",
#& ,
BoxID -> "ResourceTag-Sample Data: France Bastides-Input",
AutoDelete->True]\), "Data"], {"PValue", "TestConclusion"}] // Column
Out[11]=

Gosia Konwerska, "Sample Data: France Bastides" from the Wolfram Data Repository (2022)  

Data Resource History

Source Metadata

Publisher Information