All examples below use data from two samples tablets (X202071 and X202470).
Get the number of observations:
In[]:=
Length[obsData]
Out[]=
236
The types of observations:
In[]:=
obsData[Counts/*ReverseSort,"Type"]
Out[]=
Weather
86
RelativePosition
65
ZodiacalPosition
18
MarketRates
9
FirstAppearenceEast
6
Other
6
HistoricalEvent
5
LastAppearenceWest
5
NotVisible
4
SunsetToMoonrise
4
Visible
3
MoonriseToSunrise
3
SunriseToMoonset
3
MoonriseToSunset
3
MoonsetToSunrise
2
LastAppearenceEast
2
RiverLevel
2
RiverLevelChange
2
FirstStationaryPoint
2
WindDirection
2
rows 1–20 of
24
An example observation:
In[]:=
Normal[obsData[22]]
Out[]=
TypeRelativePosition,ContentObject
Moon
,Reference
Spica
,Displacements{{Cubits4,FingersMissing[Unmentioned],TotalCubits4,IdealDegrees8,RealDegrees9.08,Below},Missing[Destroyed]},DateJulianDateYear-207,Month4,Day15,BabylonianDateYear{SE,104},MonthI,Day13,TimeBeginningOfTheNight,DateObject
Day: Thu 15 Apr -208(Juliancalendar)
,ProvenanceCreatorChristopher Wolfram,TabletIDX202071,LineNumber8,NotesMissing[]
Get the number of observations by month (for the two sample tablets):
In[]:=
obsData[Counts/*(BarChart[#,ChartLabelsAutomatic]&),"Date","BabylonianDate","Month"]
Out[]=
Use JPL ephemerides to calculate the real positions of synodic objects and references in observations of relative position:
In[]:=
Module[{obs,lines,paths},​​obs=Normal[obsData[Select[#Type==="RelativePosition"&],{Query["Content",{"Object","Reference"}],Query["Date","DateObject"]}]];​​lines={DiaryAstronomicalPosition[#1["Object"],#2],DiaryAstronomicalPosition[#1["Reference"],#2]}&@@@obs;​​paths=Arrow@*GeoPath/@Map[GeoPosition@*Reverse,DeleteMissing[lines,1,2],{2}];​​GeoGraphics[{Arrowheads[Small],paths},GeoBackgroundWhite,GeoRange{{-15,15},{0,360}},FrameTrue]​​]
Out[]=
Observation type by time of day:
In[]:=
obsData[GroupBy[Query["Date","Time"]]/*KeySelect[Not@*MissingQ],Counts/*ReverseSort,"Type"]
Out[]=
Night
Weather
22
RelativePosition
5
SunsetToMoonrise
4
MoonriseToSunset
3
5 total ›
Afternoon
Weather
3
WindDirection
1
Day
Weather
22
ZodiacalPosition
5
FirstAppearenceEast
3
MoonriseToSunrise
3
15 total ›
BeginningOfTheNight
RelativePosition
26
Weather
3
LastPartOfTheNight
RelativePosition
23
Weather
3
Mean observations per day (only for dated observations):
In[]:=
N@obsData[Counts/*Mean,"Date","BabylonianDate"]
Out[]=
2.18519
More examples are available at the end of the Format Discussion.