dataset : <\/p>\n\n\n\n
photos with timestamps<\/p>\n\n\n\n
UBX from M2 and from base station<\/p>\n\n\n\n
<\/p>\n\n\n\n
Processing:<\/p>\n\n\n\n
1- Post-Process base station and get corrected coordinates<\/p>\n\n\n\n
2- with Emlid Studio (or RTKLib) calculate events (triggers) positions from UBX<\/p>\n\n\n\n
as output one gets an …._events.pos file that gives corrected coordinates for each event together with a timestamp<\/p>\n\n\n\n
3- extract timestamps from exif of pictures<\/p>\n\n\n\n
use exiftool with this command<\/p>\n\n\n\n
exiftool -n -csv -filemodifydate -api TimeZone=UTC *.JPG > metadata.csv<\/p>\n\n\n\n
the output file “metadata.csv” gives the name of the picture and the timestamp in UTC format<\/p>\n\n\n\n
4- Merge …_events.pos and pictures events<\/p>\n\n\n\n
For now I’ve got this ugly Python script:<\/p>\n\n\n\n
#! \/usr\/bin\/env python\n\n\"\"\"\nUpdate Emlid Reach Survey points with PPK position output from RTKLIB\nDavid Shean\ndshean@gmail.com\nEdited to fix Pandas datetime\/Timestamp tz issues, and a few key changes likely based on Emlid updates\n\"\"\"\n\nimport os\nimport argparse\nimport numpy as np\nimport pandas as pd\n\ndef getparser():\n parser = argparse.ArgumentParser(description='Update Emlid Reach Survey points with \\\n PPK positions from RTKLIB')\n parser.add_argument('survey_pts_csv_fn', type=str, help='Survey point csv filename')\n parser.add_argument('ppk_pos_fn', type=str, help='PPK pos filename')\n return parser\n\n\ndef main():\n parser = getparser()\n args = parser.parse_args()\n\n survey_pts_csv_fn = args.survey_pts_csv_fn\n ppk_pos_fn = args.ppk_pos_fn\n\n print('Loading: %s' % survey_pts_csv_fn)\n survey_pts = pd.read_csv(survey_pts_csv_fn, parse_dates=[1], header=0)\n survey_pts['date']=pd.to_datetime(survey_pts['FileModifyDate'],format=\"%Y:%m:%d %H:%M:%S+00:00\")\n survey_pts.sort_values('date', inplace=True)\n survey_pts.index=survey_pts['date']\n print(survey_pts.dtypes)\n print(survey_pts)\n header = 'Date UTC latitude(deg) longitude(deg) height(m) Q ns sdn(m) sde(m) sdu(m) sdne(m) sdeu(m) sdun(m) age(s) ratio'\n print('Loading: %s' % ppk_pos_fn)\n ppk_pos = pd.read_csv(ppk_pos_fn, comment='%', delim_whitespace=True, names=header.split(), parse_dates=[[0,1]])\n ppk_pos['date']=pd.to_datetime(ppk_pos['Date_UTC'])\n ppk_pos.index=ppk_pos['Date_UTC']\n print(ppk_pos.dtypes)\n print(ppk_pos)\n\n # Applying merge_asof on data and store it\n # in a variable\n merged_dataframe = pd.merge_asof(ppk_pos, survey_pts, right_index=True,left_index=True,direction='nearest',tolerance=pd.Timedelta(\"1s\"))\n\n print(merged_dataframe)\n\n #Write out new file\n out_fn = os.path.splitext(survey_pts_csv_fn)[0]+'_merged.csv'\n print(\"Writing out: %s\" % out_fn)\n merged_dataframe.to_csv(out_fn) \n\nif __name__ == \"__main__\":\n main()<\/code><\/pre>\n\n\n\n5- In Metashape, import the coordinates of the cameras from the new file<\/p>\n","protected":false},"excerpt":{"rendered":"
dataset : photos with timestamps UBX from M2 and from base station Processing: 1- Post-Process base station and get corrected coordinates 2- with Emlid Studio (or RTKLib) calculate events (triggers) positions from UBX as output one gets an …._events.pos file that gives corrected coordinates for each event together with a timestamp 3- extract timestamps from … Continue reading “Use Reach M2 with Phantom 4 Pro”<\/span><\/a><\/p>\n","protected":false},"author":3,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[24,8,1],"tags":[],"class_list":["post-711","post","type-post","status-publish","format-standard","hentry","category-gps","category-photogrammetrie","category-sigeo"],"jetpack_featured_media_url":"","_links":{"self":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/711","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=711"}],"version-history":[{"count":4,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/711\/revisions"}],"predecessor-version":[{"id":719,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=\/wp\/v2\/posts\/711\/revisions\/719"}],"wp:attachment":[{"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=711"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=711"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/sigeo.cerege.fr\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=711"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}