Analyzing Data from the Cuba Meteorite Fall of Feb 1, 2019

Source: https://community.wolfram.com/groups/-/m/t/1605478

meteoritesite=Entity["AdministrativeDivision",{"Vinales","PinarDelRio","Cuba"}]["Position"]

finalposition=GeoPosition[Join[meteoritesite[[1]],{Quantity[26000.,"Feet"],DateObject[{2019,2,1,18,21},TimeZone"UTC"]}]]

url="https://www.amsmeteors.org/members/imo_kml/view_trajectory_kml?event_id=513&event_year=2019&avg_start_height=80000&avg_end_height=41274.569406897&impact_location_lon=-83.947826518033&impact_location_lat=22.6421042042&end_visible_location_lon=-83.216019070606&end_visible_location_lat=23.60842985294&start_visible_location_lon=-82.505386608356&start_visible_location_lat=24.502650658255";alldata=Quiet[Import[url,"KML"]]

kmldata=Quiet[Import[url,{"KML","Data"}]];

"LayerName"/.kmldata

trajectorydata=Select[kmldata,!FreeQ[#,"Trajectory"]&]

pointsofinterest=Rule@@@Transpose[({"PlacemarkNames","Geometry"}/.trajectorydata)[[1]]]

geometricimpactpoint=("Geometric Impact Point"/.pointsofinterest)[[1]];

path=SortBy[Append[Cases[pointsofinterest,_Point,Infinity][[All,1]],geometricimpactpoint],-#[[1,2]]&];

estimatedtrajectory=GeoPosition[DeleteDuplicates[SortBy[("Estimated Trajectory"/.pointsofinterest)[[1,1]],-#[[2]]&]]];

GeoGraphics[{GeoMarker[finalposition],Blue,Line[path],GeoMarker[GeoPosition[geometricimpactpoint],"Color"Blue],Red,Arrow[estimatedtrajectory]},GeoRangeQuantity[200,"Miles"]]

GeoGraphics[{Red,Arrow[estimatedtrajectory],GeoMarker[finalposition],Orange,GeoVisibleRegion[finalposition]},GeoRangeQuantity[300,"Miles"]]

texture=GeoGraphics[{meteoritesite},GeoRangeQuantity[200,"Kilometers"],GeoGridLinesQuantity[20,"ArcMinutes"],ImageSize1920];Graphics3D[{ParametricPlot3D[{x,y,0},{x,-200,200},{y,-200,200},MeshNone,PlotStyleTexture[texture]][[1]],Red,PointSize[.01],Point[{0,0,QuantityMagnitude[Quantity[26000.,"Feet"],"Kilometers"]}],Line[{{0,0,QuantityMagnitude[Quantity[26000.,"Feet"],"Kilometers"]},{0,0,0}}]},Lighting"Neutral",BoxedFalse,AxesFalse,SphericalRegionTrue,ViewAnglePi/100,ImageSize550,ViewPoint{0.94,-3.14,0.83},ViewVertical{0,0,1}]

GeoGraphics[{RGBColor[1,0,0],GeoMarker[finalposition],Arrowheads[Medium],Arrow[estimatedtrajectory]},GeoRange{{25,28},{-100,-70}},GeoProjection{"VerticalPerspective","Centering"{30,-40}},GeoGridLinesQuantity[5,"AngularDegrees"],GeoGridLinesStyleDirective[Opacity[0.5`],GrayLevel[1]],GeoRangePaddingFull,GeoBackgroundGeoStyling["ReliefMap"],BackgroundHue[0.6`,1,0.2`],PlotRangePadding{{0.05`,0},{0.03`,0}}]